Graphical Conversion between Compliance and Modulus, Permittivity and Electric Modulus, and Impedance and Admittance

Spectrometries probing relaxation and retardation phenomena, such as dielectric,mechanical, and impedance spectroscopies, often require the analyses with both susceptibilities spectra z and its reciprocals 1/z (e.g., complex permittivity and electric modulus, mechanical compliance and mechanical modulus, and impedance and admittance). In the present paper, the geometric relation between z and 1/z is derived and the procedure to convert z into 1/z on a Cole-Cole diagram is proposed.This method helps us to relate them intuitively and yields clearer understanding on their interrelations. Moreover, it opens the new route for the geometric approach to derive many mathematical properties of spectra. The relation between peak position of z spectrum and that of 1/z spectrum and the shape of spectra are discussed on the basis of this method.


Introduction
In spectrometries focusing on relaxation or retardation phenomena, such as dielectric spectroscopy [1], dynamic mechanical spectroscopy [2], and impedance spectroscopy [3], not only frequency dependent susceptibilities but also their reciprocals are often subject to consideration.In dielectric spectroscopy, for example, the complex permittivity is usually employed to discuss the dynamic nature of localized charges, such as the rotational diffusion of dipoles, while its reciprocal, electric modulus, is also employed to assess the dynamics of delocalized charges, such as the migration of ions [4][5][6].This is also the case for the mechanical spectroscopy.Actually, both of the mechanical compliance and the mechanical modulus which is the reciprocal of the compliance are measured and discussed depending of the phenomena of interests [7].
Generally speaking, extremely broad frequency range can be probed in these techniques: it often exceeds 10 orders of magnitude in the frequency [1][2][3].Thus, obtained broad spectra often contain a variety of the dynamic processes.Some of them are compatible with the analysis in susceptibility whereas the others may be better analyzed in its reciprocal representation.In this case, one needs to switch the quantity of consideration from one to the other depending on the purpose.Therefore, the straightforward method, which enables us to change from one to the other, helps to understand the experimental data and makes their interrelations clearer.
Graphical methods can be a powerful tool to satisfy this demand.For example, the Smith chart which is indispensably used in the field of microwave engineering [8][9][10] graphically relates the complex reflection coefficient and the load impedance and so forth.Of course, since these quantities are related by mathematically rigorous formulae, one can algebraically calculate any quantities once another quantity is given.Nevertheless, the intuitive manner provides clear perspectives of the analysis.Furthermore, it can also provide another route to analyze the mathematical properties of spectra, that is, geometric method, in addition to the conventional algebraic approaches.
A Cole-Cole diagram [11,12] is one of the most frequently used representations for relaxation spectra.Briefly speaking, it is a plot of a complex quantity  (or its complex conjugate, depending on the definition) on the complex plane, where complex permittivity, compliance, and impedance are often chosen as .Plotting  varying the frequency yields characteristic trajectory reflecting the phenomenon.Remarkably, some elementally models, such as single exponential relaxation function in dielectric spectra and the behavior of a parallel RC circuit in impedance spectra, describe circular trajectories as frequency is varied.An analogous diagram is also employed to plot transfer functions and is called Nyquist diagram.The parameters of the transfer function can be determined and the stability of a system can be analyzed thereby with it.
In the present paper, I propose a straightforward method to convert the spectrum of arbitrary complex quantity  into its reciprocal 1/ based on the graphical procedure on Cole-Cole diagram.Thereby, interrelation between  (permittivity, compliance, and impedance) spectra and 1/ (moduli, admittance) can be graphically recognized.Moreover, it also opens the route to the geometric method to assess the mathematical properties of spectra.

Theory
2.1.Definition.We consider a pair of reciprocally related complex quantities,  and 1/.Plotting  or its complex conjugate in the complex plane yields a Cole-Cole diagram.Loci of  trace a characteristic trajectory on the diagram as an internal parameter is varied (typically, frequency), reflecting the phenomenon represented by .
From the viewpoint of the similarity of the typical trajectories on a Cole-Cole diagram, it is convenient to define generalized compliance  *  () representing permittivity, mechanical compliance, and impedance, whereas electric modulus, mechanical modulus, and admittance are represented by generalized modulus  * (), as summarized in Table 1.The generalized modulus is the reciprocal of the generalized compliance and vice versa.Although both  *  () and  * () can be plotted onto the Cole-Cole diagram as the complex quantity , only the former situation is to be described in the following discussion to avoid redundancy.

