METHOD AND SYSTEM FOR ESTIMATING A DIFFUSION POTENTIAL OF A DIFFUSIVE PROPERTY
Methods, systems, and techniques for estimating a transient diffusion potential of a diffusive property involve modeling, as a circuit, diffusive behavior of a diffusion region and then simulating operation of the circuit to estimate the transient diffusion potential at a location in the diffusion region by determining circuit potential at a node in the circuit that corresponds to the location in the diffusion region. The circuit has steady-state and transient portions that model the steady-state and transient behavior of the diffusion region, respectively. The transient behavior is modeled using a capacitive circuit element. The diffusive property diffuses linearly within the diffusion region and generation of the diffusive property is distributed within the diffusion region.
Pursuant to 35 U.S.C. §119(e), this application claims the benefit of U.S. Patent Application No. 61/643,107, filed May 4, 2012 and entitled “Method and System for Estimating a Diffusion Potential of a Diffusive Property”, which is hereby incorporated by reference in its entirety.
TECHNICAL FIELDThe present disclosure is directed at methods, systems, and techniques for estimating a diffusion potential, such as a thermal potential (colloquially referred to as “temperature”), of a diffusive property, such as thermal energy (colloquially referred to as “heat”). More particularly, the present disclosure is directed at estimating a transient diffusion potential obeying linear diffusion; that is, a diffusion potential that is varying with time and whose gradient is a force that is linearly proportional to the flux of the underlying diffusive property.
BACKGROUNDLinear diffusion under various names is used to describe the flux of different diffusive properties. For example, Fick's law describes diffusion of a number of molecules of one gas species within a mixture of several species (the potential is called the concentration or partial pressure); Ohm's law describes the diffusion of electric charge (the electric flux is called electric current, the electric potential is called the voltage); and Fourier's law describes the diffusion of heat in conductors (the thermal potential is called the temperature). The constant of proportionality between the gradient and the flux is called the conductance, e.g. thermal conductance in Fourier's law and electrical conductance in Ohm's law, or diffusivity in Fick's law.
In particular, the effects of heat and its diffusion are becoming increasingly problematic when manufacturing and using integrated circuits (ICs). The dimensions of ICs manufactured using complementary metal-oxide semiconductor (CMOS) technology continue to shrink, which increases their power density. As ICs shrink, their power density tends to increase for two reasons. First, the ICs generally shrink at a rate that is faster than the rate at which their supply voltages decrease. Second, the frequencies at which ICs are operated tend to increase as they shrink, resulting in increasing power losses related to high frequency switching.
Research and development accordingly continue into techniques to address the problems that heat and its diffusion pose to ICs. More generally, research and development continue into ways to more generally apply theories describing linear diffusion to solve various practical problems.
SUMMARYAccording to a first aspect, there is provided a method for estimating a diffusion potential of a diffusive property. The method comprises modeling, as a circuit, diffusive behavior of a resistor, wherein a portion of the circuit models the steady-state diffusive behavior of the resistor and another portion of the circuit models transient diffusive behavior of the resistor using a capacitive circuit element representing the diffusive property capacity of the resistor, and wherein the diffusive property diffuses linearly within the resistor and generation or loss of the diffusive property is distributed within the resistor; and estimating the diffusion potential at the endpoints of the resistor by simulating operation of the circuit and determining the circuit potential at nodes in the circuit that correspond to the endpoints of the resistor.
Estimating the diffusion potential at a location within the resistor may be done using the diffusion potential at the endpoints of the resistor.
The circuit may be an electric circuit and the circuit potential may be voltage.
The diffusive property may be transient heat and the diffusive property capacity may be heat capacity.
The capacitive circuit element may be selected such that the heat capacity of the resistor is scaled by a factor representing a quasistatic temperature trajectory between the two nodes of the resistor.
The resistor may be a via and the factor may be approximately 0.5.
The resistor may have Glat≠0, a temperature distribution
and the factor may be
According to another aspect, there is provided a method for estimating a transient diffusion potential of a diffusive property that comprises modeling, as a circuit, diffusive behavior of a diffusion region, wherein one portion of the circuit models the steady-state diffusive behavior of the diffusion region and another portion of the circuit models transient diffusive behavior of the diffusion region using a capacitive circuit element representing the diffusive property capacity of the diffusion region, and wherein the diffusive property diffuses linearly within the diffusion region and generation or loss of the diffusive property is distributed within the diffusion region; and simulating operation of the circuit to estimate the transient diffusion potential at a location in the diffusion region by determining circuit potential at a node in the circuit that corresponds to the location in the diffusion region.
The circuit used to model the diffusive behavior may be an electric circuit.
The portion of the circuit that models the steady-state behavior of the diffusion region may comprise a pi-network.
The capacitive circuit element may comprise two capacitors and one of the capacitors may be connected in parallel across each of the ports of the pi-network.
The diffusion region may be a resistor, the diffusion potential may be temperature, voltage of the circuit may model the temperature, charge of the circuit may model heat, and the diffusive property capacity may be heat capacity.
The diffusion region may be an IC interconnect.
The resistor may be adjacent an environment having a uniform reference temperature, and lateral thermal conductance of the resistor relative to the environment may be zero.
The resistor may be adjacent an environment having a uniform reference temperature, lateral thermal conductance of the resistor relative to the environment may be non-zero, and modeling the temperature of the resistor may comprise assuming that temperature distribution at any moment within the resistor is quasistatic.
The resistor may comprise a single electrical net.
Temperature along the length of the resistor may be given by
and the average temperature of the resistor may be given by
where ξ=√(Glat/Glong), Glat is lateral conductance of the resistor, Glong is longitudinal conductance of the resistor, L is length of the resistor, T0 is temperature at one end of the resistor, TL is temperature at another end of the resistor, and Tref is environment reference temperature.
Alternatively, temperature along the resistor may be given by
and the average temperature of the resistor may be given by
where ξ=√(Glat/Glong), T∞=Φ/Glat+Tref, Glat is lateral conductance of the resistor, Glong is longitudinal conductance of the resistor, L is length of the resistor, T0 is temperature at one end of the resistor, TL is temperature at another end of the resistor, and Tref is environment reference temperature.
According to another aspect, there is provided a system for estimating a diffusion potential of a diffusive property that comprises a controller; and a non-transitory computer readable medium communicatively coupled to the controller and having encoded thereon statements and instructions to cause the controller to perform any of the foregoing aspects of methods.
According to another aspect, there is provided a non-transitory computer readable medium having encoded thereon statements and instructions to cause a controller to perform any of the foregoing aspects of methods.
This summary does not necessarily describe the entire scope of all aspects. Other aspects, features and advantages will be apparent to those of ordinary skill in the art upon review of the following description of specific embodiments.
In the accompanying drawings, which illustrate one or more exemplary embodiments:
Directional terms such as “top”, “bottom”, “left”, “right”, “front”, and “rear” are used in the following description for the purpose of providing relative reference only, and are not intended to suggest any limitations on how any article is to be positioned during use, or to be mounted in an assembly or relative to an environment.
Electromigration in ICs can prejudice their reliability. Electromigration refers to the migration and resulting degradation of portions of metallic conductors over time in response to both current flowing through the conductors and to relatively high temperatures to which the conductors are exposed. Since current flowing through a conductor generates heat in a distributed fashion throughout the conductor, high temperatures and current are closely linked in practice. By a “distributed fashion”, in this example it is meant that heat is generated throughout the conductor, and not simply, for example, at the endpoints of the conductor.
The degree to which ICs are susceptible to electromigration depends on their layout. Accurately predicting the effect an IC's layout will have on its operating temperature can be difficult. Unfortunately, if a problem with an IC's layout is found only after tape-out, the costs in terms of both time and money of subsequently changing the layout to rectify the problem are relatively high.
Some of the embodiments discussed herein are directed at a method and system for estimating thermal potential, also referred to herein as “temperature”, of thermal energy, also referred to herein as “heat”.
In some embodiments directed at modeling temperature of electric circuits, the transient (i.e. non steady-state, or time varying) temperature at various points in a resistive element (hereinafter interchangeably referred to as a “resistor”) is determined at various times. In these embodiments, a “resistor” is not limited to a discrete, two-terminal electrical component that ideally has only a real impedance, but rather to any electrical component whose impedance has a real component. For example, and as illustrated in the examples below, a “resistor” includes vias (in both ICs and printed circuit boards [PCBs]), contacts, circuit traces on PCBs, and interconnects within ICs. These embodiments utilize a model that take into account that heat is generated in a distributed fashion along the resistor.
More generally, in alternative embodiments the diffusion potentials of different types of diffusive properties may be estimated, so long as the diffusive property is generated or lost in a distributed fashion and diffuses linearly throughout the region through which it may diffuse (“diffusion region”). For example, in one alternative embodiment a tube of length L is open at both ends and has an inner surface that is coated with a catalyst. The tube is filled with a reactive gas and sustains a chemical reaction on its interior surface that releases a gaseous product. If the tube forms part of a larger system of connected tubes, it may be convenient to model the tube as a single entity with two ends. The rate of generation of the gaseous product and its concentration, which is the diffusion potential, results from chemical activity within the tube. In this example, the underlying diffusive property is the number of particles of the gaseous product and the diffusion region is the volume within the tube. In another alternative embodiment, a long iron bar is fixed at both ends to permafrost, but is exposed to sunlight. Heat is accordingly introduced to the portion of the bar exposed to sunlight. Heat is conducted through the two ends of the bar into the bodies that anchor it, making it likely to melt the permafrost. In this embodiment, the temperature of the bar is the diffusion potential and the diffusive property is heat, and the diffusion region is the bar itself through which heat may be conducted. These embodiments utilize a model that takes into account that the diffusive property is generated or lost in a distributed fashion throughout the diffusion region, and not just at discrete locations, such as the endpoints, of the diffusion region.
Heat EquationA model for heating in a resistor 100, as shown in
The heat per unit length in a differential slice of the resistor is mCpT(x)Δx, where the length of the differential slice is Δx, T(x) is the temperature along the resistor's length, m is the mass of the resistor per unit length, and Cp is the heat capacity of the resistor. The time rate of change of the heat is then given by
where F(x) is the diffusive heat flux. The heat flux is oriented along the resistor 100's length and is related to the temperature gradient by Fourier's law of heat conduction:
where Glong is the longitudinal thermal conductance of the resistor 100.
The Joule heat generated per unit length of the resistor 100 is given by φ=Irms2R; R is the resistor 100's electrical resistance per unit length, which is assumed to be constant within the resistor 100, and Irms is the time-averaged rms current flowing through the resistor 100. The function ƒ(x) in Equation 1 models the heat flux through the sides of the resistor 100 to the environment, which is at a reference temperature Tref. An exemplary environment is a substrate on which the resistor 100 is mounted; for example, when the resistor 100 is an IC interconnect segment, the environment may be the IC's substrate 202 (shown in
ƒ(x)=GlatT(x)−Tref (3)
Building on the discussion above, the following embodiment is directed at a method and system for estimating a diffusion potential of a diffusive property that is transient, or evolving, in nature (the “transient embodiment”). For example, the transient embodiment can be used to estimate temperature at a location within an electric circuit in which the temperature is varying with time because, for example, the circuit has not been operating for sufficiently long to reach a steady-state temperature, or because the circuit includes time-varying independent current and voltage sources. The following embodiment utilizes one example of the resistor 100, an interconnect segment 200, as an exemplary diffusion region, but is applicable to other types of diffusion regions as well.
Consider a macroscopic (i.e. not a differential slice), one-dimensional body, such as the interconnect segment 200 shown in
-
- (i) a mass per unit length (kg/μm) of “m”;
- (ii) a heat capacity (J/kg·K) of “Cp”;
- (iii) an electrical resistance per unit length (Ω/μm) or “R”;
- (iv) a longitudinal thermal conductance times length L (W·μm/K) of “Glong”; and
- (v) a lateral thermal conductance to the substrate 202 at temperature Tref per unit length (W/μm·K) of “Glat”.
The values of mass “m” and heat capacity “Cp” are adjustable to better model the heat content of the segment 200 being modeled and some portion of its immediate surroundings.
Then the component heat flows into the segment 200 are as follows:
Equation (8) reflects that heat is presumed to be generated uniformly throughout the segment 200; however, in alternative embodiments (not depicted), heat (or another diffusive property) may be generated or lost non-uniformly within the segment 200. For example, resistance may vary with location along the segment 200, which would likely prevent uniform heat generation. The average temperature
When Glat≠0, a quasistatic assumption can be made such that the temperature distribution at any moment within the segment 200 follows the following steady-state form:
where ξ=√(Glat/Glong) and T∞=Φ/Glat+Tref. Generally, making a “quasistatic assumption” refers to assuming the solution to a transient problem is the solution to the corresponding steady-state (“static”) problem.
Then the average temperature at any time over the length of the segment 200 is:
When Glat=0, the temperature distribution within the segment 200 is:
And the average temperature at any time over the length of the segment 200 is:
The treatment of the heat equation on macroscopic elements, as shown in
As discussed above, when an electric circuit is used to model heat generation and subsequent diffusion, the circuit analogy for temperature is electric potential (V), and for heat is electric charge (q). The time-dependent statement of Kirchoff's Current Law (KCL) at a node j is:
The expression for the capacitive terms is:
Equation 15 is directly analogous to Equation 9 and the charging of the capacitor is akin to raising the average temperature of the segment 200. The analogy is thus
For elements obeying Equation 11:
which says that raising the average temperature of the segment 200 is tantamount to raising the temperatures at the two nodes at lengths 0 and L of the segment 200, as if there were an equivalent capacitor γ (“thermal capacitor γ”) at each of these nodes, with the value
For elements obeying Equation 13:
which leads to the equivalent thermal capacitor
A thermal capacitor is then added to the equivalent pi-network for the segment 200 as described by either Equation 10 (Glat≠0) or Equation 12 (Glat=0), as depicted in
Coefficients for the elements of the pi-network of
And for the pi-network of
Given a thermal model of the segment 200 that is completely analogous to the electrical model, solving a single net can be performed using a simulator such as SPICE. This can be done using time-varying independent sources.
ExampleAn example involving heating a single net by a current square wave (2 ns period) follows.
The interconnect 400 is heated by the current square wave with Irms=0.707 mA at 500 MHz. A continuous analysis path is chosen through the interconnect 400 that corresponds to a lumped circuit, with component values calculated from the four discrete bodies extractable from the layout. A thermal technology file based on a generic 0.18 μm CMOS semiconductor process technology is used to generate appropriate component values α, β, γ, and ψ/2. Significant parameters are provided below, which is a SPICE netlist with three nodes used for simulation:
A circuit 500 corresponding to the netlist is shown in
The thermal evolution as determined by SPICE shows the levelling-out of the temperatures to the steady-state values; this is depicted in
The temperature trajectories along the dashed path shown in
In
As alluded to above, a controller can be used to perform the time-transient method. A system 800, such as the one shown in
For example, the system 800 may be used to partially or entirely implement a method 900 for estimating a transient diffusion potential, as shown in
Once the circuit is modeled, the engineer proceeds to block 904 where he or she simulates operation of the circuit to estimate the transient diffusion potential at a location in the diffusion region by determining circuit potential at a node in the circuit that corresponds to the location in the diffusion region. An example of this is determining voltages at various nodes of the circuit 500 shown in
The depicted embodiments model diffusive behavior using electric circuits. However, in alternative embodiments (not depicted), different types of circuits may be used. For example, a hydraulic circuit may be used in which the capacitive circuit element is a hydraulic accumulator. Also, while the depicted embodiment uses heat as an example of a diffusive property, temperature as an example of a diffusion potential, and a resistor as an example of a diffusion region, as mentioned above alternative embodiments are possible in which any one or more of the diffusive property, diffusion potential, and diffusion region are different. In these alternative embodiments, the capacitive circuit element represents not heat capacity as it does in the depicted embodiments but rather the diffusive property capacity. For example, in an embodiment in which the diffusive property is the number of particles of a gaseous product, the diffusion potential is concentration of the particles, and the diffusion region is a volume within a tube, and the mass-heat capacity product m·Cp in the depicted examples is unity.
Furthermore, while the depicted embodiments focus on generation of the diffusive property (heat generation), as mentioned above in alternative embodiments (not depicted) the diffusive property may be generated or lost. For example, instead of temperature within the resistor 100 rising, the temperature within the resistor 100 may decrease when T(x)>Tref and Glat>0. As another example, when a gaseous chemical product species is produced via a catalytic reaction of a reactant species on an interior surface of a tube, while the product species is generated the reactant species is lost. In this example, if T(x) models the reactant species, Φ(x) may be negative to represent destruction of the reactant species and f(X) can represent diffusion of the reactant species through porous walls out of the tube and into the ambient.
The circuits shown in
For the sake of convenience, the example embodiments above are described as various interconnected functional blocks or distinct software modules. This is not necessary, however, and there may be cases where these functional blocks or modules are equivalently aggregated into a single logic device, program or operation with unclear boundaries. In any event, the functional blocks or software modules can be implemented by themselves, or in combination with other operations in either hardware or software.
It is contemplated that any part of any aspect or embodiment discussed in this specification can be implemented or combined with any part of any other aspect or embodiment discussed in this specification.
While particular embodiments have been described in the foregoing, it is to be understood that other embodiments are possible and are intended to be included herein. It will be clear to any person skilled in the art that modifications of and adjustments to the foregoing embodiments, not shown, are possible. The scope of the claims should not be limited by the embodiments set forth in the examples, but should be given the broadest possible interpretation consistent with the description as a whole.
Claims
1. A method for estimating a diffusion potential of a diffusive property, the method comprising:
- (a) modeling, as a circuit, diffusive behavior of a resistor, wherein a portion of the circuit models the steady-state diffusive behavior of the resistor and another portion of the circuit models transient diffusive behavior of the resistor using a capacitive circuit element representing the diffusive property capacity of the resistor, and wherein the diffusive property diffuses linearly within the resistor and generation or loss of the diffusive property is distributed within the resistor; and
- (b) estimating the diffusion potential at the endpoints of the resistor by simulating operation of the circuit and determining the circuit potential at nodes in the circuit that correspond to the endpoints of the resistor.
2. The method of claim 1 further comprising estimating the diffusion potential at a location within the resistor from the diffusion potential at the endpoints of the resistor.
3. The method of claim 2 wherein the circuit is an electric circuit and the circuit potential is voltage.
4. The method of claim 2 wherein the diffusive property is transient heat and the diffusive property capacity is heat capacity.
5. The method of claim 4 wherein the capacitive circuit element is selected such that the heat capacity of the resistor is scaled by a factor representing a quasistatic temperature trajectory between the two nodes of the resistor.
6. The method of claim 5 wherein the resistor is a via and the factor is approximately 0.5.
7. The method of claim 5 wherein the resistor has: T ( x ) = T L - T ∞ - ( T 0 - T ∞ ) cosh ( ξ L ) sinh ( ξ L ) sinh ( ξ x ) + ( T 0 - T ∞ ) cosh ( ξ x ) + T ∞; and cosh ( ξ L ) - 1 ξ L sinh ( ξ L ).
- (a) Glat≠0;
- (b) a temperature distribution
- (c) the factor is
8. A method for estimating a transient diffusion potential of a diffusive property, the method comprising:
- (a) modeling, as a circuit, diffusive behavior of a diffusion region, wherein one portion of the circuit models the steady-state diffusive behavior of the diffusion region and another portion of the circuit models transient diffusive behavior of the diffusion region using a capacitive circuit element representing the diffusive property capacity of the diffusion region, and wherein the diffusive property diffuses linearly within the diffusion region and generation or loss of the diffusive property is distributed within the diffusion region; and
- (b) simulating operation of the circuit to estimate the transient diffusion potential at a location in the diffusion region by determining circuit potential at a node in the circuit that corresponds to the location in the diffusion region.
9. The method of claim 8 wherein the circuit used to model the diffusive behavior is an electric circuit.
10. The method of claim 9 wherein the portion of the circuit that models the steady-state behavior of the diffusion region comprises a pi-network.
11. The method of claim 10 wherein the capacitive circuit element comprises two capacitors and wherein one of the capacitors is connected in parallel across each of the ports of the pi-network.
12. The method of claim 11 wherein the diffusion region is a resistor, the diffusion potential is temperature, voltage of the circuit models the temperature, charge of the circuit models heat, and the diffusive property capacity is heat capacity.
13. The method of claim 12 wherein the diffusion region is an IC interconnect.
14. The method of claim 12 wherein the resistor is adjacent an environment having a uniform reference temperature, and wherein lateral thermal conductance of the resistor relative to the environment is zero.
15. The method of claim 12 wherein the resistor is adjacent an environment having a uniform reference temperature, wherein lateral thermal conductance of the resistor relative to the environment is non-zero, and wherein modeling the temperature of the resistor comprises assuming that temperature distribution at any moment within the resistor is quasistatic.
16. The method of claim 12 wherein the resistor comprises a single electrical net.
17. The method of claim 12 wherein temperature along the length of the resistor is given by T ( x ) = T L - T ∞ - ( T 0 - T ∞ ) cosh ( ξ L ) sinh ( ξ L ) sinh ( ξ x ) + ( T 0 - T ∞ ) cosh ( ξ x ) + T ∞ and average temperature is given by T _ = ( T 0 + T L ) cosh ( ξ L ) - 1 ξ L sinh ( ξ L ) + T ∞ ( 1 - 2 cosh ( ξ L ) - 1 ξ L sinh ( ξ L ) ), wherein ξ=√(Glat/Glong), Glat is lateral conductance of the resistor, Glong is longitudinal conductance of the resistor, L is length of the resistor, T0 is temperature at one end of the resistor, TL is temperature at another end of the resistor, and Tref is environment reference temperature.
18. The method of claim 12 wherein temperature along the resistor is given by T ( x ) = I 2 Rx 2 2 G long + ( T L - T 0 L + 1 2 · Φ L G long ) x + T 0 and average temperature is given by T _ = T L + T 0 2 + Φ L 2 12 G long, wherein ξ=√(Glat/Glong), T∞=Φ/Glat+Tref, Glat is lateral conductance of the resistor, Glong is longitudinal conductance of the resistor, L is length of the resistor, T0 is temperature at one end of the resistor, TL is temperature at another end of the resistor, and Tref is environment reference temperature.
19. A system for estimating a diffusion potential of a diffusive property, the system comprising:
- (a) a controller; and
- (b) a non-transitory computer readable medium communicatively coupled to the controller and having encoded thereon statements and instructions to cause the controller to perform a method comprising: (i) modeling, as a circuit, diffusive behavior of a diffusion region, wherein one portion of the circuit models the steady-state diffusive behavior of the diffusion region and another portion of the circuit models transient diffusive behavior of the diffusion region using a capacitive circuit element representing the diffusive property capacity of the diffusion region, wherein the diffusive property diffuses linearly within the diffusion region and generation of the diffusive property is distributed within the diffusion region; and (ii) simulating operation of the circuit to estimate the transient diffusion potential at a location in the diffusion region by determining circuit potential at a node in the circuit that corresponds to the location in the diffusion region.
20. A non-transitory computer readable medium having encoded thereon statements and instructions to cause a controller to perform a method comprising:
- (a) modeling, as a circuit, diffusive behavior of a diffusion region, wherein one portion of the circuit models the steady-state diffusive behavior of the diffusion region and another portion of the circuit models transient diffusive behavior of the diffusion region using a capacitive circuit element representing the diffusive property capacity of the diffusion region, wherein the diffusive property diffuses linearly within the diffusion region and generation of the diffusive property is distributed within the diffusion region; and
- (b) simulating operation of the circuit to estimate the transient diffusion potential at a location in the diffusion region by determining circuit potential at a node in the circuit that corresponds to the location in the diffusion region.
Type: Application
Filed: May 6, 2013
Publication Date: May 1, 2014
Patent Grant number: 8966417
Inventor: Andrew LABUN (Kelowna)
Application Number: 13/888,208