Source author record

Jürgen Horbach

Jürgen Horbach appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

22works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

Brittle yielding in supercooled liquids below the critical temperature of mode coupling theory

Molecular Dynamics (MD) computer simulations of a polydisperse soft-sphere model under shear are presented. Starting point for these simulations are deeply supercooled samples far below the critical temperature, $T_c$, of mode coupling theory. These samples are fully equilibrated with the aid of the swap Monte Carlo technique. For states below $T_c$, we identify a life time $τ_{\rm lt}$ that measures the time scale on which the system can be considered as an amorphous solid. The temperature dependence of $τ_{\rm lt}$ can be well described by an Arrhenius law. The existence of transient amorphous solid states below $T_c$ is associated with the possibility of brittle yielding, as manifested by a sharp stress drop in the stress-strain relation and shear banding. We show that brittle yielding requires on the one hand low shear rates and on the other hand, the time scale corresponding to the inverse shear rate has to be smaller or of the order of $τ_{\rm lt}$. Both conditions can be only met for large life time $τ_{\rm lt}$, i.e.~for states far below $T_c$.

preprint2022arXiv

Interatomic machine learning potentials for aluminium: application to solidification phenomena

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.

preprint2020arXiv

Flow heterogeneities in supercooled liquids and glasses under shear

Using extensive non-equilibrium molecular dynamics simulations, we investigate a glassforming binary Lennard-Jones mixture under shear. Both supercooled liquids and glasses are considered. Our focus is on the characterization of inhomogeneous flow patterns such as shear bands that appear as a transient response to the external shear. For the supercooled liquids, we analyze the crossover from Newtonian to non-Newtonian behavior with increasing shear rate $\dotγ$. Above a critical shear rate $\dotγ_c$ where a non-Newtonian response sets in, the transient dynamics are associated with the occurrence of short-lived vertical shear bands, i.e. bands of high mobility that form perpendicular to the flow direction. In the glass states, long-lived horizontal shear bands, i.e. bands of high mobility parallel to the flow direction, are observed in addition to vertical ones. The systems with shear bands are characterized in terms of mobility maps, stress-strain relations, mean-squared displacements, and (local) potential energies. The initial formation of a horizontal shear band provides an efficient stress release, corresponds to a local minimum of the potential energy, and is followed by a slow broadening of the band towards the homogeneously flowing fluid in the steady-state. Whether a horizontal or a vertical shear band forms cannot be predicted from the initial undeformed sample. Furthermore, we show that with increasing system size the probability for the occurrence of horizontal shear bands increases.

preprint2020arXiv

Wrinkles, folds and ripplocations: unusual deformation structures of confined elastic sheets at non-zero temperatures

We study the deformation of a fluctuating crystalline sheet confined between two flat rigid walls as a simple model for layered solids where bonds among atoms {\it within} the same layer are much stronger than those {\it between} layers. When subjected to sufficiently high loads in an appropriate geometry, these solids deform and fail in unconventional ways. Recent experiments suggest that configurations named {\it ripplocations}, where a layer folds backwards over itself, are involved. These structures are distinct and separated by large free energy barriers from smooth {\it ripples} of the atomic layers that are always present at any non-zero temperature. We use Monte Carlo simulation in combination with an umbrella sampling technique to obtain conditions under which such structures form and study their specific experimental signatures.

preprint2019arXiv

A nucleation theory for yielding of nearly defect-free crystals: understanding rate dependent yield points

Experiments and simulations show that when an initially defect free rigid crystal is subjected to deformation at a constant rate, irreversible plastic flow commences at the so-called {\em yield point}. The yield point is a weak function of the deformation rate, which is usually expressed as a power law with an extremely small non-universal exponent. We re-analyze a representative set of published data on nanometer sized, mostly defect free, Cu, Ni and Au crystals in the light of a recently proposed theory of yielding based on nucleation of stable stress-free regions inside the metastable rigid solid. The single relation derived here, which is {\em not} a power law, explains data covering {\em fifteen} orders of magnitude in time scales.

preprint2016arXiv

Heterogeneous dynamics during yielding of glasses: effect of aging