Modulus Grids Superimposed onto the Cole-Cole Diagram.
When the modulus is decomposed into real and imaginary parts as  * () =   +   , the real and imaginary parts of the corresponding compliance ( *  () =    −    ) are If the compliances at various frequencies are plotted on complex plane with    along the horizontal axis and    along the vertical axis, the well-known Cole-Cole diagram for compliance  *  () is yielded.Now, consider how iso-  and iso-  , the lines along which   and   are constant, appear in the Cole-Cole diagram.Sets of these iso-  and iso-  lines give a coordinate grid for the modulus  * () superimposed onto the Cole-Cole diagram for the compliance  *  ().These   and   grids enable us to convert compliance to modulus graphically once a trajectory on the Cole-Cole diagram is given.
To find the iso-  line (the trajectory along which   is constant) for given   , substitute   =   and   =  1 in (2a) and (2b), where   is just a constant and  1 is a parameter to be eliminated: Equation (3a) yields and dividing (3b) by (3a) yields Substituting (3c) for (3a  ) to eliminate the parameter  1 , we obtain a quadratic equation in    and    .Then, completing square leads to an equation for a circle: The radius of the circle is 1/2  and the center is at (   ,    ) = (1/2  , 0).In addition, one point of this circle is fixed at the origin, irrespective of the radius.As the value of the real part of the modulus,   , is decreased, the size of the circle increases.Therefore, iso-  lines on the Cole-Cole diagram of compliance appear as circles whose center is on the horizontal axis with one point fixed at the origin, as shown in Figure 1.
The iso-  line (the trajectories along which   is constant) for given   can be found similarly, by substituting  The grid is composed of iso-  and iso-  lines, along which values of   and   are constant, respectively (see Section 2.2 for the detail).Orange dash dotted line circles indicate iso-  lines, whereas green dashed circles indicate iso-  lines.Full line (i) indicates a trajectory circumscribed to an iso-  line, exhibiting maximum   at the tangent point.Full line (ii) indicates a trajectory inscribed to an iso-  line, exhibiting minimum   at the tangent point.
=  2 and   =   into (2a) and (2b).Iso-  lines are circles with radius 1/2  and center (   ,    ) = (0, 1/2  ) on the vertical axis.Again, one point of this circle lies at the origin, regardless of the radius, and the size of the circle increases as   decreases (see Figure 1).
The calculated iso-  and iso-  lines form coordinate grids superimposed on the Cole-Cole diagram of the compliance.All these circular grid lines pass through the origin and increase in size as the value on the grid line decreases.The Cauchy-Riemann theorem tells us that an arbitrary mapping via an analytic complex function is a conformal mapping in which the angle encompassed by an infinitesimally small triangle is preserved [13].Conversion from modulus to compliance is also conformal mapping since the complex function employed, () = 1/, is analytic except at the origin.It is obvious that iso-  and iso-  lines are orthogonal to one another on the complex plane of  * .According to the Cauchy-Riemann theorem, this relation is conserved under mapping by the analytic function () = 1/, and therefore iso-  and iso-  curves still cross orthogonally even on the complex plane of  *  .It can be seen in Figure 1 that two arbitrary intersecting circles always cross orthogonally.

