Dielectric Materials with Memory II: Free Energies in Non-Magnetic Materials
An isothermal theory of free energies and free enthalpies, corresponding to linear constitutive relations with memory, is presented for isotropic non-magnetic materials. This is a second paper, following recent work on a general tensor theory of isothermal dielectrics and on the form of the minimum free energy. Both papers are based on continuum thermodynamics. For a standard choice of relaxation function, the minimum and maximum free energies are given explicitly, using a method previously developed in a mechanics context. Also, a new family of intermediate free energy functionals is derived for dielectrics. All these are solutions of a constrained optimization problem.
A general theory of free energies for dielectric materials under isothermal conditions was given recently . The present work presents more detailed formulae and results for free energies in non-magnetic isotropic materials. We shall derive expressions for the minimum, maximum and other intermediate free energy functionals which also have an extremum property. Two approaches are possible: (1) To present the developments as in mechanics , based explicitly on thermodynamics; and (2) To use the method given independently [3-5] for dielectric materials with memory, but with quite diﬀerent methods, notation and terminology. In this work, as in , the fist approach will be adopted but with reference to the connection between the two methods. The intermediate free energies are analogous to quantities known in mechanics . They are new and particularly physically relevant in the context of dielectrics, in that they apply to memory models that are standard for such materials, but not usually applied to viscoelastic behaviour.
Applications of the results of this work, from a physical point of view, include the establishment of bounds on the level of dissipation in the material, using the fact that the total dissipation associated with the minimum(maximum) free energy is an upper(lower) bound on the actual physical dissipation.
On the matter of notation, vectors and tensors are denoted by lowercase and uppercase boldface characters respectively and scalars by ordinary script. The real line is denoted by , the non-negative reals by and the strictly positive reals by . Similarly, is the set of non-positive reals and the strictly negative reals. Complex quantities arise in the frequency domain so we have complex vector spaces for which the dot product involves using the complex conjugate of objects in the dual space. The magnitude squared, denoted by refers to the dot product of objects with their complex conjugates.
Consider a rigid non-magnetic isotropic dielectric material subject to a varying electric field. Let the body under consideration occupy a volume . A typical point in B is x while t is a given time. The electric field on this region is E(x,t), with electric displacement denoted by D(x,t). The magnetic field and induction contributions are neglected. The space variables are generally omitted.
Let be any free energy of the material, for the isothermal case, and D(t) the rate of dissipation. Then the first and second laws of thermodynamics can be written as 
Let be defined bya
aThis quantity would be denoted as in [3-5], while the future continuation, defined by (3.3), would be written as . The future continuation is referred to in  as the recovery field.
Equation (3.2) gives the history and also the current value E(t) of the electric field. A given future continuation is denoted by
Let us define the free enthalpy  as
In terms of this quantity, (3.1) becomes
The relative history is defined by
A relative future continuation is also defined by (3.6) for
We define the equilibrium free enthalpy to be that given for the static history .
emphasizing that the function depends only on E(t).
A. Required properties of a free enthalpy
Let us state the characteristic properties of a free enthalpy, provable within a general framework [1,6-8]:
P2: For any history and current value
P3: Condition (3.5) holds.
These will be referred to as the Graffi conditions by analogy with those for a free energy in mechanics ( for example).
A Linear Memory Model
A special case of the linear model introduced in  is described in this section. This is the material considered in [3-5], namely a passive, homogeneous isotropic non-magnetic dielectric. All kernels are scalar quantities. We have
where (4.1b) follows from (4.1a) by partial integrations. Also
with similar limits at large u holding for fixed s. The parameter in (4.1a) is defined by (4.3b). Relation (3.8) gives the linear constitutive relation
The quantity G(s) is the relaxation function of the material. In writing the final form of (4.4), we are assuming that vanishes; we furthermore assume that it goes to zero sufficiently strongly so that various integrals, introduced below, exist. The quantity is related to the relaxation function through
We deduce from (3.5), (4.1a) and the time derivative of (4.1b) that
where subscripts refer to partial diﬀerentiation with respect to the first or second argument. Equation (4.7b) follows from (4.7a) by partial integrations, using (4.1a), (4.3), and the fact that vanishes. Relation (4.1a), expressed in terms of histories, becomes
From (3.4), (4.2) and (4.8), we deduce that
The integral terms, with the negative signs included, are non-negative . The total work done by the electric field up to time t is
Under the assumption that , an integration of (3.1) over yields
is the total dissipation up to time t. We take to be equal to the vacuum permittivity [3-5], giving
where, using (4.4b), one sees that the polarization P(t) is given by the causal relationship
Then, (4.11a) becomes
where is the quantity introduced in (4.10) and is given by
is the accumulation of energy (density) transferred to the medium at the point under consideration, from the beginning of the pulse-medium interaction until time t. We can write (4.12) as
where is given by the integral term in (4.9a).
A. Minimal states
The fundamental definition of the state of a material with memory at time t is the history of the independent field variable and its current value (Et, E(t)). Also, diﬀerent histories may be members of the same minimal state . Two states , are equivalent, or in the same minimal state if from a time t onwards, we have
Where D1, D2 are defined by (4.4) for these states. Then it follows that
A fundamental distinction between materials is that for certain relaxation functions, namely those with only isolated singularities (in the frequency domain), the set of minimal states is non-singleton, while if some branch cuts are present in the relaxation function, the material has only singleton minimal states [2,9]. For relaxation functions with only isolated singularities, there is a maximum free energy that is less than the work function W(t) and also a range of related intermediate free energies, which are discussed later. On the other hand, if branch cuts are present, the maximum free energy is W(t).
In this work, we will deal exclusively with the case where the relaxation function has only isolated singularities, so that the minimal states will be non-singleton.
For such materials, the free energy functional is positive semi-definite (, page 152).
Note that the statement that and are equivalent is the same as the assertion that is equivalent to the zero state (0,0), where 0 is the zero in (and also the zero history), while
A functional of (Et, E(t)) which yields the same value for all members of the same minimal state will be referred to as a functional of the minimal state (abbreviated to FMS) or as a minimal state variable.
Let , be any equivalent states. Then, a free energy is a functional of the minimal state if
It is not necessary that a free energy have this property, though it holds for the minimum and all related free energies introduced later.
Kernels and Field Variables in the Frequency Domain
For any , we denote its Fourier transform byb
bThe quantity , defined by (5.1), would be denoted in [3-5] by or for real .
Only real valued functions will be considered so that
Where the bar denotes the complex conjugate. For complex values of ω, (5.2) becomes
Functions defined on are identified with functions on R which vanish identically on . For such functions,
Where , are respectively the Fourier cosine and sine transforms. A property of Fourier transforms which will be used later is
If is non-zero.
A. The kernel G(u)
We can write the Fourier transforms of and in (4.4) as
The quantity is the susceptibility, denoted by in [3- 5]. By partial integration, one can show that
giving, in particular, that
Using (5.2), we see that
which will be useful later. By considering periodic behavior in , we can show that a consequence of the second law is that (, page 140)
It follows from (5.11) that . We have , so that
B. The complex frequency plane and the function H(ω)
We will be considering frequency domain quantities, defined by analytic continuation from integral definitions, as functions on the complex ω plane, denoted by Ω, where
Similarly, Ω− and Ω(−) are the lower half-planes including and excluding the real axis, respectively.
The quantities , defined by (5.1), are analytic in Ω(∓) respectively (, page 547). Thus, the quantity is analytic on Ω(−). It will be assumed that is analytic on and thus on Ω−, or more precisely, on an open set containing Ω−. Also, it will be taken to be analytic at infinity. This function is defined by analytic continuation in regions of Ω+ where the Fourier integral does not converge. The quantity has singularities in both Ω(+) and Ω(−) that are mirror images of each other. It goes to zero at the origin and must also be analytic there. A quantity central to our considerations is defined by
It is a non-negative, even function of the frequency and goes to zero quadratically at the origin. The relation (see (5.5))
For the model considered later, H(ω) goes to zero at large ω so that H∞ = 0. The Fourier transforms of the history and continuation are denoted by and respectivelyc. These are particular examples of (5.1c) and (5.1d). The quantity is analytic on Ω(−) and is analytic on Ω(+). Both are assumed to be analytic on an open set including . It is further assumed that they are analytic at infinity.cReferring to footnote a, we see that would be denoted in [3-5] by .
The Fourier transforms of the relative history and continuation, defined by (3.6), have the form
where the notation ω± is defined as ω ± iα, α > 0. The parameter α is assumed to tend to zero after any integrations have been carried out (see for example , page 551). We have
The second relation follows from (3.2) and an integration by parts in the Fourier integral defining . Also, based on arguments from  and , page 146, for example, we can express the constitutive equation (4.4) in the form
The Work Function
Various forms for the work function, defined by (4.11), are given for the general tensor case in . Those required in this work for scalar kernels are summarized here. The analogous mechanics version of these relations may be found in , page 153, for example. It can be shown that
where U(t) and are defined by (4.10). We see from (6.1) that W(t) can be cast in the form (4.9) by putting
In terms of frequency domain quantities, we find that 
The quantity W(t)-D(t)•E(t) can be shown to obey the properties specified in subsection ⅢA of a free enthalpy, with zero dissipation. Because of the vanishing dissipation, it must be the maximum free energy associated with the material or greater than this quantity, an observation which follows from (4.12). For relaxation functions with only isolated singu-larities, as introduced later, there is a maximum free energy which is a functional of the minimum state and is less that W(t).
Consider the scalar product on the function space of electric fields, defined by [2,5,9-12]
Recalling (5.1), we see that the quantity is the Fourier transform of , in the time domain. Also, the subscript F in the bracket notation indicates that the frequency domain version is being used. The related norm is
In the notation of (6.5), the work function is given by
where , in (6.6a) and , in (6.6b).
Factorization of H (ω)
We consider materials such that (or ) has only a finite number of isolated singularities in Ω(+). Thus, H(ω) has only isolated singularities in Ω(±), which are mirror images of each other in the real axis, as ascribed to before (5.13)). This means that it can be put in the form of a ratio of polynomials. We will take the singularities to be a finite number of simple poles. The quantity H(ω) has a finite number of zeros in Ω±(ω), also mirror images of each other. It is real and non-negative on , an even function of ω and therefore a function of ω2, in view of its analyticity about the origin. It vanishes quadratically at the origin. This non-negative quantity can be factorized in general as outlined in [2,13]. Thus, we have that
The factorization (7.1) is unique up to multiplication by a phase factor eiα, where α is a constant. Relation (7.2a) reduces this arbitrariness to multiplication by a factor ±1.
The factorization (7.1) is equivalent to a homogeneous Riemann-Hilbert problem on the half-plane [1,3-5].
The quantity H+(ω) (H−(ω)) has all its singularities in Ω(+) (Ω(−)) and all its zeros in Ω+ (Ω−) (, page 239). There are many other factorizations, obtained by interchanging some or all of the zeros of H+(ω) and H−(ω), while retaining the same singularity structure. The diﬀerent factorizations are labeled by the subscript or superscript f. We have
Certain types of exchanges must be excluded to ensure that (7.3b) are true , (, page 338). Each factorization generally yields a diﬀerent free energy, though there may be exceptions. All these free energies are FMSs . The factorization with no exchange of zeros, which is that given by (7.1), yields the minimum free energy Ψm(t).
Remark Ⅶ.1.Each exchange of zeros, starting from these factors, can be shown to yield a free energy which is greater than or equal to the previous quantity (, page 363).
Note that the zeros of at the origin play no part in these exchanges.
Remark Ⅶ.2. A particularly interesting choice of , which we denote by , is obtained by interchanging all the zeros. This can be identified as the maximum free energy among all those that are FMSsd. It is less than the work function, which is not a FMS for materials with only isolated singularities .dThe minimum and maximum free energies are given also in . The intermediate free energies, introduced in this work, have not been discussed previously in the context of dielectrics, but are analogous to those derived in mechanics in for example ,10] -12].
If there are N diﬀerent factorizations of H(ω), then f = N is chosen to denote the maximum free energy.
Note that there are several (indeed many, for a large number of isolated singularities) diﬀerent zero exchange pathways leading from the minimum to the maximum free energy.
The most general free energy and rate of dissipation arising from these factorizations is given by
A particular case of this linear form is the free energy proposed in .
Remark Ⅶ.3. The set of all free energies at time t associated with a given material is denoted by Φ(t). The boundary of this set is where one or more of the fundamental properties listed in subsection ⅢA is breaking down; an example, which will arise later, would be where, for a given non-zero history, the rate of dissipation is zero, so that a small shift in its kernel parameters results in it becoming negative.
All the emerge from extremum arguments. This is apparent in the case of the minimum free energy for a given state, which is obtained by finding an optimal continuation , yielding the maximum recoverable work from this state. Also, the maximum free energy is determined by finding the optimal history , which minimizes the work done to achieve the given state. It must be an equivalent state to the given history. In other words, the two histories must be in the same minimal state. The intermediate free energies are also obtained from an extremum principle, but involving an optimization of the history/continuation Et(s), , as we shall see in section Ⅷ.
The relevant theoretical developments motivating the results of section Ⅷ are presented in , chapter 15 for mechanics; they also apply to dielectrics. We will present an abbreviated version of these arguments, and also demonstrate that the functionals which emerge are in fact free energies, with the required properties. It can be shown that all these quantities, including the intermediate functionals, are on the boundary of Φ(t) for the material (, page 365), which is consistent with the fact that they emerge from an extremum principle. This property is discussed briefly below.
The Free Energy Associated with a Particular Factorization
Let us define the quantities and by the relation
where is the given history, which is zero for . The quantity is typically non-zero on and follows from (8.1). In the frequency domain, the relation becomes
where the quantities are particular cases of (5.1c) and (5.1d). Consider the norm , defined by (6.5). This can be written in a form similar to the quantity W(∞) in , from which the minimum free energy was deduced. We have, from (6.4b) and (6.5), that
Let us choose a particular factorization of H(ω), as given by (7.3). Also, we define 
where the are analytic in Ω∓, respectively. The singularities of are the same as those of , as may be perceived by closing the contour in (8.4b) on Ω(−). Singularities on the real axis are excluded by assumption. Using the factorization in (8.4), we can express in the form
Let us further define
Note that, by virtue of Cauchy's integral formula,
Where are the quantities defined by (8.4b) but with replaced by
The optimization problem which yields the free energy functional associated with a given factorization of H(ω) is given as follows (, pages 346, 350). We minimize , subject to the constraint
The resultant free energy is
where is the history/continuation that is the solution of the constrained minimization problem. The formulation of this problem is discussed at length in , chapter 15 and earlier papers. In particular, it is shown that each has all the required properties of a free energy, which will in any case be shown later for the present context.
Using (8.6), (8.7c) and (8.9) in (8.5) together with proposition 1 of , we find that
where is given by (8.8). Since depends only on the given history , the solution to the minimization problem is obtained by choosing to give that
From (8.11a) and (8.12), we obtain that
Then, (8.4a) gives
The observation after (8.4) about the singularities of and yield that the singularity structure of is determined by the zeros of .
Remark Ⅷ.1. Note that for , N, the quantity has some of its zeros in Ω− and some in Ω+. For f = 1, all the zeros of are in Ω−, so that has singularities only in Ω− and its inverse Fourier transform is thus non-zero only on . This latter quantity is the negative of the optimal continuation, yielding the minimum free energy, which is discussed for the general tensor case in . It can be shown that both the constraint (8.9) and the minimization of result in the same relation, namely (8.16) for f = 1.
For f = N, all the zeros of are in Ω+, which gives that has singularities only in Ω+. Therefore, its inverse Fourier transform is non-zero only on and is the optimal history yielding the maximum free energy. In this case, the quantity to be minimized is the work function and (8.9) is the condition that the optimum history and the given history are in the same minimal state.
For intermediate cases, is non-zero on and minimizes , subject to the constraints (8.9), for various factorizations. These constraints are transitional between the two extremes discussed above.
We will see that all of these cases result in free energies that are on the boundary of Φ(t), defined by remark Ⅶ.3. This property is apparent for the minimum and maximum free energies. For these and all intermediate cases, we will see that the corresponding rates of dissipation are positive semi-definite rather than positive definite quadratic forms. They are zero for certain non-zero histories so that small variations in their kernels may lead to negative rates of dissipation, for these histories, in contradiction to the second law, as given by (3.1).
Thus, the fundamental physical content of the constrained optimizations is that they lead to free energies on the boundary of Φ(t).
The free energy (8.10) has the form
by virtue of (8.13) and (8.15). We have the following relations ([1,10,13,14] for example)
Remark Ⅷ.2. For the material considered in section Ⅸ, the quantity H(ω) vanishes for large values of ω. Thus, Hsr in (8.18) and (8.19) is zero. Also, , (8.20)
as can be seen by closing the contour on Ω(+), so that, recalling (5.16), it is clear that (8.18c) is replaced by
The work function, given by (6.3), takes the form 
From (4.12), (8.17) and (8.22), we deduce that the total dissipation corresponding to the minimum free energy is given by
Differentiating this relation with respect to t and using (8.18a), (8.19c), gives
Since Kf(t), given by (8.21) can vanish for non-zero histories, Df(t) is a positive semi-definite rather than a positive definite quadratic form, as noted in remark Ⅷ.1.
We can re-express these results in terms of relative histories . Let us put
By closing the contour on Ω+, it emerges that
Also, the work function, given by (8.22), becomes
while (8.17b) and (8.23) are given by
Equation (8.24) is unchanged.
We now briefly present double integral forms of Ψf(t), Df(t) and D(t), which are analogous to (though simpler than) the corresponding formulae in  for the minimum free energy in the full tensor case. Using (8.17b) and (8.28a), we can write Ψf(t) in the forms
Also, Df(t), given by (8.24), can be expressed as
From (8.23) and (8.28b), we deduce that
One finds that
by integrating over ω2 for example and closing the contour on Ω−, since and have no singularity in the lower half plane. This allows us to write (8.29a) in the explicitly convergent form
The fact that the integral term is non-negative implies that 
By invoking (5.17), we see that the quantity may be replaced by in the above double integral expressions.
The free enthalpy Ff(t) corresponding to Ψf(t) may be deduced from (3.4) to be
and Fe(t) is given by (4.1c). It is easy to show that Ff(t) obeys the Graﬃ conditions listed in subsection ⅢA. Property P2 is immediately apparent. The relation (4.12) holds by virtue of (8.17), (8.22) and (8.23). The time derivative of (4.12) yields (3.1) with Df(t) given by (8.24), which is equivalent to P3. Property P1 can be proved with the aid of (8.35b), by showing that
Remark Ⅷ.3. It was shown in  (see also , page 340) that , defined by (8.4), is a function of the minimal state in the sense defined after (4.21). This result transfers to the present context without alteration. From (8.17), we deduce that Ψf(t) is a function of the minimal state as defined by (4.22).
For the minimum free energy Ψm(t), we denote Ψint(t), introduced in (4.18), as Ψrec(t), given by
where is defined by (8.4) for f = 1. This is the recoverable energy defined in . The quantity , given by (8.23) for f = 1, is that referred to as the irrecoverable energy in . Corresponding quantities and , for each f, may be similarly labeled.
The optimal history/continuation, given by the inverse transform of (8.1), can be shown to have a discontinuity at the origin which is infinite if Hsr = 0, while at t = ±∞, it is non-zero ([1,2], page 354).
A Detailed Dielectric Model
by virtue of (5.11). As in [3-5], we take the susceptibility to be modeled by a sum of Lorentz oscillatorseeThis behavior, while often used to describe memory behavior of dielectrics, is not usually applied to viscoelastic materials. The singularity structure given by (9.4) corresponds to exponential decay with sinusoidal behavior in the time domain, while viscoelastic materials are generally modeled by simple exponential decay, which in the frequency domain yields simple poles on the positive imaginary axis.
in terms of the oscillator strength fn, the plasma frequency ωpn, the resonant frequency ωn and the damping rate γn of each oscillator. All of these are positive parameters. It is assumed that
The singularities of χ+(ω) are of course isolated. They are simple poles at
Thus, they occur in pairs, equidistant from the positive imaginary axis, as required by (5.3). Separating the poles in (9.2), we can write
In the time domain, this translates into
which are decaying exponentials multiplying sine functions , so that there will be oscillatory behavior superimposed on the exponential decay. From (9.2) and (9.1), we have
which can be written as a ratio of two (factored) polynomials:
where the denominator of (9.8a) uses the notation of (9.4) while the numerator is factored to yield the zeros of H(ω); these must occur also in pairs, as in (9.4), so that for each l. We have explicitly included the fact that H(ω) vanishes quadratically at the origin. Note that the smaller number of zeros reflects the fact that H(ω) behaves as for large ω.
The factorization can be carried out by inspection. We obtain
The most general case of a rational function, which is considered in [2,11], can be obtained from (9.8) by allowing singularities to coalesce. As a result of this, some zeros may also coalesce. Indeed, we cannot exclude the possibility that some of the zeros in (9.8) have a power higher than unity, even if all the singularities are simple poles. For simplicity, it is assumed that this does not happen for our choice of parameters.
Let us define for l = 1, 2, . . ., 2N −2,
We can write (9.5) as
where (5.10) has been used and
In this notation, (9.6) becomes
Using (9.14b) in (4.20), we see that two electric field histories are in the same minimal state if and only if
Thus, a minimal state is defined by the quantities (e1,e1,...e2N), where
Relation (5.3) gives
Formula (5.17) for yields
Equation (4.4b) can be written as
with the aid of (9.14), so that
If a free energy or free enthalpy is to depend only on the el(t), l = 1,2,...,2N, then the most general forms of and in (4.2) and (4.9) are
where C and are non-negative hermitian matrices. Then, (4.9) gives that
Relation (4.5a) gives
Taking the complex conjugate of this relation yields
Equation (3.1) is satisfied if
The matrix with components Γij must also be hermitian and non-negative.
Let us now consider the free energies Ψf(t). We can put (9.9) in the form
where (7.3) has been used. The relationship
must hold since tends to at large ω. Therefore, we can also write (9.26b), (9.28c) and (9.28d) as
We identify also a much larger class of factorizations of H(ω), determined by interchanging particular βl in (9.26a) with in (9.28a). These diﬀerent factorizations are labeled by the subscript or superscript f. Thus,
Observe that (9.29) also holds for the , a property which has been used in writing (9.31c). There are 22N−2 diﬀerent factorizations. Referring to the discussion in sections Ⅶ and Ⅷ, we note that if all the are equal to one, then Ψf(t) is the minimum free energy Ψm(t), while if all are zero, we obtain the maximum free energy ΨM(t). All other possibilities yield functionals that are intermediate between these two extremes. These observations follow from remark Ⅶ.1. Observe that
Now, we have the relation
which follows by applying an integration over to (9.1) and completing the contour over Ω(−). The integral over vanishes. Thus, we have, on integrating over Ω+,
The quantity , defined by (8.4), may be evaluated by closing the contour on Ω(−), giving
From (8.4a), we have
which has singularities at those of in Ω(+) but none in Ω(−). These explicit relations allow analytic continuation of to the whole complex plane, excluding singular points. From (8.26a), (9.18), (9.31) and (9.35), we see that the quantity is given by
The optimal history/continuation in the frequency domain, which is a special case of (8.16), is given by formulae which are generalizations of those in , page 366. We deduce from (8.17b) and (9.35) that
Also, from (8.28a) and (9.37),
We write (9.38a) and (9.40a) in the form (9.22), giving
Comparing (9.34b) with (9.12), we see that
where (9.29), for the , has been used. The relation (4.5a), which yields that , where the latter quantity is given by (9.14), can be confirmed with the aid of (9.41) and (9.43). Relation (9.23), which is equivalent to (4.5a), can be confirmed for Cf.
Recall from (9.16) that the quantities e1, i = 1,2,...,2N represent a particular minimal state and all the Ψf(t) are manifestly FMSs.
From (8.18d), (8.21), (9.31b), and (9.31c), we have
This is in agreement with (9.25). Relations (9.38) and (9.45) can be derived also from (8.29) and (8.30).
Note that Df(t) vanishes if
Solutions to this equation will exist for non-zero values of ei(t). Therefore, (9.45b) is a positive semi-definite rather than a positive definite quadratic form, so that the associated matrix will have some zero eigenvalues. If one of these zero eigenvalues were to become slightly negative, then the second law would no longer hold. Thus, the free energy Ψf(t) is on the boundary of Φ(t), defined in remark Ⅶ.3. This observation is of course a special case of that after (8.24), and relates to remark Ⅷ.1.
A. Sinusoidal histories for non-magnetic materials
Let us consider the formulae (9.44) and (9.38) for histories
giving that el(t), defined by (9.16), has the form
The relation (9.38a) becomes
Also, from (9.44), Kf(t), which determines the rate of dissipation, has the form
Remark Ⅸ.1. General formulae for the minimum free energy and the associated rate of dissipation were given for sinusoidal histories in . These generalize immediately to apply to Ψf(t) and Df(t) by simply replacing H±(ω) (or the tensor version of these quantities used in ) with
Thus, the formulae for Ψf(t) and Kf(t) are as follows
where ∆f(ω) is given by (8.34). It will now be shown that relations (9.50) and (9.51) can also be obtained from these, but after somewhat difficult manipulations. We see from (9.26), (9.28) and remark Ⅸ.1 that
Referring to (9.34), we have
which agrees with the coefficient of in (9.50). From (9.55a), we see that
Also, using (9.30), we see that the quantity ∆f(ω0), defined by (8.34), is given by
Now, by virtue of (9.29), for the , it follows that
We deduce the relation
Using the symmetry of ∆f(ω0) with respect to ω0, we find that
which agrees with the constant term (proportional to |A|2) in (9.50).
The main results are listed below.
1. An isothermal theory of free energies, based on continuum thermodynamics and corresponding to linear constitutive relations with memory, is presented for isotropic non-magnetic materials. For a standard choice of relaxation function, the minimum and maximum free energies are given explicitly.
2. The central new result of this work is the determination of a family of intermediate free energy functionals, for the same relaxation function. These are analogous, though not identical, to those known within a mechanics framework. However, the assumed singularity structure of the frequency domain version of the relaxation function (the susceptibility, for dielectrics) is more relevant, in a physical sense, to dielectrics than to mechanics, as discussed in footnote 5, and this singularity structure is the basis of the new results.
3. An important aspect of these developments, from a physical point of view, is that all of these quantities are solutions of a constrained minimization problem, leading to free energy functionals that are on the boundary of the set of all free energies. Going outside of this boundary would yield functionals that do not obey properties imposed by the laws of thermodynamics.
- S Glasgow, JM Golden, Dielectric materials with memory I: Minimum and general free energies.
- G Amendola, M Fabrizio, M Golden (2012) Thermodynamics of Materials with Memory: Theory and Applications, Springer, New York.
- S Glasgow, M Meilstrup, J Peatross, et al. (2007) Real-time recoverable and irrecoverable energy in dispersive-dissipative dielectrics, Phys Rev E Stat Nonlin Soft Matter Phys 75.
- S Glasgow, M Ware (2009) Real-time dissipation of optical pulses in passive dielectrics, Physical Review.
- SA Glasgow, J Corson, C Verhaaren (2010) Dispersive dielectrics and time reversal: Free energies, orthogonal spectra and parity in dissipative media. Phys Rev E Stat Nonlin Soft Matter Phys 82.
- BD Coleman, EH Dill (1971) On the thermodynamics of electromagnetic fields in materials with memory. Archive for Rational Mechanics and Analysis 41: 132-162.
- BD Coleman, EH Dill (1971) Thermodynamic restrictions on the constitutive equations of electromagnetic theory, Zeitschrift für angewandte Mathematik und Physik 22: 691-702.
- V Berti, G Gentili (1999) The minimum free energy for isothermal dielectrics with memory. Journal of Non-Equilibrium Thermodynamics 24: 154.
- L Deseri, M Fabrizio, M Golden (2006) The concept of a minimal state in viscoelasticity: New free energies and applications to PDEs. Archive for Rational Mechanics and Analysis 181: 43-96.
- M Fabrizio, J Golden (2002) Maximum and minimum free energies for a linear viscoelastic material, Quaterly of Applied Mathematics 60: 341-381.
- J Golden (2005) A proposal concerning the physical dissipation of materials with memory. Quaterly of Applied Mathematics 63: 117-155.
- J Golden (2007) A proposal concerning the physical dissipation of materials with memory: The non-isothermal case. Mathematics and Mechanics of Solids 12: 403-449.
- J Golden (2000) Free energies in the frequency domain: The scalar case. Quaterly of Applied Mathematics 58: 127-150.
- L Deseri, G Gentili, JM Golden (1999) An explicit formula for the minimum free energy in linear viscoelasticity. Journal of Elasticity 54: 141-185.
- M Fabrizio, A Morro (2003) Electromagnetism of continuous media, Oxford University Press.
S Glasgow, Department of Mathematics, Brigham Young University, Provo, Utah 84602, USA.
© 2017 Glasgow S, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.