Source author record

Jean-Louis Barrat

Jean-Louis Barrat 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

43works
9topics
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

43 published item(s)

preprint2026arXiv

Mechanisms of Spatiotemporal Damage Evolution in Double Polymer Networks

Double polymer networks exhibit a striking enhancement of toughness compared to single networks, yet the microscopic mechanisms governing stress redistribution, damage evolution, and fracture remain incompletely understood. Using large-scale coarse-grained molecular dynamics simulations under uniaxial deformation, we resolve bond scission statistics, local stress redistribution following individual bond-breaking events, and the spatiotemporal evolution of damage in single- and double-network architectures. We show that while the early mechanical response is dominated by the pre-stretched sacrificial network, damage evolution in double networks follows a qualitatively distinct pathway. In contrast to single networks, where anisotropic stress redistribution promotes rapid localization and catastrophic fracture, the presence of a soft matrix in double networks induces a screening of stress redistribution generated by sacrificial bond scission. This screening suppresses correlated rupture events and stabilizes multiple damage zones, leading to a strongly delocalized damage landscape over a broad deformation range. At larger strains, when the matrix becomes load-bearing, damage progressively localizes, ultimately triggering fracture. By isolating the dynamics of individual damage zones, we further demonstrate that matrix-mediated stress screening stabilizes defects and delays localization. Together, these results identify stress-screening-induced damage delocalization as a central microscopic mechanism underlying toughness enhancement in multiple-network elastomers.

preprint2022arXiv

Temperature Dependence of Thermodynamic, Dynamical and Dielectric Properties of Water Models

We investigate the temperature dependence of thermodynamic (density, isobaric heat capacity), dynamical (self-diffusion coefficient, shear viscosity), and dielectric properties of several water models, the commonly employed TIP3P water model, well-established 4-point water model TIP4P-2005, and recently developed 4-point water model TIP4P-D. We focus on the temperature range of interest for the field of computational biophysics and soft matter (280-350 K). The 4-point water models lead to a spectacularly improved agreement with experimental data, strongly suggesting that the use of more modern parametrizations should be favored compared to the more traditional TIP3P for modeling temperature-dependent phenomena in biomolecular systems.

preprint2019arXiv

Elastic avalanches reveal marginal behaviour in amorphous solids

Mechanical deformation of amorphous solids can be described as consisting of an elastic part in which the stress increases linearly with strain, up to a yield point at which the solid either fractures or starts deforming plastically. It is well established, however, that the apparent linearity of stress with strain is actually a proxy for a much more complex behavior, with a microscopic plasticity that is reflected in diverging nonlinear elastic coefficients. Very generally, the complex structure of the energy landscape is expected to induce a singular response to small perturbations. In the athermal quasistatic regime, this response manifests itself in the form of a scale free plastic activity. The distribution of the corresponding avalanches should reflect, according to theoretical mean field calculations (Franz and Spigler, Phys. Rev. E., 2017, 95, 022139), the geometry of phase space in the vicinity of a typical local minimum. In this work, we characterize this distribution for simple models of glass forming systems, and we find that its scaling is compatible with the mean field predictions for systems above the jamming transition. These systems exhibit marginal stability, and scaling relations that hold in the stationary state are examined and confirmed in the elastic regime. By studying the respective influence of system size and age, we suggest that marginal stability is systematic in the thermodynamic limit.

preprint2016arXiv

Cut-off nonlinearities in the low-temperature vibrations of glasses and crystals

We present a computer simulation study of glassy and crystalline states using the standard Lennard-Jones interaction potential that is truncated at a finite cut-off distance, as is typical of many computer simulations. We demonstrate that the discontinuity at the cut-off distance in the first derivative of the potential (corresponding to the interparticle force) leads to the appearance of cut-off nonlinearities. These cut-off nonlinearities persist into the very-low-temperature regime thereby affecting low-temperature thermal vibrations, which leads to a breakdown of the harmonic approximation for many eigen modes, particularly for low-frequency vibrational modes. Furthermore, while expansion nonlinearities which are due to higher order terms in the Taylor expansion of the interaction potential are usually ignored at low temperatures and show up as the temperature increases, cut-off nonlinearities can become most significant at the lowest temperatures. Anharmonic effects readily show up in the elastic moduli which not only depend on the eigen frequencies, but are crucially sensitive to the eigen vectors of the normal modes. Whereas, those observables that rely mainly on static structural information or just the eigen frequencies, such as the vibrational density of states, total potential energy, and specific heat, show negligible dependence on the presence of the cut-off. Similar aspects of nonlinear behavior have recently been reported in model granular materials, where the constituent particles interact through finite-range, purely-repulsive potentials. These nonlinearities have been ascribed to the nature of the sudden cut-off at contact in the force-law, thus we demonstrate that cut-off nonlinearities emerge as a general feature of ordered and disordered solid state systems interacting through truncated potentials.