Molecular dynamics computer simulations of a binary Lennard-Jones glass under shear are presented. The mechanical response of glassy states having different thermal histories is investigated by imposing a wide range of external shear rates, at different temperatures. The stress-strain relations exhibit an overshoot at a strain of around 0.1, marking the yielding of the glass sample and the onset of plastic flow. The amplitude of the overshoot shows a logarithmic behavior with respect to a dimensionless variable, given by the age of the sample times the shear rate. Dynamical heterogeneities having finite lifetimes, in the form of shear bands, are observed as the glass deforms under shear. By quantifying the spatial fluctuations of particle mobility, we demonstrate that such shearbanding occurs only under specific combinations of imposed shear-rate, age of glass and ambient temperature.

preprint2015arXiv

Crystal growth kinetics in Lennard-Jones and Weeks-Chandler-Andersen systems along the solid-liquid coexistence line

Kinetics of crystal-growth is investigated along the solid-liquid coexistence line for the (100), (110) and (111) orientations of the Lennard-Jones and Weeks-Chandler-Andersen fcc crystal-liquid interface, using non-equilibrium molecular dynamics simulations. A slowing down of the growth kinetics along the coexistence line is observed, which is mostly a temperature effect, with other quantities such as the melting pressure and liquid self-diffusion coefficient having a negligible impact. The growth kinetics of the two potentials become similar at large values of the melting temperature and pressure, when both resemble a purely repulsive soft-sphere potential. Classical models of crystallization from the melt are in reasonable qualitative agreement with our simulation data. Finally, several one-phase empirical melting/freezing rules are studied with respect to their validity along the coexistence line.

preprint2015arXiv

Exploring structural inhomogeneities in glasses during cavitation

Using large-scale molecular dynamics simulations for a system of $10^6$ particles, the response of a dense amorphous solid to the continuous expansion of its volume is investigated. We find that the spatially uniform glassy state becomes unstable via the formation of cavities, which eventually leads to failure. By scanning through a wide range of densities and temperatures, we determine the state points at which the instability occurs and thereby provide estimates of the co-existence density of the resultant glass phase. Evidence for long-lived, inhomogeneous configurations with a negative pressure is found, where the frozen-in glass structure contains spherical cavities or a network of void space. Furthermore, we demonstrate the occurrence of hysteretic effects when the cavitated solid is compressed to regain the dense glassy state. As a result, a new glass state is obtained, the pressure of which is different from the initial one due to small density inhomogeneities that are generated by the dilation-compression cycle.

preprint2015arXiv

Formation and growth of shear bands in glasses: existence of an underlying directed percolation transition

The response of glasses to mechanical loading often leads to the formation of inhomogeneous flow patterns that strongly affect materials properties. Among them, shear bands are ubiquitous in a wide variety of materials, ranging from soft matter systems to metallic alloys. Shear banding is associated with strain localization, i.e. the deformation of the sheared glassy solid is localized in space in form of band-like structures. These structures are often precursors to catastrophic failure, implying that a proper understanding of the underlying mechanisms could lead to the design of smarter materials. However, despite its importance in material science, the microscopic origin of shear banding in glassy materials is only poorly understood. Here, the formation of shear banding in glassy systems is revealed by non-equilibrium molecular dynamics simulations (NEMD) of a binary Lennard-Jones mixture, subject to a constant strain rate. In its glass state, this system exhibits for all considered strain rates the formation of a percolating cluster of mobile regions at a critical strain. We show that this percolation transition belongs to the universality class of directed percolation. Only at low shear rates, where the steady-state stress is close to the yielding threshold, the percolating cluster evolves into a transient (but long-lived) shear band with a diffusive growth of its width.

preprint2015arXiv

Free energy cost of forming an interface between a crystal and its frozen version

Using a thermodynamic integration scheme, we compute the free energy cost per unit area, $γ$, of forming an interface between a crystal and a frozen structured wall, formed by particles frozen into the same equilibrium structure as the crystal. Even though the structure and potential energy of the crystalline phase in the vicinity of the wall is same as in the bulk, $γ$ has a non-zero value and increases with increasing density of the crystal and the wall. Investigating the effect of several interaction potentials between the particles, we observe a positive $γ$ at all crystalline densities if the potential is purely repulsive. For models with attractive interactions, such as the Lennard-Jones potential, a negative value for $γ$ is obtained at low densities.