Examples for Applications
With the aid of the circular modulus grids obtained above, one can obtain graphically the appearance of the modulus spectrum from the trajectory on the Cole-Cole diagram.In this section, I will present several applications for the modulus grids.
One of the major applications for the modulus grid is to find a peak in the modulus spectrum.As is clear from the circular shape of the modulus grids, maxima or minima of   and   are at points where a circular modulus grid line is tangent to the trajectory of the compliance (lines (i) and (ii) in Figure 1).Rigorously speaking, frequency derivatives of   or   are zero at the point where an   -grid line or an  grid line is tangent to the compliance trajectory, respectively.Moreover, if the trajectory is circumscribed to the circular grid line (line (i) in Figure 1), the frequency dependence of   or   is concave down, whereas if the trajectory is inscribed to the circular grid (line (ii) in Figure 1), the frequency dependence of   or   is concave up.By using these properties of the modulus grids, it is possible to find International Journal of Spectroscopy The trajectory of the Debye function describes a semicircle on the Cole-Cole diagram given by Apparently, the peak in the compliance spectrum is at the top of the semicircle for this trajectory.To find the peak in the modulus spectrum, a circular   -grid is drawn keeping its center along the vertical axis and a point fixed to the origin (see the orange dash dotted line in Figure 2).When the grid line is tangent to the semicircle of the trajectory of the Debye relaxation, the peak in the modulus spectrum is found at the tangent point.Then, the value of the peak height (radius of the circular   -grid) can be straightforwardly calculated based on a geometrical analysis.As is clear from Figure 2, the peak in the modulus spectrum is always on the left-hand side of the semicircle of the trajectory.Taking into account that the locus on the Debye semicircle moves from right to left with increasing frequency, we can graphically confirm the well-known statement that the peak is at higher frequency in the modulus spectrum than the compliance spectrum.Considering   -grid line which is tangent to the trajectory of the Debye relaxation, two situations are possible: one circumscribes to it and the other inscribes (see green dashed lines in Figure 2).The tangent point of the former is at infinite frequency and indicates that there is a maximum in   at this point, while that of the latter is at zero frequency and it is seen that a minimum in   shows up there.As a second example, let us find an   -peak for the complex compliance composed of the constants    and the    inversely proportional to the frequency (Figure 3): where  ∞ and  are parameters independent of frequency.If the generalized compliance represents complex permittivity, this represents a parallel circuit of a capacitor and resistor with constant permittivity and dc conductivity.If the compliance represents mechanical compliance, this represents a Maxwell model [2] of a spring and a dashpot connected in series.As obvious from the straight line of the trajectory drawn on the Cole-Cole diagram shown in Figure 3, the compliance spectrum exhibits neither maxima nor minima within finite frequency ranges.On the other hand, the existence of a circular   -grid tangent to the straight line of the trajectory indicates that a peak shows up in the modulus spectrum.From the geometric restriction on the   -grid, it is graphically proved that the peak of the modulus spectrum of ( 8) occurs at the frequency where    =    .Before proceeding with further examples, I will state a theorem that is useful for finding modulus peaks.When the trajectory on the Cole-Cole diagram has a circular shape, the peak position of the modulus spectrum can be found via a routine procedure, irrespective of the position of the center of the circle.Such circular trajectories are found in a variety of cases: not only the Debye function but also the Cole-Cole function (CC function) [14] and Van Vleck-Weisskopf function (VW function) [15][16][17].The statement of this theorem is as follows.

Theorem 1. When the trajectory on the Cole-Cole diagram is circular, the peak of 𝑀 󸀠󸀠 appears at the intersection of the circle and the straight line running from the top of the circle to the origin, irrespective of the position of the circle center.
Proof.The proof, which is based on elementary geometry, is straightforward.Consider a circle centered at arbitrary position  as drawn in Figure 4, and suppose that another circle with center   and passing through the origin is tangent to the first circle at .The line segment    cuts both circles at the tangent point .Since   ‖ , alternate-interior angles are equal, ∠   = ∠.Therefore, the isosceles triangles R   and R are homologous to one another.It follows that ∠   = ∠.Since vertically opposite angles, ∠   and ∠, are equal, line segments  and  turn out to be a continuous straight line .Therefore, the straight line running from  to  crosses circle with center  at the common tangent point.

Two applications of this theorem are presented as follows.
Although it is less important, an analogous theorem for the   -grid can also be proved."When the trajectory on the Cole-Cole diagram is circular, a peak of   appears at the intersection of the circle and the straight line running from the extreme right-hand point on the circle to the origin, irrespective of the position of the circle center." The CC function [14] is given by where  is a parameter that controls the broadness of the spectrum.This spectrum function also has a circular trajectory on the Cole-Cole diagram although the center of the circle is in the lower half plane.The trajectory on the Cole-Cole diagram is Figure 5 shows an example of such a trajectory (labeled "CC").Invoking the theorem, the peak position of the modulus spectrum can be found straightforwardly.As shown in Figure 5, by drawing a straight line from the top of the circle to the origin, the peak in modulus spectrum is found to be at the crossing point.It should be stressed that the peak was found without drawing an   -grid.