preprint2016arXiv

Effects of inertia on the steady-shear rheology of disordered solids

We study the finite-shear-rate rheology of disordered solids by means of molecular dynamics simulations in two dimensions. By systematically varying the damping magnitude $ζ$ in the low-temperature limit, we identify two well defined flow regimes, separated by a thin (temperature-dependent) crossover region. In the overdamped regime, the athermal rheology is governed by the competition between elastic forces and viscous forces, whose ratio gives the Weissenberg number $Wi= ζ\dotγ$ (up to elastic parameters); the macroscopic stress $Σ$ follows the frequently encountered Herschel-Bulkley law $Σ= Σ\_0 + k \sqrt{Wi}$, with yield stress $Σ\_0\textgreater{}0$. In the underdamped (inertial) regime, dramatic changes in the rheology are observed for low damping: the flow curve becomes non-monotonic. This change is not caused by longer-lived correlations in the particle dynamics at lower damping; instead, for weak dissipation, the sample heats up considerably due to, and in proportion to, the driving. By suitably thermostatting more or less underdamped systems, we show that their rheology only depends on their kinetic temperature and the shear rate, rescaled with Einstein's vibration frequency.

preprint2016arXiv

Plastic response and correlations in athermally sheared amorphous solids

The onset of irreversible deformation in low-temperature amorphous solids is due to the accumulation of elementary events, consisting of spacially and temporally localized atomic rearrangements involving only a few tens of atoms. Recently, numerical and experimental work addressed the issue of spatio-temporal correlations between these plastic events. Here, we provide further insight into these correlations by investigating, via molecular dynamics (MD) simulations, the plastic response of a two-dimensional amorphous solid to artificially triggered local shear transformations. We show that while the plastic response is virtually absent in as-quenched configurations, it becomes apparent if a shear strain was previously imposed on the system. Plastic response has a four-fold symmetry which is characteristic of the shear stress redistribution following the local transformation. At high shear rate we report evidence for a fluctuation-dissipation relation, connecting plastic response and correlation, which seems to break down if lower shear rates are considered.

preprint2016arXiv

Relation of vibrational excitations and thermal conductivity to elastic heterogeneities in disordered solids

In crystals, molecules thermally vibrate around the periodic lattice sites. Vibrational motions are well understood in terms of phonons, which carry heat and control heat transport. The situation is notably different in disordered solids, where vibrational excitations are not phonons and can be even localized. Recent numerical work has established the concept of elastic heterogeneity: disordered solids show inhomogeneous local mechanical response. Clearly, the heterogeneous nature of elastic properties strongly influences vibrational and thermal properties, and it is expected to be the origin of anomalous features, including boson peak, vibrational localization, and temperature dependence of thermal conductivity. These are all crucial long-standing problems in materials physics, which we address in the present work. We have considered a toy model able to stabilize different states of matter, by introducing an increasing amount of size disorder. The phase diagram generated by molecular dynamics simulations encompasses the perfect crystalline state with a spatially homogeneous elastic moduli distribution, multiple defective phases with increasing moduli heterogeneities, and eventually a series of amorphous states. We have established clear correlations among the heterogeneous local mechanical response, vibrational states, and thermal conductivity. We provide evidence that elastic heterogeneity controls both vibrational and thermal properties, and is a key concept to understand the anomalous puzzling features of disordered solids.

preprint2015arXiv

Beating the amorphous limit in thermal conductivity by superlattices design

The value measured in the amorphous structure with the same chemical composition is often considered as a lower bound for the thermal conductivity of any material: the heat carriers are strongly scattered by disorder, and their lifetimes reach the minimum time scale of thermal vibrations. An appropriate design at the nano-scale, however, may allow one to reduce the thermal conductivity even below the amorphous limit. In the present contribution, using molecular-dynamics simulation and the Green-Kubo formulation, we study systematically the thermal conductivity of layered phononic materials (superlattices), by tuning different parameters that can characterize such structures. We discover that the key to reach a lower-than-amorphous thermal conductivity is to block almost completely the propagation of the heat carriers, the superlattice phonons. We demonstrate that a large mass difference in the two intercalated layers, or weakened interactions across the interface between layers result in materials with very low thermal conductivity, below the values of the corresponding amorphous counterparts.

preprint2015arXiv

Driving rate dependence of avalanche statistics and shapes at the yielding transition

We study stress time series caused by plastic avalanches in athermally sheared disordered materials. Using particle-based simulations and a mesoscopic elasto-plastic model, we analyze size and shear-rate dependence of the stress-drop durations and size distributions together with their average temporal shape. We find critical exponents different from mean-field predictions, and a clear asymmetry for individual avalanches. We probe scaling relations for the rate dependency of the dynamics and we report a crossover towards mean-field results for strong driving.

