Showing 1–40 of 40 results for author: Gendelman, O

Unified physical framework for stretched-exponential, compressed-exponential, and logarithmic relaxation phenomena in glassy polymers

Authors: Valeriy V. Ginzburg, Oleg V. Gendelman, Alessio Zaccone

Abstract: We develop a simple yet comprehensive nonlinear model to describe relaxation phenomena in amorphous glass-formers near the glass transition temperature. The model is based on the two-state, two-(time)scale (TS2) framework, and describes the isothermal relaxation of specific volume, enthalpy, or shear stress via a simple first-order nonlinear differential equation (the Trachenko-Zaccone [TZ] equati… ▽ More We develop a simple yet comprehensive nonlinear model to describe relaxation phenomena in amorphous glass-formers near the glass transition temperature. The model is based on the two-state, two-(time)scale (TS2) framework, and describes the isothermal relaxation of specific volume, enthalpy, or shear stress via a simple first-order nonlinear differential equation (the Trachenko-Zaccone [TZ] equation) for local cooperative events. These nonlinear dynamics of cooperatively rearranging regions (CRR) naturally arise from the TS2 framework. We demonstrate that the solutions of the TZ equation comprehensively encompass the Debye exponential relaxation, the Kohlrausch-Williams-Watts (KWW) stretched and compressed relaxations, and the Guiu-Pratt logarithmic relaxation. Furthermore, for the case of stress relaxation modeling, our model recovers, as one of its limits, the Eyring law for plastic flow, where the Eyring activation volume is related to thermodynamic parameters of the material. Using the example of polystyrene (PS), we demonstrate how our model successfully describes the Kovacs "asymmetry of approach" specific volume and enthalpy experiments, as well as the stress relaxation. Other potential applications of the model, including the dielectric relaxation, are also discussed. The presented approach disentangles the physical origins of different relaxation laws within a single general framework based on the underlying physics. △ Less

Submitted 1 December, 2023; v1 submitted 15 November, 2023; originally announced November 2023.

Comments: 34 pages, 5 figures; submitted to Macromolecules

Hysteresis resulting from Lennard-Jones interactions

Authors: Dmitrii Rachinskii, Andrei Zagvozdkin, Oleg Gendelman

Abstract: The fundamental mechanism of hysteresis in the quasistatic limit of multi-stable systems is associated with transitions of the system from one local minimum of the potential energy to another. In this scenario, as system parameters are (quasistatically) varied, the transition is prompted when a saddle-node bifurcation eliminates the minimum where the system resides in. The objective of the present… ▽ More The fundamental mechanism of hysteresis in the quasistatic limit of multi-stable systems is associated with transitions of the system from one local minimum of the potential energy to another. In this scenario, as system parameters are (quasistatically) varied, the transition is prompted when a saddle-node bifurcation eliminates the minimum where the system resides in. The objective of the present work is to specify this generic mechanism for systems of interacting particles assuming a natural single-well (Lennard-Jones) interaction potential for each pair of particles. We show multi-stability and present details of hysteresis scenarios with the associated bifurcations and transitions in a case study of constrained four-degrees-of-freedom four particle systems on the plane. △ Less

Submitted 17 September, 2023; originally announced September 2023.

MSC Class: 34C55; 70F40

Effects of Asymmetric Cooling and Surface Wettability on the Orientation of the Freezing Tip

Authors: Anton Starostin, Vladimir Strelnikov, Leonid A. Dombrovsky, Shraga Shoval, Oleg Gendelman, Edward Bormashenko

Abstract: Freezing of water droplets placed on the bare and superhydrophobic surfaces of polymer wedges are studied both experimentally and computationally. Two-dimensional numerical calculations of the transient temperature field in a chilled polymer wedge show that the direction of heat flux from the droplet through the thermal contact region with the wedge differs significantly from the normal to the wed… ▽ More Freezing of water droplets placed on the bare and superhydrophobic surfaces of polymer wedges are studied both experimentally and computationally. Two-dimensional numerical calculations of the transient temperature field in a chilled polymer wedge show that the direction of heat flux from the droplet through the thermal contact region with the wedge differs significantly from the normal to the wedge surface. This is the physical cause of the recently observed asymmetric cooling of the droplet. A novel approximate computational model is proposed that takes into account the variable area of the water freezing front in the droplet. This model gives a quantitative estimate of the faster freezing of the droplet on the bare surface. The obtained numerical results agree with the data of laboratory experiments. The velocity of the crystallization front and the droplet deformation including the so-called freezing tip formation are monitored in the experiment. The direction of the freezing cone axis appears to be noticeably different for the cases of bare and superhydrophobic wedge surfaces. This deviation is explained by the fact that the direction of the freezing cone axis is controlled by the local direction of the heat flux. For a hydrophobic wedge surface, the deviation of the freezing tip from the vertical is smaller, because the reduced thermal contact area reduces the influence of the heat flux direction at the wedge surface. △ Less

Submitted 10 January, 2023; originally announced February 2023.

Comments: 13 pages, 7 figures

Kapitza resistance at a domain boundary in linear and nonlinear chains

Authors: Jithu Paul, O. V. Gendelman

Abstract: We explore Kapitza thermal resistance on the boundary between two homogeneous chain fragments with different characteristics. For a linear model, an exact expression for the resistance is derived, and well-defined in the thermodynamic limit. However, the resistance in this case depends on the thermostat properties and therefore is not a local property of the considered domain boundary. If the doma… ▽ More We explore Kapitza thermal resistance on the boundary between two homogeneous chain fragments with different characteristics. For a linear model, an exact expression for the resistance is derived, and well-defined in the thermodynamic limit. However, the resistance in this case depends on the thermostat properties and therefore is not a local property of the considered domain boundary. If the domains are nonlinear, but integrable - Toda lattice, elastically colliding particles - the anomalies are similar to the case of the linear chain, besides well-articulated thermal dependence of the resistance. For the case of elastically colliding particles, this dependence follows a simple scaling law - the resistance is proportional to the inverse square root of the temperature. For Fermi-Pasta-Ulam domains, both the temperature drop and the heat flux decrease with the chain length, but with different exponents, so the resistance vanishes in the thermodynamic limit. For the domains comprised of rotators, the thermal resistance exhibits the expected normal behavior. △ Less

Submitted 22 August, 2021; originally announced August 2021.

Comments: 31 pages, 23 figures

Kapitza thermal resistance in linear and nonlinear chain models: isotopic defect

Authors: O. V. Gendelman, Jithu Paul

Abstract: Kapitza resistance in the chain models with internal defects is considered. For the case of the linear chain, the exact analytic solution for the boundary resistance is derived for arbitrary linear time-independent conservative inclusion or defect. A simple case of isolated isotopic defects is explored in more detail. Contrary to the bulk conductivity in the linear chain, the Kapitza resistance is… ▽ More Kapitza resistance in the chain models with internal defects is considered. For the case of the linear chain, the exact analytic solution for the boundary resistance is derived for arbitrary linear time-independent conservative inclusion or defect. A simple case of isolated isotopic defects is explored in more detail. Contrary to the bulk conductivity in the linear chain, the Kapitza resistance is finite. However, the universal thermodynamic limit does not exist in this case. In other terms, the exact value of the resistance is not uniquely defined and depends on the way of approaching the infinite lengths of the chain fragments. For this reason, and also due to the explicit dependence on the parameters of the thermostats, the resistance cannot be considered as a local property of the defect. Asymptotic scaling behavior of the heat flux in the case of very heavy defect is explored and compared to the nonlinear counterparts; similarities in the scaling behavior are revealed. For the lightweight isotopic defect in the linear chain, one encounters a typical dip of the temperature profile, related to weak excitation of the localized mode in the attenuation zone. If the nonlinear interactions are included, this dip can still appear at a relatively short time scale, with subsequent elimination due to the nonlinear interactions. This observation implies that even in the nonlinear chains, the linear dynamics can predict the main features of the short-time evolution of the thermal profile if the temperature is low enough. △ Less

Submitted 13 May, 2021; v1 submitted 4 February, 2021; originally announced February 2021.

Journal ref: Phys. Rev. E 103, 052113 (2021)

Giant Amplification of Small Perturbations in Frictional Amorphous Solid

Authors: Harish Charan, Oleg Gendelman, Itamar Procaccia, Yarden Sheffer

Abstract: Catastrophic events in Nature can be often triggered by small perturbations, with "remote triggering" of earthquakes being an important example. Here we present a mechanism for the giant amplification of small perturbations that is expected to be generic in systems whose dynamics is not derivable from a Hamiltonian. We offer a general discussion of the typical instabilities involved (being oscilla… ▽ More Catastrophic events in Nature can be often triggered by small perturbations, with "remote triggering" of earthquakes being an important example. Here we present a mechanism for the giant amplification of small perturbations that is expected to be generic in systems whose dynamics is not derivable from a Hamiltonian. We offer a general discussion of the typical instabilities involved (being oscillatory with an exponential increase of noise) and examine in detail the normal forms that determine the relevant dynamics. The high sensitivity to external perturbations is explained for systems with and without dissipation. Numerical examples are provided using the dynamics of frictional granular matter. Finally we point out the relationship of the presently discussed phenomenon to the highly topical issue of "exceptional points" in quantum models with non-Hermitian Hamiltonians. △ Less

Submitted 19 January, 2020; originally announced January 2020.

Comments: 9 pages, o table, and 8 figures

Journal ref: Phys. Rev. E 101, 062902 (2020)

Breather arrest in a chain of damped oscillators with Hertzian contact

Authors: M. Strozzi, O. V. Gendelman

Abstract: We explore breather propagation in the damped oscillatory chain with essentially nonlinear (non-linearizable) nearest-neighbour coupling. Combination of the damping and the substantially nonlinear coupling leads to rather unusual two-stage pattern of the breather propagation. The first stage occurs at finite fragment of the chain and is characterized by power-law decay of the breather amplitude. T… ▽ More We explore breather propagation in the damped oscillatory chain with essentially nonlinear (non-linearizable) nearest-neighbour coupling. Combination of the damping and the substantially nonlinear coupling leads to rather unusual two-stage pattern of the breather propagation. The first stage occurs at finite fragment of the chain and is characterized by power-law decay of the breather amplitude. The second stage is characterized by extremely small breather amplitudes that decay hyper-exponentially with the site number. Thus, practically, one can speak about finite penetration depth of the breather. This phenomenon is referred to as breather arrest (BA). As particular example, we explore the chain with Hertzian contacts. Dependencies of the breather penetration depth on the initial excitation and on the damping coefficient on the breather penetration depth obey power laws. The results are rationalized by considering beating responses in a system of two damped linear oscillators with strongly nonlinear (non-linearizable) coupling. Initial excitation of one of these oscillators leads to strictly finite number of beating cycles. Then, the beating cycle in this simplified system is associated with the passage of the discrete breather between the neighbouring sites in the chain. Somewhat surprisingly, this simplified model reliably predicts main quantitative features of the breather arrest in the chain, including the exponents in numerically observed power laws. △ Less

Submitted 29 July, 2019; originally announced July 2019.

Kapitza resistance in basic chain models with isolated defects

Authors: Jithu Paul, O. V. Gendelman

Abstract: Kapitza thermal resistance is a common feature of material interfaces. It is defined as the ratio of the thermal drop at the interface to the heat flux flowing across the interface. One expects that this resistance will depend on the structure of the interface and on the temperature. We address the heat conduction in one-dimensional chain models with isotopic and/or coupling defects and explore th… ▽ More Kapitza thermal resistance is a common feature of material interfaces. It is defined as the ratio of the thermal drop at the interface to the heat flux flowing across the interface. One expects that this resistance will depend on the structure of the interface and on the temperature. We address the heat conduction in one-dimensional chain models with isotopic and/or coupling defects and explore the relationship between the interaction potentials and simulated properties of the Kapitza resistance. It is revealed that in linear models the Kapitza resistance is well-defined and size-independent (contrary to the bulk heat conduction coefficient), but depends on the parameters of thermostats used in the simulation. For $β$-FPU model one also encounters the dependence on the thermostats; in addition, the simulated boundary resistance strongly depends on the total system size. Finally, in the models characterized by convergent bulk heat conductivity (chain of rotators, Frenkel-Kontorova model) the boundary resistance is thermostat- and size-independent, as one expects. In linear chains, the Kapitza resistance is temperature-independent; thus, its temperature dependence allows one to judge on significance of the nonlinear interactions in the phonon scattering processes at the interface. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Nucleation and propagation of excitation fronts in self-excited systems

Authors: I. B. Shiroky, A. Papangelo, N. Hoffmann, O. V. Gendelman

Abstract: Frictional interfaces exhibits a very rich dynamical behavior due to frictional interactions. This work focuses on the transition between spatially localized and propagating stick-slip motion. To overcome the difficulties related to the non-smoothness of the friction law we have studied two models: (model I) a friction-excited chain of weakly coupled oscillators excited by a frictional moving belt… ▽ More Frictional interfaces exhibits a very rich dynamical behavior due to frictional interactions. This work focuses on the transition between spatially localized and propagating stick-slip motion. To overcome the difficulties related to the non-smoothness of the friction law we have studied two models: (model I) a friction-excited chain of weakly coupled oscillators excited by a frictional moving belt and (model II) a chain of Van der Pol oscillators. The two models show very similar dynamical features. Both exhibit discrete breathers (i.e. spatially localized periodic solutions) solutions for low coupling and show propagation along the chain of high amplitude limit cycles for stronger elastic coupling between the unit cells. In both models the transition from discrete breathers to propagating limit cycles happens through a nucleation process and, above a certain critical coupling, the velocity of propagation scales linearly with the elastic coupling coefficient. We suggest that dynamical features similar to what was demonstrated in this work are expected for models that are different quantitatively but preserve the same single-unit topology. △ Less

Submitted 16 April, 2019; originally announced April 2019.

Noise Amplification in Frictional Systems: Oscillatory Instabilities

Authors: Joyjit Chattoraj, Oleg Gendelman, Massimo Pica Ciamarra, Itamar Procaccia

Abstract: It was discovered recently that frictional granular materials can exhibit an important mechanism for instabilities, i.e the appearance of pairs of complex eigenvalues in their stability matrix. The consequence is an oscillatory exponential growth of small perturbations which are tamed by dynamical nonlinearities. The amplification can be giant, many orders of magnitude, and it ends with a catastro… ▽ More It was discovered recently that frictional granular materials can exhibit an important mechanism for instabilities, i.e the appearance of pairs of complex eigenvalues in their stability matrix. The consequence is an oscillatory exponential growth of small perturbations which are tamed by dynamical nonlinearities. The amplification can be giant, many orders of magnitude, and it ends with a catastrophic system-spanning plastic event. Here we follow up on this discovery, explore the scaling laws characterizing the onset of the instability, the scenarios of the development of the instability with and without damping, and the nature of the eventual system spanning events. The possible relevance to earthquake physics and to the transition from static to dynamic friction is discussed. △ Less

Submitted 25 March, 2019; originally announced March 2019.

Comments: arXiv admin note: text overlap with arXiv:1901.02376

Journal ref: Phys. Rev. E 100, 042901 (2019)

Oscillatory Instabilities in Frictional Granular Matter

Authors: Joyjit Chattoraj, Oleg Gendelman, Massimo Pica Ciamarra, Itamar Procaccia

Abstract: Frictional granular matter is shown to be fundamentally different in its plastic responses to external strains from generic glasses and amorphous solids without friction. While regular glasses exhibit plastic instabilities due to a vanishing of a real eigenvalue of the Hessian matrix, frictional granular materials can exhibit a previously unnoticed additional mechanism for instabilities, i.e. the… ▽ More Frictional granular matter is shown to be fundamentally different in its plastic responses to external strains from generic glasses and amorphous solids without friction. While regular glasses exhibit plastic instabilities due to a vanishing of a real eigenvalue of the Hessian matrix, frictional granular materials can exhibit a previously unnoticed additional mechanism for instabilities, i.e. the appearance of a pair of complex eigenvalues leading to oscillatory exponential growth of perturbations which are tamed by dynamical nonlinearities. This fundamental difference appears crucial for the understanding of plasticity and failure in frictional granular materials. The possible relevance to earthquake physics is discussed. △ Less

Submitted 9 January, 2019; v1 submitted 8 January, 2019; originally announced January 2019.

Comments: Supplementary materials follows the main text Animation of the instability: https://vimeo.com/310132278

Journal ref: Phys. Rev. Lett. 123, 098003 (2019)

Universal Scaling Laws for Shear Induced Dilation in Frictional Granular Media

Authors: Mahesh M. Bandi, Prasenjit Das, Oleg Gendelman, H. George E. Hentschel, Itamar Procaccia

Abstract: Compressed frictional granular matter cannot flow without dilation. Upon forced shearing to generate flow, the amount of dilation may depend on the initial preparation and a host of material variables. On the basis of both experiments and numerical simulations we show that as a result of training by repeated compression-decompression cycles the amount of dilation induced by shearing the system dep… ▽ More Compressed frictional granular matter cannot flow without dilation. Upon forced shearing to generate flow, the amount of dilation may depend on the initial preparation and a host of material variables. On the basis of both experiments and numerical simulations we show that as a result of training by repeated compression-decompression cycles the amount of dilation induced by shearing the system depends only on the shear rate and on the (pre-shearing) packing fraction. Relating the rheological response to structural properties allows us to derive a scaling law for the amount of dilation after $n$ cycles of compression-decompression. The resulting scaling law has a universal exponent that for trained systems is independent of the inter-granules force laws, friction parameters and strain rate. The amplitude of the scaling law is analytically computable, and it depends only on the shear rate and the asymptotic packing fraction. △ Less

Submitted 7 May, 2019; v1 submitted 30 August, 2018; originally announced August 2018.

Comments: 8 pages, 10 figures, Published Version

Journal ref: Gran. Mat. 40, 21 (2019)

Front propagation in lattices with on-site bistable non-degenerate potential: multiplicity, bifurcations and route to chaos

Authors: I. B. Shiroky, O. V. Gendelman

Abstract: Propagation of transition fronts in models of coupled oscillators with non-degenerate on-site potential is usually considered in terms of travelling waves. We show that the system dynamics can be reformulated as an implicit map structure, and the travelling waves correspond to stable fixed points. Therefore, the loss of stability of such waves should follow well-known generic bifurcation scenarios… ▽ More Propagation of transition fronts in models of coupled oscillators with non-degenerate on-site potential is usually considered in terms of travelling waves. We show that the system dynamics can be reformulated as an implicit map structure, and the travelling waves correspond to stable fixed points. Therefore, the loss of stability of such waves should follow well-known generic bifurcation scenarios. Then, one can expect a plethora of qualitatively different propagating-front solutions - multistable, multi-periodic, quasiperiodic and chaotic. △ Less

Submitted 25 October, 2017; originally announced October 2017.

Journal ref: Phys. Rev. E 98, 012220 (2018)

Reply to Comment on `What Determines the Static Force Chains in Stressed Granular Media?'

Authors: Oleg Gendelman, Yoav G. Pollack, Itamar Procaccia, Shiladitya Sengupta, Jacques Zylberg

Abstract: The formalism proposed in Phys. Rev. Lett. 116, 078001 (2016) for determination of the normal and tangential inter-particle forces in frictional disks from visual data was criticized in a comment to that paper. However the theory developed in the original publication is aimed at finding the forces when provided with accurate measurements while the comment addresses a different problem, being the s… ▽ More The formalism proposed in Phys. Rev. Lett. 116, 078001 (2016) for determination of the normal and tangential inter-particle forces in frictional disks from visual data was criticized in a comment to that paper. However the theory developed in the original publication is aimed at finding the forces when provided with accurate measurements while the comment addresses a different problem, being the statistics of forces in a configuration with uncertainty of the relevant parameters. We address this difference in perspectives and point out that even with experimental errors one can improve the results using an iterative procedure proposed in Phys. Rev. E. Rapid Communication, 93, 060601(R) (2016). △ Less

Submitted 26 June, 2017; originally announced June 2017.

Journal ref: Phys. Rev. Lett. 117, 159802 (2016)

Flat bands and compactons in mechanical lattices

Authors: Nathan Perchikov, O. V. Gendelman

Abstract: Local configurational symmetry in lattice structures may give rise to stationary, compact solutions, even in the absence of disorder and nonlinearity. These compact solutions are related to the existence of flat dispersion curves (bands). Nonlinearity can destabilize such compactons. One common flat-band-generating system is the 1D cross-stitch model, in which compactons were shown to exist for th… ▽ More Local configurational symmetry in lattice structures may give rise to stationary, compact solutions, even in the absence of disorder and nonlinearity. These compact solutions are related to the existence of flat dispersion curves (bands). Nonlinearity can destabilize such compactons. One common flat-band-generating system is the 1D cross-stitch model, in which compactons were shown to exist for the photonic lattice with Kerr nonlinearity. The compactons exist there already in the linear regime and are not generally destructed by that nonlinearity. Smooth nonlinearity of this kind does not allow doing a complete stability analysis for this chain. We consider a discrete mechanical system with flat dispersion bands, in which the nonlinearity exists due to impact constraints. There, one can use the concept of the saltation matrix for the analytic construction of the monodromy matrix. Besides, we consider a smooth nonlinear lattice with linearly connected massless boxes, each containing two symmetric anharmonic oscillators. In this model, the flat bands and discrete compactons also readily emerge. This system also permits performing comprehensive stability analysis, at least in the anti-continuum limit, due to the reduced number of degrees-of-freedom. In both systems, there exist two types of localization. The first one is the complete localization, and the second one is the more common exponential localization. The latter type is associated with discrete breathers (DBs). Two principal mechanisms for the loss of stability are revealed. The first one is the possible internal instability of the symmetric and/or antisymmetric solution in the individual unit cell of the chain. One can interpret this instability pattern as internal resonance between the compacton and the DB. The other mechanism is global instability related to resonance of the stationary solution with the propagation frequencies. △ Less

Submitted 9 November, 2017; v1 submitted 20 May, 2017; originally announced May 2017.

Comments: 26 pages, 27 figures

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

Propagation of transition fronts in nonlinear chains with non-degenerate on-site potentials

Authors: I. B. Shiroky, O. V. Gendelman

Abstract: We address the problem of a front propagation in chains with a bi-stable nondegenerate on-site potential and a nonlinear gradient coupling. For a generic nonlinear coupling, one encounters a special regime of transitions, characterized by extremely narrow fronts, far supersonic velocities of propagation and long waves in the oscillatory tail. This regime can be qualitatively associated with a shoc… ▽ More We address the problem of a front propagation in chains with a bi-stable nondegenerate on-site potential and a nonlinear gradient coupling. For a generic nonlinear coupling, one encounters a special regime of transitions, characterized by extremely narrow fronts, far supersonic velocities of propagation and long waves in the oscillatory tail. This regime can be qualitatively associated with a shock wave. The front propagation can be described with the help of a simple reduced-order model; the latter delivers a kinetic law, which is almost not sensitive to fine details of the on-site potential. Besides, it is possible to predict all main characteristics of the transition front, including its shape and frequency and amplitude of the oscillatory tail. The numerical results are in a good agreement with the analytical predictions. The suggested approach allows one to consider the effects of an external pre-load and on-site damping. When the damping is moderate, the analysis remains in the frame of the reduced-order model. It is possible to consider the solution for the front propagating in the damped chain as a perturbation of the undamped dynamics. This approach yield reasonable predictions. When the damping is high, the transition front enters a completely different asymptotic regime. The gradient nonlinearity generically turns negligible, and the propagating front converges to the exact solution obtained from a simple linear continuous model. △ Less

Submitted 16 May, 2017; originally announced May 2017.

Compact localized states and flatband generators in one dimension

Authors: Wulayimu Maimaiti, Alexei Andreanov, Hee Chul Park, Oleg Gendelman, Sergej Flach

Abstract: Flat bands (FB) are strictly dispersionless bands in the Bloch spectrum of a periodic lattice Hamiltonian, recently observed in a variety of photonic and dissipative condensate networks. FB Hamiltonians are finetuned networks, still lacking a comprehensive generating principle. We introduce a FB generator based on local network properties. We classify FB networks through the properties of compact… ▽ More Flat bands (FB) are strictly dispersionless bands in the Bloch spectrum of a periodic lattice Hamiltonian, recently observed in a variety of photonic and dissipative condensate networks. FB Hamiltonians are finetuned networks, still lacking a comprehensive generating principle. We introduce a FB generator based on local network properties. We classify FB networks through the properties of compact localized states (CLS) which are exact FB eigenstates and occupy $U$ unit cells. We obtain the complete two-parameter FB family of two-band $d=1$ networks with nearest unit cell interaction and $U=2$. We discover a novel high symmetry sawtooth chain with identical hoppings in a transverse dc field, easily accessible in experiments. Our results pave the way towards a complete description of FBs in networks with more bands and in higher dimensions. △ Less

Submitted 25 December, 2016; v1 submitted 5 October, 2016; originally announced October 2016.

Comments: 11 pages, 4 figures

Journal ref: Phys. Rev. B 95, 115135 (2017)

Heat conduction in a chain of colliding particles with stiff repulsive potential

Authors: O. V. Gendelman, A. V. Savin

Abstract: One-dimensional billiard, i.e. a chain of colliding particles with equal masses, is well-known example of completely integrable system. Billiards with different particles are generically not integrable, but still exhibit divergence of a heat conduction coefficient (HCC) in thermodynamic limit. Traditional billiard models imply instantaneous (zero-time) collisions between the particles. We lift thi… ▽ More One-dimensional billiard, i.e. a chain of colliding particles with equal masses, is well-known example of completely integrable system. Billiards with different particles are generically not integrable, but still exhibit divergence of a heat conduction coefficient (HCC) in thermodynamic limit. Traditional billiard models imply instantaneous (zero-time) collisions between the particles. We lift this condition and consider the heat transport in a chain of stiff colliding particles with power-law potential of the nearest-neighbor interaction. The instantaneous collisions correspond to the limit of infinite power in the interaction potential; for finite powers, the interactions take nonzero time. This modification of the model leads to profound physical consequence -- probability of multiple, in particular, triple particle collisions becomes nonzero. Contrary to the integrable billiard of equal particles, the modified model exhibits saturation of the heat conduction coefficient for large system size. Moreover, identification of scattering events with the triple particle collisions leads to simple definition of characteristic mean free path and kinetic description of the heat transport. This approach allows prediction both of temperature and density dependencies for the HCC limit values. The latter dependence is quite counterintuitive - the HCC is inversely proportional to the particle density in the chain. Both predictions are confirmed by direct numeric simulations. △ Less

Submitted 2 September, 2016; originally announced September 2016.

Comments: 6 pages, 7 figures

Determining the Inter-Particle Force-Laws in Amorphous Solids from a Visual Image

Authors: Oleg Gendelman, Yoav G. Pollack, Itamar Procaccia

Abstract: We consider the problem of how to determine the force laws in an amorphous system of interacting particles. Given the positions of the centers of mass of the constituent particles we propose a new algorithm to determine the inter-particle force-laws. Having $n$ different types of constituents we determine the coefficients in the Laurent polynomials for the $n(n+1)/2$ possibly different force-laws.… ▽ More We consider the problem of how to determine the force laws in an amorphous system of interacting particles. Given the positions of the centers of mass of the constituent particles we propose a new algorithm to determine the inter-particle force-laws. Having $n$ different types of constituents we determine the coefficients in the Laurent polynomials for the $n(n+1)/2$ possibly different force-laws. A visual providing the particle positions in addition to a measurement of the pressure is all that is required. The algorithm proposed includes a part that can correct for experimental errors in the positions of the particles. Such a correction of unavoidable measurement errors is expected to benefit many experiments in the field. △ Less

Submitted 3 May, 2016; v1 submitted 22 March, 2016; originally announced March 2016.

Comments: Very slight changes from original version. Mostly visual improvement and clarifications

Journal ref: Phys. Rev. E 93, 060601 (2016)

Normal heat conductivity in two-dimensional scalar lattices

Authors: A. V. Savin, V. Zolotarevskiy, O. V. Gendelman

Abstract: The paper revisits recent counterintuitive results on divergence of heat conduction coefficient in two-dimensional lattices. It was reported that in certain lattices with on-site potential, for which one-dimensional chain has convergent conductivity, for the 2D case it turns out to diverge. We demonstrate that this conclusion is an artifact caused by insufficient size of the simulated system. To o… ▽ More The paper revisits recent counterintuitive results on divergence of heat conduction coefficient in two-dimensional lattices. It was reported that in certain lattices with on-site potential, for which one-dimensional chain has convergent conductivity, for the 2D case it turns out to diverge. We demonstrate that this conclusion is an artifact caused by insufficient size of the simulated system. To overcome computational restrictions, a ribbon of relatively small width is simulated instead of the square specimen. It is further demonstrated that the heat conduction coefficient in the "long" direction of the ribbon ceases to depend on the width, as the latter achieves only 10 to 20 interparticle distances. So, one can consider the dynamics of much longer systems, than in the traditional setting, and still can gain a reliable information regarding the 2D lattice. It turns out that for all considered models, for which the conductivity is convergent in the 1D case, it is indeed convergent in the 2D case. In the same time, however, the length of the system, necessary to reveal the convergence in the 2D case, may be much bigger than in its 1D counterpart. △ Less

Submitted 14 November, 2015; originally announced November 2015.

Comments: 6 pages, 6 figures

A Generalized Electrowetting Equation: Its Derivation and Consequences

Authors: Edward Bormashenko, Oleg Gendelman

Abstract: The thermodynamics of electrowetting is treated. A general equation of electrowetting is derived from the first principles. It is demonstrated that the well-known Lippmann Equation describes a particular case of electrowetting when the radial derivative of the capacitance of the double layer is constant. The apparent contact angle of electrowetting depends on the gradient of capacity of a double l… ▽ More The thermodynamics of electrowetting is treated. A general equation of electrowetting is derived from the first principles. It is demonstrated that the well-known Lippmann Equation describes a particular case of electrowetting when the radial derivative of the capacitance of the double layer is constant. The apparent contact angle of electrowetting depends on the gradient of capacity of a double layer in the vicinity of the triple line. The role of the area adjacent to the triple line in constituting the equilibrium apparent contact angle of electrowetting is emphasized. △ Less

Submitted 5 November, 2015; originally announced November 2015.

Comments: 11 pages, 3 Figures

What determines static force chains in stressed granular media?

Authors: Oleg Gendelman, Yoav G. Pollack, Itamar Procaccia, Shiladitya Sengupta, Jacques Zylberg

Abstract: The determination of the normal and transverse (frictional) inter-particle forces within a granular medium is a long standing, taunting, and yet unresolved problem. We present a new formalism which employs the knowledge of the external forces and the orientations of contacts between particles (of any given sizes), to compute all the inter-particle forces. Having solved this problem we exemplify th… ▽ More The determination of the normal and transverse (frictional) inter-particle forces within a granular medium is a long standing, taunting, and yet unresolved problem. We present a new formalism which employs the knowledge of the external forces and the orientations of contacts between particles (of any given sizes), to compute all the inter-particle forces. Having solved this problem we exemplify the efficacy of the formalism showing that the force chains in such systems are determined by an expansion in the eigenfunctions of a newly defined operator. △ Less

Submitted 21 February, 2016; v1 submitted 25 May, 2015; originally announced May 2015.

Comments: v2: Supplementary Material added with examples of the M,T and Q matrices v3&v4: Figures improved and minor modifications

Journal ref: Phys. Rev. Lett. 116, 078001 (2016)

Heat conduction in diatomic chains with correlated disorder

Authors: Alexander V. Savin, Vadim Zolotarevskiy, Oleg V. Gendelman

Abstract: The paper considers heat transport in diatomic one-dimensional lattices, containing equal amounts of particles with different masses. Ordering of the particles in the chain is governed by single correlation parameter -- the probability for two neighboring particles to have the same mass. As this parameter grows from zero to unity, the structure of the chain varies from regular staggering chain to… ▽ More The paper considers heat transport in diatomic one-dimensional lattices, containing equal amounts of particles with different masses. Ordering of the particles in the chain is governed by single correlation parameter -- the probability for two neighboring particles to have the same mass. As this parameter grows from zero to unity, the structure of the chain varies from regular staggering chain to completely random configuration, and then -- to very long clusters of particles with equal masses. Therefore, this correlation parameter allows a control of typical cluster size in the chain. In order to explore different regimes of the heat transport, two interatomic potentials are considered. The first one is an infinite potential wall, corresponding to instantaneous elastic collisions between the neighboring particles. In homogeneous chains such interaction leads to an anomalous heat transport. The other one is classical Lennard-Jones interatomic potential, which leads to a normal heat transport. The simulations demonstrate that the correlated disorder of the particle arrangement does not change the convergence properties of the heat conduction coefficient, but essentially modifies its value. For the collision potential, one observes essential growth of the coefficient for fixed chain length as the limit of large homogeneous clusters is approached. The thermal transport in these models remains superdiffusive. In the Lennard-Jones chain the effect of correlation appears to be not monotonous in the limit of low temperatures. This behavior stems from the competition between formation of long clusters mentioned above, and Anderson localization close to the staggering ordered state. △ Less

Submitted 12 August, 2015; v1 submitted 23 March, 2015; originally announced March 2015.

Comments: 9 pages, 10 figures

Elastic Properties of Liquid Surfaces Coated with Colloidal Particles

Authors: Edward Bormashenko, Gene Whyman, Oleg Gendelman

Abstract: The physical mechanism of elasticity of liquid surfaces coated with colloidal particles is proposed. It is suggested that particles are separated by water clearings and the capillary interaction between them is negligible. The case is treated when the colloidal layer is deformed normally to its surface. The elasticity arises as an interfacial effect. The effective Young modulus of a surface depend… ▽ More The physical mechanism of elasticity of liquid surfaces coated with colloidal particles is proposed. It is suggested that particles are separated by water clearings and the capillary interaction between them is negligible. The case is treated when the colloidal layer is deformed normally to its surface. The elasticity arises as an interfacial effect. The effective Young modulus of a surface depends on the interfacial tension, equilibrium contact angle, radius of colloidal particles and their surface density. For the nanometrically scaled particles the line tension becomes essential and has an influence on the effective Young modulus. △ Less

Submitted 10 March, 2015; originally announced March 2015.

Comments: 20 pages, 7 figures

Heat conduction in a chain of dissociating particles: effect of dimensionality

Authors: V. Zolotarevskiy, A. V. Savin, O. V. Gendelman

Abstract: The paper considers heat conduction in a model chain of composite particles with hard core and elastic external shell. Such model mimics three main features of realistic interatomic potentials - hard repulsive core, quasilinear behavior in a ground state and possibility of dissociation. It has become clear recently, that this latter feature has crucial effect on convergence of the heat conduction… ▽ More The paper considers heat conduction in a model chain of composite particles with hard core and elastic external shell. Such model mimics three main features of realistic interatomic potentials - hard repulsive core, quasilinear behavior in a ground state and possibility of dissociation. It has become clear recently, that this latter feature has crucial effect on convergence of the heat conduction coefficient in thermodynamic limit. We demonstrate that in one-dimensional chain of elastic particles with hard core the heat conduction coefficient also converges, as one could expect. Then we explore effect of dimensionality on the heat transport in this model. For this sake, longitudinal and transversal motions of the particles are allowed in a long narrow channel. With varying width of the channel, we observe sharp transition from "one-dimensional" to "two-dimensional" behavior. Namely, the heat conduction coefficient drops by about order of magnitude for relatively small widening of the channel. This transition is not unique for the considered system. Similar phenomenon of transition to quasi-1D behavior with growth of aspect ratio of the channel is observed also in a gas of densely packed hard (billiard) particles, both for two- and three-dimensional cases. It is the case despite the fact that the character of transition in these two systems is not similar, due to different convergence properties of the heat conductivity. In the billiard model, the divergence of the heat conduction coefficient smoothly changes from logarithmic to power-like law with increase of the length. △ Less

Submitted 23 December, 2014; originally announced December 2014.

Comments: 11 pages, 15 figures

Journal ref: Physical Review E 91, 032127 (2015)

Shear Transformation Zones: State Determined or Protocol Dependent?

Authors: Oleg Gendelman, Prabhat K. Jaiswal, Itamar Procaccia, Bhaskar Sen Gupta, Jacques Zylberg

Abstract: The concept of a Shear Transformation Zone (STZ) refers to a region in an amorphous solid that undergoes a plastic event when the material is put under an external mechanical load. An important question that had accompanied the development of the theory of plasticity in amorphous solids for many years now is whether an STZ is a {\em region} existing in the material (which can be predicted by analy… ▽ More The concept of a Shear Transformation Zone (STZ) refers to a region in an amorphous solid that undergoes a plastic event when the material is put under an external mechanical load. An important question that had accompanied the development of the theory of plasticity in amorphous solids for many years now is whether an STZ is a {\em region} existing in the material (which can be predicted by analyzing the unloaded material), or is it an {\em event} that depends on the loading protocol (i.e., the event cannot be predicted without following the protocol itself). In this Letter we present strong evidence that the latter is the case. Infinitesimal changes of protocol result in macroscopically big jumps in the positions of plastic events, meaning that these can never be predicted from considering the unloaded material. △ Less

Submitted 18 August, 2014; originally announced August 2014.

Comments: 4 pages, 5 figures

Journal ref: Europhys. Lett. 109, 16002 (2015)

Elasticity and Plasticity in Stiff and Flexible Oligomeric Glasses

Authors: Oleg Gendelman, H. George E. Hentschel, Pankaj K. Mishra, Itamar Procaccia, Jacques Zylberg

Abstract: In this paper we focus on the mechanical properties of oligomeric glasses (waxes), employing a microscopic model that provides, via numerical simulations, information about the shear modulus of such materials, the failure mechanism via plastic instabilities and about the geometric responses of the oligomers themselves to a mechanical load. We present a microscopic theory that explains the numerica… ▽ More In this paper we focus on the mechanical properties of oligomeric glasses (waxes), employing a microscopic model that provides, via numerical simulations, information about the shear modulus of such materials, the failure mechanism via plastic instabilities and about the geometric responses of the oligomers themselves to a mechanical load. We present a microscopic theory that explains the numerically observed phenomena, including an exact theory of the shear modulus and of the plastic instabilities, both local and system spanning. In addition we present a model to explain the geometric changes in the oligomeric chains under increasing strains. △ Less

Submitted 7 July, 2014; originally announced July 2014.

Comments: 15 pages, 13 figures

Low Voltage Reversible Electrowetting Exploiting Lubricated Polymer Honeycomb Substrates

Authors: Edward Bormashenko, Roman Pogreb, Yelena Bormashenko, Roman Grynyov, Oleg Gendelman

Abstract: Low-voltage electrowetting-on-dielectric scheme realized with lubricated honeycomb polymer surfaces is reported. Polycarbonate honeycomb reliefs manufactured with the breath-figures self-assembly were impregnated with silicone and castor oils. The onset of the reversible electrowetting for silicone oil impregnated substrates occurred at 35 V, whereas for castor oil impregnated ones it took place a… ▽ More Low-voltage electrowetting-on-dielectric scheme realized with lubricated honeycomb polymer surfaces is reported. Polycarbonate honeycomb reliefs manufactured with the breath-figures self-assembly were impregnated with silicone and castor oils. The onset of the reversible electrowetting for silicone oil impregnated substrates occurred at 35 V, whereas for castor oil impregnated ones it took place at 80 V. The semi-quantitative analysis of electrowetting of impregnated surfaces is proposed. △ Less

Submitted 16 June, 2014; originally announced June 2014.

Comments: 14 pages, 5 figures

Journal ref: Applied Physics Letters 104 (2014) 171601

Normal heat conductivity in chains capable of dissociation

Authors: O. V. Gendelman, A. V. Savin

Abstract: The paper suggests a resolution for recent controversy over convergence of heat conductivity in one-dimensional chains with asymmetric nearest-neighbor potential. We conjecture that the convergence is promoted not by the mere asymmetry of the potential, but due to ability of the chain to dissociate. In other terms, the attractive part of the potential function should approach a finite value as the… ▽ More The paper suggests a resolution for recent controversy over convergence of heat conductivity in one-dimensional chains with asymmetric nearest-neighbor potential. We conjecture that the convergence is promoted not by the mere asymmetry of the potential, but due to ability of the chain to dissociate. In other terms, the attractive part of the potential function should approach a finite value as the distance between the neighbors grows. To clarify this point, we study the simplest model of this sort -- a chain of linearly elastic disks with finite diameter. If the distance between the disk centers exceeds their diameter, the disks cease to interact. Formation of gaps between the disks is the only possible mechanism for scattering of the oscillatory waves. Heat conduction in this system turns out to be convergent. Moreover, an asymptotic behavior of the heat conduction coefficient for the case of large densities and relatively low temperatures obeys simple Arrhenius-type law. In the limit of low densities, the heat conduction coefficient converges due to triple disk collisions. Numeric observations in both limits are grounded by analytic arguments. In a chain with Lennard-Jones nearest-neighbor potential the heat conductivity also saturates in a thermodynamic limit and the coefficient also scales according to the Arrhenius law for low temperatures. This finding points on a universal role played by the possibility of dissociation, as convergence of the heat conduction coefficient is considered. △ Less

Submitted 26 December, 2013; originally announced December 2013.

Percolating Plastic Failure as a Mechanism for Shear Softening in Amorphous Solids

Authors: Vijayakumar Chikkadi, Oleg Gendelman, Valery Ilyin, J Ashwin, Itamar Procaccia

Abstract: ``Shear softening" refers to the observed reduction in shear modulus when the stress on an amorphous solid is increased beyond the initial linear region. Careful numerical quasi-static simulations reveal an intimate relation between plastic failure and shear softening. The attaintment of the steady-state value of the shear modulus associated with plastic flow is identified with a percolation of th… ▽ More ``Shear softening" refers to the observed reduction in shear modulus when the stress on an amorphous solid is increased beyond the initial linear region. Careful numerical quasi-static simulations reveal an intimate relation between plastic failure and shear softening. The attaintment of the steady-state value of the shear modulus associated with plastic flow is identified with a percolation of the regions that underwent a plastic event. We present an elementary ``two-state" model that interpolates between failed and virgin regions and provides a simple and effective characterization of the shear softening. △ Less

Submitted 15 December, 2013; originally announced December 2013.

On the Effect of Micro-alloying on the Mechanical Properties of Metallic Glasses

Authors: Oleg Gendelman, J. Ashwin, Pankaj Mishra, Itamar Procaccia, Konrad Samwer

Abstract: "Micro-alloying", referring to the addition of small concentration of a foreign metal to a given metallic glass, was used extensively in recent years to attempt to improve the mechanical properties of the latter. The results are haphazard and nonsystematic. In this paper we provide a microscopic theory of the effect of micro-alloying, exposing the delicate consequences of this procedure and the la… ▽ More "Micro-alloying", referring to the addition of small concentration of a foreign metal to a given metallic glass, was used extensively in recent years to attempt to improve the mechanical properties of the latter. The results are haphazard and nonsystematic. In this paper we provide a microscopic theory of the effect of micro-alloying, exposing the delicate consequences of this procedure and the large parameter space which needs to be controlled. In particular we consider two very similar models which exhibit opposite trends for the change of the shear modulus, and explain the origins of the difference as displayed in the different microscopic structure and properties. △ Less

Submitted 19 September, 2013; originally announced September 2013.

Mechanical control of heat conductivity in microscopic models of dielectrics

Authors: Alexander V. Savin, Oleg V. Gendelman

Abstract: We discuss a possibility to control a heat conductivity in simple one-dimensional models of dielectrics by means of external mechanical loads. To illustrate such possibilities we consider first a well-studied chain with degenerate double-well potential of the interparticle interaction. Contrary to previous studies, we consider varying length of the chain with fixed number of particles. Number of p… ▽ More We discuss a possibility to control a heat conductivity in simple one-dimensional models of dielectrics by means of external mechanical loads. To illustrate such possibilities we consider first a well-studied chain with degenerate double-well potential of the interparticle interaction. Contrary to previous studies, we consider varying length of the chain with fixed number of particles. Number of possible energetically degenerate ground states strongly depends on the overall length of the chain, or, in other terms, on average length of the link between neighboring particles. These degenerate states correspond to mechanical equilibrium, therefore one can say that the transition between them mimics to some extent a process of plastic deformation. We demonstrate that such modification of the chain length can lead to quite profound (almost five-fold) reduction of the heat conduction coefficient. Even more profound effect is revealed for a model with single-well non-convex potential. It is demonstrated that in certain range of constant external forcing this model becomes "effectively"\ double-well, and has a multitude of possible states of equilibrium for the same value of the external load. Thus, the heat conduction coefficient can be reduced by two orders of magnitude. We suggest a mechanical model of a chain with periodic double-well potential, which allows control over the heat conduction. The models considered may be useful for description of heat transport in biological macromolecules and for control of the heat transport in microsystems. △ Less

Submitted 13 September, 2013; originally announced September 2013.

Comments: 10 pages, 16 figures

Shear localization in 3-Dimensional Amorphous Solids

Authors: Ratul Dasgupta, Oleg Gendelman, Pankaj Mishra, Itamar Procaccia, Carmel A. B. Z. Shor

Abstract: In this paper we extend the recent theory of shear-localization in 2-dimensional amorphous solids to 3-D. In 2-D the fundamental instability of shear-localization is related to the appearance of a line of displacement quadrupoles, that makes an angle of 45 degrees with the principal stress axis. In 3-D the fundamental plastic instability is also explained by the formation of a lattice of anisotrop… ▽ More In this paper we extend the recent theory of shear-localization in 2-dimensional amorphous solids to 3-D. In 2-D the fundamental instability of shear-localization is related to the appearance of a line of displacement quadrupoles, that makes an angle of 45 degrees with the principal stress axis. In 3-D the fundamental plastic instability is also explained by the formation of a lattice of anisotropic elastic inclusions. In the case of pure external shear stress, we demonstrate that this is a 2-dimensional triangular lattice of similar elementary events. It is shown that this lattice is arranged on a plane, that, similarly to the 2-D case, makes an angle of 45 degrees with respect to the principal stress axis. This solution is energetically favorable only if the external strain exceeds a yield-strain value, which is determined by the strain parameters of the elementary events and the Poisson ratio. The predictions of the theory are compared to numerical simulations and very good agreement is observed. △ Less

Submitted 10 July, 2013; originally announced July 2013.

The Yield-Strain and Shear-Band Direction in Amorphous Solids Under General Loading

Authors: Ashwin J., Oleg Gendelman, Itamar Procaccia, Carmel Shor

Abstract: It is well known experimentally that well-quenched amorphous solids exhibit a plastic instability in the form of a catastrophic shear localization at a well defined value of the external strain. The instability may develop to a shear-band that in some cases is followed by a fracture. It is also known that the values of the yield-strain (and yield-stress), as well as the direction of the shear band… ▽ More It is well known experimentally that well-quenched amorphous solids exhibit a plastic instability in the form of a catastrophic shear localization at a well defined value of the external strain. The instability may develop to a shear-band that in some cases is followed by a fracture. It is also known that the values of the yield-strain (and yield-stress), as well as the direction of the shear band with respect to the principal stress axis, vary considerably with variations in the external loading conditions. In this paper we present a microscopic theory of these phenomena for 2-dimensional athermal amorphous solids that are strained quasi-statically. We present analytic formulae for the yield-strains for different loading conditions, and well as for the angles of the shear bands. We explain that the external loading conditions determine the eigenvalues of the quadrupolar Eshelby inclusions which model the non-affine displacement field. These inclusions model elementary plastic events and determine both the yield-strain and the direction of the shear-band. We show that the angles of the shear bands with respect to the principal stress axis are limited theoretically between $30^o$ and $60^o$. Available experimental data conform to this prediction. △ Less

Submitted 24 April, 2013; originally announced April 2013.

Journal ref: Phys. Rev. E 88, 022310 (2013)

Atomistic theory of the shear band direction in amorphous solids

Authors: Ashwin J., Oleg Gendelman, Itamar Procaccia, Carmel Shor

Abstract: One of the major theoretical riddles in shear banding instabilities is the angle that the shear band chooses spontaneously with respect to the principal stress axis. Here we employ our recent atomistic theory to compute analytically the angle in terms of the characteristics of the Eshelby inclusion that models faithfully the eigenfunction of the Hessian matrix that goes soft at the plastic instabi… ▽ More One of the major theoretical riddles in shear banding instabilities is the angle that the shear band chooses spontaneously with respect to the principal stress axis. Here we employ our recent atomistic theory to compute analytically the angle in terms of the characteristics of the Eshelby inclusion that models faithfully the eigenfunction of the Hessian matrix that goes soft at the plastic instability. We show that loading protocols that do not conserve volume result in shear bands at angles different from 45$^o$ to the strain axis; only when the external strains preserve volume like in pure shear, the shear bands align precisely at 45$^o$ to the strain axis. We compute an analytic formula for the angle of the shear band in terms of the characteristics of the loading protocol; quantitative agreement with computer simulations is demonstrated. △ Less

Submitted 15 April, 2013; originally announced April 2013.

Comments: 4 pages

Exact solutions for discrete breathers in forced-damped chain

Authors: O. V. Gendelman

Abstract: Exact solutions for symmetric discrete breathers (DBs) are obtained in forced-damped linear chain with on-site vibro-impact constraints. The damping is related to inelastic impacts; the forcing may be chosen from broad class of periodic antisymmetric functions. Global conditions for existence and stability of the DB are established. Some unusual phenomena, like non-monotonous dependence of the sta… ▽ More Exact solutions for symmetric discrete breathers (DBs) are obtained in forced-damped linear chain with on-site vibro-impact constraints. The damping is related to inelastic impacts; the forcing may be chosen from broad class of periodic antisymmetric functions. Global conditions for existence and stability of the DB are established. Some unusual phenomena, like non-monotonous dependence of the stability boundary on the forcing amplitude, are revealed analytically for the full system and illustrated numerically for small periodic lattices. △ Less

Submitted 3 October, 2012; originally announced October 2012.

Comments: 13 pages, 10 figures

Nonstationary heat conduction in one-dimensional models with substrate potential

Authors: O. V. Gendelman, R. Shvartsman, B. Madar, A. V. Savin

Abstract: The paper investigates non-stationary heat conduction in one-dimensional models with substrate potential. In order to establish universal characteristic properties of the process, we explore three different models --- Frenkel-Kontorova (FK), phi4+ ($φ^4$+) and phi4- ($φ^4-$). Direct numeric simulations reveal in all these models a crossover from oscillatory decay of short-wave perturbations of the… ▽ More The paper investigates non-stationary heat conduction in one-dimensional models with substrate potential. In order to establish universal characteristic properties of the process, we explore three different models --- Frenkel-Kontorova (FK), phi4+ ($φ^4$+) and phi4- ($φ^4-$). Direct numeric simulations reveal in all these models a crossover from oscillatory decay of short-wave perturbations of the temperature field to smooth diffusive decay of the long-wave perturbations. Such behavior is inconsistent with parabolic Fourier equation of the heat conduction and clearly demonstrates the necessity of hyperbolic models. The crossover wavelength decreases with increase of average temperature. The decay patterns of the temperature field almost do not depend on the amplitude of the perturbations, so the use of linear evolution equations for temperature field is justified. In all model investigated, the relaxation of thermal perturbations is exponential -- contrary to linear chain, where it follows a power law. However, the most popular lowest-order hyperbolic generalization of the Fourier law, known as Cattaneo-Vernotte (CV) or telegraph equation (TE) is not valid for description of the observed behavior of the models with on-site potential. In part of the models a characteristic relaxation times exhibit peculiar scaling with respect to the temperature perturbation wavelength. Quite surprisingly, such behavior is similar to that of well-known model with divergent heat conduction (Fermi-Pasta-Ulam chain) and rather different from the model with normal heat conduction (chain of rotators). Thus, the data on the non-stationary heat conduction in the systems with on-site potential do not fit commonly accepted concept of universality classes for heat conduction in one-dimensional models. △ Less

Submitted 18 August, 2011; v1 submitted 14 August, 2011; originally announced August 2011.

Comments: 9 pages, 7 figures, corrected version

Non-stationary heat conduction in one-dimensional chains with conserved momentum

Authors: Oleg V. Gendelman, Alexander V. Savin

Abstract: The Letter addresses the relationship between hyperbolic equations of heat conduction and microscopic models of dielectrics. Effects of the non-stationary heat conduction are investigated in two one-dimensional models with conserved momentum: Fermi-Pasta-Ulam (FPU) chain and chain of rotators (CR). These models belong to different universality classes with respect to stationary heat conduction.… ▽ More The Letter addresses the relationship between hyperbolic equations of heat conduction and microscopic models of dielectrics. Effects of the non-stationary heat conduction are investigated in two one-dimensional models with conserved momentum: Fermi-Pasta-Ulam (FPU) chain and chain of rotators (CR). These models belong to different universality classes with respect to stationary heat conduction. Direct numeric simulations reveal in both models a crossover from oscillatory decay of short-wave perturbations of the temperature field to smooth diffusive decay of the long-wave perturbations. Such behavior is inconsistent with parabolic Fourier equation of the heat conduction. The crossover wavelength decreases with increase of average temperature in both models. For the FPU model the lowest order hyperbolic Cattaneo-Vernotte equation for the non-stationary heat conduction is not applicable, since no unique relaxation time can be determined. △ Less

Submitted 16 September, 2009; originally announced September 2009.

Comments: 4 pages, 5 figures

Heat Conduction in One-Dimensional chain of Hard Discs with Substrate Potential

Authors: O. V. Gendelman, A. V. Savin

Abstract: Heat conduction of one-dimensional chain of equivalent rigid particles in the field of external on-site potential is considered. Zero diameters of the particles correspond to exactly integrable case with divergent heat conduction coefficient. By means of simple analytical model it is demonstrated that for any nonzero particle size the integrability is violated and the heat conduction coefficient… ▽ More Heat conduction of one-dimensional chain of equivalent rigid particles in the field of external on-site potential is considered. Zero diameters of the particles correspond to exactly integrable case with divergent heat conduction coefficient. By means of simple analytical model it is demonstrated that for any nonzero particle size the integrability is violated and the heat conduction coefficient converges. The result of the analytical computation is verified by means of numerical simulation in a plausible diapason of parameters and good agreement is observed △ Less

Submitted 29 August, 2003; originally announced August 2003.

Comments: 14 pages, 7 figures

Heat conduction in 1D lattices with on-site potential

Authors: A. V. Savin, O. V. Gendelman

Abstract: The process of heat conduction in one-dimensional lattice with on-site potential is studied by means of numerical simulation. Using discrete Frenkel-Kontorova, $φ$--4 and sinh-Gordon we demonstrate that contrary to previously expressed opinions the sole anharmonicity of the on-site potential is insufficient to ensure the normal heat conductivity in these systems. The character of the heat conduc… ▽ More The process of heat conduction in one-dimensional lattice with on-site potential is studied by means of numerical simulation. Using discrete Frenkel-Kontorova, $φ$--4 and sinh-Gordon we demonstrate that contrary to previously expressed opinions the sole anharmonicity of the on-site potential is insufficient to ensure the normal heat conductivity in these systems. The character of the heat conduction is determined by the spectrum of nonlinear excitations peculiar for every given model and therefore depends on the concrete potential shape and temperature of the lattice. The reason is that the peculiarities of the nonlinear excitations and their interactions prescribe the energy scattering mechanism in each model. For models sin-Gordon and $φ$--4 phonons are scattered at thermalized lattice of topological solitons; for sinh-Gordon and $φ$--4 - models the phonons are scattered at localized high-frequency breathers (in the case of $φ$--4 the scattering mechanism switches with the growth of the temperature). △ Less

Submitted 28 October, 2002; v1 submitted 30 April, 2002; originally announced April 2002.

Comments: 26 pages, 18 figures