the elasto - viscoplastic deformation of solids is a complex process that couples widely separate time and length scales . in polycrystalline solids , plastic deformation
is mainly mediated by dislocations , which form complex and strongly interacting networks that give rise to a rich macroscopic phenomenology . in disordered solids
( e.g. noncrystalline / amorphous solids below their glass temperature and granular materials ) , which lack long - range structural order , plastic deformation is mediated by localized shear - driven rearrangements involving a small number of basic elements ( molecules / atoms / grains ) .
these rearrangements produce inhomogeneous mesoscopic deformation fields that affect their surroundings . in both cases ,
the small scale dynamics are highly heterogenous , intermittent and seemingly non - deterministic , yet the emergent macroscopic phenomena are quite robust . developing predictive continuum models that effectively bridge over the widely separate scales remains a major challenge of prime scientific and technological importance .
the most successful methodology for bridging over widely separated scales is equilibrium statistical thermodynamics .
however , the viscoplastic deformation of solids is an intrinsically out - of - equilibrium phenomenon , hence standard equilibrium thermodynamics is not directly applicable in this case .
nevertheless , statistical concepts and the most general form of the laws of thermodynamics might still be applicable .
our goal in this paper is to incorporate recently developed out - of - equilibrium statistical concepts into a continuum thermomechanical framework for the large elasto - viscoplastic deformation of amorphous solids .
models for the elasto - viscoplastic deformation of amorphous solids have been quite extensively discussed in recent literature , see for example @xcite .
the basic idea that has emerged in the statistical physics community in the last 15 years or so , is that the degrees of freedom of an amorphous / glassy material ( i.e. the various contributions to the internal energy and entropy ) can be approximately decomposed into two qualitatively different groups @xcite .
one group includes the vibrational motions of particles around a mechanically stable configuration ( sometimes termed an inherent structure ) .
these motions equilibrate with the thermal reservoir on a microscopic timescale ; i.e. the temperature apparent in the vibrational fluctuations quickly evolves to match that of the reservoir .
the other group includes the mechanically stable configurations themselves and the slow configurational rearrangements that allow transitions between them .
these are mainly driven by external mechanical forces , possibly assisted by thermal fluctuations .
it is the configurational degrees of freedom that fall out of internal equilibrium with the vibrational degrees when a glass is formed ; and it is these degrees of freedom that undergo structural relaxation when a glass ages and evolve when a glass deforms plastically .
numerous recent theoretical , numerical and experimental studies have provided quantitative evidence that the configurational degrees of freedom of either an aging or a deforming amorphous / glassy material attain a temperature of their own , which is markedly different from the temperature corresponding to the vibrational fluctuations of the microconstituents .
the early works demonstrating the existence of a configurational temperature , usually termed an `` effective temperature '' in the statistical physics community , were rather exclusively based on deviations from the equilibrium fluctuation - dissipation relation ( for a very recent and comprehensive review see @xcite ) .
that is , it was shown ( mainly using simple models and molecular dynamics simulations ) that the long time dynamics of an aging glassy material satisfy a generalized fluctuation - dissipation relation with a temperature that is different from the temperature associated with the vibrational motion of the microconstituents @xcite .
of special importance for our present purposes are the colloidal glass experiments of @xcite , which explicitly demonstrate that the effective temperature of a sample slowly evolves toward the thermal ( i.e. vibrational ) temperature over very long times , and eventually approaches it .
this suggests that the two temperatures interact thermodynamically through an equilibration process akin to fourier s relation to reach a common value .
generalized fluctuation - dissipation relations were also tested for steadily deforming simulated glasses , again demonstrating the existence of a well - defined effective temperature @xcite .
another line of investigation focused on `` athermal systems '' , i.e. systems where no ordinary thermal ( vibrational ) fluctuations take place either due to coupling to a zero temperature bath ( e.g. in computer simulations ) or when the basic elements are too big to experience thermal vibrations .
realistic examples include foams ( made of bubbles ) and granular materials ( made of grains ) where the large size of the basic elements induces dissipation that rapidly quenches vibrations . in @xcite , durian s foam model @xcite
was studied under slow , steady shear deformation conditions .
five different definitions commonly used to define the ordinary thermal ( vibrational ) temperature  pressure and energy fluctuation relations , the green - kubo relation for the shear viscosity , the stokes - einstein relation , and the thermodynamic definition from the rate of change of configurational energy with entropy  were tested and it was discovered that all definitions yielded a similar effective temperature , though the system itself was strictly athermal .
similar ideas were developed and tested for granular materials @xcite .
both @xcite and @xcite showed by explicit calculation that the effective temperature is equivalently given by the derivative of the configurational energy with respect to the entropy of configurational states .
this implies that the configurational subsystem behaves as a distinct thermodynamic entity , featuring a configurational temperature defined as the derivative of the configurational energy with respect to the configurational entropy , precisely as we formulate in detail below .
experimental verification of the effective temperature concept is currently limited to systems in which the basic elements are sufficiently large to be directly tracked , such as grains in granular materials and colloidal particles in colloidal glasses @xcite .
we finally note that the literature on the effective ( configurational ) temperature is immense and can not be exhaustively reviewed here .
the interested reader is referred to @xcite for additional information and relevant references .
two - temperature thermodynamic frameworks for spatially homogeneous , non - plastically deforming , glassy materials have been investigated by several authors @xcite , and extended models coupling the evolution of internal state variables to the new thermodynamic variables have begun to emerge .
applications to amorphous plasticity  recent shear - transformation - zone ( stz ) models @xcite  , to soft glassy rheology @xcite , to dislocation - mediated crystal plasticity @xcite , to glassy memory effects @xcite and to hard - sphere materials @xcite have been proposed .
very recently , this framework has been used to explain a ductile - to - brittle transition in bulk metallic glasses @xcite .    with this background ,
we present the three major goals of this paper :    * to determine the general features of a rational thermomechanical framework which permits , at the outset , the existence of a second primitive temperature .
we view this as a question of basic , fundamental interest . +
* to show that the essential two - temperature principles developed in the statistical physics community can be fit within a consistent thermomechanical framework . +
* to highlight certain non - trivial phenomena captured by the resulting constitutive framework , when choosing forms for various constitutive functions .
we proceed through several steps . in the first step
, we describe the decomposition of the standard continuum laws of energy balance and entropy imbalance ( cf .
@xcite ) into sub - balances for the two subsystems . at this point
, we demonstrate basic thermomechanical consequences of these hypotheses ( including items 1 , 2 and 3 on the list below ) by presenting a simple example of a two - temperature plasticity model . to develop a micro - statistically motivated model we next revert back to our basic hypotheses and , borrowing certain notions from @xcite , introduce a minimal set of internal state variables , and permit the energy and entropy of the configuration to depend on them .
we assume defect - like functional forms for the entropy and energy of the internal variables , and assume constitutive dependences for the other variables consistent with the subsystems decomposition .
restrictions based on the second - law of thermodynamics then lead to a set of `` fourier - type relations '' for heat transfer between the subsystems and through space , which give rise to a pair of heat equations , one for each temperature .
finally , we develop a specific plastic flow rule based on an stz - micromechanism @xcite .
the resulting model is then shown to capture a number of non - trivial phenomena observed in glassy deformation , including    1 .
structural rejuvenation ( i.e. increase in structural disorder ) under shearing .
2 .   a finite configurational temperature in slow , steady shearing .
slow aging of the structural degrees of freedom near the glass transition .
newtonian flow near the glass transition .
emergence of a flow stress for steady flows at low temperatures .
our goal in this section is to mathematically express the basic hypotheses discussed in the introduction and to systematically explore the implications of them for thermodynamics .
this will set the stage for the thermomechanical analysis to be presented later .      by `` nonequilibrium '' processes we refer to thermodynamic processes that do not immediately maximize the entropy . as a motivational example , consider a closed system containing two identical blocks of a highly heat - conductive material , which make contact on either side of a thermally conducting layer .
the layer conducts heat between the blocks ( labeled a and b ) at a rate @xmath0 , where @xmath1 and @xmath2 are the blocks temperatures . if @xmath3 , the system equilibrates on a microscopic timescale , i.e. a temperature difference can not exist for a finite time .
if @xmath4 , we have a completely decoupled nonequilibrium system , for which the system is composed of two separate equilibrium subsystems that never approach a mutual equilibrium temperature ( if @xmath5 initially ) . for finite @xmath6
, we obtain a non - trivial , nonequilibrium system : any temperature difference between the blocks gradually approaches an equilibrium temperature and the joint system gradually approaches a maximum entropy state .    as discussed above
, this latter case is analogous to observations made in amorphous materials . instead of two blocks separated in space
, the material is viewed as having two weakly coupled subsystems appearing within the internal energy and entropy functions .
consider an idealized amorphous system containing randomly and densely packed interacting particles .
focus first on the _ kinetic / vibrational _ ( k / v ) degrees of freedom , which describe random thermal vibrations around a mechanically stable microscopic configuration . in the amorphous setting , these degrees of freedom manifest as the part of the internal energy that arises due to the particles rattling about some mean positions within their cages of nearest neighbors .
this subsystem rapidly equilibrates with an external bath . in the vicinity of the glass transition temperature and below it ,
the timescale of these vibrations is much faster than the slow timescale characterizing structural relaxation and typical rates of deformation .
indeed , the subsystems separation is essentially a timescales separation . by averaging over the short vibrational timescale ,
the energy stored within the microscopic structure itself can be measured , which is classified as energy in the _ configurational _ subsystem . moreover , as the structure is deformed into new configurations , we should also account in this subsystem for the energy due to non - affine motion of the particles ( averaged over the vibrational timescale ) , resulting from the intermittent collective particle - level rearrangements .
under most circumstances , the kinetic energy of this movement comprises a very small contribution to the configurational energy .
the _ internally nonequilibrium _ nature of the material is characterized by the fact that for a given total internal energy , the distribution of energy into configurational and k / v subsystems does not generally maximize the total entropy .
rather , a finite energy exchange rate @xmath7 can exist between the two subsystems , passing energy from the higher to the lower temperature subsystem , which plays the role of @xmath8 in the previous illustrative example . indeed , the internal heat flow @xmath7 accounts for entropy production in the material , together with more conventional sources of entropy production , such as spatial heat fluxes .
the finiteness of @xmath7 in the amorphous setting can be distinguished from the example cases of weakly sheared newtonian fluids or elastically deformed crystals . in a
weakly sheared fluid , thermal vibrations and stress - relaxing hopping motions occur on timescales that are not many orders of magnitude separated , and all of the degrees of freedom can be characterized by a single temperature . hence , equipartition is satisfied , i.e. the total configurational energy remains in equilibrium with the kinetic degrees of freedom , analogous to @xmath9 . in an elastically deformed crystal , both thermal vibrations and the external driving forces
are , by the definition of the elastic regime , too weak to induce structural changes ( e.g. dislocations motion ) , and hence the configurational degrees of freedom are essentially decoupled from the vibrational ones , analogous to @xmath10 .
nevertheless , the microscopic kinetic and potential energies due to the crystal vibrations do maintain equilibrium following the equipartition theorem , which is completely analogous to our caricature of the k / v subsystem as single `` block '' in the thermodynamic analysis .
in the case of shearing an amorphous solid , as was discussed in detail above ( see , for example , @xcite ) , there exists ample evidence that the configurational degrees of freedom also form a single equilibrated `` block '' , which is distinguished from the k / v `` block '' , and hence two distinct temperatures may be sufficient to approximately describe the internal thermodynamics of the system . this physical picture
is schematically sketched in fig .
[ fig1 ] .
let us denote by @xmath11 the temperature of the configurational subsystem and by @xmath12 the temperature of the kinetic / vibrational subsystem , which will both be formalized in detail in the following sections .
the configurational temperature @xmath11 is a generalization of the relatively old notion of a fictive temperature @xcite .
it is typically _ higher _ than the thermal temperature @xmath12 , because the structure typically retains memory of configurations characteristic of the higher temperature prior to a rapid quench and is also constantly being rejuvenated by the external driving forces .
following our two - block analogy , one would expect a fourier - type relation of the form @xmath13 .
we note in passing that thermo - elastic effects ( e.g. thermal expansion , temperature - sensitive elastic moduli ) provide another form of interaction between the two subsystems , which is not mediated by the heat flux @xmath7 .
_ we neglect these phenomena in the present work _ to focus more heavily on the deviatoric ( non - volumetric ) plastic response
.    an important subtlety of the ensuing analysis pertains to the functional dependence of the conductivity @xmath14 . in the absence of plastic deformation
, only ordinary thermal fluctuations can pull the system out of a local configuration and induce flow from the configurational subsystem to the kinetic / vibrational one .
as the glass transition is approached , thermal fluctuations become increasingly ineffective in spontaneously inducing structural relaxation , and consequently @xmath14 becomes extremely small in the absence of imposed shearing ( it strictly vanishes if there exists an ideal glass transition ) . in the presence of plastic deformation
, rearrangements are activated and energy can again flow between the two subsystems , i.e. @xmath15 , at a rate that is affected by the rate of plastic deformation itself . moreover , the rearrangements themselves can mediate diffusion of @xmath11 , when the latter is spatially inhomogeneous .
this occurs microscopically through the long range eshelby - like fields that are generated by the localized shear rearrangements ( that can be viewed as eshelby - like inhomogeneities @xcite ) .
the direct dependence of the inter - subsystem conductivity and intra - configurational subsystem diffusivity on the rate of plastic deformation is a unique feature of these systems and hence of the corresponding material models .
we now turn to formulate these ideas mathematically .
consider a deformation process for which @xmath16 is the reference coordinate , @xmath17 is the deformed coordinate and @xmath18 , where @xmath19 is the motion function . the spatial velocity field is then defined by @xmath20 , where @xmath21 denotes the lagrangian time derivative .
the local first law of thermodynamics reads @xmath22 where @xmath23 is the mass density , @xmath24 is the internal energy per unit mass , @xmath25 is the cauchy stress tensor , @xmath26 is the spatial velocity gradient ( i.e @xmath27 ) , @xmath28 is the spatial heat flux and @xmath29 is the spatial divergence operator .
in the above we have assumed no internal heat sources or sinks , which would otherwise add a scalar heating term to the right hand side .
the second law of thermodynamics reads @xmath30 where @xmath31 is the internal entropy per unit mass and @xmath32 the entropy flux .    in the spirit of the previous discussion
, we now apply the assumption that the thermodynamical quantities can be approximately decomposed into their kinetic / vibrational and configurational contributions , denoted hereafter by the subscripts @xmath33 and @xmath34 respectively .
first , we separate the internal energy and entropy into configurational and kinetic / vibrational contributions @xmath35 next , the heat flux is assumed to approximately decompose as @xmath36 where @xmath37 and @xmath38 are the fluxes of configurational energy ( @xmath39 ) and k / v energy ( @xmath40 , respectively . to gain physical insight into this decomposition , note that the heat flux @xmath28 into ( out of ) a part of a body represents the flux of energy due to the flow of randomized movement of the microscopic constituents  movements that collectively have zero mean velocity relative to the continuum - scale motion . in our context
, we can envision two different manifestations of the passage of randomized movement .
the kinetic heat flux @xmath38 may be seen as the `` standard '' flux describing the passage of random vibrational motion into ( out of ) a part .
the configurational contribution @xmath37 can be understood as the energy flux due to rearrangements of the structure , which pass energy across space by disturbing the neighboring structural assembly in a random fashion , and which does not count any k / v vibrational mechanisms in so doing .
the entropy flux is similarly decomposed into two additive contributions @xmath41 in the classical clausius - duhem inequality , the temperature @xmath42 is positted as a scalar , non - negative field , which relates the heat and entropy fluxes by @xmath43 .
we introduce the notion of a second primitive temperature  the most noteworthy departure from common thermomechanics in this paper  by assuming analogously that each subsystem contribution to the entropy flux is proportional to its corresponding heat flux via an as - yet - undefined temperature field , giving @xmath44 and @xmath45 . as a result @xmath46
the case of @xmath47 corresponds to internal nonequilibrium , whereas @xmath48 in internal equilibrium and we recover the standard entropy flux @xmath43 .
we also assume each temperature is necessarily non - negative .
we should also mention that our choice of @xmath44 and @xmath45 is somewhat analogous to multi - temperature thermomechanical models for immiscible mixtures @xcite ; instead of @xmath34 and @xmath33 subsystems , each phase of the mixture is permitted its own temperature and the relation between heat and entropy fluxes within each phase , and their additive appearance in the second law , appears in direct analogy to ours .    adopting the conventional krner - lee decomposition of the total deformation gradient tensor @xmath49 @xcite , @xmath50 we can write the stress - power as @xmath51 with @xmath52 , @xmath53 , @xmath54 , @xmath55 , @xmath56 and @xmath57 @xcite . we note that @xmath58 maps the reference ( undeformed ) space into the so - called intermediate space and that @xmath59 maps the intermediate space into the deformed space .
also note that @xmath60 and hence @xmath61 .    in what follows
we assume that plastic shear deformation involves no volumetric changes , i.e. we set @xmath62 . in doing so
we neither imply that there exists no coupling between plastic shear deformation and dilation / contraction nor that this coupling is of no importance .
such coupling is well documented for amorphous solids and might have consequences for their mechanics @xcite .
we have developed a theory that includes pressure - sensitivity and shear dilation ; however , the added complexities that arise suit reporting in a subsequent publication .
to proceed , we need separate energy balance equations for the two subsystems , which sum together to give the first law of thermodynamics in eq .
( [ 1st ] ) .
we propose that these take the following form @xmath63    let us discuss in detail the physical meaning of the different terms on the right - hand - sides of these equations .
the spatial contributions to the heat flow _ within _ each subsystem are @xmath64 and @xmath65 for the kinetic and configurational subsystems respectively .
@xmath7 is the scalar internal heat flow that transfers energy _ between _ the configurational subsystem and the kinetic one .
the rest of the terms represent the way the total stress - power @xmath66 is distributed among the subsystems .
first , we note that elasticity resides in the configurational subsystem and hence reversible elastic energy appears in the @xmath34 subsystem energy balance ; recall that we neglect thermo - elastic and volumetric effects , which could conceivably admit an elastic term in the @xmath33 subsystem .
the remaining terms give the plastic stress - power @xmath67 .
to remain general , we permit the plastic stress - power to affect both subsystems , by splitting @xmath68 into the sum of two contributions @xmath69 .
the _ stored plastic power _ describes the rate of stored energy of cold work , and thus contributes in the @xmath34 energy balance .
plastic dissipation _ is the rate of energy dissipated as ordinary heat , which is in the form of vibrations and thereby affects the @xmath33 subsystem balance .
applying the chain rule to eq .
( [ 2nd ] ) and substituting eqs .
( [ firstk])-([firstc ] ) into the result to eliminate @xmath70 and @xmath71 , we obtain @xmath72 where @xmath73 is the spatial gradient operator . letting @xmath74
be the mass density in the intermediate space , we redefine quantities relative to volumes , areas , and lengths in this space as @xmath75 multiplying eq . ( [ int_space_1b ] ) by @xmath76 and substituting the above quantities , we obtain an intermediate frame entropy imbalance ( for ease of notation , the @xmath77 are omitted hereafter ) @xmath78 similarly , the energy balances become @xmath79 where @xmath80 for any vector field @xmath81 .
as was stated in the abstract , internal variables that carry information about the structural state of the deforming solid , and determine its subsequent response , play a central role in a theory of elasto - plasticity .
yet , we present here first a simple two - temperature model for rate - dependent elasto - plasticity , which does not include internal variables .
our goal is to illustrate general features of the two - temperature thermomechanical framework .
we start by defining the configurational and kinetic helmholtz free - energy per unit intermediate volume @xmath82 let @xmath83 be the plastic rate - of - deformation tensor , @xmath84 the scalar plastic flow - rate , and @xmath85 the direction of plastic flow .
let us assume plastic incompressibility , @xmath86 .
we write the following functional dependences ( note that here , and throughout the paper , we use the  @xmath87  symbol to denote constitutive functions of designated independent variables ) @xmath88 we use @xmath89 to denote the deviator , and note that the forms for @xmath90 and @xmath91 assume @xmath92
. these constitutive dependences are not unique .
rather , they are informed by past experience with constitutive theories of thermal plasticity , and are properly extended for consistency with the two - subsystem decomposition .
in addition , the @xmath93 and @xmath94 dependences are also motivated by the statistical - thermodynamics - based model to be described in the upcoming sections .
the above dependences are frame - indifferent ( see appendix [ frame_indifference ] ) , and unlike standard heat transfer theories , we allow all heat flow variables associated with the configurational subsystem to explicitly depend on the plastic rate of deformation . the latter is a key point , as plastic deformation directly influences the time - scale for configurational energy transfer .
finally , in the spirit of the subsystems separation , the above constitutive dependences aim at coupling configurational / kinetic dependent variables with configurational / kinetic independent variables .    substituting the free - energies into the inequality in ( [ int_space_2 ] ) and the expanding time - derivatives , we obtain @xmath95
as the first three terms are independent of @xmath96 , @xmath97 and @xmath98 , they must vanish in order to ensure non - violation of the inequality @xcite .
hence , @xmath99 which represent an elasticity relation and helmholtz relations for each subsystem , respectively .
since @xmath100 is assumed independent of @xmath101 , the configurational helmholtz relation implies that @xmath102 .    to ensure that the four last terms in ( [ simple_ineq ] ) do not violate the inequality , we assume  as a sufficient , but not necessary condition
that each is non - negative separately . for the last three terms ,
non - negativity is assured by the following set of fourier - type relations @xmath103 the term proportional to @xmath104 can be made non - negative under various selections of the functions for @xmath90 and @xmath91 . to obtain a simple , isotropic , incompressible bingham - like rheology
, we choose @xmath105 which is valid only for @xmath92 .
together , the plastic stress obeys @xmath106 , which makes it evident that @xmath107 is an effective plastic flow viscosity and @xmath108 acts as a yield strength .    finally , we insert these results into eqs .
( [ firstki ] ) and ( [ firstci ] ) to obtain heat equations for the two subsystems .
applying the chain rule and the definition of the free - energies , we can write @xmath109 and @xmath110 , to obtain @xmath111 \ , \label{conventional_heat_k } \\
c_c \dot\theta_c = & \ y^p_{store } d^p -{\kappa}_{ck}\ ( \theta_c-\theta_k ) + \text{div}_i \left[{\kappa}_c\ ,   { { \bm{g}}}_c \right ] \ .
\label{conventional_heat_c}\end{aligned}\ ] ]    to appreciate the implications of the two - temperature framework in general , and of the configurational heat equation in particular , we need to be more specific about how @xmath11 affects the flow rule .
that is , we need to specify the @xmath11-dependence of @xmath112 . in the subsequent sections , where we will discuss an internal variables model within a statistical thermodynamic framework
, we will show that @xmath104 is proportional to a boltzmann - like factor with @xmath11 appearing as the relevant temperature .
for the present discussion this implies that @xmath113 where @xmath114 is boltzmann s constant and @xmath115 is a constant whose physical meaning will be explained later on .
the most important thing for our purposes here is that eq .
( [ mu_p ] ) implies that the plastic rate of deformation @xmath104 in eq .
( [ conventional_flow ] ) is very ( i.e. exponentially ) sensitive to variations in @xmath11 .
we therefore focus here on the behavior of @xmath11 , according to eq .
( [ conventional_heat_c ] ) , under various thermomechanical conditions :    * _ thermodynamic aging _  let us consider a situation in which a homogeneous sample is held at a temperature @xmath12 near its glass temperature with no imposed deformation , @xmath116 .
the evolution of the system under these conditions corresponds to `` structural relaxation '' or `` aging '' .
structural relaxation is a central glassy phenomenon in which physical properties ( e.g. volume , enthalpy , elastic moduli etc . ) vary spontaneously ( i.e. in the absence of external driving forces ) and slowly with no characteristic timescale @xcite .
in many cases , the time evolution is logarithmic .
aging is one of the most basic manifestations of the nonequilibrium nature of the glassy state . under the stated conditions eq .
( [ conventional_heat_c ] ) predicts that @xmath117 hence , our model predicts that @xmath11 evolves from its initial value to @xmath12 ; the two subsystems gradually equilibrate , and the @xmath11 relaxation dynamics correspond to structural aging in our nonequilibrium thermodynamic framework ( typically the initial condition satisfies @xmath118 and consequently @xmath11 decreases during the aging process ) .
since @xmath119 typically varies slowly and nonlinearly with @xmath11 , the relaxation of @xmath11 to equilibrium at @xmath12 will be slow and strongly non - exponential .
+ the thermodynamic picture of aging is partially supported by various generalized fluctuation - dissipation analyses of aging systems @xcite . as was mentioned in the introduction
, the recent colloidal glass experiments of @xcite explicitly demonstrated that @xmath11  estimated through a generalized fluctuation - dissipation relation  equilibrates with @xmath12 for sufficiently long aging times ( cf .
3b there ) . while these observations provide strong support for the existence and importance of @xmath11 in aging dynamics , it is important to note , however , that even if the function representing @xmath120 is chosen so as to quantitatively match experimental aging data , this single - variable description of aging remains incomplete , cf .
_ low - temperature homogeneous steady state deformation _  consider then homogeneous steady state deformation at a low temperature @xmath12 . under these conditions
we can set @xmath121 in eq .
( [ conventional_heat_c ] ) ( due to the steady state conditions ) , neglect @xmath12 as compared to @xmath11 ( low @xmath12 ) and omit the spatial term ( homogeneous deformation ) . with these we can isolate @xmath11 in eq .
( [ conventional_heat_c ] ) , obtaining @xmath122 this result shows that the configurational temperature @xmath11 attains a finite value ( obtained by solving the above equation if @xmath123 and @xmath14 depend on @xmath11 ) due to a competition between a `` rejuvenation '' term @xmath124 ( acting to increase @xmath11 ) and a `` relaxation '' term @xmath125 ( acting to decrease @xmath11 ) . as discussed above , at the low temperatures @xmath12 of interest here
, it is sensible to assume that @xmath126 .
this implies that the steady state value of @xmath11 in eq .
( [ conventional_thetac ] ) is _ independent _ of the deformation rate ( if @xmath123 is independent of @xmath127 ) , as is indeed widely observed ( see , for example , @xcite ) .
* _ configurational temperature diffusion _
 the configurational heat diffusion term @xmath128 $ ] in eq .