preprint2015arXiv

Edwards thermodynamics for a driven athermal system with dry friction

We obtain, using semi-analytical transfer operator techniques, the Edwards thermodynamics of a one-dimensional model of blocks connected by harmonic springs and subjected to dry friction. The theory is able to reproduce the linear divergence of the correlation length as a function of energy density observed in direct numerical simulations of the model under tapping dynamics. We further characterize analytically this divergence using a Gaussian approximation for the distribution of mechanically stable configurations, and show that it is related to the existence of a peculiar infinite temperature critical point.

preprint2015arXiv

Elastic consequences of a single plastic event: towards a realistic account of structural disorder and shear wave propagation in models of flowing amorphous solids

Shear transformations (i.e., localised rearrangements of particles resulting in the shear deformation of a small region of the sample) are the building blocks of mesoscale models for the flow of disordered solids. In order to compute the time-dependent response of the solid material to such a shear transformation, with a proper account of elastic heterogeneity and shear wave propagation, we propose and implement a very simple Finite-Element (FE) -based method. Molecular Dynamics (MD) simulations of a binary Lennard-Jones glass are used as a benchmark for comparison, and information about the microscopic viscosity and the local elastic constants is directly extracted from the MD system and used as input in FE. We find very good agreement between FE and MD regarding the temporal evolution of the disorder-averaged displacement field induced by a shear transformation, which turns out to coincide with the response of a uniform elastic medium. However, fluctuations are relatively large, and their magnitude is satisfactorily captured by the FE simulations of an elastically heterogeneous system. Besides, accounting for elastic anisotropy on the mesoscale is not crucial in this respect. The proposed method thus paves the way for models of the rheology of amorphous solids which are both computationally efficient and realistic, in that structural disorder and inertial effects are accounted for.

preprint2015arXiv

On the relevance of disorder in athermal amorphous materials under shear

We show that, at least at a mean-field level, the effect of structural disorder in sheared amorphous media is very dissimilar depending on the thermal or athermal nature of their underlying dynamics. We first introduce a toy model, including explicitly two types of noise (thermal versus athermal). Within this interpretation framework, we argue that mean-field athermal dynamics can be accounted for by the so-called H{é}braud-Lequeux (HL) model, in which the mechanical noise stems explicitly from the plastic activity in the sheared medium. Then, we show that the inclusion of structural disorder, by means of a distribution of yield energy barriers, has no qualitative effect in the HL model, while such a disorder is known to be one of the key ingredients leading kinematically to a finite macroscopic yield stress in other mean-field descriptions, such as the Soft-Glassy-Rheology model. We conclude that the statistical mechanisms at play in the emergence of a macroscopic yield stress, and a complex stationary dynamics at low shear rate, are different in thermal and athermal amorphous systems.

preprint2014arXiv

Elastic heterogeneity, vibrational states, and thermal conductivity across an amorphisation transition

Disordered solids exhibit unusual properties of their vibrational states and thermal conductivities. Recent progresses have well established the concept of "elastic heterogeneity", i.e., disordered materials show spatially inhomogeneous elastic moduli. In this study, by using molecular-dynamics simulations, we gradually introduce "disorder" into a numerical system to control its modulus heterogeneity. The system starts from a perfect crystalline state, progressively transforms into an increasingly disordered crystalline state, and finally undergoes structural amorphisation. We monitor independently the elastic heterogeneity, the vibrational states, and the thermal conductivity across this transition, and show that the heterogeneity in elastic moduli is well correlated to vibrational and thermal anomalies of the disordered system.

preprint2014arXiv

Generic transport coefficients of a confined electrolyte solution

Physical parameters characterising electrokinetic transport in a confined electrolyte solution are reconstructed from the generic transport coefficients obtained within the classical non-equilibrium statistical thermodynamic framework. The electro-osmotic flow, the diffusio-osmotic flow, the osmotic current, as well as the pressure-driven Poiseuille-type flow, the electric conduction, and the ion diffusion, are described by this set of transport coefficients. The reconstruction is demonstrated for an aqueous NaCl solution between two parallel charged surfaces with a nanoscale gap, by using the molecular dynamic (MD) simulations. A Green-Kubo approach is employed to evaluate the transport coefficients in the linear-response regime, and the fluxes induced by the pressure, electric, and chemical potential fields are compared with the results of non-equilibrium MD simulations. Using this numerical scheme, the influence of the salt concentration on the transport coefficients is investigated. Anomalous reversal of diffusio-osmotic current, as well as that of electro-osmotic flow, is observed at high surface charge densities and high added-salt concentrations.

preprint2014arXiv

Molecular dynamics simulation of electrokinetic flow of an aqueous electrolyte solution in nanochannels

