Showing 1–50 of 66 results for author: Voigtmann, T

Feature Selection for High-Dimensional Neural Network Potentials with the Adaptive Group Lasso

Authors: Johannes Sandberg, Thomas Voigtmann, Emilie Devijver, Noel Jakse

Abstract: Neural network potentials are a powerful tool for atomistic simulations, allowing to accurately reproduce \textit{ab initio} potential energy surfaces with computational performance approaching classical force fields. A central component of such potentials is the transformation of atomic positions into a set of atomic features in a most efficient and informative way.In this work, a feature selecti… ▽ More Neural network potentials are a powerful tool for atomistic simulations, allowing to accurately reproduce \textit{ab initio} potential energy surfaces with computational performance approaching classical force fields. A central component of such potentials is the transformation of atomic positions into a set of atomic features in a most efficient and informative way.In this work, a feature selection method is introduced for high dimensional neural network potentials, based on the Adaptive Group Lasso (AGL) approach. It is shown that the use of an embedded method, taking into account the interplay between features and their action in the estimator, is necessary to optimize the number of features. The method's efficiency is tested on three different monoatomic systems, including Lennard-Jones as a simple test case, Aluminium as a system characterized by predominantly radial interactions, and Boron as representative of a system with strongly directional interactions. The AGL is compared with unsupervised filter methods and found to perform consistently better in reducing the number of features needed to reproduce the reference simulation data. {In particular, our results show the importance of taking into account model predictions in feature selection for interatomic potentials. △ Less

Submitted 26 December, 2023; originally announced December 2023.

Comments: 30 pages, 9 figures

Microscopic theory for nonequilibrium correlation functions in dense active fluids

Authors: Vincent E. Debets, Lila Sarfati, Thomas Voigtmann, Liesbeth M. C. Janssen

Abstract: One of the key hallmarks of dense active matter in the liquid, supercooled, and solid phases is so-called equal-time velocity correlations. Crucially, these correlations can emerge spontaneously, i.e., they require no explicit alignment interactions, and therefore represent a generic feature of dense active matter. This indicates that for a meaningful comparison or possible mapping between active… ▽ More One of the key hallmarks of dense active matter in the liquid, supercooled, and solid phases is so-called equal-time velocity correlations. Crucially, these correlations can emerge spontaneously, i.e., they require no explicit alignment interactions, and therefore represent a generic feature of dense active matter. This indicates that for a meaningful comparison or possible mapping between active and passive liquids one not only needs to understand their structural properties, but also the impact of these velocity correlations. This has already prompted several simulation and theoretical studies, though they are mostly focused on athermal systems and thus overlook the effect of translational diffusion. Here we present a fully microscopic method to calculate nonequilibrium correlations in systems of thermal active Brownian particles (ABPs). We use the integration through transients (ITT) formalism together with (active) mode-coupling theory (MCT) and analytically calculate qualitatively consistent static structure factors and active velocity correlations. We complement our theoretical results with simulations of both thermal and athermal ABPs which exemplify the disruptive role that thermal noise has on velocity correlations. △ Less

Submitted 23 October, 2023; originally announced October 2023.

Residual stresses couple microscopic and macroscopic scales

Authors: Sebastian Steinhäuser, Timm Treskatis, Stefan Turek, Thomas Voigtmann

Abstract: We show how residual stresses emerge in a visco-elastic material as a signature of its past flow history, through an interplay between flow-modified microscopic relaxation and macroscopic features of the flow. Long-lasting temporal-history dependence of the microscopic dynamics and nonlinear rheology are incorporated through the mode-coupling theory of the glass transition (MCT). The theory's inte… ▽ More We show how residual stresses emerge in a visco-elastic material as a signature of its past flow history, through an interplay between flow-modified microscopic relaxation and macroscopic features of the flow. Long-lasting temporal-history dependence of the microscopic dynamics and nonlinear rheology are incorporated through the mode-coupling theory of the glass transition (MCT). The theory's integral constitutive equation (ICE) is coupled to continuum mechanics in a finite-element method (FEM) scheme that tracks the flow history through the Finger tensor. The method is suitable for a calculation of residual stresses from a "first-principles" starting point following well-understood approximations. As an example, we calculate within a schematic version of MCT the stress-induced optical birefringence pattern of an amorphous solid cast into the shape of a slab with a cylindrical obstacle and demonstrate how FEM-MCT can predict the dependence of material properties on the material's processing history. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 5 pages, 3 figures

Spreading of droplets under various gravitational accelerations

Authors: Olfa D'Angelo, Felix Kuthe, Kasper van Nieuwland, Clint Ederveen Janssen, Thomas Voigtmann, Maziyar Jalaal

Abstract: We describe a setup to perform systematic studies on the spreading of droplets of complex fluids under microgravity conditions. Tweaking the gravitational acceleration under which droplets are deposited provides access to different regimes of the spreading dynamics, quantified through the Bond number. In particular, microgravity allows to form large droplets while remaining in the regime where sur… ▽ More We describe a setup to perform systematic studies on the spreading of droplets of complex fluids under microgravity conditions. Tweaking the gravitational acceleration under which droplets are deposited provides access to different regimes of the spreading dynamics, quantified through the Bond number. In particular, microgravity allows to form large droplets while remaining in the regime where surface tension effects and internal driving stresses are predominant over hydrostatic forces. The VIP-DROP2 experimental module provides a versatile platform to study a wide range of complex fluids through the deposition of axisymmetric droplets. The module offers the possibility to deposit droplets on a precursor layer, which can be composed of the same or of a different fluid. Besides, it allows to deposit four droplets simultaneously, while conducting shadowgraphy on all of them, and observing either the flow field (through particle image velocimetry), or the stress distribution inside the droplet in the case of stress birefringent fluids. Developed for a drop tower catapult system, it is designed to withstand a vertical acceleration of up to 30 times Earth's gravitational acceleration in the downwards direction, and can operate remotely, under microgravity conditions. We provide a detailed description of the module, and exemplary data analysis for droplets spreading on-ground and in microgravity. △ Less

Submitted 27 October, 2022; v1 submitted 15 April, 2022; originally announced April 2022.

Comments: Review of Scientific Instruments (2022)

From sub-aging to hyper-aging in structural glasses

Authors: Luis F. Elizondo-Aguilera, Tommaso Rizzo, Thomas Voigtmann