International Journal of Spectroscopy
A similar analysis is also possible for the VW function which represents a vibrational mode [15][16][17].The VW function is Although the trajectory of a VW function on the Cole-Cole diagram slightly deviates from a perfect circle in proximity to the horizontal axis, it is approximately circular when the vibrational frequency  0 is much faster than the damping rate,  0  ≫ 1, and | ln(/  )| ≪ 1, where   is the peak angular frequency in compliance, given by    = √( 0 ) 2 + 1.Then the approximated circle is centered in the upper half plane on the Cole-Cole diagram representing the characteristic sigmoidal shape of the    spectrum around the vibrational frequency.Again, the peak in the modulus spectrum can be found by drawing a straight line from top of the circular trajectory to the origin as shown in Figure 5 (labeled "VW").Then, the peak is found at the point where the straight line crosses the circular trajectory.Furthermore, it can be clearly seen from the Cole-Cole diagram that the real part of the modulus,   , also exhibits maxima and minima that are similar to those of the    spectrum, since it is possible to draw circular   -grid lines circumscribed and inscribed to the VW trajectory (these circles are not shown in Figure 5).This is not the case for the CC function, as   -grid lines cannot be tangent to the function's trajectory within the upper half of the Cole-Cole diagram.
The circular modulus grids tell us graphically that the peak of   always occurs at a higher frequency than that of    since circular   -grid lines that are tangent to the trajectory always lie on the left-hand side of the trajectory if it has tangent point.This is one of the well-known relations between peaks in the modulus and the compliance.Additionally, the position of the tangent point relative to the center of the circular trajectory depends on the values of  ∞ and Δ.With increasing  ∞ (moving the circle to the right while keeping its size unchanged), the tangent point approaches the peak of    , meaning that the frequency of   peak approaches that of the    peak.As Δ is increased (making the circle bigger while keeping its left-hand end point fixed), on the other hand, the tangent point approaches the point where (   ,    ) = ( ∞ ,  ∞ ), resulting in infinite frequency.Furthermore, information about loss tangent can be also associated with this graphical method.On the Cole-Cole diagram the loss tangent, tan  =    /   , corresponds to the slope of a line from the origin to an arbitrary point on the trajectory of compliance.The maximum of the loss tangent is at the point where a straight line running from the origin is tangent to the trajectory circle.Since the straight line used to find the modulus peak ( in Figure 4) has a lower slope than the tangent to the trajectory circle, the peak of the loss tangent lies between the peaks of the modulus and the compliance ( and  in Figure 4, resp.).Thus, we have reproduced several well-known relations between the peaks of modulus and compliance based on graphical considerations.It should be noted that if the modulus is plotted as Cole-Cole diagram instead of the compliance, typical trajectories run in the opposite direction: from left to right with increasing frequency.Then, a tangent point on the left-hand side of the trajectory means that peak of the compliance appears at a lower frequency than that of modulus peak, which is consistent with the previous results.
Although only circular trajectories have been considered above (straight line is infinitely large circle), some noncircular functions are also widely employed to reproduce experimental results: for example, Cole-Davidson (CD) function [18] and Havriliak-Negami (HN) function [19].Reflecting their asymmetric distribution of the relaxation times, trajectories of CD and HN functions on Cole-Cole diagram draw asymmetrically skewed arcs.In these situations, the theorem to find modulus peak described above is no longer applicable because it is limited only to circular trajectories.Hence, the present graphical method enables us to make only qualitative analysis for noncircular trajectories.By drawing  grid lines, qualitative relation between peaks of compliance, modulus, and loss tangent can be straightforwardly found.More generally, one can predict the qualitative shape of the modulus from a complicated spectrum as explained in what follows.In Figure 6, superposition of two HN functions is plotted as Cole-Cole diagram (its frequency dependence is also shown in the inset).As is clear from the diagram, the HN process at higher frequency exhibits only a shoulder instead of a peak in the compliance representation.By drawing two   -grid lines which tangent to the trajectory (orange dash dotted lines in Figure 6), appearance of two peaks in the modulus spectrum can be recognized.Since tangent line at lower frequency is nearly parallel to the trajectory for large range, it is estimated that the corresponding peak in the modulus is not well resolved (it is peak but close to shoulder).Such qualitative discussion is always possible irrespectively of the shape of spectra and provides clear understanding on the interrelation between processes in the modulus and those of compliance.