( [ conventional_heat_c ] ) is expected to play a role in the presence of inhomogeneous flows and localization . in particular , @xmath129 is inevitably related to a _ length - scale _ and is expected to give rise to , and potentially explain , shear - band diffusive widening . in this regard , it is possible that @xmath11 is related to the so - called `` fluidity '' state variable recently invoked in nonlocal flow theories @xcite .
like @xmath11 , the fluidity is hypothesized to obey diffusional dynamics , which enable it to describe shear band widths in various amorphous media .
+ in more general terms , as the configurational heat conductivity @xmath129 introduces a length - scale , the two - temperature framework might provide a new mechanism for size effects in constitutive models .
unlike other size - dependent plasticity models where scale - sensitive plastic flow properties are obtained by introducing gradient - based pairings in the internal / external power statements @xcite , the nonlocal behavior here emerges from the splitting of the usual first law into two pieces .
the diffusion of the new temperature @xmath11 then affects the flow through @xmath112 .
this coupling gives rise to a physically nonlocal constitutive mechanism , on time / space scales distinct from the usual ( i.e. k / v ) heat diffusion .
this nonlocality bears a qualitative similarity to the implicit gradient model of @xcite , where an additional diffusive variable ( the `` fictitious plastic strain '' ) is introduced , and nonlocal plastic flow behaviors are induced by dependence of the flow rule on this variable .
it is our goal now to go beyond the illustrative simple model discussed in the previous section and to derive a new constitutive model that integrates coarse - grained internal variables within the proposed framework . in
so doing we hope to develop a more basic statistical thermodynamic approach , leading to models with improved predictive powers .
internal variables implicitly carry , through their evolution , information about the history of deformation , and , as descriptors of the material state , directly influence the deformation behavior @xcite . the evolution of these variables ,
we posit , is driven by the exchange of energy and entropy with other parts of the deforming system and hence the internal variables  which are naturally associated with the configurational subsystem  must be explicitly associated with energy and entropy of their own @xcite .
this is a rather significant deviation from conventional approaches , which traditionally do not associate entropy with internal variables @xcite .    in an ideal scenario
, one should be able to define macroscopic internal variables by systematically coarse - graining over microscopic dynamics . in practice , this is a daunting task and one usually resorts to various phenomenological descriptions , which are sometimes rationalized based on physical intuition , experimental observations , symmetries , computer simulations , etc .
an alternative approach would be to first identify mesoscopic objects that are abstracted from the microscopic physics and only then to coarse - grain to obtain the internal variables and their evolution , as in shear - transformation - zone ( stz)-based theories @xcite to be discussed later .
before we specify such a microscopic mechanism ( see section [ stz_sec ] ) , our goal in this section is to develop a broad theoretical framework that permits a minimal set of internal variables , to which a more specific micro - mechanism may later be attributed .    in amorphous solids , different initial states of disorder would result from different cooling rates through which the glass transition was approached .
the different initial structural states may determine whether the stress - strain curve exhibits a stress peak or not , and whether the material shear localizes or not ( see the review papers of @xcite ) .
in addition , we know that solids ( both crystalline and amorphous ) carry orientational memory of previous plastic deformation , which affects subsequent deformation
. this is demonstrated , for example , in the classical bauschinger effect in metals and glassy polymers , where specimens undergo a decrease in flow strength in compression after experiencing a certain amount of plastic extension . to address these observations in a minimal sense , in what follows
we consider one scalar and one ( symmetric , deviatoric ) tensorial internal variable , denoted respectively by @xmath130 and @xmath131 .
these variable types are common to isotropic and kinematic hardening models @xcite . in our setting
, we assume that the internal variables @xmath130 and @xmath131 arise from dilute , weakly interacting flow defects ( whose precise nature is not specified yet ) that mediate plastic deformation . the internal variable @xmath130 stands for their density ( number per unit intermediate volume ) and the deviatoric tensor @xmath131 describes their ( dimensionless ) average orientation in the intermediate space .
defining @xmath132 , we require @xmath133 where @xmath134 represents an isotropic orientational state and @xmath135 a completely anisotropic state where all defects align in the same direction .
we use the set @xmath136 , its time derivatives , and the set @xmath137 as the independent variables , and assume the following constitutive dependences @xmath138 we stress that we use entropies and internal energies , rather than free - energies , as the basic quantities because we insist on directly associating entropy and energy with the internal variables . in our view
, this constitutes an important deviation from conventional modelling .
the above constitutive choices reflect the separation between subsystems by associating each entropy only to the variables naturally arising in that subsystem , and are also frame - indifferent ( see appendix [ frame_indifference ] ) .
we now follow , nearly verbatim , the steps described in section [ sec : conventional ] . expanding time - derivatives in eq .
( [ int_space_2 ] ) , we obtain @xmath139    in spite of the fact that we do not yet specify explicit constitutive dependences for the other dependent variables in ( [ second2 ] ) , let us assume they are all independent of @xmath140 . then , since the three first terms in the above inequality appear linear in the independent variables @xmath141 , @xmath142 , and @xmath143 respectively , the corresponding coefficients must vanish to ensure non - violation of the inequality @xcite .
consequently , @xmath144 which are the counterparts of eq .
( [ conventional_maxwell ] ) , just expressed in terms of entropies and internal energies , rather than free - energies .
now we go beyond the simple model of section [ sec : conventional ] , and similar to @xcite we express the entropy as @xmath145 here , @xmath146 represents the part of the configurational energy corresponding to the energies of formation of the flow defects and @xmath147 is the energy from bulk elasticity . hence , @xmath148 and @xmath149 are the energy and entropy of the large background structure _ surrounding _ the defects , excluding the elastic energy .
the @xmath150 degrees of freedom satisfy @xmath151 , and thus equivalently behave as the `` bath '' for @xmath152 .
also , under persistent driving , @xmath152 is expected to equilibrate with @xmath150 relatively quickly compared to the ( extremely ) slow equilibration between the @xmath33 and @xmath34 subsystems .
we choose @xmath153 such that eqs .
( [ good_cn])-([ini_eta_split ] ) give rise to the following elasticity relation @xmath154 consequently , @xmath155 is _ symmetric _ , and @xmath156 and @xmath157 are assumed symmetric as well .
we propose that the energy and entropy of the flow defects take the following forms @xmath158 \ , \ ] ] where @xmath115 is the formation energy of an isolated defect and @xmath159 is the characteristic volume of a defect .
the entropy density of defects , @xmath160 , can be obtained in the dilute limit by counting the number of ways defects can be placed on available sites and then given an orientation .
the dimensionless entropic factor @xmath161 can be modeled as follows .
clearly @xmath162 as there is only one way to align defects to obtain @xmath135 , modulo rigid rotation .
it should also approach a maximum value at @xmath134 , the isotropic state , with @xmath163 at @xmath134 due to symmetry .
we therefore assume @xmath164 where @xmath165 is the orientational entropy per flow defect when @xmath134 .
we can now write @xmath166-\frac{e_z}{\theta_c}\ , \quad\quad \frac{{\partial}\hat\eta_c}{{\partial}{{\bm{m } } } } = -2 k_b \lambda s_0 { { { \bm{m}}}}\ .\ ] ] using these with eqs .
( [ good_cn ] ) in the second law ( [ second2 ] ) yields the following reduced second - law inequality @xmath167\right){\dot{\lambda}}\ge 0 \ .\nonumber\end{aligned}\ ] ]    let us now switch to using @xmath168 and @xmath12 , instead of @xmath169 and @xmath170 , as independent variables for the remaining quantities and assume the following constitutive dependences @xmath171 which should be compared to eqs .
( [ conventional_cons2])-([conventional_cons3 ] ) .    for consistency with the new independent variables , we invert eqs .
( [ ini_const_deps ] ) and ( [ ini_eta_split ] ) to obtain @xmath172 the above gives important perspective on how the different variables influence different contributions to the total internal energy .
the configurational energy decomposes into the energy of the structure outside the defects , which depends on the state of structural disorder through @xmath11 , the defect energy , governed by @xmath130 and @xmath131 , and the reversible elastic energy , governed by @xmath173 .
our strategy to prevent a violation of the second law inequality in ( [ second3 ] ) is the same as the one discussed in section [ sec : conventional ] .
we constrain the last four terms to be non - negative under all circumstances and guarantee non - negativity of the other terms together ; it is a sufficient strategy for thermodynamic compatibility though not a unique one .
non - negativity of the second - row terms is guaranteed by adopting the `` fourier - type relations '' in eqs .
( [ fourier_conventional ] ) , together with @xmath174\ , & \ \ \ \ \kappa_{cd}=\hat{\kappa}_{cd}(\theta_c,\theta_k , d^p ) \ge 0 , \quad\quad \lambda_0=\hat{\lambda}_0(m ) = \frac{e^{(1-m^2 ) s_0}}{v_0 } \ .\label{kcd}\end{aligned}\ ] ] note that at fixed @xmath175 , the condition for @xmath176 is an equilibrium number of defects _ determined by the configurational temperature _ through a boltzmann - like factor @xmath177 .
this statistical thermodynamics result , alluded to in section [ sec : conventional ] , goes beyond conventional modelling and is one of the most important consequences of the proposed framework .
the various modes of energy transfer in this model are visualized in fig .
[ fig2]a .
\a )   b )
our goal in this section is to complete the thermodynamic formulation by deriving dissipation inequalities and heat equations , similarly to the procedure discussed in section [ sec : conventional ] . to further exploit the second law inequality to derive additional constraints on the constitutive law , i.e. dissipation inequalities
, we note that since all of the terms in the second row of the inequality in ( [ second3 ] ) were already made non - negative above , it is sufficient to satisfy the following reduced inequality @xmath178 where we have replaced @xmath179 with its symmetric part @xmath180 because it is contracted with symmetric tensors . to proceed , we need to discuss the evolution equation for @xmath131 , which represents the average orientation of the flow defects population .
we propose a minimal evolution rule of the form @xmath181 where @xmath182 and @xmath183 are positive scalar functions .
this equation is shown in appendix [ stz_app ] to emerge on its own from a mesoscopic model of stz dynamics . essentially , the first term produces orientational memory in accordance with the direction and rate of flow , and the latter is a relaxation term that pushes the system to lose orientational memory over time . using the last equation with eq .
( [ etacderivs ] ) , we obtain @xmath184 . substituting this expression into the inequality in ( [ second4 ] ) yields @xmath185
the first term is always non - negative so it is sufficient to ensure the remaining term is non - negative @xmath186 in the next section we will choose functional forms consistent with this inequality .
the last step in the thermodynamic development is the derivation of heat equations using the first law relations in eqs .
( [ firstki])-([firstci ] ) .
we focus first on eq .
( [ firstki ] ) . recalling that @xmath187 and taking into account the fourier - type relations in eq .
( [ fourier_conventional ] ) , we obtain @xmath188 \ , \ ] ] with @xmath189 .
this equation is essentially the same as eq .
( [ conventional_heat_k ] ) , discussed in the context of the simple model in section [ sec : conventional ] .
hats are included to remind us which terms have been assigned functional dependences and which are independent variables ( in the current stage of the analysis ) .
this is an ordinary heat equation for @xmath12 in which the plastic dissipation @xmath190 and the heat flow between the kinetic / vibrational and configurational subsystems appear as source terms .
more interestingly , for the configurational subsystem we apply the chain - rule to express @xmath191 in terms of @xmath98 and simplify the result using eq .
( [ firstci ] ) to obtain @xmath192 \ , \label{heatc}\end{aligned}\ ] ] with @xmath193 .
this is an evolution equation for the temperature of the configurational degrees of freedom , @xmath11 , and hence describes the evolution of structural disorder .
this is a central result within the present approach and is closely related to eq .
( [ conventional_heat_c ] ) , discussed in the context of the simple model in section [ sec : conventional ] .
this section concludes the derivation of the thermodynamic framework .
the constitutive model developed so far , along with the relevant balance laws , takes the form : + conservation of linear and angular momentum : @xmath194 .
kinematic decomposition and definition of intermediate space variables : @xmath195 , @xmath196 .
elasticity relation :  @xmath197 .
decomposition of the plastic stress : @xmath198 .
evolution law for the defect concentration internal variable : @xmath199 $ ] with @xmath200 .
kinetic / vibrational heat equation : @xmath201 $ ] .
configurational heat equation : @xmath202 $ ] .
+ the not yet determined functions for the conduction coefficients @xmath14 , @xmath203 , @xmath129 , @xmath204 must all be non - negative . the functions for @xmath180 and @xmath156 , and the coefficients in the @xmath205 relation , have to be derived from physical considerations and shown to satisfy the constraint of eq .
( [ plast_dissip ] ) .
while our focus in the present paper is on the thermomechanical framework , and by no means on specific constitutive models , we aim here at discussing a more elaborated internal variables model that goes beyond the simple example discussed in section [ sec : conventional ] .
below we present a constitutive model based on an stz micro - mechanism , which rather naturally fits in the thermomechanical framework .
the resulting complete set of equations will allow us in the next section to present a few illustrative applications of the model ( full details of the derivation can be found in appendix [ stz_app ] ) .    based on extensive observations
, we state that irreversible flow in amorphous solids is mediated by a population of flow defects in the form of oriented stz s
sparsely distributed localized zones that are more susceptible to shear - driven transformations in a particular direction than their surrounding material ( see fig .
[ fig2]b ) .
we view stz s not as permanent flow defects , but rather as configurational fluctuations that can appear and disappear during flow and/or in response to thermal fluctuations .
we adopt the main physical ingredients of the stz theory , which has been built on the flow - defect theories of turnbull , cohen , argon , spaepen and others @xcite , and incorporate them into the presently developed framework .
we define the microscopic parameters : @xmath159 is the volume of an stz ( consistent with its prior usage ) , @xmath115 is the characteristic energy required to produce an stz , @xmath206 is the energy required to shear - transition an stz along its weak direction , @xmath207 is the half - amplitude of the shear strain of an stz transition , and @xmath208 is an attempt frequency for stz transitions .
one version of the two - dimensional stz model leads the following equations ( see appendix [ stz_app ] ) : + evolution of the orientation internal variable : @xmath209 which is consistent with eq .
( [ assumed_dotm ] ) and allows a unique determination of @xmath182 and @xmath183 there . + flow rule : @xmath210 these two constitutive equations are then substituted into the second law inequality of ( [ plast_dissip ] ) , which can be satisfied by choosing the following form for @xmath156
. + stored plastic stress : @xmath211
our goal here is to discuss the properties of the internal variables constitutive model described above , focussing on features that go beyond those already discussed in the context of the simple model in section [ sec : conventional ] .
in particular , the discussion of thermodynamic aging and configurational temperature diffusion applies here as is , since these are generic features of the proposed framework .
the existence of a steady state configurational temperature is also generic , though we discuss it here again in order to gain a somewhat better physical understanding of its emergence .
near the glass transition temperature , relaxation times are large ( macroscopic ) , but still fall within the experimental measurement window and phenomena such as newtonian flow and aging can be probed @xcite .
our goal here is to explore the emergence of such phenomena within the proposed constitutive law .
we first focus on newtonian flows and hence consider spatially homogeneous , pure shear situations .
we reduce the analysis to two dimensions of space and assume the deformation is driven by a deviatoric stress of the form @xmath212 , where @xmath213 is the unit - magnitude deviatoric tensor with major principal direction @xmath214 .    then we can write @xmath215 , @xmath216 , and @xmath217 for `` @xmath218 '' indicating any of our subscripts . the viscoplastic relation in eqs
. ( [ mevolve1])-([d_1st1 ] ) gives @xmath219,\quad d^p \!=\ !
2\,v_0\gamma_s\,\lambda_s\ , \nu\ , e^{-\delta / k_b\theta_k}\left[\frac{t^p_{dis}}{\tau_0}-m\right ] \ , \end{aligned}\ ] ] where we have assumed eq .
( [ kcd ] ) remains at its fixed - point @xmath220 .
also , recall that @xmath221 .    a crucial step in analyzing these equations
is specifying the thermodynamic coupling coefficient @xmath222 that is responsible for the creation and annihilation of stz s .
as was discussed in detail above , this ( and other ) thermodynamic coupling coefficient depends explicitly on the plastic rate of deformation @xmath127 .
moreover , in the absence of plastic deformation , @xmath116 , @xmath203 may be finite due to thermal fluctuations alone .
therefore , we set @xmath223 where the dimensionless function @xmath224 remains an unspecified function of the configurational temperature @xmath11 .
it is reasonable that the thermal attempt frequency @xmath208 sets the time - scale for the @xmath12-dependent term , and the flow - rate plays the analogous role for the @xmath11-dependent term .
the physical meaning of the dimensionless function @xmath225 will be discussed soon .
we are now in a position to study eqs .
( [ simple_shear ] ) in various situations .
first , consider the steady state linear response  i.e. small variations about @xmath226  to a small applied stress @xmath227 .
the newtonian ( extensional ) viscosity is defined by @xmath228 .    to calculate @xmath229 , we first substitute eq .
( [ k_cd ] ) in eq .
( [ simple_shear ] ) to obtain @xmath230 \ .\ ] ] as we are interested in the linear response regime , the second - order term proportional to @xmath231 can be omitted .
moreover , as we focus on steady state we can set @xmath232 and use @xmath233 of eq .
( [ simple_shear ] ) to obtain @xmath234 in the vicinity of the glass transition temperature and just below it , we expect spontaneous thermally - activated creation and annihilation of stz s to be strongly suppressed , i.e. @xmath225 becomes small . therefore , @xmath235 is a small parameter and to leading order in the smallness eq .
( [ m_eta ] ) reads @xmath236    by comparing the left side to the expression for @xmath233 in eq .
( [ simple_shear ] ) , we immediately obtain @xmath237 , which tells us that the material flows persistently only because @xmath225 is finite .
the physics here is clear : when @xmath225 vanishes , stz s can not be created and annihilated by thermal fluctuations ( recall that mechanically - induced stz creation and annihilation do not contribute at linear order ) and hence existing ( properly oriented ) stz s can undergo transitions in the direction of the applied stress , but once they all do so , plastic flow can not persist . to complete the calculation of @xmath229 we use eq .
( [ tp_store_result1 ] ) for @xmath238 , which together with eq .
( [ m_eta ] ) implies @xmath239 and hence @xmath240 .
therefore , we obtain @xmath241 where we used @xmath242 in the linear response regime and the definition of @xmath243 .
we observe that @xmath244 , which shows that the newtonian viscosity increases as @xmath225 decreases with decreasing temperature @xmath12 and strictly diverges if @xmath225 vanishes . indeed
, the essence of the glass transition puzzle , which is commonly probed by measuring the dramatic increase of @xmath229 with decreasing temperature @xmath12 , is encapsulated in @xmath225 , which at present can not be calculated from first principles . in the liquid state ( at high temperatures ) @xmath245 , while @xmath225 decreases dramatically in the supercooled and glassy regime ( for `` fragile '' glasses @xcite ) .
we also note that @xmath229 strictly diverges only if there exists an ideal glass transition , @xmath246 , and that due to the appearance of the arrhenius factor @xmath247 in eq .
( [ stz_eta ] ) , @xmath225 can be identified as the non - arrhenius temperature dependence of the newtonian viscosity , similarly to @xcite .
we would like now to explore the nonlinear response in the framework of the proposed model , and in particular the emergence of a yield stress .
using eqs .
( [ tp_store_result1])-([simple_shear ] ) , we obtain @xmath248 m \equiv t^p - \bar\tau m \ .\ ] ] we focus on low bath temperatures @xmath12 for which @xmath225 is extremely small and can be neglected . in this case , eq .
( [ k_cd ] ) can be approximated as @xmath249 and the first equation in ( [ simple_shear ] ) takes the form @xmath250\left[\left(\lambda\,v_0\gamma_s\,\hat{f}(\theta_c ) \right)^{-1 } - m\   \text{sign}\left(t^p - \bar\tau m\right ) \right]\ , \ ] ] where @xmath251 and eq .
( [ dp_yielding ] ) were used .    to understand how a flow resistance naturally emerges within the model in the nonlinear response regime
, we need to analyze the fixed - points of eq .
( [ stability_exchange ] ) and their stability ( see fig .
[ fig3 ] ) . the first bracketed term in eq . ( [ stability_exchange ] ) vanishes only at the fixed point @xmath252 . in the second bracketed term
, @xmath130 reaches its fixed - point at @xmath253 , which depends itself on @xmath254 .
the @xmath254-dependence of @xmath255 is gaussian ( cf .
( [ kcd ] ) ) , so we expect the second bracketed term to have zero or two positive solutions .
we assume that the yet to be determined function @xmath224 is such that only one solution exists in the physical range @xmath256 . while this condition is not necessary on general grounds , it will be shown below to give rise to a physically appealing picture of yielding
. this single physical solution is denoted as @xmath257 and satisfies @xmath258 while @xmath11 is indeed a dynamic variable in the above expression , it is expected to vary slowly compared to @xmath254 , so we treat it as a constant for the ensuing stability analysis .    for stresses
@xmath227 such that @xmath259 , @xmath260 is the only stable fixed - point , see fig .
( [ dp_yielding ] ) implies the system approaches @xmath261 .
this is the sub - yielding regime in which after transient plastic deformation the flow arrests .
as the applied stress @xmath227 further increases , we reach the point where @xmath262 . at this point
an exchange of stability occurs , and @xmath263 becomes the stable fixed - point , see fig .
[ fig3 ] . in this case , eq .
( [ dp_yielding ] ) shows that @xmath264 , i.e. that the material flows persistently .
this happens at a threshold stress @xmath265 for @xmath266 which we interpret as a dynamic yield stress , responsible for the transition from non - flowing to flowing steady states .
it is apparently a function of the state of disorder through @xmath11 , ultimately through the details of @xmath224 , which can not be determined through thermodynamics alone .
it is clear , that eqs .
( [ stability_exchange])-([yield_stress ] ) make physical sense if @xmath224 leads to a solution with @xmath267 . in summary
, we have shown that a steady - state flow resistance in the model emerges as a dynamic exchange of stability in the stz orientation variable , and need not be stipulated a priori as done , for example , in the simple model discussed in section [ sec : conventional ] .
the steady state of configurational disorder , i.e. of @xmath11 , was discussed in section [ sec : conventional ] in the context of the simple phenomenological model , demonstrating its generic features .
we briefly discuss this issue again here in the context of a the internal variables model . to that aim
, we focus again on the configurational heat equation in ( [ heatc ] ) and observe that under rather generic situations in which @xmath268 , @xmath269 and @xmath118 , the only potentially positive term that can be responsible for rejuvenation is @xmath270 . when rejuvenation takes place , it can not persist forever and @xmath11 reaches a steady state under persistent plastic deformation , as is observed in numerous computer simulations of glass forming systems , e.g. those of @xcite . to see how it emerges
, we specialize to pure shear isothermal ( in @xmath12 ) situations as above , which allows us to omit all time and space derivatives in eq .
( [ heatc ] ) to obtain @xmath271 @xmath179 was replaced by its symmetric part @xmath180 because it is contracted with a symmetric tensor and we have simplified @xmath272 to @xmath273 , since steady pure shearing under the prescribed driving guarantees positivity of @xmath238 and @xmath233 .    to actually calculate the steady - state value
, we recall that ( as discussed several times above ) the thermodynamic coupling coefficient @xmath14 includes a deformation - driven part that depends on @xmath127 and a thermally - driven part that depends on @xmath12 . at sufficiently small @xmath12
, we have @xmath274 , where @xmath275 is a positive dimensionless function .
this leads to @xmath276 which is an implicit equation for the steady state value of @xmath11 .
the most interesting aspect of this result , beyond the existence of a steady state ( which is expected on general grounds ) , is that the steady state value of @xmath11 is independent of the plastic deformation rate @xmath127 and the externally applied stress . for sufficiently low temperatures
, @xmath12 can be also neglected .
these results are consistent with many numerical simulations as long as @xmath277 is much larger than molecular timescales , see for example @xcite .
in this paper we have developed systematically a finite - deformation , nonequilibrium , continuum thermomechanical framework for deforming solids wherein the configurational degrees of freedom evolve slowly compared to kinetic / vibrational degrees , and may be characterized by a temperature of their own .
this description of a solid as composed of weakly interacting configurational and kinetic / vibrational subsystems is inspired by recent progress in statistical physics @xcite . within the proposed framework ,
internal variables play an essential role in characterizing the out - of - equilibrium state of a deforming solid and are explicitly associated with distinct energy and entropy contributions .
an important deviation from standard thermodynamics is that coupling and transport coefficients associated with the configurational subsystem feature an explicit dependence on the plastic rate of deformation .
we derived a configurational heat equation , in addition to an ordinary heat equation , and showed how a strain - rate independent steady state of the configurational temperature naturally emerges under persistent deformation .
first , this general framework was applied to a simple phenomenological two - temperature plasticity model .
then , building on the concept of shear - transformation - zones ( stz s ) , and invoking a minimal scenario involving only one scalar and one tensorial internal variables , an internal variables based constitutive law was developed .
these constitutive models were shown to capture various amorphous plasticity phenomena , including isotropic and kinematic strain softening / hardening , newtonian viscosity and thermodynamic aging near the glass temperature , a sharp yield stress at low temperatures , flow driven structural rejuvenation and memory effects .
in addition , as was discussed in section [ sec : conventional ] , the diffusion of the configurational temperature  with the accompanying configurational heat conductivity @xmath129  might offer a new route to describing size effects in constitutive modelling .
we hope that the present development can serve as a physically - motivated and self - consistent thermomechanical framework for formulating constitutive models for the elasto - viscoplastic deformation of solids that lack internal equilibrium .
this should eventually allow to address a wide range of physical phenomena related to the shape deformations and failure modes of real amorphous materials ( e.g. bulk metallic glasses ) , which find an enormous range of modern engineering applications , but whose faithful description in the context of complex realistic problems is currently missing . in order to make progress in achieving this ambitious goal
, we plan to follow up the present investigation in several directions .
first , as was mentioned above , we plan to address pressure effects and the coupling between shear and volumetric deformation , which were not discussed here .
this is an important aspect of the deformation of amorphous solids and hence should be systematically incorporated into the proposed thermomechanical framework and constitutive modeling .
in addition , we need to perform a quantitative analysis of available experimental data to determine specific material parameters and functions , and to test the extent to which the proposed model can quantitatively account for various experimental protocols ( quasi - static loading , dynamic loading , strain - rate controlled experiment , stress controlled experiment , jumps in the applied strain - rate , creep , strain relaxation ) and control parameters ( strain rate , deformation temperature , thermal history etc . ) .
we then plan to consider spatially inhomogeneous deformation situations in 2d and to specifically focus on shear band widening dynamics , which involve a lengthscale associated with the diffusion of the configurational temperature .
finally , we plan to implement the proposed model in a finite element method code to allow testing and parameters calibration in various 3d problems .
we hope to be able to report on some progress along these lines in the not too far future .    * acknowledgements * + k.k
. acknowledges funding from nsf cbet grant 1253228 .
e.b . acknowledges support from the minerva foundation with funding from the federal german ministry for education and research , the israel science foundation ( grant no .
712/12 ) , the harold perlman family foundation and the william z. and eda bess novick young scientist fund .
* appendix *
our goal in this section is to augment our thermomechanical model with ideas from recent shear - transformation - zone ( stz ) theories @xcite . by doing this ,
we provide additional information to fill - in remaining details of the constitutive relations , culminating in a shear flow constitutive law . to motivate the connection
, we first discuss the stz micromechanism and its ability to replicate the @xmath278 and @xmath205 functional forms deduced in the previous sections .
we start by postulating the existence of a population of oriented stz s  sparsely distributed localized zones that are more susceptible to shear - driven transformations in a given direction than their surrounding material ( see fig . [ fig2]b ) .
we view stz s not as permanent flow defects , but rather as configurational fluctuations that can appear and disappear due to thermal and mechanically - generated noise ( i.e. the mesoscale stress fields that are generated by structural rearrangements ) .
therefore , the density of stz s essentially describes the probability to find such a configurational fluctuation .
the existence of stz s has been demonstrated in numerous computer simulations and experimental systems , see for example @xcite .
below we will adopt the main physical ingredients of the stz theory , which has been built on the flow - defect theories of turnbull , cohen , argon , spaepen and others @xcite , and incorporate them into the presently developed framework .    in two dimensions , each stz  during its lifetime  is described by its orientation @xmath279 with respect to a fixed reference direction . by definition , an stz oriented along @xmath279 is predisposed to undergo a shear ( volume - conserving ) transformation in the direction of the purely deviatoric , symmetric , unit - magnitude tensor @xmath280 whose major principal vector @xmath281 is oriented along the @xmath279 direction .
the components of the tensor take the form @xmath282 .
the population density of stz s at position @xmath17 at time @xmath283 oriented along @xmath279 is defined as @xmath284 , but we will generally write this without the @xmath17 and @xmath283 for notational simplicity . note that azimuthal periodicity implies @xmath285 .
we define @xmath130 to be the ( intermediate ) volume density of stz s , expressed as an integral over all orientations @xmath286 similarly , we define @xmath131 to be a normalized tensorial moment of the stz distribution in the sense that @xmath287 where the minus sign is used for consistency with previous analyses @xcite .
note also that by construction @xmath288 and that @xmath131 bears some resemblance to orientational order parameters in magnetic and liquid crystalline systems .
the stz s are assumed to be two - state systems that can be driven to `` flip '' from their original orientation @xmath279 to a perpendicular orientation @xmath289 during their lifetime .
while this is definitely an oversimplification as the stz s are multi - atom objects that can feature nontrivial internal structures , this assumption was shown to be sufficient to capture a range of glassy deformation phenomena @xcite .
one might characterize an stz oriented along @xmath279 as an elliptical region whose semi - major axis aligns with @xmath281 .
a `` flip '' is a pure shear transition , that deforms the stz inward along the @xmath281 direction and outward along @xmath290 , leaving the stz , post - transition , oriented in the @xmath290 direction
. the rate at which stz s oriented along @xmath279 undergo a shear transformation into a state oriented along @xmath291 depends on the resolved stress , @xmath292 , and thermal fluctuations , as specified by a rate - factor @xmath293 . for simplicity
we omit the arguments @xmath12 and @xmath11 in writing the rate - factor hereafter .
note that we assume that stz transitions are driven by the dissipative stress @xmath157 .
the evolution law for @xmath294 must account for the change in the stz s orientation when it undergoes shear transformations , as well as for the creation and annihilation of stz s due thermal and mechanically - generated structural fluctuations .
therefore , we suggest that the evolution of @xmath294 satisfies the following master equation @xmath295 term by term , this equation represents the rate - of - change of @xmath294 in terms of the rate of stz s transitions _ into _ a state @xmath279 , minus those transitions _ out of _ a state @xmath279 , plus a production / annihiation term that tends to equilibrate the stz s population with @xmath11 in an isotropic manner .
@xmath203 is a thermodynamic coupling coefficient that plays the role of a noise source that drives the creation and annihilation of stz s .
we must emphasize that this form is consistent with the previous thermomechanical analysis ; it agrees exactly with our second - law constraint , eq .
( [ kcd ] ) , i.e. @xmath296 the integral simplifies as above because the transition terms do not affect the total number of stz s    macroscopic plastic deformation arises in this model from the accumulated effect of all stz s transitions occurring at a given time .
hence , we express the plastic shear ( deviatoric ) flow rate @xmath180 as @xmath297 where @xmath207 is ( half ) the shear strain magnitude of an stz transition .    taking the time derivative of @xmath131 defined in eq .
( [ mdef ] ) and using eqs .
( [ lambdadef])-([int_flow ] ) , we obtain the following closed - form equation for @xmath205 @xmath298 it takes the form assumed in eq .
( [ assumed_dotm ] ) , allowing us to identify @xmath299    equation ( [ int_flow ] ) becomes an explicit flow rule when the function @xmath300 is specified . in the most general case
this would lead to an integral representation of the shear flow rate , that is not necessarily reducible to a simple closed form relation in terms of @xmath130 and @xmath131 alone . in order to explore the properties of the emerging constitutive law , we proceed by choosing a prototypical form for the @xmath294 distribution compatible with some given @xmath130 and @xmath131 values :
@xmath301    it is important at this point to reflect on the status of our analysis so far .
we have at hand a flow rule ( eq . ( [ int_flow ] ) , together with eqs .
( [ levolve ] ) , ( [ mevolve ] ) and ( [ n_approx ] ) ) , which is not fully determined until the rate factor @xmath300 is specified .
in addition , we still have to satisfy the second law inequality in ( [ plast_dissip ] ) , where @xmath302 is not yet specified ( note that @xmath303 is known , once @xmath302 is specified , through the relation @xmath304 ) . there are ( at least ) two possible routes to proceed ; either we model @xmath302 based on physical , but not purely thermodynamical , considerations and use the second law inequality in ( [ plast_dissip ] ) to constrain the rate factor @xmath300 or we model the rate factor @xmath300 based on physical considerations and use the second law inequality in ( [ plast_dissip ] ) to constrain @xmath302 .
the two routes are expected to have different consequences and it is not a priori clear which one is more physically faithful .
below we follow the second route by adopting a simple rate - factor @xmath300 and exploring the consequences through several approximations .
we consider a simple thermally activated process for stz transitions , similar to the form for atomic jumps in @xcite , for which @xmath305 \ , \ ] ] where @xmath208 is a molecular attempt rate and @xmath206 is an stz transition barrier .
using this simple rate factor in the flow rule in eq .
( [ int_flow ] ) , one can perform the integration analytically to obtain the following closed - form expression @xmath306,\ ] ] where @xmath307 and @xmath308 are the first- and second- order modified bessel functions of the first kind and @xmath309 and @xmath310 .
expanding eq .
( [ simple_flow ] ) to first - order in @xmath311 , we obtain @xmath312    this simplified flow rule has been used in the body of the manuscript for illustrative purposes without insisting on satisfying the condition @xmath313 .
next , we need to make sure that the flow rule in eq .
( [ d_1st ] ) does not violate the second law constraint in ( [ plast_dissip ] ) , which as explained previously , amounts to constraining @xmath302 .
we then assume the following constitutive structure @xmath314 , i.e. that @xmath156 is linear in the orientation tensor @xmath131 , which is then substituted together with @xmath180 of eq .
( [ d_1st ] ) in the inequality in ( [ plast_dissip ] ) .
the resulting inequality is satisfied if we set @xmath315    the expression for @xmath156 in eq .
( [ tp_store_result ] ) is interesting and deserves some discussion .
the prefactor relating @xmath156 and @xmath131 is proportional to @xmath316 , which implies that @xmath156 in this analysis is of _ entropic _ nature .
more specifically , it suggests that configurational fluctuations  whose magnitude is determined by @xmath317  generate an entropic stress proportional to the orientation tensor @xmath131 . recalling the meaning of @xmath318 ( defined in eq .
( [ m_entropy ] ) ) , a lower bound value for @xmath318 is @xmath319 , obtainable by analyzing the minimal case where flow defects have two possible internal orientations .
hence , @xmath320 , which gives positivity of the prefactor in eq .
( [ tp_store_result ] ) .
here we show the constitutive dependences discussed in section [ constit_sys ] are objective or frame - indifferent .
we suppose a rigid - body change in frame of the observer , @xmath321 where @xmath322 is an arbitrary time - varying rotation about some origin @xmath323 and @xmath324 is an arbitrary translation vector .
hence , @xmath325 is the motion as observed by the new observer ; we will continue to use @xmath326 to describe quantities as seen in the new frame .
frame - indifference is upheld if satisfaction of the constitutive functions in the original frame implies satisfaction in the new frame . to show this , we first need to transform the field variables into their corresponding forms in the @xmath326 frame and then check that the same constitutive functions correctly relate the transformed variables .
we start with the fundamental transformations : @xmath327 under the krner - lee kinematic decomposition , the @xmath49 transformation above is satisfied by letting @xmath328 ( cf .
@xcite ) . consequently
, @xmath329 and @xmath330 combining the latter expression with the transformation for @xmath25 , the elastic and plastic stresses obey @xmath331 the rate of stored plastic power @xmath332 must , as a material scalar , be frame - invariant ( i.e. @xmath333 ) . since @xmath179 has just been shown to be frame - invariant as well , we deduce @xmath334 , and thus @xmath335 .
all scalar constitutive quantities ( e.g. @xmath100 , @xmath338 , @xmath39 , @xmath339 , @xmath11 , @xmath12 , @xmath130 , @xmath7 ) are inherently frame - invariant . combining this fact with the transformation rules above , we have now shown that the dependent and independent variables related under the constitutive functions @xmath340 , @xmath341 , @xmath342 , @xmath343 , @xmath344 , @xmath345 , @xmath346 proposed in sections [ sec : conventional ] and [ constit_sys ] , are all frame - invariant . likewise , the functional dependences as proposed necessarily satisfy frame - indifference .
lastly , we check that the heat - flux constitutive functions are also objective .
let us temporarily bring back the sub- and superscript @xmath347 to denote intermediate space variables ( recall they were introduced and then dropped in section [ math ] ) .
the heat flux vector in the current configuration @xmath38 transforms under the rotating new frame as @xmath348 .
pulled back to the intermediate space , we have @xmath349 similarly the temperature gradient @xmath350 transforms as @xmath351 .
likewise @xmath352 we thus have that @xmath353 and @xmath354 are both frame - invariant . the same argument
gives that @xmath355 and @xmath356 are also frame - invariant .
consequently , all dependent and independent variables of the functions @xmath357 and @xmath358 in section [ constit_sys ] are frame - invariant , implying that the functional dependences are themselves indifferent to frame .
anand , l. , su , c. , 2005 . a theory for amorphous viscoplastic materials un- dergoing finite deformations , with application to metallic glasses .
journal of the mechanics and physics of solids 53 ( 6 ) , 13621396 .
anand , l. , aslan , o. , chester , s.  a. , 2012 . a large - deformation gradient theory for elastic - plastic materials : strain softening and regularization of shear bands
. international journal of plasticity 30 - 31 , 116143 .
bouchbinder , e. , langer , j.  s. , procaccia , i. , 2007 .
athermal shear - transformation - zone theory of amorphous plastic deformation .
ii . analysis of simulated amorphous silicon .
physical review e 75  ( 3 ) , 036108 .                                                                  gurtin , m.  e. , anand , l. , 2005 .
a theory of strain - gradient plasticity for isotropic , plastically irrotational materials .
part i : small deformations .
journal of the mechanics and physics of solids 53  ( 7 ) , 16241649 .
henann , d. , anand , l. , 2008 . a constitutive theory for the mechanical response of amorphous metals at high temperatures spanning the glass transition temperature : application to microscale thermoplastic forming .
acta materialia 56 , 3290-3305 .
henann , d. , kamrin , k. , 2013 .
a predicitive , size - sensitive continuum model for dense granular flows .
proceedings of the national academy of sciences of the united states of america 110(17 ) , 67206735 .
lu , j. , ravichandran , g. , johnson , w. , 2003 .
deformation behavior of the zr41.2ti13.8cu12.5ni10be22.5 bulk metallic glass over a wide range of strain - rates and temperatures .
acta materialia 51  ( 12 ) , 3429  3443 .
pan , d. , inoue , a. , sakurai , t. , chen , m.  w. , 2008 . experimental characterization of shear transformation zones for plastic flow of bulk metallic glasses .
proceedings of the national academy of sciences of the united states of america 105  ( 39 ) , 1476914772 .
picard , g. , ajdari , a. , lequeux , f. , bocquet , l. , 2004 .
elastic consequences of a single plastic event : a step towards the microscopic modeling of the flow of yield stress fluids .
european physical journal e 15  ( 4 ) , 371381 .
shi , y. , katz , m.  b. , li , h. , falk , m.  l. , 2007 .
evaluation of the disorder temperature and free - volume formalisms via simulations of shear banding in amorphous solids .
physical review letters 98  ( 18 ) .        song , c.  m. , wang , p. , makse , h.  a. , 2005 .
experimental measurement of an effective temperature for jammed granular materials .
proceedings of the national academy of sciences of the united states of america 102  ( 7 ) , 22992304 .