preprint2014arXiv

Crystal-liquid interfacial free energy of hard spheres via a novel thermodynamic integration scheme

The hard sphere crystal-liquid interfacial free energy, ($γ_{\rm cl}$), is determined from molecular dynamics simulations using a novel thermodynamic integration (TI) scheme. The advantage of this TI scheme compared to previous methods is to successfully circumvent hysteresis effects due to the movement of the crystal-liquid interface. This is accomplished by the use of extremely short-ranged and impenetrable Gaussian flat walls which prevent the drift of the interface while imposing a negligible free-energy penalty. We find that it is crucial to analyze finite-size effects in order to obtain reliable estimates of $γ_{\rm cl}$ in the thermodynamic limit.

preprint2014arXiv

Crystal-liquid interfacial free energy via thermodynamic integration

A novel thermodynamic integration (TI) scheme is presented to compute the crystal-liquid interfacial free energy ($γ_{\rm cl}$) from molecular dynamics simulation. The scheme is applied to a Lennard-Jones system. By using extremely short-ranged and impenetrable Gaussian flat walls to confine the liquid and crystal phases, we overcome hysteresis problems of previous TI schemes that stem from the translational movement of the crystal-liquid interface. Our technique is applied to compute $γ_{\rm cl}$ for the (100), (110) and (111) orientation of the crystalline phase at three temperatures under coexistence conditions. For one case, namely the (100) interface at the temperature $T=1.0$ (in reduced units), we demonstrate that finite-size scaling in the framework of capillary wave theory can be used to estimate $γ_{\rm cl}$ in the thermodynamic limit. Thereby, we show that our TI scheme is not associated with the suppression of capillary wave fluctuations.

preprint2014arXiv

Excess free energy of supercooled liquids at disordered walls

Using a novel thermodynamic integration scheme, we compute the excess free energy, $γ$, of a glass-forming, binary Lennard-Jones liquid in contact with a frozen amorphous wall, formed by particles frozen into a similar structure as the liquid. We find that $γ$ is non-zero, becoming negative at low temperature. This indicates that the thermodynamics of the system is perturbed by the effect of the amorphous wall.

preprint2014arXiv

Onset of flow in a confined colloidal glass under an imposed shear stress

A confined colloidal glass, under the imposition of a uniform shear stress, is investigated using numerical simulations. Both at macro- and microscales, the consequent dynamics during the onset of flow is studied. When the imposed stress is gradually decreased, the time scale for the onset of steady flow diverges, associated with long-lived spatial heterogeneities. Near this yield-stress regime, persistent creep in the form of shear-banded structures is observed.

preprint2014arXiv

Poiseuille flow of soft glasses in narrow channels: From quiescence to steady state

Using numerical simulations, the onset of Poiseuille flow in a confined soft glass is investigated. Starting from the quiescent state, steady flow sets in at a time scale which increases with a decrease in applied forcing. At this onset time scale, a rapid transition occurs via the simultaneous fluidization of regions having different local stresses. In the absence of steady flow at long times, creep is observed even in regions where the local stress is larger than the bulk yielding threshold. Finally, we show that the time scale to attain steady flow depends strongly on the history of the initial state.

preprint2014arXiv

Rounding of the localization transition in model porous media

The generic mechanisms of anomalous transport in porous media are investigated by computer simulations of two-dimensional model systems. In order to bridge the gap between the strongly idealized Lorentz model and realistic models of porous media, two models of increasing complexity are considered: a cherry-pit model with hard-core correlations as well as a soft-potential model. An ideal gas of tracer particles inserted into these structures is found to exhibit anomalous transport which extends up to several decades in time. Also, the self-diffusion of the tracers becomes suppressed upon increasing the density of the systems. These phenomena are attributed to an underlying percolation transition. In the soft potential model the transition is rounded, since each tracer encounters its own critical density according to its energy. Therefore, the rounding of the transition is a generic occurrence in realistic, soft systems.

preprint2013arXiv

Lennard-Jones systems near solid walls: Computing interfacial free energies from molecular simulation methods