Electrokinetic flows of an aqueous NaCl solution in nanochannels with negatively charged surfaces are studied using molecular dynamics (MD) simulations. The four transport coefficients that characterise the response to weak electric and pressure fields, namely the coefficients for the electrical current in response to the electric field ($M^{jj}$) and the pressure field ($M^{jm}$), and those for the mass flow in response to the same fields ($M^{mj}$ and $M^{mm}$), are obtained in the linear regime using a Green--Kubo approach. Nonequilibrium simulations with explicit external fields are also carried out, and the current and mass flows are directly obtained. The two methods exhibit good agreement even for large external field strengths, and Onsager's reciprocal relation ($M^{jm} = M^{mj}$) is numerically confirmed in both approaches. The influence of the surface charge density on the flow is also considered. The values of the transport coefficients are found to be smaller for larger surface charge density, because the counter-ions strongly bound near the channel surface interfere with the charge and mass flows. A reversal of the streaming current and of the reciprocal electro-osmotic flow, with a change of sign of $M^{mj}$ due to the excess co-ions, takes places for very high surface charge density.

preprint2014arXiv

Relaxation in yield stress systems through elastically interacting activated events

We study consequences of long-range elasticity in thermally assisted dynamics of yield stress materials. Within a two-dimensinal mesoscopic model we calculate the mean-square displacement and the dynamical structure factor for tracer particle trajectories. The ballistic regime at short time scales is associated with a compressed exponential decay in the dynamical structure factor, followed by a subdiffusive crossover prior to the onset of diffusion. We relate this crossover to spatiotemporal correlations and thus go beyond established mean field predictions.

preprint2014arXiv

Rheology of athermal amorphous solids: Revisiting simplified scenarios and the concept of mechanical noise temperature

We study the rheology of amorphous solids in the limit of negligible thermal fluctuations. On the basis of general arguments, the flow curve is shown to result from an interplay between the time scales of the macroscopic driving and the (cascades of) local particle rearrangements. Such rearrangements are known to induce a redistribution of the elastic stress in the system. Although mechanical noise, i.e., the local stress fluctuations arising from this redistribution, is widely believed to activate new particle rearrangements, we provide evidence that casts severe doubt on the analogy with thermal fluctuations: mechanical and thermal fluctuations lead to asymptotically different statistics for barrier crossing. These ideas are illustrated and supported by a simple coarse-grained model whose ingredients are directly connected with the physical processes relevant for the flow.

preprint2014arXiv

Spatiotemporal correlations between plastic events in the shear flow of athermal amorphous solids

The slow flow of amorphous solids exhibits striking heterogeneities: swift localised particle rearrangements take place in the midst of a more or less homogeneously deforming medium. Recently, experimental as well as numerical work has revealed spatial correlations between these flow heterogeneities. Here, we use molecular dynamics (MD) simulations to characterise the rearrangements and systematically probe their correlations both in time and in space. In particular, these correlations display a four-fold azimuthal symmetry characteristic of shear stress redistribution in an elastic medium and we unambiguously detect their increase in range with time. With increasing shear rate, correlations become shorter-ranged and more isotropic. In addition, we study a coarse-grained model motivated by the observed flow characteristics and challenge its predictions directly with the MD simulations. While the model captures both macroscopic and local properties rather satisfactorily, the agreement with respect to the spatiotemporal correlations is at most qualitative. The discrepancies provide important insight into relevant physics that is missing in all related coarse-grained models that have been developed for the flow of amorphous materials so far, namely the finite shear wave velocity and the impact of elastic heterogeneities on stress redistribution.

preprint2014arXiv

Time dependent elastic response to a local shear transformation in amorphous solids

The elastic response of a two-dimensional amorphous solid to induced local shear transformations, which mimic the elementary plastic events occurring in deformed glasses, is investigated via Molecular Dynamics simulations. We show that for different spatial realizations of the transformation, despite relative fluctuations of order one, the long time equilibrium response averages out to the prediction of the Eshelby inclusion problem for a continuum elastic medium. We characterize the effects of the underlying dynamics on the propagation of the elastic signal. A crossover from a propagative transmission in the case of weakly-damped dynamics to a diffusive transmission for strong damping is evidenced. In the latter case, the full time dependent elastic response is in agreement with the theoretical prediction, obtained by solving the diffusion equation for the displacement field in an elastic medium.

preprint2014arXiv

Universal and non-universal features in coarse-grained models of flow in disordered solids

