-
Sample-Efficient Generation of Novel Photo-acid Generator Molecules using a Deep Generative Model
Authors:
Samuel C. Hoffman,
Vijil Chenthamarakshan,
Dmitry Yu. Zubarev,
Daniel P. Sanders,
Payel Das
Abstract:
Photo-acid generators (PAGs) are compounds that release acids ($H^+$ ions) when exposed to light. These compounds are critical components of the photolithography processes that are used in the manufacture of semiconductor logic and memory chips. The exponential increase in the demand for semiconductors has highlighted the need for discovering novel photo-acid generators. While de novo molecule des…
▽ More
Photo-acid generators (PAGs) are compounds that release acids ($H^+$ ions) when exposed to light. These compounds are critical components of the photolithography processes that are used in the manufacture of semiconductor logic and memory chips. The exponential increase in the demand for semiconductors has highlighted the need for discovering novel photo-acid generators. While de novo molecule design using deep generative models has been widely employed for drug discovery and material design, its application to the creation of novel photo-acid generators poses several unique challenges, such as lack of property labels. In this paper, we highlight these challenges and propose a generative modeling approach that utilizes conditional generation from a pre-trained deep autoencoder and expert-in-the-loop techniques. The validity of the proposed approach was evaluated with the help of subject matter experts, indicating the promise of such an approach for applications beyond the creation of novel photo-acid generators.
△ Less
Submitted 2 December, 2021;
originally announced December 2021.
-
Molecule Generation Experience: An Open Platform of Material Design for Public Users
Authors:
Seiji Takeda,
Toshiyuki Hama,
Hsiang-Han Hsu,
Akihiro Kishimoto,
Makoto Kogoh,
Takumi Hongo,
Kumiko Fujieda,
Hideaki Nakashika,
Dmitry Zubarev,
Daniel P. Sanders,
Jed W. Pitera,
Junta Fuchiwaki,
Daiju Nakano
Abstract:
Artificial Intelligence (AI)-driven material design has been attracting great attentions as a groundbreaking technology across a wide spectrum of industries. Molecular design is particularly important owing to its broad application domains and boundless creativity attributed to progresses in generative models. The recent maturity of molecular generative models has stimulated expectations for pract…
▽ More
Artificial Intelligence (AI)-driven material design has been attracting great attentions as a groundbreaking technology across a wide spectrum of industries. Molecular design is particularly important owing to its broad application domains and boundless creativity attributed to progresses in generative models. The recent maturity of molecular generative models has stimulated expectations for practical use among potential users, who are not necessarily familiar with coding or scripting, such as experimental engineers and students in chemical domains. However, most of the existing molecular generative models are Python libraries on GitHub, that are accessible for only IT-savvy users. To fill this gap, we newly developed a graphical user interface (GUI)-based web application of molecular generative models, Molecule Generation Experience, that is open to the general public. This is the first web application of molecular generative models enabling users to work with built-in datasets to carry out molecular design. In this paper, we describe the background technology extended from our previous work. Our new online evaluation and structural filtering algorithms significantly improved the generation speed by 30 to 1,000 times with a wider structural variety, satisfying chemical stability and synthetic reality. We also describe in detail our Kubernetes-based scalable cloud architecture and user-oriented GUI that are necessary components to achieve a public service. Finally, we present actual use cases in industrial research to design new photoacid generators (PAGs) as well as release cases in educational events.
△ Less
Submitted 6 August, 2021;
originally announced August 2021.
-
Mean encounter times for multiple random walkers on networks
Authors:
Alejandro P. Riascos,
David P. Sanders
Abstract:
We introduce a general approach for the study of the collective dynamics of non-interacting random walkers on connected networks. We analyze the movement of $R$ independent (Markovian) walkers, each defined by its own transition matrix. By using the eigenvalues and eigenvectors of the $R$ independent transition matrices, we deduce analytical expressions for the collective stationary distribution a…
▽ More
We introduce a general approach for the study of the collective dynamics of non-interacting random walkers on connected networks. We analyze the movement of $R$ independent (Markovian) walkers, each defined by its own transition matrix. By using the eigenvalues and eigenvectors of the $R$ independent transition matrices, we deduce analytical expressions for the collective stationary distribution and the average number of steps needed by the random walkers to start in a particular configuration and reach specific nodes the first time (mean first-passage times), as well as global times that characterize the global activity. We apply these results to the study of mean first-encounter times for local and non-local random walk strategies on different types of networks, with both synchronous and asynchronous motion.
△ Less
Submitted 22 February, 2021; v1 submitted 28 August, 2020;
originally announced August 2020.
-
Molecular Inverse-Design Platform for Material Industries
Authors:
Seiji Takeda,
Toshiyuki Hama,
Hsiang-Han Hsu,
Victoria A. Piunova,
Dmitry Zubarev,
Daniel P. Sanders,
Jed W. Pitera,
Makoto Kogoh,
Takumi Hongo,
Yenwei Cheng,
Wolf Bocanett,
Hideaki Nakashika,
Akihiro Fujita,
Yuta Tsuchiya,
Katsuhiko Hino,
Kentaro Yano,
Shuichi Hirose,
Hiroki Toda,
Yasumitsu Orii,
Daiju Nakano
Abstract:
The discovery of new materials has been the essential force which brings a discontinuous improvement to industrial products' performance. However, the extra-vast combinatorial design space of material structures exceeds human experts' capability to explore all, thereby hampering material development. In this paper, we present a material industry-oriented web platform of an AI-driven molecular inve…
▽ More
The discovery of new materials has been the essential force which brings a discontinuous improvement to industrial products' performance. However, the extra-vast combinatorial design space of material structures exceeds human experts' capability to explore all, thereby hampering material development. In this paper, we present a material industry-oriented web platform of an AI-driven molecular inverse-design system, which automatically designs brand new molecular structures rapidly and diversely. Different from existing inverse-design solutions, in this system, the combination of substructure-based feature encoding and molecular graph generation algorithms allows a user to gain high-speed, interpretable, and customizable design process. Also, a hierarchical data structure and user-oriented UI provide a flexible and intuitive workflow. The system is deployed on IBM's and our client's cloud servers and has been used by 5 partner companies. To illustrate actual industrial use cases, we exhibit inverse-design of sugar and dye molecules, that were carried out by experimental chemists in those client companies. Compared to general human chemist's standard performance, the molecular design speed was accelerated more than 10 times, and greatly increased variety was observed in the inverse-designed molecules without loss of chemical realism.
△ Less
Submitted 16 May, 2020; v1 submitted 23 April, 2020;
originally announced April 2020.
-
AI-driven Inverse Design System for Organic Molecules
Authors:
Seiji Takeda,
Toshiyuki Hama,
Hsiang-Han Hsu,
Toshiyuki Yamane,
Koji Masuda,
Victoria A. Piunova,
Dmitry Zubarev,
Jed Pitera,
Daniel P. Sanders,
Daiju Nakano
Abstract:
Designing novel materials that possess desired properties is a central need across many manufacturing industries. Driven by that industrial need, a variety of algorithms and tools have been developed that combine AI (machine learning and analytics) with domain knowledge in physics, chemistry, and materials science. AI-driven materials design can be divided to mainly two stages; the first one is th…
▽ More
Designing novel materials that possess desired properties is a central need across many manufacturing industries. Driven by that industrial need, a variety of algorithms and tools have been developed that combine AI (machine learning and analytics) with domain knowledge in physics, chemistry, and materials science. AI-driven materials design can be divided to mainly two stages; the first one is the modeling stage, where the goal is to build an accurate regression or classification model to predict material properties (e.g. glass transition temperature) or attributes (e.g. toxic/non-toxic). The next stage is design, where the goal is to assemble or tune material structures so that they can achieve user-demanded target property values based on a prediction model that is trained in the modeling stage. For maximum benefit, these two stages should be architected to form a coherent workflow. Today there are several emerging services and tools for AI-driven material design, however, most of them provide only partial technical components (e.g. data analyzer, regression model, structure generator, etc.), that are useful for specific purposes, but for comprehensive material design, those components need to be orchestrated appropriately. Our material design system provides an end-to-end solution to this problem, with a workflow that consists of data input, feature encoding, prediction modeling, solution search, and structure generation. The system builds a regression model to predict properties, solves an inverse problem on the trained model, and generates novel chemical structure candidates that satisfy the target properties. In this paper we will introduce the methodology of our system, and demonstrate a simple example of inverse design generating new chemical structures that satisfy targeted physical property values.
△ Less
Submitted 20 January, 2020;
originally announced January 2020.
-
Exact hopping and collision times for two hard discs in a box
Authors:
Wilhelm P. K. Zapfe,
David P. Sanders,
Rosa Rodríguez-Mota
Abstract:
We study the molecular dynamics of two discs undergoing Newtonian ("inertial") dynamics, with elastic collisions in a rectangular box. Using a mapping to a billiard model and a key result from ergodic theory, we obtain exact, analytical expressions for the mean times between the following events: hops, i.e.~horizontal or vertical interchanges of the particles; wall collisions; and disc collisions.…
▽ More
We study the molecular dynamics of two discs undergoing Newtonian ("inertial") dynamics, with elastic collisions in a rectangular box. Using a mapping to a billiard model and a key result from ergodic theory, we obtain exact, analytical expressions for the mean times between the following events: hops, i.e.~horizontal or vertical interchanges of the particles; wall collisions; and disc collisions. To do so, we calculate volumes and cross-sectional areas in the four-dimensional configuration space. We compare the analytical results against Monte Carlo and molecular dynamics simulations, with excellent agreement.
△ Less
Submitted 13 August, 2019;
originally announced August 2019.
-
Monte Carlo sampling in diffusive dynamical systems
Authors:
Diego Tapias,
David P. Sanders,
Eduardo G. Altmann
Abstract:
We introduce a Monte Carlo algorithm to efficiently compute transport properties of chaotic dynamical systems. Our method exploits the importance sampling technique that favors trajectories in the tail of the distribution of displacements, where deviations from a diffusive process are most prominent. We search for initial conditions using a proposal that correlates states in the Markov chain const…
▽ More
We introduce a Monte Carlo algorithm to efficiently compute transport properties of chaotic dynamical systems. Our method exploits the importance sampling technique that favors trajectories in the tail of the distribution of displacements, where deviations from a diffusive process are most prominent. We search for initial conditions using a proposal that correlates states in the Markov chain constructed via a Metropolis-Hastings algorithm. We show that our method outperforms the direct sampling method and also Metropolis-Hastings methods with alternative proposals. We test our general method through numerical simulations in 1D (box-map) and 2D (Lorentz gas) systems.
△ Less
Submitted 4 April, 2018;
originally announced April 2018.
-
Ergodicity of one-dimensional systems coupled to the logistic thermostat
Authors:
Diego Tapias,
Alessandro Bravetti,
David P. Sanders
Abstract:
We analyze the ergodicity of three one-dimensional Hamiltonian systems, with harmonic, quartic and Mexican-hat potentials, coupled to the logistic thermostat. As criteria for ergodicity we employ: the independence of the Lyapunov spectrum with respect to initial conditions; the absence of visual "holes" in two-dimensional Poincaré sections; the agreement between the histograms in each variable and…
▽ More
We analyze the ergodicity of three one-dimensional Hamiltonian systems, with harmonic, quartic and Mexican-hat potentials, coupled to the logistic thermostat. As criteria for ergodicity we employ: the independence of the Lyapunov spectrum with respect to initial conditions; the absence of visual "holes" in two-dimensional Poincaré sections; the agreement between the histograms in each variable and the theoretical marginal distributions; and the convergence of the global joint distribution to the theoretical one, as measured by the Hellinger distance. Taking a large number of random initial conditions, for certain parameter values of the thermostat we find no indication of regular trajectories and show that the time distribution converges to the ensemble one for an arbitrarily long trajectory for all the systems considered. Our results thus provide a robust numerical indication that the logistic thermostat can serve as a single one-parameter thermostat for stiff one-dimensional systems.
△ Less
Submitted 15 November, 2016;
originally announced November 2016.
-
Geometric integrator for simulations in the canonical ensemble
Authors:
Diego Tapias,
David P. Sanders,
Alessandro Bravetti
Abstract:
In this work we introduce a geometric integrator for molecular dynamics simulations of physical systems in the canonical ensemble. In particular, we consider the equations arising from the so-called density dynamics algorithm with any possible type of thermostat and provide an integrator that preserves the invariant distribution. Our integrator thus constitutes a unified framework that allows the…
▽ More
In this work we introduce a geometric integrator for molecular dynamics simulations of physical systems in the canonical ensemble. In particular, we consider the equations arising from the so-called density dynamics algorithm with any possible type of thermostat and provide an integrator that preserves the invariant distribution. Our integrator thus constitutes a unified framework that allows the study and comparison of different thermostats and of their influence on the equilibrium and non-equilibrium (thermo-)dynamic properties of the system. To show the validity and the generality of the integrator, we implement it with a second-order, time-reversible method and apply it to the simulation of a Lennard-Jones system with three different thermostats, obtaining good conservation of the geometrical properties and recovering the expected thermodynamic results.
△ Less
Submitted 5 May, 2016;
originally announced May 2016.
-
Horizons and free path distributions in quasiperiodic Lorentz gases
Authors:
Atahualpa S. Kraemer,
Michael Schmiedeberg,
David P. Sanders
Abstract:
We study the structure of quasiperiodic Lorentz gases, i.e., particles bouncing elastically off fixed obstacles arranged in quasiperiodic lattices. By employing a construction to embed such structures into a higher dimensional periodic hyperlattice, we give a simple and efficient algorithm for numerical simulation of the dynamics of these systems. This same construction shows that quasiperiodic Lo…
▽ More
We study the structure of quasiperiodic Lorentz gases, i.e., particles bouncing elastically off fixed obstacles arranged in quasiperiodic lattices. By employing a construction to embed such structures into a higher dimensional periodic hyperlattice, we give a simple and efficient algorithm for numerical simulation of the dynamics of these systems. This same construction shows that quasiperiodic Lorentz gases generically exhibit a regime with infinite horizon, that is, empty channels through which the particles move without colliding, when the obstacles are small enough; in this case, the distribution of free paths is asymptotically a power law with exponent -3, as expected from infinite-horizon periodic Lorentz gases. For the critical radius at which these channels disappear, however, a new regime with locally-finite horizon arises, where this distribution has an unexpected exponent of -5, previously observed only in a Lorentz gas formed by superposing three incommensurable periodic lattices in the Boltzmann-Grad limit where the radius of the obstacles tends to zero.
△ Less
Submitted 1 November, 2015;
originally announced November 2015.
-
Efficient algorithms for general periodic Lorentz gases in two and three dimensions
Authors:
Atahualpa S. Kraemer,
Nikolay Kryukov,
David P. Sanders
Abstract:
We present efficient algorithms to calculate trajectories for periodic Lorentz gases consisting of square lattices of circular obstacles in two dimensions, and simple cubic lattices of spheres in three dimensions; these become increasingly efficient as the radius of the obstacles tends to 0, the so-called Boltzmann-Grad limit. The 2D algorithm applies continued fractions to obtain the exact disc w…
▽ More
We present efficient algorithms to calculate trajectories for periodic Lorentz gases consisting of square lattices of circular obstacles in two dimensions, and simple cubic lattices of spheres in three dimensions; these become increasingly efficient as the radius of the obstacles tends to 0, the so-called Boltzmann-Grad limit. The 2D algorithm applies continued fractions to obtain the exact disc with which a particle will collide at each step, instead of using periodic boundary conditions as in the classical algorithm. The 3D version incorporates the 2D algorithm by projecting to the three coordinate planes. As an application, we calculate distributions of free path lengths close to the Boltzmann-Grad limit for certain Lorentz gases. We also show how the algorithms may be applied to deal with general crystal lattices.
△ Less
Submitted 1 November, 2015;
originally announced November 2015.
-
Lévy walks on lattices as multi-state processes
Authors:
Giampaolo Cristadoro,
Thomas Gilbert,
Marco Lenci,
David P. Sanders
Abstract:
Continuous-time random walks combining diffusive scattering and ballistic propagation on lattices model a class of Lévy walks. The assumption that transitions in the scattering phase occur with exponentially-distributed waiting times leads to a description of the process in terms of multiple states, whose distributions evolve according to a set of delay differential equations, amenable to analytic…
▽ More
Continuous-time random walks combining diffusive scattering and ballistic propagation on lattices model a class of Lévy walks. The assumption that transitions in the scattering phase occur with exponentially-distributed waiting times leads to a description of the process in terms of multiple states, whose distributions evolve according to a set of delay differential equations, amenable to analytic treatment. We obtain an exact expression of the mean squared displacement associated with such processes and discuss the emergence of asymptotic scaling laws in regimes of diffusive and superdiffusive (subballistic) transport, emphasizing, in the latter case, the effect of initial conditions on the transport coefficients. Of particular interest is the case of rare ballistic propagation, in which case a regime of superdiffusion may lurk underneath one of normal diffusion.
△ Less
Submitted 21 January, 2015;
originally announced January 2015.
-
Three-edge-colouring doublecross cubic graphs
Authors:
Katherine Edwards,
Daniel P. Sanders,
Paul Seymour,
Robin Thomas
Abstract:
A graph is apex if there is a vertex whose deletion makes the graph planar, and doublecross if it can be drawn in the plane with only two crossings, both incident with the infinite region in the natural sense. In 1966, Tutte conjectured that every two-edge-connected cubic graph with no Petersen graph minor is three-edge-colourable. With Neil Robertson, two of us showed that this is true in general…
▽ More
A graph is apex if there is a vertex whose deletion makes the graph planar, and doublecross if it can be drawn in the plane with only two crossings, both incident with the infinite region in the natural sense. In 1966, Tutte conjectured that every two-edge-connected cubic graph with no Petersen graph minor is three-edge-colourable. With Neil Robertson, two of us showed that this is true in general if it is true for apex graphs and doublecross graphs. In another paper, two of us solved the apex case, but the doublecross case remained open. Here we solve the doublecross case; that is, we prove that every two-edge-connected doublecross cubic graph is three-edge-colourable. The proof method is a variant on the proof of the four-colour theorem.
△ Less
Submitted 30 November, 2015; v1 submitted 16 November, 2014;
originally announced November 2014.
-
Machta-Zwanzig regime of anomalous diffusion in infinite-horizon billiards
Authors:
Giampaolo Cristadoro,
Thomas Gilbert,
Marco Lenci,
David P. Sanders
Abstract:
We study diffusion on a periodic billiard table with infinite horizon in the limit of narrow corridors. An effective trapping mechanism emerges according to which the process can be modeled by a Lévy walk combining exponentially-distributed trapping times with free propagation along paths whose precise probabilities we compute. This description yields an approximation of the mean squared displacem…
▽ More
We study diffusion on a periodic billiard table with infinite horizon in the limit of narrow corridors. An effective trapping mechanism emerges according to which the process can be modeled by a Lévy walk combining exponentially-distributed trapping times with free propagation along paths whose precise probabilities we compute. This description yields an approximation of the mean squared displacement of infinite-horizon billiards in terms of two transport coefficients which generalizes to this anomalous regime the Machta-Zwanzig approximation of normal diffusion in finite-horizon billiards [Phys. Rev. Lett. 50, 1959 (1983)].
△ Less
Submitted 2 August, 2014;
originally announced August 2014.
-
Transport properties of Lévy walks: an analysis in terms of multistate processes
Authors:
Giampaolo Cristadoro,
Thomas Gilbert,
Marco Lenci,
David P. Sanders
Abstract:
Continuous time random walks combining diffusive and ballistic regimes are introduced to describe a class of Lévy walks on lattices. By including exponentially-distributed waiting times separating the successive jump events of a walker, we are led to a description of such Lévy walks in terms of multistate processes whose time-evolution is shown to obey a set of coupled delay differential equations…
▽ More
Continuous time random walks combining diffusive and ballistic regimes are introduced to describe a class of Lévy walks on lattices. By including exponentially-distributed waiting times separating the successive jump events of a walker, we are led to a description of such Lévy walks in terms of multistate processes whose time-evolution is shown to obey a set of coupled delay differential equations. Using simple arguments, we obtain asymptotic solutions to these equations and rederive the scaling laws for the mean squared displacement of such processes. Our calculation includes the computation of all relevant transport coefficients in terms of the parameters of the models.
△ Less
Submitted 10 November, 2014; v1 submitted 1 July, 2014;
originally announced July 2014.
-
Zero density of open paths in the Lorentz mirror model for arbitrary mirror probability
Authors:
Atahualpa S. Kraemer,
David P. Sanders
Abstract:
We show, incorporating results obtained from numerical simulations, that in the Lorentz mirror model, the density of open paths in any finite box tends to 0 as the box size tends to infinity, for any mirror probability.
We show, incorporating results obtained from numerical simulations, that in the Lorentz mirror model, the density of open paths in any finite box tends to 0 as the box size tends to infinity, for any mirror probability.
△ Less
Submitted 18 June, 2014;
originally announced June 2014.
-
Measuring logarithmic corrections to normal diffusion in infinite-horizon billiards
Authors:
Giampaolo Cristadoro,
Thomas Gilbert,
Marco Lenci,
David P. Sanders
Abstract:
We perform numerical measurements of the moments of the position of a tracer particle in a two-dimensional periodic billiard model (Lorentz gas) with infinite corridors. This model is known to exhibit a weak form of super-diffusion, in the sense that there is a logarithmic correction to the linear growth in time of the mean-squared displacement. We show numerically that this expected asymptotic be…
▽ More
We perform numerical measurements of the moments of the position of a tracer particle in a two-dimensional periodic billiard model (Lorentz gas) with infinite corridors. This model is known to exhibit a weak form of super-diffusion, in the sense that there is a logarithmic correction to the linear growth in time of the mean-squared displacement. We show numerically that this expected asymptotic behavior is easily overwhelmed by the subleading linear growth throughout the time-range accessible to numerical simulations. We compare our simulations to the known analytical results for the variance of the anomalously-rescaled limiting normal distributions.
△ Less
Submitted 6 May, 2014; v1 submitted 5 May, 2014;
originally announced May 2014.
-
Discharging cartwheels
Authors:
Neil Robertson,
Daniel P. Sanders,
Paul Seymour,
Robin Thomas
Abstract:
In [J. Combin. Theory Ser. B 70 (1997), 2-44] we gave a simplified proof of the Four-Color Theorem. The proof is computer-assisted in the sense that for two lemmas in the article we did not give proofs, and instead asserted that we have verified those statements using a computer. Here we give additional details for one of those lemmas, and we include the original computer programs and data as "anc…
▽ More
In [J. Combin. Theory Ser. B 70 (1997), 2-44] we gave a simplified proof of the Four-Color Theorem. The proof is computer-assisted in the sense that for two lemmas in the article we did not give proofs, and instead asserted that we have verified those statements using a computer. Here we give additional details for one of those lemmas, and we include the original computer programs and data as "ancillary files" accompanying this submission.
△ Less
Submitted 24 January, 2014;
originally announced January 2014.
-
Reducibility in the Four-Color Theorem
Authors:
Neil Robertson,
Daniel P. Sanders,
Paul Seymour,
Robin Thomas
Abstract:
In [J. Combin. Theory Ser. B 70 (1997), 2-44] we gave a simplified proof of the Four-Color Theorem. The proof is computer-assisted in the sense that for two lemmas in the article we did not give proofs, and instead asserted that we have verified those statements using a computer. Here we give additional details for one of those lemmas, and we include the original computer programs and data as "anc…
▽ More
In [J. Combin. Theory Ser. B 70 (1997), 2-44] we gave a simplified proof of the Four-Color Theorem. The proof is computer-assisted in the sense that for two lemmas in the article we did not give proofs, and instead asserted that we have verified those statements using a computer. Here we give additional details for one of those lemmas, and we include the original computer programs and data as "ancillary files" accompanying this submission.
△ Less
Submitted 24 January, 2014;
originally announced January 2014.
-
Emergence of a sustained current by coupling equilibrating systems: Making a NESS out of equilibrium
Authors:
Hernán Larralde,
David P. Sanders
Abstract:
We show that coupling together two closed thermodynamic systems that independently attain equilibrium may give rise to a nonequilibrium stationary state (NESS) with a persistent, non-vanishing current. We study a simple example that is exactly soluble, consisting of random walkers with different biases towards a reflecting boundary, modelling, for example, Brownian particles with different charge…
▽ More
We show that coupling together two closed thermodynamic systems that independently attain equilibrium may give rise to a nonequilibrium stationary state (NESS) with a persistent, non-vanishing current. We study a simple example that is exactly soluble, consisting of random walkers with different biases towards a reflecting boundary, modelling, for example, Brownian particles with different charge states in an electric field. We obtain exact analytical expressions for the (generating functions of) concentrations and currents in the NESS for this model, and exhibit the main features by numerical simulation.
△ Less
Submitted 21 January, 2014;
originally announced January 2014.
-
Encounter times in overlapping domains: application to epidemic spread in a population of territorial animals
Authors:
Luca Giuggioli,
Sebastian Pérez-Becker,
David P. Sanders
Abstract:
We develop an analytical method to calculate encounter times of two random walkers in one dimension when each individual is segregated in its own spatial domain and shares with its neighbor only a fraction of the available space, finding very good agreement with numerically-exact calculations. We model a population of susceptible and infected territorial individuals with this spatial arrangement,…
▽ More
We develop an analytical method to calculate encounter times of two random walkers in one dimension when each individual is segregated in its own spatial domain and shares with its neighbor only a fraction of the available space, finding very good agreement with numerically-exact calculations. We model a population of susceptible and infected territorial individuals with this spatial arrangement, and which may transmit an epidemic when they meet. We apply the results on encounter times to determine analytically the macroscopic propagation speed of the epidemic as a function of the microscopic characteristics: the confining geometry, the animal diffusion constant, and the infection transmission probability.
△ Less
Submitted 21 November, 2012; v1 submitted 10 July, 2012;
originally announced July 2012.
-
Periodizing quasicrystals: Anomalous diffusion in quasiperiodic systems
Authors:
Atahualpa S. Kraemer,
David P. Sanders
Abstract:
We introduce a construction to embed a quasiperiodic lattice of obstacles into a single unit cell of a higher-dimensional space, with periodic boundary conditions. This construction transparently shows the existence of channels in these systems,in which particles may travel without colliding, up to a critical obstacle radius. It provides a simple and efficient algorithm for numerical simulation of…
▽ More
We introduce a construction to embed a quasiperiodic lattice of obstacles into a single unit cell of a higher-dimensional space, with periodic boundary conditions. This construction transparently shows the existence of channels in these systems,in which particles may travel without colliding, up to a critical obstacle radius. It provides a simple and efficient algorithm for numerical simulation of dynamics in quasiperiodic structures, as well as giving a natural notion of uniform distribution (measure) and averages. As an application, we simulate diffusion in a two-dimensional quasicrystal, finding three different regimes, in particular atypical weak super-diffusion in the presence of channels, and sub-diffusion when obstacles overlap.
△ Less
Submitted 5 June, 2012;
originally announced June 2012.
-
Structure and evolution of strange attractors in non-elastic triangular billiards
Authors:
Aubin Arroyo,
Roberto Markarian,
David P. Sanders
Abstract:
We study pinball billiard dynamics in an equilateral triangular table. In such dynamics, collisions with the walls are non-elastic: the outgoing angle with the normal vector to the boundary is a uniform factor $λ< 1$ smaller than the incoming angle. This leads to contraction in phase space for the discrete-time dynamics between consecutive collisions, and hence to attractors of zero Lebesgue measu…
▽ More
We study pinball billiard dynamics in an equilateral triangular table. In such dynamics, collisions with the walls are non-elastic: the outgoing angle with the normal vector to the boundary is a uniform factor $λ< 1$ smaller than the incoming angle. This leads to contraction in phase space for the discrete-time dynamics between consecutive collisions, and hence to attractors of zero Lebesgue measure, which are almost always fractal strange attractors with chaotic dynamics, due to the presence of an expansion mechanism. We study the structure of these strange attractors and their evolution as the contraction parameter $λ$ is varied. For $λ$ in the interval (0, 1/3), we prove rigorously that the attractor has the structure of a Cantor set times an interval, whereas for larger values of $λ$ the billiard dynamics gives rise to nonaccessible regions in phase space. For $λ$ close to 1, the attractor splits into three transitive components, the basins of attraction of which have fractal basin boundaries.
△ Less
Submitted 6 December, 2011;
originally announced December 2011.
-
Comment on "Capturing correlations in chaotic diffusion by approximation methods"
Authors:
Thomas Gilbert,
David P. Sanders
Abstract:
This is a comment on [G. Knight and R. Klages, Phys. Rev. E 84, 041135 (2011); also available at arXiv:1107.5293v2 [math-ph]].
This is a comment on [G. Knight and R. Klages, Phys. Rev. E 84, 041135 (2011); also available at arXiv:1107.5293v2 [math-ph]].
△ Less
Submitted 27 November, 2011;
originally announced November 2011.
-
Diffusive properties of persistent walks on cubic lattices with application to periodic Lorentz gases
Authors:
Thomas Gilbert,
Huu Chuong Nguyen,
David P Sanders
Abstract:
We calculate the diffusion coefficients of persistent random walks on cubic and hypercubic lattices, where the direction of a walker at a given step depends on the memory of one or two previous steps. These results are then applied to study a billiard model, namely a three-dimensional periodic Lorentz gas. The geometry of the model is studied in order to find the regimes in which it exhibits norma…
▽ More
We calculate the diffusion coefficients of persistent random walks on cubic and hypercubic lattices, where the direction of a walker at a given step depends on the memory of one or two previous steps. These results are then applied to study a billiard model, namely a three-dimensional periodic Lorentz gas. The geometry of the model is studied in order to find the regimes in which it exhibits normal diffusion. In this regime, we calculate numerically the transition probabilities between cells to compare the persistent random-walk approximation with simulation results for the diffusion coefficient.
△ Less
Submitted 20 September, 2010;
originally announced September 2010.
-
Stable and unstable regimes in higher-dimensional convex billiards with cylindrical shape
Authors:
Thomas Gilbert,
David P. Sanders
Abstract:
We introduce a class of convex, higher-dimensional billiard models which generalise stadium billiards. These models correspond to the free motion of a point-particle in a region bounded by cylinders cut by planes. They are motivated by models of particles interacting via a string-type mechanism, and confined by hard walls. The combination of these elements may give rise to a defocusing mechanism,…
▽ More
We introduce a class of convex, higher-dimensional billiard models which generalise stadium billiards. These models correspond to the free motion of a point-particle in a region bounded by cylinders cut by planes. They are motivated by models of particles interacting via a string-type mechanism, and confined by hard walls. The combination of these elements may give rise to a defocusing mechanism, similar to that in two dimensions, which allows large chaotic regions in phase space. The remaining part of phase space is associated with marginally stable behaviour. In fact periodic orbits in these systems generically come in continuous parametric families, sociated with a pair of parabolic eigen-directions: the periodic orbits are unstable in the presence of a defocusing mechanism, but marginally stable otherwise. By performing the stability analysis of families of periodic orbits at a nonlinear level, we establish the conditions under which families are nonlinearly stable or unstable. As a result, we identify regions in the parameter space of the models which admit non-linearly stable oscillations in the form of whispering gallery modes. Where no families of periodic orbits are stable, the billiards are completely chaotic, i.e.\ the Lyapunov exponents of the billiard map are non-zero.
△ Less
Submitted 2 September, 2010;
originally announced September 2010.
-
Chaos and stability in a two-parameter family of convex billiard tables
Authors:
Péter Bálint,
Miklós Halász,
Jorge Hernández-Tahuilán,
David P. Sanders
Abstract:
We study, by numerical simulations and semi-rigorous arguments, a two-parameter family of convex, two-dimensional billiard tables, generalizing the one-parameter class of oval billiards of Benettin--Strelcyn [Phys. Rev. A 17, 773 (1978)]. We observe interesting dynamical phenomena when the billiard tables are continuously deformed from the integrable circular billiard to different versions of comp…
▽ More
We study, by numerical simulations and semi-rigorous arguments, a two-parameter family of convex, two-dimensional billiard tables, generalizing the one-parameter class of oval billiards of Benettin--Strelcyn [Phys. Rev. A 17, 773 (1978)]. We observe interesting dynamical phenomena when the billiard tables are continuously deformed from the integrable circular billiard to different versions of completely-chaotic stadia. In particular, we conjecture that a new class of ergodic billiard tables is obtained in certain regions of the two-dimensional parameter space, when the billiards are close to skewed stadia. We provide heuristic arguments supporting this conjecture, and give numerical confirmation using the powerful method of Lyapunov-weighted dynamics.
△ Less
Submitted 15 June, 2010;
originally announced June 2010.
-
Chaos in cylindrical stadium billiards via a generic nonlinear mechanism
Authors:
Thomas Gilbert,
David P. Sanders
Abstract:
We describe conditions under which higher-dimensional billiard models in bounded, convex regions are fully chaotic, generalizing the Bunimovich stadium to dimensions above two. An example is a three-dimensional stadium bounded by a cylinder and several planes; the combination of these elements may give rise to defocusing, allowing large chaotic regions in phase space. By studying families of margi…
▽ More
We describe conditions under which higher-dimensional billiard models in bounded, convex regions are fully chaotic, generalizing the Bunimovich stadium to dimensions above two. An example is a three-dimensional stadium bounded by a cylinder and several planes; the combination of these elements may give rise to defocusing, allowing large chaotic regions in phase space. By studying families of marginally-stable periodic orbits that populate the residual part of phase space, we identify conditions under which a nonlinear instability mechanism arises in their vicinity. For particular geometries, this mechanism rather induces stable nonlinear oscillations, including in the form of whispering-gallery modes.
△ Less
Submitted 3 September, 2010; v1 submitted 28 August, 2009;
originally announced August 2009.
-
Diffusion coefficients for multi-step persistent random walks on lattices
Authors:
Thomas Gilbert,
David P. Sanders
Abstract:
We calculate the diffusion coefficients of persistent random walks on lattices, where the direction of a walker at a given step depends on the memory of a certain number of previous steps. In particular, we describe a simple method which enables us to obtain explicit expressions for the diffusion coefficients of walks with two-step memory on different classes of one-, two- and higher-dimensional…
▽ More
We calculate the diffusion coefficients of persistent random walks on lattices, where the direction of a walker at a given step depends on the memory of a certain number of previous steps. In particular, we describe a simple method which enables us to obtain explicit expressions for the diffusion coefficients of walks with two-step memory on different classes of one-, two- and higher-dimensional lattices.
△ Less
Submitted 10 August, 2009;
originally announced August 2009.
-
Persistence effects in deterministic diffusion
Authors:
Thomas Gilbert,
David P. Sanders
Abstract:
In systems which exhibit deterministic diffusion, the gross parameter dependence of the diffusion coefficient can often be understood in terms of random walk models. Provided the decay of correlations is fast enough, one can ignore memory effects and approximate the diffusion coefficient according to dimensional arguments. By successively including the effects of one and two steps of memory on t…
▽ More
In systems which exhibit deterministic diffusion, the gross parameter dependence of the diffusion coefficient can often be understood in terms of random walk models. Provided the decay of correlations is fast enough, one can ignore memory effects and approximate the diffusion coefficient according to dimensional arguments. By successively including the effects of one and two steps of memory on this approximation, we examine the effects of ``persistence'' on the diffusion coefficients of extended two-dimensional billiard tables and show how to properly account for these effects, using walks in which a particle undergoes jumps in different directions with probabilities that depend on where they came from.
△ Less
Submitted 5 August, 2009;
originally announced August 2009.
-
Exact encounter times for many random walkers on regular and complex networks
Authors:
David P. Sanders
Abstract:
The exact mean time between encounters of a given particle in a system consisting of many particles undergoing random walks in discrete time is calculated, on both regular and complex networks. Analytical results are obtained both for independent walkers, where any number of walkers can occupy the same site, and for walkers with an exclusion interaction, when no site can contain more than one wa…
▽ More
The exact mean time between encounters of a given particle in a system consisting of many particles undergoing random walks in discrete time is calculated, on both regular and complex networks. Analytical results are obtained both for independent walkers, where any number of walkers can occupy the same site, and for walkers with an exclusion interaction, when no site can contain more than one walker. These analytical results are then compared with numerical simulations, showing very good agreement.
△ Less
Submitted 3 June, 2009;
originally announced June 2009.
-
Long-range correlations in a simple stochastic model of coupled transport
Authors:
Hernán Larralde,
David P. Sanders
Abstract:
We study coupled transport in the nonequilibrium stationary state of a model consisting of independent random walkers, moving along a one-dimensional channel, which carry a conserved energy-like quantity, with density and temperature gradients imposed by reservoirs at the ends of the channel. In our model, walkers interact with other walkers at the same site by sharing energy at each time step,…
▽ More
We study coupled transport in the nonequilibrium stationary state of a model consisting of independent random walkers, moving along a one-dimensional channel, which carry a conserved energy-like quantity, with density and temperature gradients imposed by reservoirs at the ends of the channel. In our model, walkers interact with other walkers at the same site by sharing energy at each time step, but the amount of energy carried does not affect the motion of the walkers. We find that already in this simple model long-range correlations arise in the nonequilibrium stationary state which are similar to those observed in more realistic models of coupled transport.
△ Less
Submitted 11 August, 2009; v1 submitted 18 March, 2009;
originally announced March 2009.
-
Bifurcations of periodic and chaotic attractors in pinball billiards with focusing boundaries
Authors:
Aubin Arroyo,
Roberto Markarian,
David P. Sanders
Abstract:
We study the dynamics of billiard models with a modified collision rule: the outgoing angle from a collision is a uniform contraction, by a factor lambda, of the incident angle. These pinball billiards interpolate between a one-dimensional map when lambda=0 and the classical Hamiltonian case of elastic collisions when lambda=1. For all lambda<1, the dynamics is dissipative, and thus gives rise t…
▽ More
We study the dynamics of billiard models with a modified collision rule: the outgoing angle from a collision is a uniform contraction, by a factor lambda, of the incident angle. These pinball billiards interpolate between a one-dimensional map when lambda=0 and the classical Hamiltonian case of elastic collisions when lambda=1. For all lambda<1, the dynamics is dissipative, and thus gives rise to attractors, which may be periodic or chaotic. Motivated by recent rigorous results of Markarian, Pujals and Sambarino, we numerically investigate and characterise the bifurcations of the resulting attractors as the contraction parameter is varied. Some billiards exhibit only periodic attractors, some only chaotic attractors, and others have coexistence of the two types.
△ Less
Submitted 11 June, 2009; v1 submitted 9 February, 2009;
originally announced February 2009.
-
Deterministic Diffusion in Periodic Billiard Models
Authors:
David P. Sanders
Abstract:
We investigate statistical properties of several classes of periodic billiard models which are diffusive. An introductory chapter gives motivation, and then a review of statistical properties of dynamical systems is given in chapter 2. In chapter 3, we study the geometry dependence of diffusion coefficients in a two-parameter 2D periodic Lorentz gas model, including a discussion of how to estima…
▽ More
We investigate statistical properties of several classes of periodic billiard models which are diffusive. An introductory chapter gives motivation, and then a review of statistical properties of dynamical systems is given in chapter 2. In chapter 3, we study the geometry dependence of diffusion coefficients in a two-parameter 2D periodic Lorentz gas model, including a discussion of how to estimate them from data. In chapter 4, we study the shape of position and displacement distributions, which occur in the central limit theorem. We show that there is an oscillatory fine structure and what its origin is. This allows us to conjecture a refinement of the central limit theorem in these systems. A non-Maxwellian velocity distribution is shown to lead to a non-Gaussian limit distribution. Chapter 5 treats polygonal billiard channels, developing a picture of when normal and anomalous diffusion occur, the latter being due to parallel scatterers in the billiard causing a channelling effect. We also characterize the crossover from normal to anomalous diffusion. In chapter 6, we extend our methods to a 3D periodic Lorentz gas model, showing that normal diffusion occurs under certain conditions. In particular, we construct an explicit finite-horizon model, and we discuss the effect that holes in configuration space have on the diffusive properties of the system. We finish with conclusions and directions for future research.
△ Less
Submitted 16 August, 2008;
originally announced August 2008.
-
Normal diffusion in crystal structures and higher-dimensional billiard models with gaps
Authors:
David P. Sanders
Abstract:
We show, both heuristically and numerically, that three-dimensional periodic Lorentz gases -- clouds of particles scattering off crystalline arrays of hard spheres -- often exhibit normal diffusion, even when there are gaps through which particles can travel without ever colliding, i.e., when the system has an infinite horizon. This is the case provided that these gaps are not "too big", as meas…
▽ More
We show, both heuristically and numerically, that three-dimensional periodic Lorentz gases -- clouds of particles scattering off crystalline arrays of hard spheres -- often exhibit normal diffusion, even when there are gaps through which particles can travel without ever colliding, i.e., when the system has an infinite horizon. This is the case provided that these gaps are not "too big", as measured by their dimension. The results are illustrated with simulations of a simple three-dimensional model having different types of diffusive regime, and are then extended to higher-dimensional billiard models, which include hard-sphere fluids.
△ Less
Submitted 26 January, 2009; v1 submitted 16 August, 2008;
originally announced August 2008.
-
How rare are diffusive rare events?
Authors:
David P. Sanders,
Hernán Larralde
Abstract:
We study the time until first occurrence, the first-passage time, of rare density fluctuations in diffusive systems. We approach the problem using a model consisting of many independent random walkers on a lattice. The existence of spatial correlations makes this problem analytically intractable. However, for a mean-field approximation in which the walkers can jump anywhere in the system, we obt…
▽ More
We study the time until first occurrence, the first-passage time, of rare density fluctuations in diffusive systems. We approach the problem using a model consisting of many independent random walkers on a lattice. The existence of spatial correlations makes this problem analytically intractable. However, for a mean-field approximation in which the walkers can jump anywhere in the system, we obtain a simple asymptotic form for the mean first-passage time to have a given number k of particles at a distinguished site. We show numerically, and argue heuristically, that for large enough k, the mean-field results give a good approximation for first-passage times for systems with nearest-neighbour dynamics, especially for two and higher spatial dimensions. Finally, we show how the results change when density fluctuations anywhere in the system, rather than at a specific distinguished site, are considered.
△ Less
Submitted 7 April, 2008;
originally announced April 2008.
-
Competitive nucleation and the Ostwald rule in a generalized Potts model with multiple metastable phases
Authors:
David P. Sanders,
Hernán Larralde,
François Leyvraz
Abstract:
We introduce a simple nearest-neighbor spin model with multiple metastable phases, the number and decay pathways of which are explicitly controlled by the parameters of the system. With this model we can construct, for example, a system which evolves through an arbitrarily long succession of metastable phases. We also construct systems in which different phases may nucleate competitively from a…
▽ More
We introduce a simple nearest-neighbor spin model with multiple metastable phases, the number and decay pathways of which are explicitly controlled by the parameters of the system. With this model we can construct, for example, a system which evolves through an arbitrarily long succession of metastable phases. We also construct systems in which different phases may nucleate competitively from a single initial phase. For such a system, we present a general method to extract from numerical simulations the individual nucleation rates of the nucleating phases. The results show that the Ostwald rule, which predicts which phase will nucleate, must be modified probabilistically when the new phases are almost equally stable. Finally, we show that the nucleation rate of a phase depends, among other things, on the number of other phases accessible from it.
△ Less
Submitted 9 April, 2007; v1 submitted 3 April, 2007;
originally announced April 2007.
-
Metastability in Markov processes
Authors:
Hernan Larralde,
Francois Leyvraz,
David P. Sanders
Abstract:
We present a formalism to describe slowly decaying systems in the context of finite Markov chains obeying detailed balance. We show that phase space can be partitioned into approximately decoupled regions, in which one may introduce restricted Markov chains which are close to the original process but do not leave these regions. Within this context, we identify the conditions under which the deca…
▽ More
We present a formalism to describe slowly decaying systems in the context of finite Markov chains obeying detailed balance. We show that phase space can be partitioned into approximately decoupled regions, in which one may introduce restricted Markov chains which are close to the original process but do not leave these regions. Within this context, we identify the conditions under which the decaying system can be considered to be in a metastable state. Furthermore, we show that such metastable states can be described in thermodynamic terms and define their free energy. This is accomplished showing that the probability distribution describing the metastable state is indeed proportional to the equilibrium distribution, as is commonly assumed. We test the formalism numerically in the case of the two-dimensional kinetic Ising model, using the Wang--Landau algorithm to show this proportionality explicitly, and confirm that the proportionality constant is as derived in the theory. Finally, we extend the formalism to situations in which a system can have several metastable states.
△ Less
Submitted 19 August, 2006;
originally announced August 2006.
-
Occurrence of normal and anomalous diffusion in polygonal billiard channels
Authors:
David P. Sanders,
Hernan Larralde
Abstract:
From extensive numerical simulations, we find that periodic polygonal billiard channels with angles which are irrational multiples of pi generically exhibit normal diffusion (linear growth of the mean squared displacement) when they have a finite horizon, i.e. when no particle can travel arbitrarily far without colliding. For the infinite horizon case we present numerical tests showing that the…
▽ More
From extensive numerical simulations, we find that periodic polygonal billiard channels with angles which are irrational multiples of pi generically exhibit normal diffusion (linear growth of the mean squared displacement) when they have a finite horizon, i.e. when no particle can travel arbitrarily far without colliding. For the infinite horizon case we present numerical tests showing that the mean squared displacement instead grows asymptotically as t log t. When the unit cell contains accessible parallel scatterers, however, we always find anomalous super-diffusion, i.e. power-law growth with an exponent larger than 1. This behavior cannot be accounted for quantitatively by a simple continuous-time random walk model. Instead, we argue that anomalous diffusion correlates with the existence of families of propagating periodic orbits. Finally we show that when a configuration with parallel scatterers is approached there is a crossover from normal to anomalous diffusion, with the diffusion coefficient exhibiting a power-law divergence.
△ Less
Submitted 12 February, 2006; v1 submitted 25 October, 2005;
originally announced October 2005.
-
A Definition of Metastability for Markov Processes with Detailed Balance
Authors:
Francois Leyvraz,
Hernan Larralde,
David P. Sanders
Abstract:
A definition of metastable states applicable to arbitrary finite state Markov processes satisfying detailed balance is discussed. In particular, we identify a crucial condition that distinguishes genuine metastable states from other types of slowly decaying modes and which leads to properties similar to those postulated in the restricted ensemble approach \cite{pen71}. The intuitive physical mea…
▽ More
A definition of metastable states applicable to arbitrary finite state Markov processes satisfying detailed balance is discussed. In particular, we identify a crucial condition that distinguishes genuine metastable states from other types of slowly decaying modes and which leads to properties similar to those postulated in the restricted ensemble approach \cite{pen71}. The intuitive physical meaning of this condition is simply that the total equilibrium probability of finding the system in the metastable state is negligible. As a concrete application of our formalism we present preliminary results on a 2D kinetic Ising model.
△ Less
Submitted 28 September, 2005;
originally announced September 2005.
-
Fine structure of distributions and central limit theorem in diffusive billiards
Authors:
David P. Sanders
Abstract:
We investigate deterministic diffusion in periodic billiard models, in terms of the convergence of rescaled distributions to the limiting normal distribution required by the central limit theorem; this is stronger than the usual requirement that the mean square displacement grow asymptotically linearly in time. The main model studied is a chaotic Lorentz gas where the central limit theorem has b…
▽ More
We investigate deterministic diffusion in periodic billiard models, in terms of the convergence of rescaled distributions to the limiting normal distribution required by the central limit theorem; this is stronger than the usual requirement that the mean square displacement grow asymptotically linearly in time. The main model studied is a chaotic Lorentz gas where the central limit theorem has been rigorously proved. We study one-dimensional position and displacement densities describing the time evolution of statistical ensembles in a channel geometry, using a more refined method than histograms. We find a pronounced oscillatory fine structure, and show that this has its origin in the geometry of the billiard domain. This fine structure prevents the rescaled densities from converging pointwise to gaussian densities; however, demodulating them by the fine structure gives new densities which seem to converge uniformly. We give an analytical estimate of the rate of convergence of the original distributions to the limiting normal distribution, based on the analysis of the fine structure, which agrees well with simulation results. We show that using a Maxwellian (gaussian) distribution of velocities in place of unit speed velocities does not affect the growth of the mean square displacement, but changes the limiting shape of the distributions to a non-gaussian one. Using the same methods, we give numerical evidence that a non-chaotic polygonal channel model also obeys the central limit theorem, but with a slower convergence rate.
△ Less
Submitted 4 November, 2004;
originally announced November 2004.