Different computational techniques in combination with molecular dynamics computer simulation are used to to determine the wall-liquid and the wall-crystal interfacial free energies of a modified Lennard-Jones (LJ) system in contact with a solid wall. Two different kinds of solid walls are considered: a flat structureless wall and a structured wall consisting of an ideal crystal with the particles rigidly attached to fcc lattice sites. Interfacial free energies are determined by a thermodynamic integration scheme, the anisotropy of the pressure tensor, the non-equilibrium work method based on Bennett acceptance criteria, and a method using Cahn's adsorption equations based on the interfacial thermodynamics of Gibbs. For the flat wall, interfacial free energies as a function of different densities of the LJ liquid and as a function of temperature along the coexistence curve are calculated. In case of a structured wall, the interaction strength between the wall and the LJ system and the lattice constant of the structured wall are varied. Using the values of the wall-liquid and wall-crystal interfacial energies along with the value for the crystal-liquid interfacial free energy determined previously for the same system by the "cleaving potential method", we obtain the contact angle as a function of various parameters; in particular the conditions are found under which partial wetting occurs.

preprint2012arXiv

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

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.

preprint2012arXiv

Tension and stiffness of the hard sphere crystal-fluid interface

A combination of fundamental measure density functional theory and Monte Carlo computer simulation is used to determine the orientation-resolved interfacial tension and stiffness for the equilibrium hard-sphere crystal-fluid interface. Microscopic density functional theory is in quantitative agreement with simulations and predicts a tension of 0.66 kT/σ^2 with a small anisotropy of about 0.025 kT and stiffnesses with e.g. 0.53 kT/σ^2 for the (001) orientation and 1.03 kT/σ^2 for the (111) orientation. Here kT is denoting the thermal energy and σthe hard sphere diameter. We compare our results with existing experimental findings.

preprint2012arXiv

Wall-liquid and wall-crystal interfacial free energies via thermodynamic integration: A molecular dynamics simulation study

A method is proposed to compute the interfacial free energy of a Lennard-Jones system in contact with a structured wall by molecular dynamics simulation. Both the bulk liquid and bulk face-centered-cubic crystal phase along the (111) orientation are considered. Our approach is based on a thermodynamic integration scheme where first the bulk Lennard-Jones system is reversibly transformed to a state where it interacts with a structureless flat wall. In a second step, the flat structureless wall is reversibly transformed into an atomistic wall with crystalline structure. The dependence of the interfacial free energy on various parameters such as the wall potential, the density and orientation of the wall is investigated. The conditions are indicated under which a Lennard-Jones crystal partially wets a flat wall.

preprint2010arXiv

Localization phenomena in models of ion-conducting glass formers

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.

preprint2006arXiv

Electrophoretic Properties of Highly Charged Colloids: A Hybrid MD/LB Simulation Study

Using computer simulations, the electrophoretic motion of a positively charged colloid (macroion) in an electrolyte solution is studied in the framework of the primitive model. Hydrodynamic interactions are fully taken into account by applying a hybrid simulation scheme, where the charged ions (i.e. macroion and electrolyte), propagated via molecular dynamics (MD), are coupled to a Lattice Boltzmann (LB) fluid. In a recent experiment it was shown that, for multivalent salt ions, the mobility $μ$ initially increases with charge density $σ$, reaches a maximum and then decreases with further increase of $σ$. The aim of the present work is to elucidate the behaviour of $μ$ at high values of $σ$. Even for the case of monovalent microions, we find a decrease of $μ$ with $σ$. A dynamic Stern layer is defined that includes all the counterions that move with the macroion while subject to an external electrical field. The number of counterions in the Stern layer, $q_0$, is a crucial parameter for the behavior of $μ$ at high values of $σ$. In this case, the mobility $μ$ depends primarily on the ratio $q_0/Q$ (with $Q$ the valency of the macroion). The previous contention that the increase in the distortion of the electric double layer (EDL) with increasing $σ$ leads to the lowering of $μ$ does not hold for high $σ$. In fact, we show that the deformation of the EDL decreases with increase of $σ$. The role of hydrodynamic interactions is inferred from direct comparisons to Langevin simulations where the coupling to the LB fluid is switched off. Moreover, systems with divalent counterions are considered. In this case, at high values of $σ$ the phenomenon of charge inversion is found.