We study the two-dimensional (2D) shear flow of amorphous solids within variants of an elastoplastic model, paying particular attention to spatial correlations and time fluctuations of, e.g., local stresses. The model is based on the local alternation between an elastic regime and plastic events during which the local stress is redistributed. The importance of a fully tensorial description of the stress and of the inclusion of (coarse-grained) convection in the model is investigated; scalar and tensorial models yield very similar results, while convection enhances fluctuations and breaks the spurious symmetry between the flow and velocity gradient directions, for instance when shear localisation is observed. Besides, correlation lengths measured with diverse protocols are discussed. One class of such correlation lengths simply scale with the spacing between homogeneously distributed, simultaneous plastic events. This leads to a scaling of the correlation length with the shear rate as $\dotγ^{\frac{-1}{2}}$ in 2D in the athermal regime, regardless of the details of the model. The radius of the cooperative disk, defined as the near-field region in which plastic events induce a stress redistribution that is not amenable to a mean-field treatment, notably follows this scaling. On the other hand, the cooperative volume measured from the four-point stress susceptibility and its dependence on the system size and the shear rate are model-dependent.

preprint2013arXiv

Computation of the thermal conductivity using classical and quantum molecular dynamics based methods

The thermal conductivity of a model for solid argon is investigated using nonequilibrium molecular dynamics methods, as well as the traditional Boltzmann transport equation approach with input from molecular dynamics calculations, both with classical and quantum thermostats. A surprising result is that, at low temperatures, only the classical molecular dynamics technique is in agreement with the experimental data. We argue that this agreement is due to a compensation of errors, and raise the issue of an appropriate method for calculating thermal conductivities at low (below Debye) temperatures.

preprint2013arXiv

Extension of Classical Nucleation Theory for Uniformly Sheared Systems

Nucleation is an out-of-equilibrium process, which can be strongly affected by the presence of external fields. In this letter, we report a simple extension of classical nucleation theory to systems submitted to an homogeneous shear flow. The theory involves accounting for the anisotropy of the critical nucleus formation, and introduces a shear rate dependent effective temperature. This extended theory is used to analyze the results of extensive molecular dynamics simulations, which explore a broad range of shear rates and undercoolings. At fixed temperature, a maximum in the nucleation rate is observed, when the relaxation time of the system is comparable to the inverse shear rate. In contrast to previous studies, our approach does not require a modification of the thermodynamic description, as the effect of shear is mainly embodied into a modification of the kinetic prefactor and of the temperature.

preprint2013arXiv

Measuring Spatial Distribution of Local Elastic Modulus in Glasses

Glasses exhibit spatially inhomogeneous elastic properties, which can be investigated by measuring their elastic moduli at a local scale. Various methods to evaluate the local elastic modulus have been proposed in the literature. A first possibility is to measure the local stress-local strain curve and to obtain the local elastic modulus from the slope of the curve, or equivalently to use a local fluctuation formula. Another possible route is to assume an affine strain and to use the applied global strain instead of the local strain for the calculation of the local modulus. Most recently a third technique has been introduced, which is easy to be implemented and has the advantage of low computational cost. In this contribution, we compare these three approaches by using the same model glass and reveal the differences among them caused by the non-affine deformations.

preprint2013arXiv

Nano-scale buckling in lamellar block polymers: a molecular dynamics simulation approach

Oriented block copolymers exhibit a buckling instability when submitted to a tensile test perpendicular to the lamellae direction. In this paper we study this behavior using a coarse grained molecular dynamics simulation approach. Coarse grained models of lamellar copolymers with alternate glassy rubbery layers are generated using the Radical Like Polymerization method, and their mechanical response is studied. For large enough systems, uniaxial tensile tests perpendicular to the direction of the lamellae reveal the occurrence of the buckling instability at low strain. The results that emerge from molecular simulation are compared to an elastic theory of the buckling instability introduced by Read and coworkers. At high strain rates, significant differences are observed between elastic theory and simulation results for the buckling strain and the buckling wavelength. We explain this difference by the strain rate dependence of the mechanical response. A simple model that takes into account the influence of the strain rate in the mechanical response is presented to rationalize the results at low and moderate strain rates. At very high strain rates, cavitation takes place in the rubbery phase of the sample and limits the validity of the approach.

preprint2013arXiv

Numerical study of a slip-link model for polymer melts and nanocomposites

We present a numerical study of the slip link model introduced by Likhtman for describing the dy- namics of dense polymer melts. After reviewing the technical aspects associated with the implemen- tation of the model, we extend previous work in several directions. The dependence of the relaxation modulus with the slip link density and the slip link stiffness is reported. Then the nonlinear rheolog- ical properties of the model, for a particular set of parameters, are explored. Finally, we introduce excluded volume interactions in a mean field such as manner in order to describe inhomogeneous systems, and we apply this description to a simple nanocomposite model. With this extension, the slip link model appears as a simple and generic model of a polymer melt, that can be used as an alternative to molecular dynamics for coarse grained simulations of complex polymeric systems.

preprint2013arXiv

On the Green-Kubo relationship of the liquid-solid friction coefficient