Abstract: We demonstrate non-equilibrium scaling laws for the aging dynamics in glass formers that emerge from combining a recent application of Onsager's theory of irreversible processes with the equilibrium scaling laws of glassy dynamics. Different scaling regimes are predicted for the evolution of the system's structural relaxation time $τ$ with age (waiting time $t_w$), depending on the depth of the qu… ▽ More We demonstrate non-equilibrium scaling laws for the aging dynamics in glass formers that emerge from combining a recent application of Onsager's theory of irreversible processes with the equilibrium scaling laws of glassy dynamics. Different scaling regimes are predicted for the evolution of the system's structural relaxation time $τ$ with age (waiting time $t_w$), depending on the depth of the quench from the liquid into the glass: \emph{simple aging} ($τ\sim t_w$) applies for quenches close to the critical point of mode-coupling theory (MCT) and implies \emph{sub-aging} ($τ\approx\ t_w^δ$ with $δ<1$) as a broad cross-over for quenches to nearly-arrested equilibrium states; \emph{hyper-aging} (or \emph{super-aging}, $τ\sim t_w^{δ'}$ with $δ'>1$) emerges for quenches deep into the glass. The latter is cut off by non-mean-field fluctuations that we account for within a recent extension of MCT, the stochastic $β$-relaxation theory (SBR). We exemplify the scaling laws by a schematic model that allows to quantitatively fit recent simulation results for density-quenched hard-sphere-like particles. △ Less

Submitted 27 February, 2022; originally announced February 2022.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 129, 238003 (2022)

Interatomic machine learning potentials for aluminium: application to solidification phenomena

Authors: Noel Jakse, Johannes Sandberg, Leon F. Granz, Anthony Saliou, Philippe Jarry, Emilie Devijver, Thomas Voigtmann, Jürgen Horbach, Andreas Meyer

Abstract: In studying solidification process by simulations on the atomic scale, the modeling of crystal nucleation or amorphisation requires the construction of interatomic interactions that are able to reproduce the properties of both the solid and the liquid states. Taking into account rare nucleation events or structural relaxation under deep undercooling conditions requires much larger length scales an… ▽ More In studying solidification process by simulations on the atomic scale, the modeling of crystal nucleation or amorphisation requires the construction of interatomic interactions that are able to reproduce the properties of both the solid and the liquid states. Taking into account rare nucleation events or structural relaxation under deep undercooling conditions requires much larger length scales and longer time scales than those achievable by \textit{ab initio} molecular dynamics (AIMD). This problem is addressed by means of classical MD simulations using a well established high dimensional neural network potential trained on a relevant set of configurations generated by AIMD. Our dataset contains various crystalline structures and liquid states at different pressures, including their time fluctuations in a wide range of temperatures considering only their energy labels. Applied to elemental aluminium, the resulting potential is shown to be efficient to reproduce the basic structural, dynamics and thermodynamic quantities in the liquid and undercooled states without the need to include neither explicitly the forces nor all kind of configurations in the training procedure. The early stage of crystallization is further investigated on a much larger scale with one million atoms, allowing us to unravel features of the homogeneous nucleation mechanisms in the fcc phase at ambient pressure as well as in the bcc phase at high pressure with unprecedented accuracy close to the \textit{ab initio} one. In both case, a single step nucleation process is observed. △ Less

Submitted 5 August, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: 40 pages and 10 figures + supplement file

Non-monotonic dynamic correlations beneath the surface of glass-forming liquids

Authors: Hailong Peng, Huashan Liu, Thomas Voigtmann

Abstract: Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate the emergence of distinct static and dynamic length scales probed near the free surface in fully equilibrated glass-forming liquid films, and their connection to the bulk properties of the system. In contrast to a monotonically growing static correlation length, the dynamic correl… ▽ More Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate the emergence of distinct static and dynamic length scales probed near the free surface in fully equilibrated glass-forming liquid films, and their connection to the bulk properties of the system. In contrast to a monotonically growing static correlation length, the dynamic correlation length that measures the extent of surface-dynamics acceleration into the bulk, displays a striking non-monotonic temperature evolution that is robust also against changes in detailed interatomic interaction. The maximum of dynamic correlations defines a cross-over temperature $T_*$ that we show to agree with a morphological change of cooperative rearrangement regions (CRR) of fast particles in the bulk liquids. The cross-over occurs at a temperature larger than the critical temperature Tc of mode-coupling theory (MCT). We link it to the point where fast-particle motion decouples from structural relaxation that can be defined rigorously within a recent extension of MCT, the stochastic $β$-relaxation theory (SBR). △ Less

Submitted 12 July, 2021; v1 submitted 2 July, 2021; originally announced July 2021.

Journal ref: Phys. Rev. Lett.129, 215501 (2022)

A gravity-independent powder-based additive manufacturing process tailored for space applications

Authors: Olfa D'Angelo, Felix Kuthe, Szu-Jia Liu, Raphael Wiedey, Joe M. Bennett, Martina Meisnar, Andrew Barnes, W. Till Kranz, Thomas Voigtmann, Andreas Meyer

Abstract: The future of space exploration missions will rely on technologies increasing their endurance and self-sufficiency, including for manufacturing objects on-demand. We propose a process for handling and additively manufacturing powders that functions independently of the gravitational environment and with no restriction on feedstock powder flowability. Based on a specific sequence of boundary loads… ▽ More The future of space exploration missions will rely on technologies increasing their endurance and self-sufficiency, including for manufacturing objects on-demand. We propose a process for handling and additively manufacturing powders that functions independently of the gravitational environment and with no restriction on feedstock powder flowability. Based on a specific sequence of boundary loads applied to the granular packing, powder is transported to the printing zone, homogenized and put under compression to increase the density of the final part. The powder deposition process is validated by simulations that show the homogeneity and density of deposition to be insensitive to gravity and cohesion forces within the DEM model. We further provide an experimental proof of concept of the process by successfully 3D printing parts on-ground and in weightlessness, on parabolic flight. Powders exhibiting high and low flowability are used as model feedstock material to demonstrate the versatility of the process, opening the way for additive manufacturing of recycled material. △ Less

Submitted 8 October, 2021; v1 submitted 19 February, 2021; originally announced February 2021.

Journal ref: Additive Manufacturing Additive Manufacturing 47 (2021) 102349

Transport Coefficients in Dense Active Brownian Particle Systems: Mode-Coupling Theory and Simulation Results

Authors: Julian Reichert, Leon Granz, Thomas Voigtmann

Abstract: We discuss recent advances in developing a mode-coupling theory of the glass transition (MCT) of two-dimensional systems of active Brownian particles (ABP). We specifically discuss the case of a single ABP tracer in a glass-forming passive host suspension; a case that has recently been studied in experiments on colloidal Janus particles. We employ event-driven Brownian dynamics (ED-BD) computer si… ▽ More We discuss recent advances in developing a mode-coupling theory of the glass transition (MCT) of two-dimensional systems of active Brownian particles (ABP). We specifically discuss the case of a single ABP tracer in a glass-forming passive host suspension; a case that has recently been studied in experiments on colloidal Janus particles. We employ event-driven Brownian dynamics (ED-BD) computer simulations to test the ABP-MCT, and find good agreement between the two for the MSD. The ED-BD simulation results also compare well to experimental data, although a peculiar non-monotonic mapping of self-propulsion velocities is required. The ABP-MCT predicts a specific self-propulsion dependence of the Stokes-Einstein relation between the long-time diffusion coefficient and the host-system viscosity that matches well the results from simulation. An application of ABP-MCT within the integration-through transients (ITT) framework to calculate the density-renormalized effective swim velocity of the interacting ABP agrees qualitatively with the ED-BD simulation data at densities close to the glass transition, and quantitatively for the full density range only after the mapping of packing fractions employed for the passive system. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Tracer Dynamics in Crowded Active-Particle Suspensions

Authors: Julian Reichert, Thomas Voigtmann

Abstract: We derive equations of motion for the mean-squared displacement (MSD) of an active Brownian particle (ABP) in a crowded environment modeled by a dense system of passive Brownian particles, and of a passive tracer particle in a dense active-Brownian particle system, using a projection-operator scheme. The interaction of the tracer particle with the dense host environment gives rise to strong memory… ▽ More We derive equations of motion for the mean-squared displacement (MSD) of an active Brownian particle (ABP) in a crowded environment modeled by a dense system of passive Brownian particles, and of a passive tracer particle in a dense active-Brownian particle system, using a projection-operator scheme. The interaction of the tracer particle with the dense host environment gives rise to strong memory effects. Evaluating these approximately in the framework of a recently developed mode-coupling theory for the glass transition in active Brownian particles (ABP-MCT), we discuss the various regimes of activity-induced super-diffusive motion and density-induced sub-diffusive motion. The predictions of the theory are shown to be in good agreement with results from an event-driven Brownian dynamics simulation scheme for the dynamics of two-dimensional active Brownian hard disks. △ Less

Submitted 26 July, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Mode-Coupling Theory for Tagged-Particle Motion of Active Brownian Particles

Authors: Julian Reichert, Suvendu Mandal, Thomas Voigtmann

Abstract: We derive a mode-coupling theory (MCT) to describe the dynamics of tracer particles in dense systems of active Brownian particles (ABPs) in two spatial dimensions. The ABP undergo translational and rotational Brownian dynamics, and are equipped with a fixed self-propulsion speed along their orientational vector that describes their active motility. The resulting equations of motion for the tagged-… ▽ More We derive a mode-coupling theory (MCT) to describe the dynamics of tracer particles in dense systems of active Brownian particles (ABPs) in two spatial dimensions. The ABP undergo translational and rotational Brownian dynamics, and are equipped with a fixed self-propulsion speed along their orientational vector that describes their active motility. The resulting equations of motion for the tagged-particle density correlation functions describe the various cases of tracer dynamics close to the glass transition: that of a passive colloidal particle in a suspension of ABP, that of a single active particle in a glass-forming passive host suspensions, and that of active tracers in a bath of active particles. Numerical results are presented for these cases assuming hard-sphere interactions among the particles. The qualitative and quantitative accuracy of the theory is tested against event-driven Brownian dynamics (ED-BD) simulations of active and passive hard disks. Simulation and theory are found in quantitative agreement, provided one adjusts the overall density (as known from the passive description of glassy dynamics), and allows for a rescaling of self-propulsion velocities in the active host system. These adjustments account for the fact that ABP-MCT generally overestimates the tendency for kinetic arrest. We also confirm in the simulations a peculiar feature of the transient and stationary dynamical density correlation functions regarding their lack of symmetry under time reversal, demonstrating the non-equilibrium nature of the system and how it manifests itself in the theory. △ Less

Submitted 27 July, 2021; v1 submitted 26 October, 2020; originally announced October 2020.

Scaling equations for mode-coupling theories with multiple decay channels

Authors: Gerhard Jung, Thomas Voigtmann, Thomas Franosch

Abstract: Multiple relaxation channels often arise in the dynamics of liquids where the momentum current associated to the particle-conservation law splits into distinct contributions. Examples are strongly confined liquids for which the currents in lateral and longitudinal direction to the walls are very different, or fluids of nonspherical particles with distinct relaxation patterns for translational and… ▽ More Multiple relaxation channels often arise in the dynamics of liquids where the momentum current associated to the particle-conservation law splits into distinct contributions. Examples are strongly confined liquids for which the currents in lateral and longitudinal direction to the walls are very different, or fluids of nonspherical particles with distinct relaxation patterns for translational and rotational degrees of freedom. Here, we perform an asymptotic analysis of the slow structural relaxation close to kinetic arrest as described by mode-coupling theory (MCT) with several relaxation channels. Compared to standard MCT, the presence of multiple relaxation channels significantly changes the structure of the underlying equations of motion and leads to additional, non-trivial terms in the asymptotic solution. We show that the solution can be rescaled, and thus prove that the well-known $ β$-scaling equation of MCT remains valid even in the presence of multiple relaxation channels. The asymptotic treatment is validated using a novel schematic model. We demonstrate that the numerical solution of this schematic model can indeed be described by the derived asymptotic scaling laws close to kinetic arrest. Additionally, clear traces of the existence of two distinct decay channels are found in the low-frequency susceptibility spectrum, suggesting that clear footprints of the additional relaxation channels can in principle be detected in simulations or experiments of confined or molecular liquids. △ Less

Submitted 27 May, 2020; originally announced May 2020.

Comments: Accepted by J. Stat. Mech.: Theory and Experiments, https://iopscience.iop.org/journal/1742-5468

Multiple character of non-monotonic size-dependence for relaxation dynamics in polymer-particle and binary mixtures

Authors: Elias M. Zirdehi, Thomas Voigtmann, Fathollah Varnik

Abstract: Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules [The Journal of Chemical Physics 150 (2019) 024903]. In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in… ▽ More Adding plasticizers is a well-known procedure to reduce the glass transition temperature in polymers. It has been recently shown that this effect shows a non-monotonic dependence on the size of additive molecules [The Journal of Chemical Physics 150 (2019) 024903]. In this work, we demonstrate that, as the size of the additive molecules is changed at fixed concentration, multiple extrema emerge in the dependence of the system's relaxation time on the size ratio. The effect occurs on all relevant length scales including single monomer dynamics, decay of Rouse modes and relaxation of the chain's end-to-end vector. A qualitatively similar trend is found within mode-coupling theoretical results for a binary hard-sphere (HS) mixture. An interpretation of the effect in terms of local packing efficiency and coupling between the dynamics of minority and majority species is provided. △ Less

Submitted 10 December, 2019; originally announced December 2019.

Comments: 16 pages, 7 figures

Arrested dynamics of the dipolar hard-sphere model

Authors: Luis Fernando Elizondo-Aguilera, Ernesto Carlos Cortés-Morales, Pablo F. Zubieta Rico, Magdaleno Medina-Noyola, Ramón Castañeda-Priego, Thomas Voigtmann, Gabriel Pérez-Ángel

Abstract: We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (softcore) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down… ▽ More We report the combined results of molecular dynamics simulations and theoretical calculations concerning various dynamical arrest transitions in a model system representing a dipolar fluid, namely, N (softcore) rigid spheres interacting through a truncated dipole-dipole potential. By exploring different regimes of concentration and temperature, we find three distinct scenarios for the slowing down of the dynamics of the translational and orientational degrees of freedom: At low ({$η$} = 0.2) and intermediate ($η$ = 0.4) volume fractions, both dynamics are strongly coupled and become simultaneously arrested upon cooling. At high concentrations ({$η$} $\lt$ 0.6), the translational dynamics shows the features of an ordinary glass transition, either by compressing or cooling down the system, but with the orientations remaining ergodic, thus indicating the existence of partially arrested states. In this density regime, but at lower temperatures, the relaxation of the orientational dynamics also freezes. The physical scenario provided by the simulations is discussed and compared against results obtained with the self-consistent generalized Langevin equation theory, and both provide a consistent description of the dynamical arrest transitions in the system. Our results are summarized in an arrested states diagram which qualitatively organizes the simulation data and provides a generic picture of the glass transitions of a dipolar fluid. △ Less

Submitted 14 March, 2019; originally announced March 2019.

Solvable Models of Supercooled Liquids in Three Dimensions

Authors: Tommaso Rizzo, Thomas Voigtmann

Abstract: We introduce a supercooled liquid model and obtain parameter-free quantitative predictions that are in excellent agreement with numerical simulations, notably in the hard low-temperature region characterized by strong deviations from Mode-Coupling-Theory behavior. The model is the Fredrickson-Andersen Kinetically-Constrained-Model on the three-dimensional $M$-layer lattice. The agreement has impli… ▽ More We introduce a supercooled liquid model and obtain parameter-free quantitative predictions that are in excellent agreement with numerical simulations, notably in the hard low-temperature region characterized by strong deviations from Mode-Coupling-Theory behavior. The model is the Fredrickson-Andersen Kinetically-Constrained-Model on the three-dimensional $M$-layer lattice. The agreement has implications beyond the specific model considered because the theory is potentially valid for many more systems, including realistic models and actual supercooled liquids. △ Less

Submitted 4 May, 2020; v1 submitted 5 March, 2019; originally announced March 2019.

Comments: Letter with supplemental material text attached. Added discussion of finite M and realistic models in the renormalization group framework. To appear on Phys. Rev. Lett

Journal ref: Phys. Rev. Lett. 124, 195501 (2020)

On the dynamics of reaction coordinates in classical, time-dependent, many-body processes

Authors: Hugues Meyer, Thomas Voigtmann, Tanja Schilling

Abstract: Complex microscopic many-body processes are often interpreted in terms of so-called `reaction coordinates', i.e. in terms of the evolution of a small set of coarse-grained observables. A rigorous method to produce the equation of motion of such observables is to use projection operator techniques, which split the dynamics of the observables into a main contribution and a marginal one. The basis of… ▽ More Complex microscopic many-body processes are often interpreted in terms of so-called `reaction coordinates', i.e. in terms of the evolution of a small set of coarse-grained observables. A rigorous method to produce the equation of motion of such observables is to use projection operator techniques, which split the dynamics of the observables into a main contribution and a marginal one. The basis of any derivation in this framework is the classical (or quantum) Heisenberg equation for an observable. If the Hamiltonian of the underlying microscopic dynamics and the observable under study do not explicitly depend on time, this equation is obtained by a straight-forward derivation. However, the problem is more complicated if one considers Hamiltonians which depend on time explicitly as e.g. in systems under external driving, or if the observable of interest has an explicit dependence on time. We use an analogy to fluid dynamics to derive the classical Heisenberg picture and then apply a projection operator formalism to derive the non-stationary generalized Langevin equation for a coarse-grained variable. We show, in particular, that the results presented for time-independent Hamiltonians and observables in J. Chem. Phys. 147, 214110 (2017) can be generalized to the time-dependent case. △ Less

Submitted 15 March, 2019; v1 submitted 28 January, 2019; originally announced January 2019.

Glassy dynamics in asymmetric binary mixtures of hard-spheres

Authors: Edilio Lázaro-Lázaro, Jorge Adrián Perera-Burgos, Patrick Laermann, Tatjana Sentjabrskaja, Gabriel Pérez-Ángel, Marco Laurati, Stefan U. Egelhaaf, Magdaleno Medina-Noyola, Thomas Voigtmann, Ramón Castañeda-Priego, Luis Fernando Elizondo-Aguilera

Abstract: The binary hard-sphere mixture is one of the simplest representations of a many-body system with competing time and length scales. This model is relevant to fundamentally understand both the structural and dynamical properties of materials, such as metallic melts, colloids, polymers and bio-based composites. It also allows us to study how different scales influence the physical behavior of a multi… ▽ More The binary hard-sphere mixture is one of the simplest representations of a many-body system with competing time and length scales. This model is relevant to fundamentally understand both the structural and dynamical properties of materials, such as metallic melts, colloids, polymers and bio-based composites. It also allows us to study how different scales influence the physical behavior of a multicomponent glass-forming liquid; a question that still awaits a unified description. In this contribution, we report on distinct dynamical arrest transitions in highly asymmetric binary colloidal mixtures, namely, a single glass of big particles, in which the small species remains ergodic, and a double glass with the simultaneous arrest of both components. When the mixture approaches any glass transition, the relaxation of the collective dynamics of both species becomes coupled. In the single glass domain, spatial modulations occur due to the structure of the large spheres, a feature not observed in the two-glass domain. The relaxation of the \emph{self} dynamics of small and large particles, in contrast, become decoupled at the boundaries of both transitions; the large species always displays dynamical arrest, whereas the small ones appear arrested only in the double glass. Thus, in order to obtain a complete picture of the distinct glassy states, one needs to take into account the dynamics of both species. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Journal ref: Phys. Rev. E 99, 042603 (2019)

Static Structure of Active Brownian Hard Disks

Authors: Nadine de Macedo Biniossek, Hartmut Löwen, Thomas Voigtmann, Frank Smallenburg

Abstract: We explore the changes in static structure of a two-dimensional system of active Brownian particles (ABP) with hard-disk interactions, using event-driven Brownian dynamics simulations. In particular, the effect of the self-propulsion velocity and the rotational diffusivity on the orientationally-averaged fluid structure factor is discussed. Typically activity increases structural ordering and gene… ▽ More We explore the changes in static structure of a two-dimensional system of active Brownian particles (ABP) with hard-disk interactions, using event-driven Brownian dynamics simulations. In particular, the effect of the self-propulsion velocity and the rotational diffusivity on the orientationally-averaged fluid structure factor is discussed. Typically activity increases structural ordering and generates a structure factor peak at zero wave vector which is a precursor of motility-induced phase separation. Our results provide reference data to test future statistical theories for the fluid structure of active Brownian systems. △ Less

Submitted 24 December, 2017; originally announced December 2017.

Comments: Published in Journal of Physics: Condensed Matter. This is the pre-review version

How Glassy Relaxation Slows Down by Increasing Mobility

Authors: Suvendu Mandal, Thomas Franosch, Thomas Voigtmann

Abstract: We investigate how structural relaxation in mixtures with strong dynamical asymmetry is affected by the microscopic dynamics. Brownian and Newtonian dynamics simulations of dense mixtures of fast and slow hard spheres reveal a striking trend reversal. Below a critical density, increasing the mobility of the fast particles fluidizes the system, yet, above that critical density, the same increase in… ▽ More We investigate how structural relaxation in mixtures with strong dynamical asymmetry is affected by the microscopic dynamics. Brownian and Newtonian dynamics simulations of dense mixtures of fast and slow hard spheres reveal a striking trend reversal. Below a critical density, increasing the mobility of the fast particles fluidizes the system, yet, above that critical density, the same increase in mobility strongly hinders the relaxation of the slow particles. The critical density itself does not depend on the dynamical asymmetry and can be identified with the glass-transition density of the mode-coupling theory. The asymptotic dynamics close to the critical density is universal, but strong pre-asymptotic effects prevail in mixtures with additional size asymmetry. This observation reconciles earlier findings of a strong dependence on kinetic parameters of glassy dynamics in colloid--polymer mixtures with the paradigm that the glass transition is determined by the properties of configuration space alone. △ Less

Submitted 28 September, 2017; originally announced September 2017.

Mode-Coupling Theory for Active Brownian Particles

Authors: Alexander Liluashvili, Jonathan Onody, Thomas Voigtmann

Abstract: We present a mode-coupling theory (MCT) for the high-density dynamics of two-dimensional spherical active Brownian particles (ABP). The theory is based on the integration-through-transients (ITT) formalism and hence provides a starting point for the calculation of non-equilibrium averages in active-Brownian particle systems. The ABP are characterized by a self-propulsion velocity $v_0$, and by the… ▽ More We present a mode-coupling theory (MCT) for the high-density dynamics of two-dimensional spherical active Brownian particles (ABP). The theory is based on the integration-through-transients (ITT) formalism and hence provides a starting point for the calculation of non-equilibrium averages in active-Brownian particle systems. The ABP are characterized by a self-propulsion velocity $v_0$, and by their translational and rotational diffusion coefficients, $D_t$ and $D_r$. The theory treats both the translational and the orientational degrees of freedom of ABP explicitly. This allows to study the effect of self-propulsion of both weak and strong persistence of the swimming direction, also at high densities where the persistence length $\ell_p=v_0/D_r$ is large compared to the typical interaction length scale. While the low-density dynamics of ABP is characterized by a single Péclet number, $Pe=v_0^2/D_rD_t$, close to the glass transition the dynamics is found to depend on $Pe$ and $\ell_p$ separately. At fixed density, increasing the self-propulsion velocity causes structural relaxatino to speed up, while decreasing the persistence length slows down the relaxation. The theory predicts a non-trivial idealized-glass-transition diagram in the three-dimensional parameter space of density, self-propulsion velocity and rotational diffusivity. The active-MCT glass is a nonergodic state where correlations of initial density fluctuations never fully decay, but also an infinite memory of initial orientational fluctuations is retained in the positions. △ Less

Submitted 23 July, 2017; originally announced July 2017.

Journal ref: Phys. Rev. E 96, 062608 (2017)

Crystal Growth in Fluid Flow: Nonlinear Response Effects

Authors: H. L. Peng, D. M. Herlach, Th. Voigtmann

Abstract: We investigate crystal-growth kinetics in the presence of strong shear flow in the liquid, using molecular-dynamics simulations of a binary-alloy model. Close to the equilibrium melting point, shear flow always suppresses the growth of the crystal-liquid interface. For lower temperatures, we find that the growth velocity of the crystal depends non-monotonically on the shear rate. Slow enough flow… ▽ More We investigate crystal-growth kinetics in the presence of strong shear flow in the liquid, using molecular-dynamics simulations of a binary-alloy model. Close to the equilibrium melting point, shear flow always suppresses the growth of the crystal-liquid interface. For lower temperatures, we find that the growth velocity of the crystal depends non-monotonically on the shear rate. Slow enough flow enhances the crystal growth, due to an increased particle mobility in the liquid. Stronger flow causes a growth regime that is nearly temperature-independent, in striking contrast to what one expects from the thermodynamic and equilibrium kinetic properties of the system, which both depend strongly on temperature. We rationalize these effects of flow on crystal growth as resulting from the nonlinear response of the fluid to strong shearing forces. △ Less

Submitted 18 July, 2017; originally announced July 2017.

Comments: to appear in Phys. Rev. Materials

On the "generalized Generalized Langevin Equation"

Authors: Hugues Meyer, Thomas Voigtmann, Tanja Schilling

Abstract: In molecular dynamics simulations and single molecule experiments, observables are usually measured along dynamic trajectories and then averaged over an ensemble ("bundle") of trajectories. Under stationary conditions, the time-evolution of such averages is described by the generalized Langevin equation. In contrast, if the dynamics is not stationary, it is not a priori clear which form the equati… ▽ More In molecular dynamics simulations and single molecule experiments, observables are usually measured along dynamic trajectories and then averaged over an ensemble ("bundle") of trajectories. Under stationary conditions, the time-evolution of such averages is described by the generalized Langevin equation. In contrast, if the dynamics is not stationary, it is not a priori clear which form the equation of motion for an averaged observable has. We employ the formalism of time-dependent projection operator techniques to derive the equation of motion for a non-equilibrium trajectory-averaged observable as well as for its non-stationary auto-correlation function. The equation is similar in structure to the generalized Langevin equation, but exhibits a time-dependent memory kernel as well as a fluctuating force that implicitly depends on the initial conditions of the process. We also derive a relation between this memory kernel and the autocorrelation function of the fluctuating force that has a structure similar to a fluctuation-dissipation relation. In addition, we show how the choice of the projection operator allows to relate the Taylor expansion of the memory kernel to data that is accessible in MD simulations and experiments, thus allowing to construct the equation of motion. As a numerical example, the procedure is applied to Brownian motion initialized in non-equilibrium conditions, and is shown to be consistent with direct measurements from simulations. △ Less

Submitted 2 October, 2017; v1 submitted 2 June, 2017; originally announced June 2017.

Clearing out a maze: The hungry random walker and its anomalous diffusion

Authors: Tanja Schilling, Thomas Voigtmann

Abstract: We study chemotaxis in a porous medium using as a model a biased ("hungry") random walk on a percolating cluster. In close resemblance to the 1980s arcade game Pac-Man, the hungry random walker consumes food, which is initially distributed in the maze, and biases its movement towards food-filled sites. We observe that, on the percolating cluster, the mean-squared displacement of the pacman process… ▽ More We study chemotaxis in a porous medium using as a model a biased ("hungry") random walk on a percolating cluster. In close resemblance to the 1980s arcade game Pac-Man, the hungry random walker consumes food, which is initially distributed in the maze, and biases its movement towards food-filled sites. We observe that, on the percolating cluster, the mean-squared displacement of the pacman process shows anomalous dynamics, which follow a power law with a dynamical exponent different from both that of a self avoiding random walk as well as that of an unbiased random walk. The change in dynamics with the propensity to move towards food is well described by a dynamical exponent that depends continuously on this propensity, and results in slower differential growth when compared to the unbiased random walk. △ Less

Submitted 25 July, 2016; v1 submitted 5 July, 2016; originally announced July 2016.

Anomalous dynamics of intruders in a crowded environment of mobile obstacles

Authors: Tatjana Sentjabrskaja, Emanuela Zaccarelli, Cristiano De Michele, Francesco Sciortino, Piero Tartaglia, Thomas Voigtmann, Stefan U. Egelhaaf, Marco Laurati

Abstract: Many natural and industrial processes rely on constrained transport, such as proteins moving through cells, particles confined in nanocomposite materials or gels, individuals in highly dense collec- tives and vehicular traffic conditions. These are examples of motion through crowded environments, in which the host matrix may retain some glass-like dynamics. Here we investigate constrained transpor… ▽ More Many natural and industrial processes rely on constrained transport, such as proteins moving through cells, particles confined in nanocomposite materials or gels, individuals in highly dense collec- tives and vehicular traffic conditions. These are examples of motion through crowded environments, in which the host matrix may retain some glass-like dynamics. Here we investigate constrained transport in a colloidal model system, in which dilute small spheres move in a slowly rearranging, glassy matrix of large spheres. Using confocal differential dynamic microscopy and simulations, we discover a critical size asymmetry at which anomalous collective transport of the small particles appears, manifested as a logarithmic decay of the density autocorrelation functions. We demonstrate that the matrix mobility is central for the observed anomalous behaviour. These results, crucially depending on size-induced dynamic asymmetry, are of relevance for a wide range of phenomena ranging from glassy systems to cell biology. △ Less

Submitted 8 April, 2016; originally announced April 2016.

Comments: 13 pages, 7 figures

Journal ref: Nature Communications 7, 11133 (2016)

On The Nature of the Glass Crossover

Authors: Tommaso Rizzo, Thomas Voigtmann

Abstract: Stochastic Beta Relaxation (SBR) is a model for the dynamics of glass- forming liquids close to the glass transition singularity of the idealized mode- coupling theory (MCT) that has been derived from generic MCT-like theories by applying dynamical field-theory techniques. SBR displays a rich phenomenology common to most super-cooled liquids. In its simplest version it naturally explains two promi… ▽ More Stochastic Beta Relaxation (SBR) is a model for the dynamics of glass- forming liquids close to the glass transition singularity of the idealized mode- coupling theory (MCT) that has been derived from generic MCT-like theories by applying dynamical field-theory techniques. SBR displays a rich phenomenology common to most super-cooled liquids. In its simplest version it naturally explains two prominent features of the dynamical crossover: the change from a power-law to exponential increase in the structural relaxation time and the violation of the Stokes-Einstein relation between diffusion and viscosity. The solution of the model in three dimensions unveils a qualitative change at the crossover in the structure of dynamical fluctuations from a regime characterized by power-law increases of their amplitude and size to a regime dominated by strong Dynamical Heterogeneities: rare regions where dynamics is relatively much faster than in the rest of the system. While the relaxation time changes by orders of magnitude, the size of these regions does not change significantly and actually decreases below the crossover temperature. SBR cannot sustain too large fluctuations and could fail well below the crossover temperature. There it could be replaced by non-conventional activated dynamics characterized by elementary events with intrinsic time and length scales of an unusual large (but not necessarily increasing) size (mesoscopic vs. microscopic). △ Less

Submitted 23 April, 2015; originally announced April 2015.

Lattice Boltzmann simulations of a viscoelastic shear-thinning fluid

Authors: S. Papenkort, Th. Voigtmann

Abstract: We present a hybrid lattice Boltzmann algorithm for the simulation of flow glass-forming fluids, characterized by slow structural relaxation, at the level of the Navier-Stokes equation. The fluid is described in terms of a nonlinear integral constitutive equation, relating the stress tensor locally to the history of flow. As an application, we present results for an integral nonlinear Maxwell mode… ▽ More We present a hybrid lattice Boltzmann algorithm for the simulation of flow glass-forming fluids, characterized by slow structural relaxation, at the level of the Navier-Stokes equation. The fluid is described in terms of a nonlinear integral constitutive equation, relating the stress tensor locally to the history of flow. As an application, we present results for an integral nonlinear Maxwell model that combines the effects of (linear) viscoelasticity and (nonlinear) shear thinning. We discuss the transient dynamics of velocities, shear stresses, and normal stress differences in planar pressure-driven channel flow, after switching on (startup) and off (cessation) of the driving pressure. This transient dynamics depends nontrivially on the channel width due to an interplay between hydrodynamic momentum diffusion and slow structural relaxation. △ Less

Submitted 23 April, 2015; originally announced April 2015.

Dissipation in a Crystallization Process

Authors: Sven Dorosz, Thomas Voigtmann, Tanja Schilling

Abstract: We discuss the crystallization process from the supersaturated melt in terms of its non-equilibrium properties. In particular, we quantify the amount of heat that is produced irreversibly when a suspension of hard spheres crystallizes. This amount of heat can be interpreted as arising from the resistance of the system against undergoing phase transition. We identify an intrinsic compression rate t… ▽ More We discuss the crystallization process from the supersaturated melt in terms of its non-equilibrium properties. In particular, we quantify the amount of heat that is produced irreversibly when a suspension of hard spheres crystallizes. This amount of heat can be interpreted as arising from the resistance of the system against undergoing phase transition. We identify an intrinsic compression rate that separates a quasi-static regime from a regime of rapid crystallization. In the former the disspated heat grows linearly in the compression rate. In the latter the system crystallizes more easily, because new relaxation channels are opened, at the cost of forming a higher fraction of non-equilibrium crystal structures. In analogy to a shear-thinning fluid, the system shows a decreased resistance when it is driven rapidly. △ Less

Submitted 16 April, 2015; originally announced April 2015.

Diffusion and Interdiffusion in Binary Metallic Melts

Authors: P. Kuhn, J. Horbach, F. Kargl, A. Meyer, Th. Voigtmann

Abstract: We discuss the dependence of self- and interdiffusion coefficients on temperature and composition for two prototypical binary metallic melts, Al-Ni and Zr-Ni, in molecular-dynamics (MD) computer simulations and the mode-coupling theory of the glass transition (MCT). Dynamical processes that are mainly entropic in origin slow down mass transport (as expressed through self diffusion) in the mixture… ▽ More We discuss the dependence of self- and interdiffusion coefficients on temperature and composition for two prototypical binary metallic melts, Al-Ni and Zr-Ni, in molecular-dynamics (MD) computer simulations and the mode-coupling theory of the glass transition (MCT). Dynamical processes that are mainly entropic in origin slow down mass transport (as expressed through self diffusion) in the mixture as compared to the ideal-mixing contribution. Interdiffusion of chemical species is a competition of slow kinetic modes with a strong thermodynamic driving force that is caused by non-entropic interactions. The combination of both dynamic and thermodynamic effects causes qualitative differences in the concentration dependence of self-diffusion and interdiffusion coefficients. At high temperatures, the thermodynamic enhancement of interdiffusion prevails, while at low temperatures, kinetic effects dominate the concentration dependence, rationalized within MCT as the approach to its ideal-glass transition temperature $T_c$. The Darken equation relating self- and interdiffusion qualitatively reproduces the concentration-dependence in both Zr-Ni and Al-Ni, but quantitatively, the kinetic contributions to interdiffusion can be slower than the lower bound suggested by the Darken equation. As temperature is decreased, the agreement with Darken's equation improves, due to a strong coupling of all kinetic modes that is a generic feature predicted by MCT. △ Less

Submitted 9 August, 2014; originally announced August 2014.

Comments: 16 pages, 12 figures

Journal ref: Physical Review B 90, 024309 (2014)

Qualitative Features at the Glass Crossover

Authors: Tommaso Rizzo, Thomas Voigtmann

Abstract: We discuss some generic features of the dynamics of glass-forming liquids close to the glass transition singularity of the idealized mode-coupling theory (MCT). The analysis is based on a recent model by one of the authors for the intermediate-time dynamics ($β$ relaxation), derived by applying dynamical field-theory techniques to the idealized MCT. Combined with the assumption of time-temperature… ▽ More We discuss some generic features of the dynamics of glass-forming liquids close to the glass transition singularity of the idealized mode-coupling theory (MCT). The analysis is based on a recent model by one of the authors for the intermediate-time dynamics ($β$ relaxation), derived by applying dynamical field-theory techniques to the idealized MCT. Combined with the assumption of time-temperature superposition for the slow structural ($α$) relaxation, the model naturally explains three prominent features of the dynamical crossover: the change from a power-law to exponential increase in the structural relaxation time, the replacement of the Stokes-Einstein relation between diffusion and viscosity by a fractional law, and two distinct growth regimes of the thermal susceptibility that has been associated to dynamical heterogeneities. △ Less

Submitted 19 June, 2014; v1 submitted 11 March, 2014; originally announced March 2014.

Journal ref: EPL 111 (2015) 56008

Nonlinear rheology of dense colloidal systems with short-ranged attraction: A mode-coupling theory analysis

Authors: Madhu Priya, Thomas Voigtmann

Abstract: The nonlinear rheology of glass-forming colloidal suspensions with short-ranged attractions is discussed within the integration-through transients framework combined with the mode-coupling theory of the glass transition (ITT-MCT). Calculations are based on the square-well system (SWS), as a model for colloid-polymer mixtures. The high-density regime featuring reentrant melting of the glass upon in… ▽ More The nonlinear rheology of glass-forming colloidal suspensions with short-ranged attractions is discussed within the integration-through transients framework combined with the mode-coupling theory of the glass transition (ITT-MCT). Calculations are based on the square-well system (SWS), as a model for colloid-polymer mixtures. The high-density regime featuring reentrant melting of the glass upon increasing the attraction strength, and the crossover from repulsive glasses formed at weak attraction to attractive glasses formed at strong attraction, are discussed. Flow curves are found in qualitative agreement with experimental data, featuring a strong increase in the yield stress, and, for suitable interaction parameters, the crossover between two yield stresses. The yield strain, defined as the position of the stress overshoot under startup flow, is found to be proportional to the attraction range for strong attraction. At weak and intermediate attraction strength, the combined effects of hard-core caging and attraction-driven bonding result in a richer dependence on the parameters. The first normal-stress difference exhibits a weaker dependence on short-ranged attractions as the shear stress, since the latter is more sensitive the short-wavelength features of the static structure. △ Less

Submitted 29 January, 2014; originally announced January 2014.

Comments: 14 pages, 12 figures

Channel Flow of a Tensorial Shear-Thinning Maxwell Model: Lattice Boltzmann Simulations

Authors: Simon Papenkort, Thomas Voigtmann

Abstract: We introduce a nonlinear generalized tensorial Maxwell-type constitutive equation to describe shear-thinning glass-forming fluids, motivated by a recent microscopic approach to the nonlinear rheology of colloidal suspensions. The model captures a nonvanishing dynamical yield stress at the glass transition and incorporates normal-stress differences. A modified lattice-Boltzmann (LB) simulation sche… ▽ More We introduce a nonlinear generalized tensorial Maxwell-type constitutive equation to describe shear-thinning glass-forming fluids, motivated by a recent microscopic approach to the nonlinear rheology of colloidal suspensions. The model captures a nonvanishing dynamical yield stress at the glass transition and incorporates normal-stress differences. A modified lattice-Boltzmann (LB) simulation scheme is presented that includes non-Newtonian contributions to the stress tensor and deals with flow-induced pressure differences. We test this scheme in pressure-driven 2D Poiseuille flow of the nonlinear generalized Maxwell fluid. In the steady state, comparison with an analytical solution shows good agreement. The transient dynamics after startup and cessation of the pressure gradient are studied; the simulation reproduces a finite stopping time for the cessation flow of the yield-stress fluid in agreement with previous analytical estimates. △ Less

Submitted 28 October, 2013; originally announced October 2013.

Residual Stresses in Glasses

Authors: M. Ballauff, J. M. Brader, S. U. Egelhaaf, M. Fuchs, J. Horbach, N. Koumakis, M. Krüger, M. Laurati, K. J. Mutch, G. Petekidis, M. Siebenbürger, Th. Voigtmann, J. Zausch

Abstract: The history dependence of the glasses formed from flow-melted steady states by a sudden cessation of the shear rate $\dotγ$ is studied in colloidal suspensions, by molecular dynamics simulations, and mode-coupling theory. In an ideal glass, stresses relax only partially, leaving behind a finite persistent residual stress. For intermediate times, relaxation curves scale as a function of $\dotγt$, e… ▽ More The history dependence of the glasses formed from flow-melted steady states by a sudden cessation of the shear rate $\dotγ$ is studied in colloidal suspensions, by molecular dynamics simulations, and mode-coupling theory. In an ideal glass, stresses relax only partially, leaving behind a finite persistent residual stress. For intermediate times, relaxation curves scale as a function of $\dotγt$, even though no flow is present. The macroscopic stress evolution is connected to a length scale of residual liquefaction displayed by microscopic mean-squared displacements. The theory describes this history dependence of glasses sharing the same thermodynamic state variables, but differing static properties. △ Less

Submitted 15 February, 2013; originally announced February 2013.

Comments: submitted to Physical Review

Stress-strain relations in bulk metallic glasses and colloidal dispersions

Authors: C. P. Amann, M. Ballauff, S. U. Egelhaaf, S. Fritschi, M. Fuchs, M. Krüger, M. Laurati, K. J. Mutch, K. Samwer, M. Siebenbürger, Th. Voigtmann, F. Weysser

Abstract: A comparison is made between the nonlinear rheological response of bulk metallic glass formers and of colloidal dispersions. Stress-strain curves measured after switch-on of constant deformation rates are analyzed quantitatively using a schematic model of mode coupling theory generalized to homogeneous and incompressible flows. A mapping between metallic and dispersion rheology is possible when st… ▽ More A comparison is made between the nonlinear rheological response of bulk metallic glass formers and of colloidal dispersions. Stress-strain curves measured after switch-on of constant deformation rates are analyzed quantitatively using a schematic model of mode coupling theory generalized to homogeneous and incompressible flows. A mapping between metallic and dispersion rheology is possible when stresses are rescaled by an entropic scale, accumulated strains by geometrical factors, and rates by the intrinsic relaxation time. Exploiting this similarity and the possibility to directly observe individual colloidal particles, we investigate the structural distortions in the colloidal system using confocal microscopy. The distortions exhibit the (from elasticity theory) expected quadrupolar but also a strong isotropic component. △ Less

Submitted 28 September, 2017; v1 submitted 8 February, 2013; originally announced February 2013.

Comments: This paper has been withdrawn because the underlying cooperation has ended and an updated theoretical description has been published elsewhere

On the Bauschinger effect in supercoooled melts under shear: results from mode coupling theory and molecular dynamics simulations

Authors: Fabian Frahsa, Amit Kumar Bhattacharjee, Jürgen Horbach, Matthias Fuchs, Thomas Voigtmann

Abstract: We study the nonlinear rheology of a glass-forming binary mixture under the reversal of shear flow using molecular dynamics simulations and a schematic model of the mode-coupling theory of the glass transition (MCT). Memory effects lead to a history-dependent response, as exemplified by the vanishing of a stress-overshoot phenomenon in the stress--strain curves of the sheared liquid, and a change… ▽ More We study the nonlinear rheology of a glass-forming binary mixture under the reversal of shear flow using molecular dynamics simulations and a schematic model of the mode-coupling theory of the glass transition (MCT). Memory effects lead to a history-dependent response, as exemplified by the vanishing of a stress-overshoot phenomenon in the stress--strain curves of the sheared liquid, and a change in the apparent elastic coefficients around states with zero stress. We investigate the various retarded contributions to the stress response at a given time schematically within MCT. The connection of this macroscopic response to single-particle motion is demonstrated using molecular-dynamics simulation. △ Less

Submitted 29 November, 2012; originally announced November 2012.

Dynamic arrest in model porous media -- intermediate scattering functions

Authors: Markus Spanner, Simon K. Schnyder, Felix Höfling, Thomas Voigtmann, Thomas Franosch

Abstract: Heterogeneous media constitute random disordered environments where transport is drastically hindered. Employing extensive molecular dynamics simulations, we investigate the spatio-temporal dynamics of tracer particles in the Lorentz model in the vicinity of the localization transition. There transport becomes anomalous and non-gaussian due to the presence of self-similar spatial structures, and d… ▽ More Heterogeneous media constitute random disordered environments where transport is drastically hindered. Employing extensive molecular dynamics simulations, we investigate the spatio-temporal dynamics of tracer particles in the Lorentz model in the vicinity of the localization transition. There transport becomes anomalous and non-gaussian due to the presence of self-similar spatial structures, and dynamic scaling behavior is anticipated. The interplay of different time and length scales is revealed by the intermediate scattering functions, which are sensitive both to the underlying spatial fractal as well as the anomalous transport. We compare our numerical results in the transition regime to a mode-coupling approach, and find that certain aspects are surprisingly well predicted. △ Less

Submitted 19 November, 2012; originally announced November 2012.

Journal ref: Soft Matter, 2013, 9, 1604-1611

Probability densities of a forced probe particle in glass: results from mode coupling theory and simulations of active microrheology

Authors: Christian J. Harrer, Antonio M. Puertas, Thomas Voigtmann, Matthias Fuchs

Abstract: We investigate the displacements of a probe particle inside a glass, when a strong external force is applied to the probe (active nonlinear microrheology). Calculations within mode coupling theory are presented for glasses of hard spheres and compared to Langevin and Brownian dynamics simulations. Under not too strong forces where the probe remains trapped, the probe density distribution becomes a… ▽ More We investigate the displacements of a probe particle inside a glass, when a strong external force is applied to the probe (active nonlinear microrheology). Calculations within mode coupling theory are presented for glasses of hard spheres and compared to Langevin and Brownian dynamics simulations. Under not too strong forces where the probe remains trapped, the probe density distribution becomes anisotropic. It is shifted towards the direction of the force, develops an enhanced tail in that direction (signalled by a positive skewness), and exhibits different variances along and perpendicular to the force direction. A simple model of an harmonically trapped probe rationalizes the low force limit, with strong strain softening setting in at forces of the order of a few thermal energies per particle radius. △ Less

Submitted 6 July, 2012; originally announced July 2012.

Asymptotic analysis of mode-coupling theory of active nonlinear microrheology

Authors: Manuel Victor Gnann, Thomas Voigtmann

Abstract: We discuss a schematic model of mode-coupling theory for force-driven active nonlinear microrheology, where a single probe particle is pulled by a constant external force through a dense host medium. The model exhibits both a glass transition for the host, and a force-induced delocalization transition, where an initially localized probe inside the glassy host attains a nonvanishing steady-state ve… ▽ More We discuss a schematic model of mode-coupling theory for force-driven active nonlinear microrheology, where a single probe particle is pulled by a constant external force through a dense host medium. The model exhibits both a glass transition for the host, and a force-induced delocalization transition, where an initially localized probe inside the glassy host attains a nonvanishing steady-state velocity by locally melting the glass. Asymptotic expressions for the transient density correlation functions of the schematic model are derived, valid close to the transition points. There appear several nontrivial time scales relevant for the decay laws of the correlators. For the nonlinear friction coeffcient of the probe, the asymptotic expressions cause various regimes of power-law variation with the external force, and two-parameter scaling laws. △ Less

Submitted 11 April, 2012; originally announced April 2012.

Comments: 17 pages, 12 figures

Journal ref: Phys. Rev. E 86, 011406 (2012)

Schematic Mode Coupling Theory of Glass Rheology: Single and Double Step Strains

Authors: Th. Voigtmann, J. M. Brader, M. Fuchs, M. E. Cates

Abstract: Mode coupling theory (MCT) has had notable successes in addressing the rheology of hard-sphere colloidal glasses, and also soft colloidal glasses such as star-polymers. Here, we explore the properties of a recently developed MCT-based schematic constitutive equation under idealized experimental protocols involving single and then double step strains. We find strong deviations from expectations bas… ▽ More Mode coupling theory (MCT) has had notable successes in addressing the rheology of hard-sphere colloidal glasses, and also soft colloidal glasses such as star-polymers. Here, we explore the properties of a recently developed MCT-based schematic constitutive equation under idealized experimental protocols involving single and then double step strains. We find strong deviations from expectations based on factorable, BKZ-type constitutive models. Specifically, a nonvanishing stress remains long after the application of two equal and opposite step strains; this residual stress is a signature of plastic deformation. We also discuss the distinction between hypothetically instantaneous step strains and fast ramps. These are not generally equivalent in our MCT approach, with the latter more likely to capture the physics of experimental `step' strains. The distinction points to the different role played by reversible anelastic, and irreversible plastic rearrangements. △ Less

Submitted 28 March, 2012; originally announced March 2012.

Journal ref: Soft Matter, 2012, 8 (15), 4244 - 4253

Long-Wavelength Anomalies in the Asymptotic Behavior of Mode-Coupling Theory

Authors: S. K. Schnyder, F. Höfling, T. Franosch, Th. Voigtmann

Abstract: We discuss the dynamic behavior of a tagged particle close to a classical localization transition in the framework of the mode-coupling theory of the glass transition. Asymptotic results are derived for the order parameter as well as the dynamic correlation functions and the mean-squared displacement close to the transition. The influence of an infrared cutoff is discussed. We discuss the dynamic behavior of a tagged particle close to a classical localization transition in the framework of the mode-coupling theory of the glass transition. Asymptotic results are derived for the order parameter as well as the dynamic correlation functions and the mean-squared displacement close to the transition. The influence of an infrared cutoff is discussed. △ Less

Submitted 11 January, 2011; originally announced January 2011.

Comments: 15 pages, 8 figures, to appear in J Phys Condens Matter

Journal ref: Journal of Physics: Condensed Matter 23(23), 234121 (2011)

Schematic Models for Active Nonlinear Microrheology

Authors: Manuel V. Gnann, Igor Gazuz, Antonio M. Puertas, Matthias Fuchs, Thomas Voigtmann

Abstract: We analyze the nonlinear active microrheology of dense colloidal suspensions using a schematic model of mode-coupling theory. The model describes the strongly nonlinear behavior of the microscopic friction coefficient as a function of applied external force in terms of a delocalization transition. To probe this regime, we have performed Brownian dynamics simulations of a system of quasi-hard spher… ▽ More We analyze the nonlinear active microrheology of dense colloidal suspensions using a schematic model of mode-coupling theory. The model describes the strongly nonlinear behavior of the microscopic friction coefficient as a function of applied external force in terms of a delocalization transition. To probe this regime, we have performed Brownian dynamics simulations of a system of quasi-hard spheres. We also analyze experimental data on hard-sphere-like colloidal suspensions [Habdas et al., Europhys. Lett., 2004, 67, 477]. The behavior at very large forces is addressed specifically. △ Less

Submitted 14 October, 2010; originally announced October 2010.

Journal ref: Soft Matter, 7, 1390-1396 (2011)

Multiple Glasses in Asymmetric Binary Hard Spheres

Authors: Th. Voigtmann

Abstract: Multiple distinct glass states occur in binary hard-sphere mixtures with constituents of very disparate sizes according to the mode-coupling theory of the glass transition (MCT), distinguished by considering whether small particles remain mobile or not, and whether small particles contribute significantly to perturb the big-particle structure or not. In the idealized glass, the four different glas… ▽ More Multiple distinct glass states occur in binary hard-sphere mixtures with constituents of very disparate sizes according to the mode-coupling theory of the glass transition (MCT), distinguished by considering whether small particles remain mobile or not, and whether small particles contribute significantly to perturb the big-particle structure or not. In the idealized glass, the four different glasses are separated by sharp transitions that give rise to higher-order transition phenomena involving logarithmic decay laws, and to anomalous power-law-like diffusion. The phenomena are argued to be expected generally in glass-forming mixtures. △ Less

Submitted 3 October, 2010; originally announced October 2010.

Comments: 4 pages, 2 figures

Localization phenomena in models of ion-conducting glass formers

Authors: Jürgen Horbach, Thomas Voigtmann, Felix Höfling, Thomas Franosch

Abstract: The mass transport in soft-sphere mixtures of small and big particles as well as in the disordered Lorentz gas (LG) model is studied using molecular dynamics (MD) computer simulations. The soft-sphere mixture shows anomalous small-particle diffusion signifying a localization transition separate from the big-particle glass transition. Switching off small-particle excluded volume constraints slows d… ▽ More The mass transport in soft-sphere mixtures of small and big particles as well as in the disordered Lorentz gas (LG) model is studied using molecular dynamics (MD) computer simulations. The soft-sphere mixture shows anomalous small-particle diffusion signifying a localization transition separate from the big-particle glass transition. Switching off small-particle excluded volume constraints slows down the small-particle dynamics, as indicated by incoherent intermediate scattering functions. A comparison of logarithmic time derivatives of the mean-squared displacements reveals qualitative similarities between the localization transition in the soft-sphere mixture and its counterpart in the LG. Nevertheless, qualitative differences emphasize the need for further research elucidating the connection between both models. △ Less

Submitted 4 June, 2010; originally announced June 2010.

Comments: to appear in Eur. Phys. J. Special Topics

Journal ref: Eur. Phys. J. Special Topics 189, 141-145 (2010)

From Equilibrium to Steady-State Dynamics after Switch-On of Shear

Authors: Matthias Krüger, Fabian Weysser, Thomas Voigtmann

Abstract: A relation between equilibrium, steady-state, and waiting-time dependent dynamical two-time correlation functions in dense glass-forming liquids subject to homogeneous steady shear flow is discussed. The systems under study show pronounced shear thinning, i.e., a significant speedup in their steady-state slow relaxation as compared to equilibrium. An approximate relation that recovers the exact li… ▽ More A relation between equilibrium, steady-state, and waiting-time dependent dynamical two-time correlation functions in dense glass-forming liquids subject to homogeneous steady shear flow is discussed. The systems under study show pronounced shear thinning, i.e., a significant speedup in their steady-state slow relaxation as compared to equilibrium. An approximate relation that recovers the exact limit for small waiting times is derived following the integration through transients (ITT) approach for the nonequilibrium Smoluchowski dynamics, and is exemplified within a schematic model in the framework of the mode-coupling theory of the glass transition (MCT). Computer simulation results for the tagged-particle density correlation functions corresponding to wave vectors in the shear-gradient directions from both event-driven stochastic dynamics of a two-dimensional hard-disk system and from previously published Newtonian-dynamics simulations of a three-dimensional soft-sphere mixture are analyzed and compared with the predictions of the ITT-based approximation. Good qualitative and semi-quantitative agreement is found. Furthermore, for short waiting times, the theoretical description of the waiting time dependence shows excellent quantitative agreement to the simulations. This confirms the accuracy of the central approximation used earlier to derive fluctuation dissipation ratios (Phys. Rev. Lett. 102, 135701). For intermediate waiting times, the correlation functions decay faster at long times than the stationary ones. This behavior is predicted by our theory and observed in simulations. △ Less

Submitted 1 April, 2010; originally announced April 2010.

Comments: 16 pages, 12 figures, submitted to Phys Rev E

Structural relaxation of polydisperse hard spheres: comparison of the mode-coupling theory to a Langevin dynamics simulation

Authors: F. Weysser, A. M. Puertas, M. Fuchs, Th. Voigtmann

Abstract: We analyze the slow, glassy structural relaxation as measured through collective and tagged-particle density correlation functions obtained from Brownian dynamics simulations for a polydisperse system of quasi-hard spheres in the framework of the mode-coupling theory of the glass transition (MCT). Asymptotic analyses show good agreement for the collective dynamics when polydispersity effects are t… ▽ More We analyze the slow, glassy structural relaxation as measured through collective and tagged-particle density correlation functions obtained from Brownian dynamics simulations for a polydisperse system of quasi-hard spheres in the framework of the mode-coupling theory of the glass transition (MCT). Asymptotic analyses show good agreement for the collective dynamics when polydispersity effects are taken into account in a multi-component calculation, but qualitative disagreement at small $q$ when the system is treated as effectively monodisperse. The origin of the different small-$q$ behaviour is attributed to the interplay between interdiffusion processes and structural relaxation. Numerical solutions of the MCT equations are obtained taking properly binned partial static structure factors from the simulations as input. Accounting for a shift in the critical density, the collective density correlation functions are well described by the theory at all densities investigated in the simulations, with quantitative agreement best around the maxima of the static structure factor, and worst around its minima. A parameter-free comparison of the tagged-particle dynamics however reveals large quantiative errors for small wave numbers that are connected to the well-known decoupling of self-diffusion from structural relaxation and to dynamical heterogeneities. While deviations from MCT behaviour are clearly seen in the tagged-particle quantities for densities close to and on the liquid side of the MCT glass transition, no such deviations are seen in the collective dynamics. △ Less

Submitted 15 October, 2010; v1 submitted 19 March, 2010; originally announced March 2010.

Comments: 23 pages, 26 figures

Journal ref: Phys. Rev. E 82, 011504 (2010)

Glass Rheology: From mode-coupling theory to a dynamical yield criterion

Authors: J. M. Brader, Th. Voigtmann, M. Fuchs, R. G. Larson, M. E. Cates

Abstract: The mode coupling theory (MCT) of glasses, while offering an incomplete description of glass transition physics, represents the only established route to first-principles prediction of rheological behavior in nonergodic materials such as colloidal glasses. However, the constitutive equations derivable from MCT are somewhat intractable, hindering their practical use and also their interpretation.… ▽ More The mode coupling theory (MCT) of glasses, while offering an incomplete description of glass transition physics, represents the only established route to first-principles prediction of rheological behavior in nonergodic materials such as colloidal glasses. However, the constitutive equations derivable from MCT are somewhat intractable, hindering their practical use and also their interpretation. Here, we present a schematic (single-mode) MCT model which incorporates the tensorial structure of the full theory. Using it, we calculate the dynamic yield surface for a large class of flows. △ Less

Submitted 28 September, 2009; originally announced September 2009.

Journal ref: Proc.Natl.Acad.Sci.USA, 106 15186 (2009)

A Double-Transition Scenario for Anomalous Diffusion in Glass-Forming Mixtures

Authors: Th. Voigtmann, J. Horbach

Abstract: We study by molecular dynamics computer simulation a binary soft-sphere mixture that shows a pronounced decoupling of the species' long-time dynamics. Anomalous, power-law-like diffusion of small particles arises, that can be understood as a precursor of a double-transition scenario, combining a glass transition and a separate small-particle localization transition. Switching off small-particle… ▽ More We study by molecular dynamics computer simulation a binary soft-sphere mixture that shows a pronounced decoupling of the species' long-time dynamics. Anomalous, power-law-like diffusion of small particles arises, that can be understood as a precursor of a double-transition scenario, combining a glass transition and a separate small-particle localization transition. Switching off small-particle excluded-volume constraints slows down, rather than enhances, small-particle transport. The data are contrasted with results from the mode-coupling theory of the glass transition. △ Less

Submitted 18 December, 2008; originally announced December 2008.

Active and Nonlinear Microrheology in Dense Colloidal Suspensions

Authors: I. Gazuz, A. M. Puertas, Th. Voigtmann, M. Fuchs

Abstract: We present a first-principles theory for the active nonlinear microrheology of colloidal model systems: for constant external force on a spherical probe particle embedded in a dense host dispersion, neglecting hydrodynamic interactions, we derive an exact expression for the friction. Within mode-coupling theory (MCT), we discuss the threshold external force needed to delocalize the probe from a… ▽ More We present a first-principles theory for the active nonlinear microrheology of colloidal model systems: for constant external force on a spherical probe particle embedded in a dense host dispersion, neglecting hydrodynamic interactions, we derive an exact expression for the friction. Within mode-coupling theory (MCT), we discuss the threshold external force needed to delocalize the probe from a host glass, and its relation to strong nonlinear velocity-force curves in a host fluid. Experimental microrheology data and simulations, which we performed, are explained with a simplified model. △ Less

Submitted 15 October, 2008; originally announced October 2008.

Journal ref: Physical Review Letters 102, 248302 (2009)

Idealized glass transitions under pressure: dynamics versus thermodynamics

Authors: Th. Voigtmann

Abstract: The interplay of slow dynamics and thermodynamic features of dense liquids is studied by examinining how the glass transition changes depending on the presence or absence of Lennard-Jones-like attractions. Quite different thermodynamic behavior leaves the dynamics unchanged, with important consequences for high-pressure experiments on glassy liquids. Numerical results are obtained within mode-co… ▽ More The interplay of slow dynamics and thermodynamic features of dense liquids is studied by examinining how the glass transition changes depending on the presence or absence of Lennard-Jones-like attractions. Quite different thermodynamic behavior leaves the dynamics unchanged, with important consequences for high-pressure experiments on glassy liquids. Numerical results are obtained within mode-coupling theory (MCT), but the qualitative features are argued to hold more generally. A simple square-well model can be used to explain generic features found in experiment. △ Less

Submitted 5 August, 2008; originally announced August 2008.

Comments: to be published in Phys. Rev. Lett

From Equilibrium to Steady State: The Transient Dynamics of Colloidal Liquids under Shear

Authors: J. Zausch, J. Horbach, M Laurati, S. U. Egelhaaf, J. M. Brader, Th. Voigtmann, M. Fuchs

Abstract: We investigate stresses and particle motion during the start up of flow in a colloidal dispersion close to arrest into a glassy state. A combination of molecular dynamics simulation, mode coupling theory and confocal microscopy experiment is used to investigate the origins of the widely observed stress overshoot and (previously not reported) super-diffusive motion in the transient dynamics. A li… ▽ More We investigate stresses and particle motion during the start up of flow in a colloidal dispersion close to arrest into a glassy state. A combination of molecular dynamics simulation, mode coupling theory and confocal microscopy experiment is used to investigate the origins of the widely observed stress overshoot and (previously not reported) super-diffusive motion in the transient dynamics. A link between the macro-rheological stress versus strain curves and the microscopic particle motion is established. Negative correlations in the transient auto-correlation function of the potential stresses are found responsible for both phenomena, and arise even for homogeneous flows and almost Gaussian particle displacements. △ Less

Submitted 24 July, 2008; originally announced July 2008.

Comments: 24 pages, 14 figures, J. Phys.: Condens. Matter, in press

The Dynamics of Silica Melts under High Pressure: Mode-Coupling Theory Results

Authors: Th. Voigtmann, J. Horbach

Abstract: The high-pressure dynamics of a computer-modeled silica melt is studied in the framework of the mode-coupling theory of the glass transition (MCT) using static-structure input from molecular-dynamics (MD) computer simulation. The theory reproduces the experimentally known viscosity minimum (diffusivity maximum) as a function of density or pressure and explains it in terms of a corresponding mini… ▽ More The high-pressure dynamics of a computer-modeled silica melt is studied in the framework of the mode-coupling theory of the glass transition (MCT) using static-structure input from molecular-dynamics (MD) computer simulation. The theory reproduces the experimentally known viscosity minimum (diffusivity maximum) as a function of density or pressure and explains it in terms of a corresponding minimum in its critical temperature. This minimum arises from a gradual change in the equilibrium static structure which shifts from being dominated by tetrahedral ordering to showing the cageing known from high-density liquids. The theory is in qualitative agreement with computer simulation results. △ Less

Submitted 14 April, 2008; originally announced April 2008.

Comments: Presented at ESF EW Glassy Liquids under Pressure, to be published in Journal of Physics