Conclusions
General properties of modulus grids (formed from iso-  and iso-  lines) on the Cole-Cole diagram for compliance were studied in detail.Both iso-  and iso-  lines are International Journal of Spectroscopy Maxima and minima of   and   correspond to the tangent points where   and   grid lines, respectively, meet the trajectory of compliance tangentially.Therefore, the peak of the modulus can be found graphically from the trajectory.When the trajectory is a circle (Debye, Cole-Cole, or Van Vleck-Weisskopf functions), the exact position of the modulus peak can be found by drawing a straight line from top of the circle to the origin and taking the point at which the straight line intersects with the circle.With the aid of modulus grids, shapes of modulus spectra can be graphically derived once the Cole-Cole plot of the compliance has been obtained.This procedure enables us to relate compliance with modulus much easier and makes their interrelation clearer.

Figure 1 :
Figure 1: Modulus grid superimposed on a Cole-Cole diagram.The grid is composed of iso-  and iso-  lines, along which values of   and   are constant, respectively (see Section 2.2 for the detail).Orange dash dotted line circles indicate iso-  lines, whereas green dashed circles indicate iso-  lines.Full line (i) indicates a trajectory circumscribed to an iso-  line, exhibiting maximum   at the tangent point.Full line (ii) indicates a trajectory inscribed to an iso-  line, exhibiting minimum   at the tangent point.
This modulus grid enables us to predict the appearance of the modulus spectrum for a given compliance trajectory on the Cole-Cole diagram.If the modulus is plotted on Cole-Cole diagram instead of compliance, similar circular    and    gridlines are obtained.Therefore, the following discussion is relevant in both cases.

Figure 2 :
Figure 2: Cole-Cole diagram for a Debye function ( ∞ = 4 and Δ = 16) (Black full circle).An iso-  circle tangent to the trajectory of the Debye function is shown as an Orange dash dotted line circle (  = 0.1).Two iso-  circles tangent to the Debye circle are also depicted as green dashed lines (  = 0.25 and   = 0.05).The left-hand end of the trajectory corresponds to infinite frequency whereas the right-hand end corresponds to infinitesimal frequency.

Figure 3 :
Figure 3: Cole-Cole diagram for  *  () =  ∞ − / ( ∞ = 10) (black full straight line).An iso-  circle touching the trajectory of the compliance is shown as an Orange dash dotted line circle (  = 0.05).Since both    and    have to coincide with the radius of the iso-  circle, the tangent point satisfies    =    .

Figure 4 :
Figure 4: Proof that the straight line from the top () of the trajectory circle with center  to the origin () passes through the tangent point () between the iso-  circle with center   and the trajectory circle.

Figure 5 :
Figure 5: Cole-Cole diagram for Cole-Cole (CC) and Van Vleck-Weisskopf (VW) functions, which are represented as off-center circles.Dark colored lines and symbols are for the VW function, while light colored ones are for the CC function.Full curve circles stand for trajectory and dash dotted circles indicate iso-  contours.The left-hand end of a trajectory corresponds to infinite frequency whereas the right-hand end corresponds to infinitesimal frequency.

Figure 6 :
Figure 6: Cole-Cole diagram for the compliance spectrum composed of superposition of two Havriliak-Negami (HN) processes.The compliance spectrum and corresponding modulus spectrum are shown in the inset.Full bold line stands for the trajectory of the total compliance and full thin lines indicate each component (HN1 and HN2).Dash dotted circles indicate iso-  contours.

Table 1 :
Summary of the definitions of compliances and moduli.() and () are electric displacement and electric field, respectively.  () and   () are shear strain and shear stress.  () and   () are longitudinal strain and longitudinal stress.() and () are voltage and current.