In this paper, we propose a new derivation for the Green-Kubo relationship for the liquid-solid friction coefficient, characterizing hydrodynamic slippage at a wall. It is based on a general Langevin approach for the fluctuating wall velocity, involving a non-markovian memory kernel with vanishing time integral. The calculation highlights some subtleties of the wall-liquid dynamics, leading to superdiffusive motion of the fluctuating wall position.

preprint2013arXiv

Spatial cooperativity in microchannel flows of soft jammed materials: A mesoscopic approach

The flow of amorphous solids results from a combination of elastic deformation and local structural rearrangements, which induce non-local elastic deformations. These elements are incorporated into a mechanically-consistent mesoscopic model of interacting elastoplastic blocks. We investigate the specific case of channel flow with numerical simulations, paying particular attention to situations of strong confinement. We find that the simple picture of plastic events embedded in an elastic matrix successfully accounts for manifestations of spatial cooperativity. Shear rate fluctuations are observed in seemingly quiescent regions, and the velocity profiles in confined flows at high applied pressure deviate from those expected in the absence of non-local effects, in agreement with experimental data. However, we suggest a different physical origin for the large deviations observed when walls have rough surfaces, associated with "bumps" of the particles against the asperities of the walls.

preprint2012arXiv

Nanoscale buckling deformation in layered copolymer materials

In layered materials, a common mode of deformation involves buckling of the layers under tensile deformation in the direction perpendicular to the layers. The instability mechanism, which operates in elastic materials from geological to nanometer scales, involves the elastic contrast between different layers. In a regular stacking of "hard" and "soft" layers, the tensile stress is first accommodated by a large deformation of the soft layers. The inhibited Poisson contraction results in a compressive stress in the direction transverse to the tensile deformation axis. The hard layers sustain this transverse compression until buckling takes place and results in an undulated structure. Using molecular simulations, we demonstrate this scenario for a material made of triblock copolymers. The buckling deformation is observed to take place at the nanoscale, at a wavelength that depends on strain rate. In contrast to what is commonly assumed, the wavelength of the undulation is not determined by defects in the microstructure. Rather, it results from kinetic effects, with a competition between the rate of strain and the growth rate of the instability. http://www.pnas.org/content/early/2011/12/23/1111367109.abstract

preprint2012arXiv

Shear flow of non-Brownian suspensions close to jamming

The dynamical mechanisms controlling the rheology of dense suspensions close to jamming are investigated numerically, using simplified models for the relevant dissipative forces. We show that the velocity fluctuations control the dissipation rate and therefore the effective viscosity of the suspension. These fluctuations are similar in quasi-static simulations and for finite strain rate calculations with various damping schemes. We conclude that the statistical properties of grain trajectories -- in particular the critical exponent of velocity fluctuations with respect to volume fraction ϕ-- only weakly depend on the dissipation mechanism. Rather they are determined by steric effects, which are the main driving forces in the quasistatic simulations. The critical exponent of the suspension viscosity with respect to ϕcan then be deduced, and is consistent with experimental data.

preprint2012arXiv

Spontaneous formation of permanent shear bands in a mesoscopic model of flowing disordered matter

This study proposes a coherent scenario of the formation of permanent shear bands in the flow of yield stress materials. It is a well accepted point of view that flow in disordered media is occurring via local plastic events, corresponding to small size rearrangements, that yield a long range stress redistribution over the system. Within a minimalistic mesoscopic model that incorporates these local dynamics, we study the spatial organisation of the local plastic events. The most important parameter in this study is the typical restructuring time needed to regain the original structure after a local rearrangement. In agreement with a recent mean field study [Coussot \textit{et al., Eur. Phys. J. E}, 2010, \textbf{33}, 183] we observe a spontaneous formation of permanent shear bands, when this restructuring time is large compared to the typical stress release time in a rearrangement. The bands consist of a large number of plastic events within a solid region that remains elastic. This heterogeneous flow behaviour is different in nature from the transient dynamical heterogeneities that one observes in the small shear rate limit in flow without shear-banding [Martens \textit{et al., Phys. Rev. Lett.}, 2011, \textbf{106}, 156001]. We analyse in detail the dependence of the shear bands on system size, shear rate and restructuring time. Further we rationalise the scenario within a mean field version of the spatial model, that produces a non monotonous flow curve for large restructuring times. This explains the instability of the homogeneous flow below a critical shear rate, that corresponds to the minimum of the curve. Our study therefore strongly supports the idea that the characteristic time scales involved in the local dynamics are at the physical origin of permanent shear bands.

preprint2011arXiv

Connecting diffusion and dynamical heterogeneities in actively deformed amorphous systems

We present an extensive numerical study of dynamical heterogeneities and their influence on diffusion in an athermal mesoscopic model for actively deformed amorphous solids. At low strain rates the stress dynamics are governed by cooperative regions of plastic events. On the basis of scaling arguments as well as an extensive numerical study of an athermal elasto-plastic model, we show that there is a direct link between the self-diffusion coefficient and the size of cooperative regions at low strain rates. Both depend strongly on rate and on system size. A measure of the mean square displacement of passive tracers in deformed amorphous media thus gives information about the microscopic rheology, such as the geometry of the cooperative regions and their scaling with strain rate and system size.

