ISO/TTA 4:2002
(Main)Measurement of thermal conductivity of thin films on silicon substrates
Measurement of thermal conductivity of thin films on silicon substrates
TTA 4:2002 proposes a standard procedure for the three-omega method for measuring the thermal conductivity of a thin, electrically insulating film, on a substrate having a thermal conductivity significantly greater than the thermal conductivity of the film.This method is applicable to a film on a silicon substrate with the following characteristics: a) the film is electrically insulating; b) the film has a thermal conductivity that is less than one tenth the thermal conductivity of silicon; c) the film is uniform in thickness and the thickness lies in the range 0,25 to 1 micrometres. d) the maximum dimensions of the film are limited by the sizes of the preparation and measurement apparatus; e) the minimum dimensions of the film are limited by the minimum size of the circuit element that can be placed on the film surface. NOTE A specimen approximately 15 mm by 25 mm is of an appropriate size although specimens as small as 10 mm by 10 mm are usable. The method is directly applicable to films of silicon dioxide on silicon wafer substrates. The method may be applicable to insulating films on other high-thermal conductivity substrates provided that the parameters of the substrate material are substituted for the parameters of silicon used in this method and the associated computer program. The method is applicable to measurements near room temperature.
Mesurage de la conductivité thermique des films minces sur substrat de silicium
General Information
Standards Content (Sample)
TECHNOLOGY ISO/TTA 4
TRENDS
ASSESSMENT
First edition
2002-11-15
Measurement of thermal conductivity of
thin films on silicon substrates
Mesurage de la conductivité thermique des films minces sur substrat de
silicium
Reference number
ISO/TTA 4:2002(E)
©
ISO 2002
---------------------- Page: 1 ----------------------
ISO/TTA 4:2002(E)
PDF disclaimer
This PDF file may contain embedded typefaces. In accordance with Adobe's licensing policy, this file may be printed or viewed but
shall not be edited unless the typefaces which are embedded are licensed to and installed on the computer performing the editing. In
downloading this file, parties accept therein the responsibility of not infringing Adobe's licensing policy. The ISO Central Secretariat
accepts no liability in this area.
Adobe is a trademark of Adobe Systems Incorporated.
Details of the software products used to create this PDF file can be found in the General Info relative to the file; the PDF-creation
parameters were optimized for printing. Every care has been taken to ensure that the file is suitable for use by ISO member bodies. In
the unlikely event that a problem relating to it is found, please inform the Central Secretariat at the address given below.
© ISO 2002
All rights reserved. Unless otherwise specified, no part of this publication may be reproduced or utilized in any form or by any means,
electronic or mechanical, including photocopying and microfilm, without permission in writing from either ISO at the address below or
ISO's member body in the country of the requester.
ISO copyright office
Case postale 56 • CH-1211 Geneva 20
Tel. + 41 22 749 01 11
Fax + 41 22 749 09 47
E-mail copyright@iso.org
Web www.iso.org
Published in Switzerland
ii © ISO 2002 – All rights reserved
---------------------- Page: 2 ----------------------
ISO/TTA 4:2002(E)
Contents
Foreword.iv
Introduction .v
1 Scope. 1
2 Symbols . 1
3 Specimen preparation and characterization .2
4 Measurement apparatus. 5
5 Measurement procedure . 5
6 Calculations. 7
7 Uncertainty . 8
8 Test report . 8
Annex A Computer programs. 10
Annex B Various methods of measuring thin-film thermal conductivity. 14
Bibliography . 19
© ISO 2002 – All rights reserved iii
---------------------- Page: 3 ----------------------
ISO/TTA 4:2002(E)
Foreword
ISO (the International Organization for Standardization) is a worldwide federation of national standards bodies
(ISO member bodies). The work of preparing International Standards is normally carried out through ISO
technical committees. Each member body interested in a subject for which a technical committee has been
established has the right to be represented on that committee. International organizations, governmental and
non-governmental, in liaison with ISO, also take part in the work. ISO collaborates closely with the
International Electrotechnical Commission (IEC) on all matters of electrotechnical standardization.
International Standards are drafted in accordance with the rules given in the ISO/IEC Directives, Part 2.
The main task of technical committees is to prepare International Standards. Draft International Standards
adopted by the technical committees are circulated to the member bodies for voting. Publication as an
International Standard requires approval by at least 75 % of the member bodies casting a vote.
To respond to the need for global collaboration on standardization questions at early stages of technological
innovation, the ISO CounciI, following recommendations of the ISO/IEC Presidents' Advisory Board on
Technological Trends, decided to establish a new series of ISO publications named “Technology Trends
Assessments” (ISO/TTA). These publications are the results of either direct cooperation with
prestandardization organizations or ad hoc Workshops of experts concerned with standardization needs and
trends in emerging fields.
Technology Trends Assessments are thus the result of prestandardization work or research. As a condition of
publication by ISO, ISO/TTAs shall not conflict with existing International Standards or draft International
Standards (DIS), but shall contain information that would normally form the basis of standardization. ISO has
decided to publish such documents to promote the harmonization of the objectives of ongoing
prestandardization work with those of new initiatives in the Research and Development environment. It is
intended that these publications will contribute towards rationalization of technological choice prior to market
entry. Whilst ISO/TTAs are not Standards, it is intended that they will be able to be used as a basis for
standards development in the future by the various existing standards agencies.
Attention is drawn to the possibility that some of the elements of this document may be the subject of patent
rights. ISO shall not be held responsible for identifying any or all such patent rights.
VAMAS, the Versailles project on Advanced Materials and Standards, was formed in 1982 at the Economic
Summit of the G7 Heads of State as a way to support world trade in products dependent on advanced
materials technologies by providing the technical basis for harmonized measurements, testing, specifications,
and standards. Through its technical programmes, VAMAS emphasizes international collaborative pre-
standards measurement research in support of development of national and international standards by
standards development organizations such as ISO.
ISO/TTA 4 was prepared by VAMAS and published under a memorandum of understanding concluded
between ISO and VAMAS.
iv © ISO 2002 – All rights reserved
---------------------- Page: 4 ----------------------
ISO/TTA 4:2002(E)
Introduction
The purpose of this document is to propose a standard procedure for measuring the thermal conductivity of
insulating thin films on silicon substrates. Based on a recent interlaboratory comparison, a recommendation is
made for the adoption of the three-omega method as a standard measurement method. A procedure for the
three-omega method is proposed for measuring the thermal conductivity of a thin, electrically insulating film,
on a substrate having a thermal conductivity significantly greater than the thermal conductivity of the film.
Annex B contains a review of several measurement methods that have been used to measure the thermal
conductivity of such films (see reference [1]).
© ISO 2002 – All rights reserved v
---------------------- Page: 5 ----------------------
TECHNOLOGY TRENDS ASSESSMENT ISO/TTA 4:2002(E)
Measurement of thermal conductivity of thin films on silicon
substrates
1 Scope
1.1 A standard procedure for the three-omega method is proposed for measuring the thermal conductivity of
a thin, electrically insulating film, on a substrate having a thermal conductivity significantly greater than the
thermal conductivity of the film. This method is applicable to a film on a silicon substrate with the following
characteristics:
a) the film is electrically insulating;
b) the film has a thermal conductivity that is less than one tenth the thermal conductivity of silicon;
c) the film is uniform in thickness and the thickness lies in the range 0,25 µm to 1 µm;
d) the maximum dimensions of the film are limited by the sizes of the preparation and measurement
apparatus;
e) the minimum dimensions of the film are limited by the minimum size of the circuit element that can be
placed on the film surface.
NOTE A specimen approximately 15 mm by 25 mm is of an appropriate size although specimens as small as
10 mm × 10 mm are usable.
1.2 The method is directly applicable to films of silicon dioxide on silicon wafer substrates.
1.3 The method may be applicable to insulating films on other high-thermal conductivity substrates provided
that the parameters of the substrate material are substituted for the parameters of silicon used in this method
and the associated computer program.
1.4 The method is applicable to measurements near room temperature.
2 Symbols
See Figure 1.
f frequency of excitation voltage
ω angular frequency of excitation voltage = 2 πf
Λ thermal conductivity of substrate
Λ thermal conductivity of film
f
w width of specimen resistor
t film thickness
L length of specimen resistor
V excitation voltage at ω
0
V voltage at ω across the specimen
R mean resistance of the specimen
© ISO 2002 – All rights reserved 1
---------------------- Page: 6 ----------------------
ISO/TTA 4:2002(E)
∆R variation at 2ω of the specimen resistance
V voltage at ω across offset resistor
OS
R resistance of offset resistor
OS
V voltage at ω across calibration resistor
CAL
R resistance of calibration resistor = 10 Ω
CAL
V voltage at 3ω across specimen
3ω
T temperature setting
SP
T mean temperature of the specimen, also known as the measurement temperature
∆T temperature variation of specimen at 2ω
∆T bare substrate thermal signal
b
I electrical current at ω = V /(R + R + R )
0 REF CAL
2
P power dissipated in specimen resistor = I R
dR/dT temperature coefficient of resistance
Figure 1 — Nomenclature
3 Specimen preparation and characterization
3.1 Apparatus
3.1.1 Photomask, a mask for photolithography used in manufacture of integrated electronic circuits. See
Figure 2 for the recommended configuration. It is recommended that the photomask contain several
duplicates of the desired circuit elements in case a deposited element is not usable.
3.1.2 Clean room.
3.1.3 Spin coater.
3.1.4 Cleaning solution.
3.1.5 Deionized water.
2 © ISO 2002 – All rights reserved
---------------------- Page: 7 ----------------------
ISO/TTA 4:2002(E)
3.1.6 Photoresist.
3.1.7 Photoresist developer.
3.1.8 Acetone.
3.1.9 Aluminium etch.
3.1.10 Vacuum deposition chamber, equipped for deposition of aluminium and set up to hold the
specimen.
3.1.11 Ultraviolet exposure system.
3.1.12 Baking ovens.
3.1.13 Plasma etcher.
3.1.14 Assorted plastic cups and tweezers.
3.1.15 Optical comparitor for line width measurement.
3.1.16 Ellipsometer for thickness measurement.
3.1.17 Electrically conducting silver paste.
3.2 Determination
Measure the film thickness using the ellipsometer and record the value as t in Table 1.
3.3 Circuit element preparation
3.3.1 Using the cleaning solution, clean the specimen using an established cleaning procedure.
3.3.2 Place the specimen in the vacuum deposition chamber and deposit an aluminium film approximately
300 nm thick onto the film surface of the specimen. Remove the specimen from the chamber.
3.3.3 Spin coat the metallized surface of the specimen with photoresist.
3.3.4 Bake the specimen in an oven at 95 °C for 25 min.
3.3.5 Soak the specimen in water for 2 min.
3.3.6 Mount the specimen and the photomask in the ultraviolet exposure system.
3.3.7 Expose the specimen to ultraviolet radiation through the photomask for the time appropriate to the
exposure system.
3.3.8 Develop the specimen in the photoresist developer diluted with water (volume ratio of photoresist
developer to water = 1:4) for 30 to 45 s.
3.3.9 Bake the specimen in an oven at 125 °C for 25 min.
3.3.10 Plasma etch the specimen in oxygen for 45 s to remove any scum on the specimen surface.
3.3.11 Etch the specimen in the aluminium etchant at 50 °C for 30 to 45 s.
3.3.12 Soak in acetone for about 1 min to remove remaining photoresist.
3.3.13 Rinse in deionized water and dry in air.
© ISO 2002 – All rights reserved 3
---------------------- Page: 8 ----------------------
ISO/TTA 4:2002(E)
3.3.14 With the optical comparitor, measure the line width of the circuit element produced. Record this value
as w in Table 1.
Use a blank table, Table 1 displays a set of representative data.
NOTE The line width will usually differ significantly from the line width of the photomask.
3.3.15 Record in Table 1 the length of the metal strip between the two voltage pads shown in the photomask
design as L.
NOTE The value for L used in the photomask design in Figure 2 is satisfactory for the purposes of this method.
3.3.16 Using the silver paste, attach fine copper wires to the four electrical pads of the circuit element chosen
to be used in the experiment.
3.3.17 Check for electrical continuity in the circuit element with a volt-ohm-resistance meter between all pairs
of electrical pads. If an open circuit condition exists or a short circuit condition exists, choose a different circuit
element and repeat this step in the procedure with this new circuit element. If none of the circuit elements is
satisfactory, reject the specimen as unacceptable for measurement.
The resistance between the voltage pads (the two left pads of a circuit element shown in Figure 2) should be
between 10 Ω and 100 Ω.
Dimensions in millimetres
Figure 2 — Recommended photomask pattern for 3-Omega Method
4 © ISO 2002 – All rights reserved
---------------------- Page: 9 ----------------------
ISO/TTA 4:2002(E)
4 Measurement apparatus
4.1 A variable temperature environmental chamber to hold the specimen. The chamber specifications are:
4.1.1 Temperature adjustable to 20 °C and 60 °C with an uncertainty of 0,1 °C or less.
4.1.2 A thermocouple to measure the temperature at the position of the specimen and a thermocouple
readout. The uncertainty in the temperature measurement should be 0,1 °C or less.
4.1.3 Four electrical feed-throughs to accommodate two current leads and two voltage leads to the
specimen.
4.2 An electrical measurement system as shown in Figure 3 consisting of the following components:
4.2.1 A digital lock-in amplifier for performing the electrical measurements. The specifications are:
4.2.1.1 A signal generator that generates a reference signal having an output voltage adjustable to 0,1 V and
to 2,0 V at the frequency of 333 Hz.
4.2.1.2 Ability to measure single-ended input voltages and differential input voltages.
4.2.1.3 Settings capable of measuring signals at the fundamental frequency (f) and at the third harmonic
(3f ).
4.2.1.4 A voltage measurement uncertainty of 0,1 % or less.
4.2.2 A switch box for switching input A between measuring V (position y) and measuring V or V
CAL 3ω
(position x).
4.2.3 Additional components as shown in Figure 3.
5 Measurement procedure
5.1 Mount the specimen in the environmental chamber making all of the necessary electrical connections.
5.2 Set T to 20 °C and wait until the temperature stabilizes. Read the temperature on the thermocouple
SP
readout and record the value in Table 1 as T (20) to within 0,1 °C.
SP
5.3 Set the lock-in frequency to 333 Hz. Set the time constant to 3 s and the roll-off to 24 db/oct. Set the
lock-in to read the in-phase signal.
5.4 Set V to 0,1 V.
0
5.5 Set R to an intermediate position.
NULL
5.6 Wait 20 min for specimen temperature to stabilize.
5.7 Set the switch to position x. Set the lock-in to read the differential input voltage.
5.8 Adjust R until approximately zero voltage is measured. Maximize the gain on the lock-in and adjust
OS
R until a minimum voltage is obtained.
NULL
5.9 Lower the gain setting of the lock-in amplifier to the minimum setting. Switch the lock-in to single ended
input.
5.10 Set the lock-in gain to a convenient range. Use the autogain feature if the lock-in is so equipped. Read
V and record this value in Table 1, col. 1.
© ISO 2002 – All rights reserved 5
---------------------- Page: 10 ----------------------
ISO/TTA 4:2002(E)
Key
OA = two operational amplifiers (AD524), set for unity gain
G = earth (ground)
A, B = voltage inputs
R = nulling resistor consisting of a 10 turn potentiometer, resistance 10 kΩ, in series with a fixed 50 kΩ resistor
NULL
that is connected to ground
R = offset resistor consisting of a 10 turn potentiometer, resistance 5 kΩ
OS
R = calibration resistor, fixed resistance of 10,0 ± 0,1 Ω
CAL
SW = single-pole double-throw switch
Figure 3 — Electrical circuit
6 © ISO 2002 – All rights reserved
---------------------- Page: 11 ----------------------
ISO/TTA 4:2002(E)
5.11 Set the switch to position y. Measure V and record this value in Table 1, col. 1.
CAL
5.12 Set V to 2,0 V.
0
5.13 Wait 20 min for specimen temperature to stabilize.
5.14 Repeat steps 5.6 to 5.10, except enter the values of V and V in Table 1, col. 2.
CAL
5.15 Set the lock-in to read the voltage at 3ω .
5.16 Measure V and enter the value in Table 1, col. 2.
3ω
5.17 Set T to 60 °C and wait for the temperature to stabilize. Record the temperature value in Table 1 as
SP
T (60) to within 0,1 °C.
SP
5.18 Repeat steps 5.2 through 5.10 except record the indicated values in Table 1, col. 3.
5.19 Repeat steps 5.11 through 5.14 except record the indicated values in Table 1, col. 4.
5.20 Repeat steps 5.15 and 5.16 except record the respective values of V in Table 1, col. 4.
3ω
6 Calculations
6.1 For each column of Table 1, calculate I, P and R using the formulae
IV= R
CAL
CAL
R=VI
2
PI= R
Insert the calculated values in the corresponding locations in Table 1.
6.2 For the initial temperature, extrapolate R to zero input power using the formula
R −R
21
RR20=−P
()
21
PP−
21
Insert the value of R(20) into Table 1.
6.3 For the second temperature, extrapolate R to zero input power using the formula
R −R
43
RR60=−P
()
44
PP−
43
Insert the value of R(60) into Table 1.
6.4 Calculate the temperature of measurement for col. 2 using the formula
TT60 − 20
( ) ( )
SP SP
TT=+20 R−R 20
() ()
2SP 2
RR60 − 20
() ( )
and insert the value into Table 1.
6.5 Calculate the temperature coefficient of the specimen resistance from the formula below and record the
value in Table 1.
RR60 − 20
( ) ( )
dR
=
d6TT 0 −T20
() ()
SP SP
© ISO 2002 – All rights reserved 7
---------------------- Page: 12 ----------------------
ISO/TTA 4:2002(E)
6.6 Compute ∆T .
b
6.6.1 Using a text editor, generate a data input file with the name “input.dat” using the format of A.2 by
substituting the appropriate data from Table 1, col. 2.
6.6.2 Using the executable version of the program in A.1, compute the bare substrate thermal signal, ∆T ,
b
and record the value in Table 1.
NOTE The program contains formulae based on handbook data for the temperature dependence of the thermal
conductivity and thermal diffusivity of Si. If another substrate material is desired, modify the appropriate lines in the
program. The formulae are located in lines 90 and 96 of program OMEGA3.
6.7 Calculate the total thermal signal using the formula below and record the value in Table 1.
dTR
∆=TV2
3ω
dRV
−9 2 −1
6.8 Calculate the thermal resistance of the thin film (units, 10 m K·W ) using the formula below and record
the value in Table 1.
(∆−TT∆ )
b
−9
RL=−w 23×10
T
P
–9
NOTE The value 23 × 10 is the estimated thermal resistance due to interfaces in the specimen and is the value
derived in reference [1].
6.9 Calculate the thermal conductivity of the film using the formula below and record the value in Table 1.
t
Λ=
R
T
7 Uncertainty
The standard uncertainty of a measurement is estimated to be ± 10 %. A reanalysis of round robin data by the
round robin participants would permit a better determination of the uncertainty. The deviation of the sample
−1 −1
data shown in Table 1, obtained by one laboratory, gives a value 1,28 W⋅m K , which deviates from the
−1 −1 −1 −1
expected bulk value of 1,37 W⋅m K by 0,09 W⋅m K or 6,7 %.
8 Test report
The report shall contain the following information:
a) specimen identification;
b) substrate material;
c) film material;
d) film thickness;
e) temperature of measurement;
f) film thermal conductivity.
8 © ISO 2002 – All rights reserved
---------------------- Page: 13 ----------------------
ISO/TTA 4:2002(E)
Table 1 — Example of representative data
Col. 1 Col. 2 P = 0 Col. 3 Col. 4 P = 0
−3
L (m)
4,00 × 10
−6
w (m)
28,45 × 10
−6
t (m)
0,488 × 10
f (Hz)
332,6
T (°C) T (20) = 19,2 T (60) = 60,0
SP SP SP
−2 −1 −2 −1
V (V)
7,126 × 10 1,778 × 10 6,732 × 10 1,685 × 10
CAL
−1 −1 −1 −1
V (V)
1,925 × 10 4,824 × 10 2,065 × 10 5,188 × 10
−3 −2
I (A)
7,126 × 10 1,778 × 10
R = 27,01 R = 27,13 R(20) = 27,27 R = 30,67 R = 30,79 R(60) = 30,93
R (Ω)
1 2 3 4
−3 −3 −3 −3
P (W)
1,372 × 10 8,58 × 10 1,390 × 10 8,74 × 10
T = 20,8
T (°C)
2
−5
V (V)
3,59 × 10
3ω
−1 −2
dR/dT (Ω⋅K ) 8,96 × 10
−2
∆T (K) 1,46 × 10
b
−2
∆T (K)
4,51 × 10
2 −1 −7
R (m KW ) 3,82 × 10
T
−1 −1
Λ (W⋅m K ) 1,28
© ISO 2002 – All rights reserved 9
---------------------- Page: 14 ----------------------
ISO/TTA 4:2002(E)
Annex A
Computer programs
A.1 Fortran file for computation of ∆T
b
program omega3
c version 1, 01-Nov-2001
c Calculate real temperature variation at a given frequency
c for a silicon plate heated by a thin metal strip.
c The calculation can be applied to plates composed of
c material other than silicon if the thermal conductivity
c and thermal diffusivity of that material are substituted
c in the indicated lines below.
c This is used for analyzing data obtained by the 3 omega
c method
c
c library routine 'dqagi' used for integration
c --------------------------------------------------------
c A data file named 'input.par' written in ASCII is needed
c
c A sample file is given in section A2. Just substitute data
c in the appropriate lines. Double precision is used.
c
c Function Fcomplex is written for a two layer system. The
c main program defines one of these layers as having zero
c thickness and having the same material parameters as
c the other layer.
c --------------------------------------------------------
c Nomenclature (not necessarily the same as in the standard)
c D =diffusivity
c k =conductivity
c L =substrate thickness
c lnth=heater length
c b =heater half-width
c b2 =heater full width
c w =angular frequency at 2 omega
c f =fundamental frequency
c f2 =frequency at 2 omega
c q =power input to specimen
c T =measurement temperature
c tr =substrate temperature signal at 2 omega
c
c --------------------------------------------------------
c define integration parameters for dqagi
double precision, external :: fr
c
double precision, parameter :: bound=0.d0
double precision, parameter :: epsabs=1.d-6
double precision, parameter :: epsrel=1.d-12
integer, parameter :: inf=1
integer, parameter :: limit=1000
integer, parameter :: lenw=limit*4
integer, parameter :: idim=2
10 © ISO 2002 – All rights reserved
---------------------- Page: 15 ----------------------
ISO/TTA 4:2002(E)
c
double precision result, abserr, work(lenw), T
double precision f, f2, b2
integer neval, ier, last, iwork(limit)
character*30, label
c --------------------------------------------------------
c common variables
double precision w, b
double precision d(idim), k(idim), L(idim)
common w, d, k, L, b
c
double precision q, pi
double precision tr, lnth
c
1 format(A12)
2 format(A30)
3 format(6(1x,d13.6))
4 format(1x,A33,D12.4)
c --------------------------------------------------------
pi=4.d0*atan(1.d0)
c --------------------------------------------------------
open(1, 'input.par')
c
read (1,2)
read(1,2) label ; read(1,*) f
read(1,2) label ; read(1,*) q
read(1,2) label ; read(1,*) b2
read(1,2) label ; read(1,*) lnth
read(1,2) label ; read(1,*) L(1)
read(1,2) label ; read(1,*) T
c
f2=2*f
w=2.d0*pi*f2
b=b2/2.d0
c +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
c The thermal conductivity of Si as a function of
c temperature is computed in the line 90 below. If another
c material is wanted, replace line 90
90 k(1)=1.685d0-8.73d-3*T+3.62d-5*T*T-9.0d-8*T*T*T
c The thermal diffusivity of Si as a function of the
c thermal conductivity of Si is computed in line 96 below.
c If another material is wanted, replace line 96.
96 D(1)=0.093d0+0.268d0*k(1)+0.180d0*k(1)*k(1)
c ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
write(*,4)'thermal conductivity of Si, W/cm/K ', k(1)
write(*,4)'thermal diffusivity of Si, cm2/s ', D(1)
k(2)=k(1)
D(2)=D(1)
L(2)=0.d0
c
c --------------------------------------------------------
c Calculate the temperature variation
c
c calculate average temperature (real and imaginary) using
c dqagi
c
© ISO 2002 – All rights reserved 11
---------------------- Page: 16 ----------------------
ISO/TTA 4:2002(E)
call dqagi(fr, bound, inf, epsabs, epsrel,
*result, abserr, neval, ier, limit, lenw, last, iwork, work)
if (ier.gt.0) then
call list(bound, inf, epsabs, epsrel,
* result, abserr, neval, ier, limit, lenw, last, 'fr')
end if
c
tr=result*q/pi/lnth
c
write(*,4) 'substrate thermal signal ', tr
c --------------------------------------------------------
close(2)
end
c
c *******************************************************
function fr(x)
c real part of integrand for dqagi
double precision x, fr
complex*8 fcomplex
c
fr=real(fcomplex(x))
return
end function fr
c
c *******************************************************
function fcomplex(x)
c
c complex integrand for 2 layers in vacuum
c
c common variables
double precision w, b
double precision d(2), k(2), L(2)
common w, d, k, L, b
c
double precision x
complex*8 u(2), ea1, ea2
complex*8 ci, gp, gm, ex1, ex2, bp, bm, ftr, fcomplex
c
ci=(0,1.d0)
u(1)=sqrt(x*x-ci*w/d(1))
u(2)=sqrt(x*x-ci*w/d(2))
gp=u(1)*k(1)+u(2)*k(2)
gm=u(1)*k(1)-u(2)*k(2)
ea1= 2.d0*u(1)*L(1)
if (real(ea1).gt.160) then
ex1=0.d0
else
ex1= exp(-ea1)
end if
ea2= 2.d0*u(2)*L(2)
if (real(ea2).gt.160) then
ex2=0.d0
else
ex2= exp(-ea2)
end if
bp=(gp*ex2 + gm)*ex1
bm=(gm*ex2 + gp)
ftr=(bm+bp)/(bm-bp)/u(1)/k(1)
if(x.eq.0.d0) then
fcomplex=ftr
12 © ISO 2002 – All rights reserved
---------------------- Page: 17 ----------------------
ISO/TTA 4:2002(E)
else
fcomplex=ftr*sin(x*b)*sin(x*b)/x/x/b/b
end if
c print*,'end of fcomplex'
return
end function fcomplex
c
c *******************************************************
subroutine list(bound, inf, epsabs, epsrel,
* result, abserr, neval, ier, limit, lenw, last, ri)
c
double precision bound, epsabs, epsrel, result, abserr
integer inf, limit, lenw, neval, ier, last
character*2 ri
c
print*,'list parameters from dqagi call to ', ri
print*, 'epsabs=', epsabs
print*, 'epsrel=', epsrel
print*, ' limit=', limit
print*, 'result=', result
print*, ' neval=', neval
print*, ' ier=', ier
print*, ' last=', last
print*, ' '
return
end subroutine list
c *******************************************************
A.2 Data input file containing representative data
The data are written in double precision. Substitute appropriate values from Table 1.
'file input.par'
'fundamental frequency, Hz'
332.6d0
'power, W'
8.58d-3
'line full width, cm'
28.45d-4
'line length, cm'
0.4d0
'thickness Si, cm'
0.038d0
'measurement temp. T2, deg C'
20.8d0
A.3 Library routine dqagi.for
National Institute of Standards and Technology (NIST) Guide to Avail
...
Questions, Comments and Discussion
Ask us and Technical Secretary will try to provide an answer. You can facilitate discussion about the standard in here.