preprint2011arXiv

Effective temperatures of a heated Brownian particle

We investigate various possible definitions of an effective temperature for a particularly simple nonequilibrium stationary system, namely a heated Brownian particle suspended in a fluid. The effective temperature based on the fluctuation dissipation ratio depends on the time scale under consideration, so that a simple Langevin description of the heated particle is impossible. The short and long time limits of this effective temperature are shown to be consistent with the temperatures estimated from the kinetic energy and Einstein relation, respectively. The fluctuation theorem provides still another definition of the temperature, which is shown to coincide with the short time value of the fluctuation dissipation ratio.

preprint2011arXiv

Heat conduction across molecular junctions between nanoparticles

We investigate the problem of heat conduction across a molecular junction connecting two nanoparticles, both in vacuum and in a liquid environment, using classical molecular dynamics simulations. In vacuum, the well-known result of a length independent conductance is recovered; its precise value, however, is found to depend sensitively on the overlap between the vibrational spectrum of the junction and the density of states of the nanoparticles that act as thermal contacts. In a liquid environment, the conductance is constant up to a crossover length, above which a standard Fourier regime is recovered.

preprint2011arXiv

Portable implementation of a quantum thermal bath for molecular dynamics simulations

Recently, Dammak and coworkers (H. Dammak, Y. Chalopin, M. Laroche, M. Hayoun, and J.J. Greffet. Quantumthermal bath for molecular dynamics simulation. Phys. Rev. Lett., 103:190601, 2009.) proposed that the quantum statistics of vibrations in condensed systems at low temperature could be simulated by running molecular dynamics simulations in the presence of a colored noise with an appropriate power spectral density. In the present contribution, we show how this method can be implemented in a flexible manner and at a low computational cost by synthesizing the corresponding noise 'on the fly'. The proposed algorithm is tested for a simple harmonic chain as well as for a more realistic model of aluminium crystal. The energy and Debye-Waller factor are shown to be in good agreement with those obtained from harmonic approximations based on the phonon spectrum of the systems. The limitations of the method associated with anharmonic effects are also briefly discussed. Some perspectives for disordered materials and heat transfer are considered.

preprint2011arXiv

Predictors of cavitation in glassy polymers under tensile strain: a coarse grained molecular dynamics investigation

The nucleation of cavities in a homogenous polymer under tensile strain is investigated in a coarse-grained molecular dynamics simulation. In order to establish a causal relation between local microstructure and the onset of cavitation, a detailed analysis of some local properties is presented. In contrast to common assumptions, the nucleation of a cavity is neither correlated to a local loss of density nor, to the stress at the atomic scale and nor to the chain ends density in the undeformed state. Instead, a cavity in glassy polymers nucleates in regions that display a low bulk elastic modulus. This criterion allows one to predict the cavity position before the cavitation occurs. Even if the localization of a cavity is not directly predictable from the initial configuration, the elastically weak zones identified in the initial state emerge as favorite spots for cavity formation.

preprint2010arXiv

Fluctuations and correlations during the shear flow of elastic particles near the jamming transition

We present a numerical study of the flow of an assembly of frictionless soft discs at zero temperature, in the vicinity of and slightly above the jamming density. We find that some of the flow properties, such as the fluctuations in the number of contacts or the shear modulus, display a critical like behaviour that is governed by the proximity to the jamming point. Dynamical correlations during a quasistatic deformation, however, are non critical and dominated by system size. At finite strain rates, these dynamical correlations acquire a finite, strain-rate dependent amplitude, that decreases when approaching the jamming point from above.

preprint2010arXiv

Molecular dynamics simulations of glassy polymers

We review recent results from computer simulation studies of polymer glasses, from chain dynamics around the glass transition temperature Tg to the mechanical behaviour below Tg. These results clearly show that modern computer simulations are able to address and give clear answers to some important issues in the field, in spite of the obvious limitations in terms of length and time scales. In the present review we discuss the cooling rate effects, and dynamic slowing down of different relaxation processes when approaching Tg for both model and chemistry-specific polymer glasses. The impact of geometric confinement on the glass transition is discussed in detail. We also show that computer simulations are very useful tools to study structure and mechanical response of glassy polymers. The influence of large deformations on mechanical behaviour of polymer glasses in general, and strain hardening effect in particular are reviewed. Finally, we suggest some directions for future research, which we believe will be soon within the capabilities of state of the art computer simulations, and correspond to problems of fundamental interest.

preprint2010arXiv

Superdiffusive, heterogeneous, and collective particle motion near the jamming transition in athermal disordered materials

We use computer simulations to study the microscopic dynamics of an athermal assembly of soft particles near the fluid-to-solid, jamming transition. Borrowing tools developed to study dynamic heterogeneity near glass transitions, we discover a number of original signatures of the jamming transition at the particle scale. We observe superdiffusive, spatially heterogeneous, and collective particle motion over a characteristic scale which displays a surprising non-monotonic behavior across the transition. In the solid phase, the dynamics is an intermittent succession of elastic deformations and plastic relaxations, which are both characterized by scale-free spatial correlations and system size dependent dynamic susceptibilities. Our results show that dynamic heterogeneities in dense athermal systems and glass-formers are very different, and shed light on recent experimental reports of `anomalous' dynamical behavior near the jamming transition of granular and colloidal assemblies.

preprint2009arXiv

Local elasticity map and plasticity in a model Lennard-Jones glass

In this work we calculate the local elastic moduli in a weakly polydisperse 2DLennard-Jones glass undergoing a quasistatic shear deformation at zero temperature. The numerical method uses coarse grained microscopic expressions for the strain, displacement and stress fields. This method allows us to calculate the local elasticity tensor and to quantify the deviation from linear elasticity (local Hooke's law) at different coarse-graining scales. From the results a clear picture emerges of an amorphous material with strongly spatially heterogeneous elastic moduli that simultaneously satisfies Hooke's law at scales larger than a characteristic length scale of the order of five interatomic distances. At this scale the glass appears as a composite material composed of a rigid scaffoldingand of soft zones. Only recently calculated in non homogeneous materials, the local elastic structure plays a crucial role in the elasto-plastic response of the amorphous material. For a small macroscopic shear strain the structures associated with the non-affine displacement field appear directly related to the spatial structure of the elastic moduli. Moreover for a larger macroscopic shear strain we show that zones of low shear modulus concentrate most of the strain in form of plastic rearrangements. The spatio-temporal evolution of this local elasticity map and its connection with long term dynamical heterogeneity as well as with the plasticity in the material is quantified. The possibility to use this local parameter as a predictor of subsequent local plastic activity is also discussed.

preprint2009arXiv

The jamming transition as probed by quasistatic shear flow

We study the rheology of amorphous packings of soft, frictionless particles close to jamming. Implementing a quasistatic simulation method we generate a well defined ensemble of states that directly samples the system at its yield-stress. A continuous jamming transition from a freely-flowing state to a yield stress situation takes place at a well defined packing fraction, where the scaling laws characteristic of isostatic solids are observed. We propose that long-range correlations observed below the transition are dominated by this isostatic point, while those that are observed above the transition are characteristic of dense, disordered elastic media.

preprint2007arXiv

Driven activation versus thermal activation

Activated dynamics in a glassy system undergoing steady shear deformation is studied by numerical simulations. Our results show that the external driving force has a strong influence on the barrier crossing rate, even though the reaction coordinate is only weakly coupled to the nonequilibrium system. This "driven activation" can be quantified by introducing in the Arrhenius expression an effective temperature, which is close to the one determined from the fluctuation-dissipation relation. This conclusion is supported by analytical results for a simplified model system.

preprint2005arXiv

Diffusion in pores and its dependence on boundary conditions

We study the influence of the boundary conditions at the solid liquid interface on diffusion in a confined fluid. Using an hydrodynamic approach, we compute numerical estimates for the diffusion of a particle confined between two planes. Partial slip is shown to significantly influence the diffusion coefficient near a wall. Analytical expressions are derived in the low and high confinement limits, and are in good agreement with numerical results. These calculations indicate that diffusion of tagged particles could be used as a sensitive probe of the solid-liquid boundary conditions.

preprint1998arXiv

Influence of wetting properties on the hydrodynamic boundary condition at a fluid-solid interface

It is well known that, at a macroscopic level, the boundary condition for a viscous fluid at a solid wall is one of "no-slip". The liquid velocity field vanishes at a fixed solid boundary. In this paper, we consider the special case of a liquid that partially wets the solid, i.e. a drop of liquid in equilibrium with its vapor on the solid substrate has a finite contact angle. Using extensive Non-Equilibrium Molecular Dynamics (NEMD) simulations, we show that when the contact angle is large enough, the boundary condition can drastically differ (at a microscopic level) from a "no-slip" condition. Slipping lengths exceeding 30 molecular diameters are obtained for a contact angle of 140 degrees, characteristic of Mercury on glass. On the basis of a Kubo expression for $δ$, we derive an expression for the slipping length in terms of equilibrium quantities of the system. The predicted behaviour is in very good agreement with the numerical results for the slipping length obtained in the NEMD simulations. The existence of large slipping lentgh may have important implications for the transport properties in nanoporous media under such "nonwetting" conditions.