Arxiv:1707.04482V3 [Physics.Ed-Ph] 19 Apr 2019

Determination of Boltzmann constant by equipartition theorem for capacitors

T. M. Mishonov∗ Laboratory for Measurements of Fundamental Constants, Faculty of Physics, St. Clement of Ohrid University at Soﬁa, 5 James Bourchier Blvd., BG-1164 Soﬁa, Bulgaria

V. N. Gourev Department of Atomic Physics, Faculty of Physics, St. Clement of Ohrid University at Soﬁa, 5 James Bourchier Blvd., BG-1164 Soﬁa, Bulgaria

I. M. Dimitrova Faculty of Chemical Technologies, University of Chemical Technology and Metallurgy, 8 Kliment Ohridski Blvd., BG-1756 Soﬁa, Bulgaria

N. S. Serafimov Department of Telecommunications, Technical University Sofia, 8 Kliment Ohridski Blvd., BG-1000 Sofia, Bulgaria

A. A. Stefanov Faculty of Mathematics, St. Clement of Ohrid University at Sofia, 5 James Bourchier Blvd., BG-1164 Sofia, Bulgaria and Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. Georgi Bonchev Str., Block 8, 1113 Sofia, Bulgaria

E. G. Petkov, A. M. Varonov† Faculty of Physics, St. Clement of Ohrid University at Sofia, 5 James Bourchier Blvd., BG-1164 Sofia, Bulgaria (Dated: 25 Dec 2018) A new experimental set-up for Boltzmann constant measurement is described. Statistically averaged square of voltage U 2 is measured for different capacitances C. Boltzmann constant is determined by the equipartition 2 theorem C U = kB T . For fixed capacitance, voltages could be measured for different temperatures. The set-up consists of low-noise high frequency operational amplifiers ADA4898-2. An instrumental amplifier is followed by an inverting amplifier, square of the voltage is created by an analog multiplier AD633 and finally the averaged signal is measured by a multimeter. More than 10 high-school students were able to measure the Boltzmann constant with the experimental set-up in the 5th Experimental Physics Olympiad with excellent accuracy compared to the price, conditions and available time for the experiment. A new derivation of the important for the statistical physics theorems by Nyquist and Callen-Welton is given in an appendix in the level of introductory courses in physics listened by teachers. For the understanding of the work of the experimental set-up, it is necessary to know only the equipartition theorem.

I. INTRODUCTION mass by C and velocity by U. Now the equipartition theorem is included in all high school textbooks. Perhaps the most The equipartition theorem1,2 describing the relation be- famous example is the mean kinetic energy of a single atom tween the mean thermal energy of a quadratic degree of free- 1 3 dom mv2 = k T, v2 = v2 + v2 + v2. (1) arXiv:1707.04482v3 [physics.ed-ph] 19 Apr 2019 2 2 B x y z 1 1 1 D E 1 1 mv2 = kx2 = Iω˜ 2 = CU 2 =...= k T 2 x 2 2 2 2 B Here it is implicitly assumed that hvi = 0; no wind in the room. Analogously, the mean voltage of the connected resis- is one of the ﬁrst quantitative results of statistical physics. The tor and capacitor shown in Fig. 1 is zero hUi = 0.

coefﬁcient kB in front of the temperature is the Boltzmann In the beginning of the development of statistical physics constant. The physical nature of the variables is irrelevant: m 3 Albert Einstein suggested that Boltzmann constant kB can be can be the mass of a molecule and vx is the x-component of determined using the equipartition theorem for and also time the velocity, k can be elastic constant of a spring and x can averaged thermally averaged square of the voltage U 2 of a be the deformation, for a torsion magnetometer x could be capacitor C the angle, I˜ can be the moment of inertia of a molecule and 2 ω is the angular velocity. It is not a thoughtcrime to denote C U = kB T. (2) 2

mechanics. In the DPG conference in Halle 1891 Boltzmann presented his famous formula for the statistical interpretation of the entropy S and the number of quantum states (in contemporary terminology) with one and the same energy. Asked provocatively “Do you consider that ...” , unexpectedly for the auditorium, Boltzmann replied “I do not exclude that energy FIG. 1. Parallely connected capacitor C, resistor R and a voltmeter can have atomic structure”.11 In such a way, the birthyear of V1 measuring the time dependent thermal fluctuations of the voltage the quantum mechanics is 1891. U(t). A DC voltmeter should show hUi = 0 because electric noise The manuscript by Waterston2 was rejected from publica- does not introduce DC-voltages. tion by the polite referee report by Sir John William Lubbock: “the paper is nothing but nonsense, ...”12 Even nowadays, the The first attempt for the realization of this idea was made by mathematical physicists feel obliged to convince the physi- Habicht brothers4 in 1910 but unfortunately their electrostatic cists that they did not lead Boltzmann to suicide. In such a amplifier with mechanically rotating parts suffered from float- way, the Einstein idea to determine the Boltzmann constant ing off the zero. It is strange that for more than a century by the equipartition theorem is in a good company. The idea this simple experiment has not yet been performed. Up to was not realised in the 20th century and had to wait for the ap- now, only Johnson and Schottky spectral noise experiments pearance of low-noise high frequency operational amplifiers are conducted in student laboratories.5–9 to be realised in the students education. The building of an ex- Imagine that we have a set of capacitors with capacities perimental set-up for Boltzmann constant determination from 2 given in Table I with V2 ∝ U . scratch was a modus for the students in the University of Sofia to pass the exam in statistical physics. The purpose of the present work is to present a simple self- C [nF] V [V] 2 made set-up for determination of Boltzmann constant in high- 6.89 0.510 school and university laboratories. This is a good method- 17.20 0.363 ological illustration of the equipartition theorem. First of all, 27.80 0.327 6 63.40 0.289 the signal should be amplified million times Y ≈ 10 using 70.00 0.286 low-noise, high-frequency operational amplifiers. The prob- 93.20 0.280 lem of the floating off the zero is solved with large fast capac- 132.00 0.272 itors sequentially connected to the gain resistors of the circuit depicted in Fig. 2. Later on the signal has to be squared by TABLE I. Experimental results for different capacitors – capacity an analog multiplier and this squared signal has to be aver- and corresponding DC voltage measured by an ordinary cheap All- aged by a low-pass filter. Finally the time averaged signal V Sun DT-830B multimeter. is measured by an ordinary multimeter. The simplicity of the experiment and set-up, as well as the availability of the electronic elements nowadays mean that this experimental set-up According to the equipartition theorem at temperature 24◦C is suitable for high-school physics labs. In such a way, our for the 6.89 nF and 132 nF we have for the voltage fluctu- ultimate goal is to turn over the deteriorating trend in physics p 2 p ations δU = hU i = kB T/C correspondingly 771 nV and engineering education worldwide, and this experimental and 176 nV. If those voltages are amplified 1 million times, set-up is a valuable tool in this endeavor. In short, the set-up 771 mV and 176 mV can be measured even by cheap multi- can be used by students to measure a fundamental constant meters. After several time constants τRC = RC, the random kB , teachers can give explanations not only to physical prin- voltage U(t) is already independent and in order to measure ciples but also to explain how the set-up works. Enthusiastic the mean square teachers together with his or her students can make the set-up in one day. 1 Z T U 2 = U 2(t)dt, (3) T 0 it is necessary the averaging time T significantly to exceed p II. THEORY OF THE EXPERIMENTAL SET-UP τRC ; the accuracy is τRC /T 1. All those theorems are actually simple illustrations of the kinetic principle of the detailed balance applied to the most As already stated in the previous section, the voltage fluc- simple physical system – the harmonic oscillator and can be tuations of nF capacitors in room temperature is of the order considered as a comment to the Planck work10 from the 19th of hundred nV. These nV should be amplified 106 times and century. In order to explain the spectral density of the black squared after that in order to be able to measure them with a body radiation, Planck postulated the equidistant spectrum relatively cheap and commercially available multimeters. The En = ~ωn; it was the simplest realization of the Boltzmann circuit of the experimental set-up performing these actions is idea for atomic structure of the energy. In some sense, the dis- shown in Fig. 2. The capacitor C, whose voltage fluctuations crete energy spectrum is associated with the quantum mechan- are to be measured is connected parallely to a resistor R. This ics and it is interesting to trace the birthday of the quantum parallel circuit is connected to the two inputs of a buffer with 3

2 2 ∗ FIG. 2. A device for measurement of the voltage fluctuations δU = (U − hUi) = V2/U of parallely connected capacitor C and resistor R. In the circuit we can recognize an instrumental amplifier,13 an inverting amplifier,13 a multiplier with inputs connected in parallel,14 and ∗ an averaging low-frequency filter. The expressions of the constants y1 and U by the circuit parameters is given in the text, Eq. (4) and 15 Eq. (14) correspondingly. The voltage V2 is measured with multimeters. The circuit comprises two dual ADA4898-2 low-noise operational 16 amplifiers. The gain capacitors CG are 10 µF MKS2 manufactured by WIMA. In short, the voltmeter schematically represented in Fig. 1 is depicted in detail. The amplified, squared and averaged voltage is measured by an ordinary multimeter. As a whole the circuit can be described as a low-noise pre-amplifier followed by a true RMS-meter. amplification amplifier according to the manufacturer instructions14 and finally there is an averaging low pass filter consisting of a resis- RF tor Rav and a WIMA MKS2 type capacitor Cav. The values of y1 = 1 + 2 , (4) rG all passive elements and voltage supplies for the ADA4898-2 operational amplifiers and AD633 analog multiplier are given which is the first step of the amplifier. Next a MKS2 WIMA in Table II. type capacitor CG is connected to each of the two outputs of the buffer to stop the voltage offset of the low-noise ADA4898-2 operational amplifiers. The second step of the Circuit element Value amplifier, which is an difference amplifier (DA) with amplifi- R 510 Ω cation rG 20 Ω RF 1 kΩ 0 C 10 pF RF F y2 = − , (5) CG 10 µF RG RG 100 Ω 0 RF 10 kΩ is connected to both CG capacitors, forming together with the 0 buffer an instrumental amplifier. The output of the DA is con- CF 10 pF R1 2 kΩ nected again to a MKS2 WIMA CG capacitor with the same R2 18 kΩ reason to stop the voltage offset and after CG the third and last Rav 1.5 MΩ step, which is an inverting amplifier with amplification Cav 10 µF R Ω17 0 V 1 M R 14,15 y = − F , (6) VCC +9 V 3 14,15 RG VEE -9 V is connected. Finally, the total amplification of the amplifier TABLE II. Table of the numerical values of the circuit elements from Fig. 2. 0 2 RF RF Y = y1y2y3 = 1 + 2 . (7) rG RG

The Y times ampliﬁed signal is squared by an AD633 ana- Lastly the resistor Rav forms another voltage divider with the 17 log multiplier by connecting the ampliﬁer output to both in- internal resistance of the used DT-830B multimeter RV. For puts of the multiplier. A voltage divider is connected after the small values of Rav this voltage division is negligible, but in 4

18 our case Rav and RV are in the same order of magnitude and multiplier can be substituted with more precise AD835 (to- we have to take it into account in the final equation for the tal error 0.1% of full scale, small signal bandwidth 250 MHz). amplification. But even in this case, the combination of a commercial low- The time dependent thermal fluctuations of the voltage U(t) noise pre-amplifier and a true RMS-meter costs one or even of the parallely connected C and R is amplified Y times two orders of magnitude more than our experimental set-up that can be easily reproduced in any school physics classroom. Uamp = Y U. (8) The even more detailed derivation of all those formulae requires only Ohm law applied to a voltage divider and a se- Then the amplified voltage Uamp is squared by an analog mul- quential chain of simple problems for high-school students is tiplier. We use the circuit depicted in Ref. 14, Fig. 17, where described in great detail in Ref. 19 and it is remarkable that U (t) ≡ U X2 = Y2 = 0, X1 = Y1 = amp and W , i.e. we have even a high-school student was able to solve the university electronic problem. U 2 R + R 2 U = amp 1 2 , The substitution of U from the equipartition theorem (9) Um R1 Eq. (2) into the property of the circuit we derived Eq. (12) gives a linear dependence between the DC voltage V2 and the 14 where the constant factor for this multiplies AD633 reciprocal capacitance 1/C Um = 10 V. Actually we have a voltage divider for which Z = WR /(R + R ) and S = 0. The squared voltage is av- 1 Zk T 2 1 2 V = q + v , q ≡ B , Z ≈ 1. (15) eraged by an averaging low pass filter with large time con- 2 0 C 0 0 U ∗ stant τav = RavCav of order of a quarter a minute. The large R resistance is comparable with the internal resistance of the The irrelevant for the experiment constant v0 describes the av internal noise of the circuit determined mainly by the volt- voltmeter RV and for time averaged DC voltage shown by the voltmeter we have to take into account a second voltage di- age noise of the first dual ADA4898-2 amplifier in the buffer vider (double non-inverting amplifier with a virtual common point). In order v0 to be small, we need to use low-noise operational R amplifiers. V = hU i V , 2 (10) RV + Rav The slope q0 of the linear regression V2 versus 1/C while hV1(t)i = 0 the fluctuations of the electric noise create dV2 q0 ≡ −1 (16) non-zero square. Combined together equation for amplifying dC regr Eq. (8), squaring Eq. (9) and averaging Eq. (10) we arrive at the transfer function of the amplifier determines the Boltzmann constant ∗ 2 q0U Y R1 + R2 RV 2 k = . (17) V2 = U (t) , (11) B ZT Um R1 RV + Rav For high-school students the correction multiplier Z = where brackets mean time averaging at steady spectral density 1 − ε is indistinguishable from one, but for university stu- of the noise. This can also be written in the form dents at undergraduate level the several percent correction ε can be calculated analyzing the frequency dependence of the U 2 = U ∗V , (12) 2 amplifier as it is described in great detail in Appendix A. where the time averaged square of the investigated voltage U 2 and the constant voltage measured by a voltmeter V is III. FREQUENCY DEPENDENT CONSIDERATIONS related by a constant with dimension of voltage

2 1 Y R1 + R2 RV Up to now, the frequency dependence have been neglected. ∗ ≡ . (13) For high school and even undergraduate level this is perfectly U Um R1 RV + Rav acceptable but for a detailed analysis of the amplifier the fre- If we substitute here the expressions for Y from Eq. (7), quency dependence has to be included. This analysis is far Eq. (4), Eq. (5) and Eq. (6) we obtain the final expression for more complicated and is performed in Sec. A. the voltage constant describing our set-up In this section we present a short analysis of time constants and frequencies of the presented experimental set-up. 2 0 4 1 1 RF RF R1 +R2 RV The ADA4898-2 operational amplifier time constant τ ∗ = 1+2 . (14) U Um r RG R1 RV +Rav is calculated from the crossover frequency f0 with τ = G 15 1/2πf0. For amplification maximally close to the frequency As a whole, the circuit can be characterized as a low noise independent one (the considered scenario in the last section), pre-amplifier (instrumental amplifier followed by an inverting τ τF , where τF = RF CF is the time constant of the buffer amplifier all of them based on ADA4898-215) and a sequential feedback. In other words, the operational amplifier should not true RMS-meter. If necessary, the cheap AD63314 (maximum “feel” that there is a capacitance in its feedback. The amplified total error 2% of full scale, small signal bandwidth 1 MHz) by the buffer signal now has a higher time constant τA = τy1 5

3,4 or fA = f0/y1 or lower frequency due to the lower amplifi- zero. Even nowadays, there are no low-noise auto-zero op- cation of higher frequencies by the operational amplifier and erational amplifiers commercially available. In order to re-

τA τF because of the requirement for the feedback capac- move this unpleasant ﬂoating off the zero and the low-noise itance. After the buffer this signal passes through a high-pass 1/f noise in frequencies over 100 Hz, large CG = 10 µF met-

filter consisting of the resistor RG and the capacitor CG with allized polyester capacitors, whose price is comparable to that a large time constant τG = RG CG , whose function is to filter of the batteries and operational amplifiers, are put on the path 0 out the voltage offsets of the ADA4898-2 operational ampli- of the signal. 2) The small ceramic CF and CF connected fiers of the buffer. For a proper operation of the amplifier up parallely to the feedback resistors reduce the amplification in 0 0 to now, the time constants should be ordered high frequencies (> 1/RFCF, > 1/RFCF) and in this way stabilise the amplifies in terms of ringing.

τ τF τA τG . (18) In short, the removal of the offset and the ringing of the amplifier causes the addition of capacitors. If we addition- There is another capacitor C0 in the feedback of the IA and F ally assume the operational amplifiers to be ideal τ = 0, in analogously to the buffer, its presence should be barely “felt”, 0 the approximation CF = C = 0 and CG = ∞, for the therefore the time constant of the IA feedback τ 0 = R0 C0 < F F F F amplification of the amplifier we have the frequency indepen- τ . And finally the averaging low pass filter should have the A dent approximation Y (ω) ≈ y1y2y3, for which we have put largest time constant τav = RavCav in order to reliably aver- lots of effort to be a very good approximation in engineer- age the squared voltage. ing the amplifier. Infinite capacitance CG = ∞ means short Now let us return to the input signal, which is the time de- circuit and signal is going directly to the gain resistor, while pendent thermal fluctuations of the voltage of the parallely 0 CF = CF = 0 means absence of the parallel connected ca- connected capacitor C and resistor R, whose time constant is pacitors. τRC = RC. This time constant should be much larger than τA for a maximum amplification but much smaller than τG for a maximum transfer between the three steps of the whole ampli- IV. EXPERIMENT fier. Therefore, the requirement for the optimal amplification of our experimental set-up is A photograph of the realized circuit on a PCB is given in 0 Fig. 3, which was produced in 200 copies and given to high τ τF τ < τA τRC τG τav. (19) F school students in the 5th Experimental Physics Olympiad19 A set of the calculated time constants, frequencies and ad- (EPO5) is shown in Fig. 3. The set of the parameters from ditional calculated parameters of the experimental set-up is given in Table III,

Calculated parameter Value 15 f0 65 MHz τ 3.94 ns

τA 398 ns

fA 396 kHz

τRC (7.65–66.3) µs

τF 10 ns

τG 1 ms

fG 159 Hz 0 τF 100 ns τav 15 s Y 1.01 × 106 14 Um 10 V U ∗ 2.45 pV

TABLE III. Table of the calculated parameters necessary for the analysis of the circuit (time constants, frequencies, voltages and linear ampliﬁcation). f0 is the -3dB bandwidth of the ADA4898,

τ ≡ 1/2πf0, fA = f0/y1, τA = τy1, τRC = RC, τF = RF CF , 0 0 0 FIG. 3. A realization of the circuit on a printable circuit board given τF = RF CF , τG = RG CG , τav = RavCav, and 0 in EPO5.19 τ τF , τF , τA τRC τG τav.

the circuit and the integral schemes is listed in Table II. The experimental data measured with an ordinary cheap All-Sun Why do we need capacitors? 1) Even during the first DT-830B multimeter corresponding to Eq. (17) is presented attempt in measuring the Boltzmann constant, the Habicht in Table I and graphically shown in Fig. 4. The obtained value brothers noted that their amplifier suffers from floating off the of the slope q0 = (1718 ± 18) pC and the achieved accuracy 6

temperature is related with motion and ﬂuctuations. The vac- uum technique will remain expensive and even in good uni- 0.75 versities Maxwell velocity distribution is not experimentally demonstrated. However, for the last century electronics made a signiﬁcant progress and prices are reduced thousand times. 0.70 That is why the described set-up can be industrially produced

V [ ] for illustration of the thermal ﬂuctuations in high schools even

0.65 in the framework of the existing school programs. But even without commercially available set-ups the experimental set- up contains 3 integral circuits: 2 operational ampliﬁers, a mul- 0.60 tiplier, 3 large fast capacitors and 9 V batteries. Such a set-up

1.0×10 7 1.5×10 7 2.0×10 7 2.5×10 7 3.0×10 7 3.5×10 7 can be performed within one day in every high-school and 1/C[1/F] complete the teaching of thermal fluctuations and electronics. The problem of Boltzmann constant measurement with the de- 19 FIG. 4. Correlation between the voltage V measured by the right- scribed experimental setup was given in EPO5 to more than hand side multimeter from Fig. 2 and the reciprocal capacitance of a hundred high school students and more than 10% of them the capacitor C at the beginning of the same circuit. The linear were able to perform the experiment and to process their data regression with high correlation coefficient ρ = 0.9997 is a con- to obtain a final answer without any special training only in sequence of the equipartition theorem Eq. (2); see also Eq. (12). the framework of the standard education. We can certainly The slope of the straight line q0 determines the Boltzmann constant conclude that the whole experiment is appropriate for high- Eq. (17). school physics laboratories all over the world, where in addition a hardware or software oscilloscope can be used to vi- sualize the 1 million times amplified thermal noise. As a by for determination of Boltzmann constant is rather good product of EPO5 the authors of the present article can send free of charge to the physics laboratory one set-up to the first −23 kB = (1.40 ± 0.08) × 10 J/K (20) 137 teachers which will write us. When the set-up PCB ver- sion is outsourced in China, the measurement of Boltzmann for a $50 set-up which can be further elaborated. The ran- constant can reach any student from the 1-st, 2-nd and per- dom thermal voltage can be observed with an oscilloscope, haps 3-rd world before WWIII. In short, we have arrived to the temperature T of the resistor R can be varied from the the conclusion that it is time to introduce in the high school freezing to boiling point of the water but we are presenting education in physics simple experimental set-ups for determi- only the simplest experiment with varying only the capaci- nation of fundamental constants not only kB , but also electron 21 22 23 tor in room temperature which can be realized in every high charge qe , speed of light c and Planck’s constant ~ as school. It is intriguing to measure a fundamental constant by well. The touching to the fundamental determines the think- a set-up which can be created from scratch within a week by ing of the next generation. a novice.

ACKNOWLEDGMENTS

V. CONCLUSIONS The authors are grateful to Vasil Yordanov for his contribu- tion at the early stages of the present research,24 to Alexander Some administrators related to education in the third world Petkov for making the ﬁrst measurements, to Gary White for countries can conclude that very few high schools or intro- stimulating comments, to Nikolay Zografov for introducing ductory college courses would have the background, inter- order in the lab, to Andreana Andreeva for animation of the est and/or resources to adopt this experiment. However, the spirit in the lab, to Riste Popeski-Dimovski, Marina Poposka, equipartition theorem is already in high school education for Sladjana Nikolic, Slavoljub Mitic and Stojan Manolev for the more than a century. More than a century it is taught that the invaluable help and unforgettable moments during EPO.

∗ E-mail: [email protected] only), Proc. R. Soc. Lond., 5, p. 604, (1851); † E-mail: [email protected]ﬁa.bg J. J. Waterston, Lord Rayleigh, Sec. R. S., 1 J. Herapath, “On the physical properties of gases”, Annals Philos- Phil. Trans. R. Soc. Lond. A, 183, pp. 1–79, (1892), ophy, 8, pp. 56–60, (1816); doi:10.1098/rsta.1892.0001. “On the Causes, Laws and Phenomena of Heat, Gases, Gravita- 3 A. Einstein, “Uber¨ die Gultigkeitsgrenze¨ des Satzes vom ther- tion”, The Annals of Philosophy, 1, pp. 273–293, (1821). modynamischen Gleichgewicht und uber¨ die Moglichkeit¨ einer 2 J. J. Waterston, “On the Physics of media that are composed of neuen Bestimmung der Elementarquanta (On the boundaries free and perfectly elastic molecules in a state of motion” (abstract of the applicability of the law of thermodynamic equilibrium 7

and on the opportunity of a new definition of the elementary 4898-2.pdf. quanta),” Ann. Phys. 22, 569-572 (1907) (in German), http: 16 Metallized Polyester (PET) Capacitors in PCM 5 mm. Ca- //myweb.rz.uni-augsburg.de/˜eckern/adp/ pacitances from 0.01 mF to 10 mF. Rated Voltages from 50 history/einstein-papers/1907_22_569-572.pdf. VDC to 630 VDC, MKS2, 10 µF, https://www.wima.de/ 4 C. Habicht and P. Habicht, “Elektrostatischer Potentialmultiplika- wp-content/uploads/media/e_WIMA_MKS_2.pdf. tor nach A. Einstein,” Phys. Ztschr. 11, 532-555 (1910) (in Ger- 17 DT830 Series 3 1/2 Digital Multimeter, man); Amplification 360000 (111 dB); δU = 300 µV. www.all-sun.com/manual/Dt830_en.pdf. 5 Y. Kraftmakher, “Two student experiments on electrical fluctua- 18 “Technical Specification AD835: 250 MHz, Voltage Out- tions”, Am. J. Phys. 63, 932–935 (1995). put, 4-Quadrant Multiplier”, http://www.analog.com/ 6 D. R. Spiegel and R. J. Helmer, “Shot-noise measurements of static/imported-files/data_sheets/AD835.pdf. the electron charge: An undergraduate experiment”, Am. J. 19 T. M. Mishonov, E. G. Petkov, A. A. Stefanov, A. P. Petkov, Phys. 63, 554–560 (1995), http://scitation.aip.org/ I. M. Dimitrova, S. G. Manolev, S. I. Ilieva, A. M. Varonov, “Mea- content/aapt/journal/ajp/63/6/10.1119/1. surement of the Boltzmann constant by Einstein. Problem of the 17867;jsessionid=IlKxP8DoJOtP8pB_dgF7AYT8. 5-th Experimental Physics Olympiad. Sofia 9 December 2017”, x-aip-live-06. arXiv:1801.00022 [physics.ed-ph]. 7 “Johnson Noise and Shot Noise: The Determination of the Boltz- 20 T. M. Mishonov, V. I. Danchev, E. G. Petkov, V. N. Gourev, mann Constant, Absolute Zero Temperature and the Charge of the I. M. Dimitrova, A. M. Varonov, “Manhattan equation for the op- Electron”, MIT Department of Physics (2013), erational amplifier”, arXiv:1802.09342 [eess.SP]. “Johnson Noise and Shot Noise”, MIT Junior Lab Regular 21 T. M. Mishonov, E, G. Petkov, N. Zh. Mihailova, A. A. Stefanov, Experiment, (2013), I. M. Dimitrova, V. N. Gourev, N. S. Serafimov, V. I. Danchev http://web.mit.edu/8.13/www/43.shtml, and A. M. Varonov, “Simple do-it-yourself experimental set-up http://web.mit.edu/8.13/www/JLExperiments/ for electron charge qe measurement”, Eur. J. Phys. 39, 065202 JLExp43.pdf, Used equipment: Low-noise voltage pream- (2018). plifier, SR560 – DC to 1 MHz, http://thinksrs. 22 T. M. Mishonov, A. M. Varonov, D. D. Maksimovski, com/downloads/PDFs/Catalog/SR560c.pdf, e = S. G. Manolev, V. N. Gourev and V. G. Yordanov, An undergrad- √ N 4 nV/ Hz, $2595; Krohn-Hite 3BS8TB-1k/50kg, 8-pole But- uate laboratory experiment for measuring ε0, µ0 and speed of terworth, cutoff accuracy 2%, http://www.krohn-hite. light c with do-it-yourself catastrophe machines: electrostatic and com/htm/filters/PDF/ModulesData.pdf. magnetostatic pendula, Eur. J. Phys. 38, 025203, (2017). 8 D. B. Pengra, R. Van Dyck and J. Stoltenberg, “Fun- 23 D. S. Damyanov, I. N. Pavlova, S. I. Ilieva, V. N. Gourev, damental noise and fundamental constants”, University of V. G. Yordanov and T. M. Mishonov Planck’s constant mea- Washington Modern Physics Laboratory – Condensed Matter surement by Landauer quantization for student laboratories, (2012), http://courses.washington.edu/phys431/ Eur. J. Phys. 36, 055047, (2015). noise/new_noise_old_box.pdf. 24 T. M. Mishonov, V. G. Yordanov, A M. Varonov, “Measurement 9 J. C. Rodriguez-Luna and J. de Urquijo, “A simple, sensitive cir- of electron charge qe and Boltzmann’s constant kB by a cheap cuit to measure Boltzmann’s constant form Johnson’s noise”, Eur. do-it-yourself undergraduate experiment,” arXiv:1703.05224 J. Phys. 31, 675–679 (2010). [physics.ed-ph], (2017). 10 M. Planck, “Zur Theorie des Gesetzes der Energieverteilung 25 M. Lee, T. Schibli, “Electronics for the Physical Sciences,” im Normalspektrum”, Verhandlungen der Dtsch. Phys. Ges., 2, Physics 3330 course, University of Colorado Boulder, (2012), pp. 237–245, (1900) (in German); http://www.colorado.edu/physics/phys3330/ “On the Theory of the Energy Distribution Law of the Normal phys3330_sp12/phys3330_sp12/Home.html; “Ex- Spectrum”, in The Old Quantum Theory edited by D. ter Haar, periment #4: Operational Amplifiers and Negative Feed- (Pergamon Press, Oxford, New York, 1967), p. 82. back,” Eq. (2), http://www.colorado.edu/physics/ 11 D. Flamm, “Ludwig Boltzmann – A Pioneer of Modern Physics”, phys3330/phys3330_sp12/phys3330_sp12/Lab_ arXiv:physics/9710007 [physics.hist-ph], (1997); Manual_files/Exp_4_Spring12.pdf. W. Ostwald, “Lebenslinien Eine Selbstbiographie”, Klasing, 26 C. Kittel and H. Kroemer, Thermal Physics (2nd Ed.), Berlin, 2 p. 187 and 188, (1927) (in German), citing the reply (W. H. Freeman and Company, San Francisco, 1980). by Boltzmann. 27 L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics, 12 S. G. Brush, “Foundations of statistical mechanics 1845– Vol. 9, E. M. Lifshitz and L. P. Pitaevskii Statistical Physics, Part 1915”, Arch. Hist. Exact Sci., 4, pp. 145–183, (1967), DOI: 2, (Pergamon Press, New York, 1980), Sec. 78 “Current fluctua- 10.1007/BF00412958; tions in linear circuits”, Eqs. (78.1-3). T. L. Hankins, “Nineteenth-Century Physics”, Science 196, Issue 28 C. W. McCombie, “Nyquist theorem and its generalisations”, 4291, pp. 783–784, (1977), DOI: 10.1126/science.196.4291.783; Chap. 24 in Problems in Thermodynamics and Statistical Physics, 13 “Technical Specification ADA4817: Low Noise, 1 GHz Fast- (PION, London, 1971), edited by P. T. Landsberg. FET Op Amps,” http://www.analog.com/static/ 29 L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics imported-files/data_sheets/ADA4817-1_ Vol. 8, Electrodynamics of Continuous Media (1st. Ed.), (Fiz- 4817-2.pdf. matlit, Moscow, 1959), (in Russian). 14 “Technical Specification AD633: Low cost Analog Four- 30 C. W. Gardiner, Handbook of Stochastic Methods for Physics, Quadrant Multiplier”, Fig. 17, http://www.analog.com/ Chemistry and the Natural Sciences (2nd. Ed.), (Springer-Verlag, static/imported-files/data_sheets/AD633.pdf. Berlin, 1985). 15 Technical Specification ADA4898: High Voltage, Low Noise, Low 31 H. Nyquist, “Thermal agitation of electric charge Distortion, Unity-Gain Stable, High Speed Op Amp, Table 1, in conductors”, Phys. Rev. 32, 110–113 (1928), http://www.analog.com/static/ http://journals.aps.org/pr/pdf/10.1103/PhysRev.32.110. imported-files/data_sheets/ADA4898-1_ 32 H. B. Callen and T. A. Welton, “Irreversibility and Generalized Noise”, Phys. Rev. 83, 34-40 (1951). 8

33 L. D. Landau and E. M. Lifshitz, Course of Theoretical Physics, 2. Buffer Vol. 5, Statistical Physics, Part 1, (Pergamon Press, New York, 1980), Eqs. (123.1) and (124.10). For the buffer, consisting of two non-inverting ampliﬁers 34 W. Schottky, “Zur Berechnung und Beurteilung des Schroteffek- tes”, Ann. d. Physik 373 (10), (1922), (in German). (NIA) depicted in Fig. 5, the input voltages U1 and U2 are ap- 35 W. Sutherland, “A Dynamical Theory for Non-Electrolytes and plied directly to the (+) inputs of each OpAmp U+ = U1, and 0 the Molecular Mass of Albumin”, Phil. Mag. 9, 781 (1905). U+ = U2, here primed notations refer to the lower NIA. Since 36 A. Einstein, “Uber¨ die von molekularkinetischen¨ Theorie der Warme¨ geforderte Bewegung von in ruhenden F¨lussigkeiten sus- pendierter Teilchen”, Ann. d. Physik 17, 549–560 (1905), (in Ger- man). 37 A. Pais, Subtle is the Lord..., (Oxford University Press, Oxford, 1982).

Appendix A: Analog electronics in a nutshell

This appendix is addressed to colleagues involved with con- struction of similar devices and modifying the scheme. This recall of the standard electronics notion is not addressed to the students.

FIG. 5. Circuit of a buffer with capacitance in the feedback. 1. Operational amplifier master equation no currents enter in the OpAmp, the only current I starts from In the beginning was the approximate master equation of the output U0, passes through the upper feedback ZF, gain operational amplifiers resistor rG , lower feedback ZF and terminates in the output 0 U0. From this closed loop, it is straightforward to obtain an −1 1 d expression for the current Gˆ U0(t) = + τ U0 = U+ − U−, (A1) G0 dt 0 (U0 − U0) ≡ −∆U0 1 1 I = , = + jωCF. (A6) 2ZF + r ZF(ω) RF giving the relation between the output voltage U0(t) and dif- G ference of (+) and (-) inputs of operational amplifier, cf. The voltage difference between the (-) inputs of both OpAmp Ref. 25. For harmonic signals, introducing j-imaginary unit of the buffer is jωt st U ∝ e = e , j = −i, s ≡ jω, (A2) 0 rG U− − U− = −IrG = ∆U0 . (A7) 2ZF + rG we have The master equations for both operational amplifiers are −1 1 −1 G (ω)U0 = + jωτ U0 = U+ − U−, (A3) U G (ω) = U − U (A8) G 0 1 − 0 0 −1 0 2 3 U0G (ω) = U2 − U−, (A9) −1 1 1 + a1s + a2s + a3s + ... G = +sτ 2 3 , (A4) G0 1 + b1s + b2s + b3s + ... and after subtracting them, we obtain a single equation

−1 0 where in the second row we present a Pade´ approximant for ∆U0G = (U1 − U2) − (U− − U−). (A10) the frequency dependent open loop gain. For the used by us operational ampliﬁer ADA489815 the static open loop gain is Substituting Eq. (A7) into the last equation, the frequency de- approximately 100 dB pendent ampliﬁcation of the buffer

∆U0 1 5 1 G0 = 10 , f0 ≡ = 65 MHz (A5) ΥNIA(ω) ≡ = −1 −1 , (A11) 2πτ U1 − U2 G (ω) + y (ω) ZF(ω) and the time τ constant is parameterized by the crossover fre- y(ω) ≡ + 1 (A12) rG quency f0. In the next subsections we will recall how the frequency in agreement with Ref. 13, Eq. (4). Using complex numbers 2 dependent open loop gain determines the frequency dependent in programming we can simply calculate |ΥNIA(ω)| , how- transmission function of different ampliﬁers. ever using only real numbers we have to apply some efforts in 9

−1 complex algebra and after a straightforward calculation from This expression is equal to G (ω)U∆ according to Eq. (A1) A11 we obtain (here U∆ is the output voltage) and therefore after a little re- 2 arrangement of the terms |ΥNIA(ω)| (A13) 2 0 N (ω) ZG R U G−1 + = −(U − U 0 ) F . (A20) = 2 , ∆ 0 0 0 0 ZG + RF ZG + RF −1 −1 2 2 2 2 τ G0 N + (y1 + ω τs ) + (ωτs) N + I0 τs The frequency dependent ampliﬁcation of the difference am- where pliﬁer therefore is

2 2 −1 CFRF U∆ −1 N (ω) = 1 + ω τs ,I0 ≡ 1 − y1 , τs ≡ . (A14) Υ∆(ω) ≡ 0 = −1 , (A21) y1 U0 − U0 Λ(ω) + G (ω)[1 + Λ(ω)]

0 The calculation of the pass-bandwidth requires the square of where Λ(ω) ≡ ZG(ω)/RF. The calculation of the pass- the modulus of the complex amplifications for all steps of the bandwidth requires the square of the modulus of the complex amplifier which we are presenting in full detail. amplification, which after a straightforward calculation from A21

2 3. Difference Amplifier |Υ∆(ω)| (A22) (ωτ )2 = g , The difference amplifier shown in Fig. 6 has a capacitor CG 2 −1 2 2 2 τ −1 at each of its inputs. The current IF flows from the first input 1+G0 −Mω ττg + (ωτg) +Λ0 +G0 M τg

RG 0 Λ0 ≡ 0 ,M ≡ 1 + Λ0, τg ≡ CGRF. RF This complicated expression is necessary if we wish to use real numbers in the numerical integration necessary to be pre- formed for the pass bandwidth. The last step of the ampliﬁer is the inverting one described in the next subsection.

4. Inverting Amplifier FIG. 6. Circuit of a difference amplifier with capacitance in the input. 0 Setting U0 = 0 in the difference amplifier from the last U0 of the difference amplifier to the output U∆ and subsection, we obtain an inverting amplifier (IA). In our case U − U 1 there is one more difference, which is a capacitor CF con- I = ∆ 0 ,Z (ω) = R + , (A15) nected parallely to the feedback resistor R0 and therefore the F Z + R0 G G jωC F G F G reciprocal of the impedance feedback 0 while the current I0 flows from the second input U0 to the 1 1 common point and 0 = 0 + jωCF. (A23) ZF(ω) RF −U 0 0 Λ(ω) Γ(ω) ≡ Z (ω)/Z0 I0 = 0 . (A16) Re-denoting to G F, the frequency de- ZG + RF pendent amplification of the IA The voltage drops at the inputs of the OpAmp 1 ΥIA(ω) = − (A24) U − U Γ(ω) + G−1(ω)[1 + Γ(ω)] U = I Z + U = ∆ 0 Z + U , (A17) − F G 0 Z + R0 G 0 G F in agreement with Ref. 13, Eq. (7). This agreement gives 0 0 −U0 0 an implicit proof of the applicability of the time-dependent U+ = I0ZG + U0 = 0 ZG + U0. (A18) ZG + RF Eq. A1 and its frequency transformation Eq. A3 for circuits with low-noise operational amplifiers ADA4898, ADA4817, Subtracting the last two equations, we obtain an expression AD711 etc. The calculation of the pass-bandwidth requires for the voltage difference of the OpAmp inputs the square of the modulus of the complex amplification, which after a straightforward but a bit more complicated calculation 0 ZG U∆ZG U+−U− = −(U0−U0) 1 − 0 − 0 . (A19) from A24 ZG +RF ZG +RF 10

2 2 (ωτg) |ΥIA(ω)| = , (A25) 0 2 −1 2 2 2 02 2 τ −1 τf −1 2 (1+G0 )(1−ω τfgτg)−Mω ττg −ω ττf + (ωτg) +Λ0 +G0 M + (G0 +1)−ω ττfg τg τg

where In the approximation of frequency independent amplification Υ(ω) ≈ Y = y y y we have 0 0 1 2 3 τf ≡ CFRF, τfg ≡ CFRG. (A26) Z ∞ 4k TR dω 4k TR 0 U 2 (t) ≈ Y 2 B = Y 2 B . For CF = 0, τf = τfg = 0 and the modulus of the amplifica- amp 2 0 1+(ωRC) 2π 4RC tion of the IA Eq. (A25) becomes equal to the modulus of the (A33) amplification of the difference amplifier Eq. (A22). That is why for simplicity we will represent the exact results Finally, for the whole amplifier we have as the approximation formula

Υ(ω) = ΥNIA(ω)Υ∆(ω)ΥIA(ω). (A27) k T U 2 (t) = Y 2 B Z, (A34) amp C 5. Application of the Nyquist theorem for our circuit corrected by the close to unity coefﬁcient

∞ Z In order to calculate the mean square of the ampliﬁed sig- 2 dω 2 |Υ(ω)| 2 nal Uamp(t) we have to apply the Nyquist theorem for the 1 + (ωRC) spectral density of the thermal noise of the parallely connected 0 Z ≡ 1 − ε ≡ ∞ . (A35) capacitor C and resistor R at the input of the circuit. Accord- Z dω Y 2 ing this theorem the spectral density of the noise is given by 1 + (ωRC)2 the real part of the impedance 0 2 The small correction ε(C) as a function of the input capaci- (U )f = 4kB TR(ω),R(ω) = <(Z(ω)), ~ω kB T. (A28) tances C is graphically presented in Fig. 7.

A new pedagogical re-derivation of Nyquist theorem and its extensions is given in Appendix B. Although the quantum physics was created from the Planck explanation of the spectral density of the electromagnetic ﬂuc- tuations and 3 Noble prizes were given for the black body radiation, there is no deﬁnite notation for the spectral density of the voltage noise in electronics. For instance, in 2 2 Ref. 26, Chap. 4, Ref. 27, Eq. 78.3, where (E )ω = (E )f /2 2 or Ref. 28, Chap. 24, where GV = (E )f /2π, or Ref. 29, 2 2 Eq. 89.3, where (E )ω = (E )f /4π, or Ref. 30, Eq. (1.4.47), 2 where S(ω) = (E )f /4π. Applying Nyquist theorem to a paralelly connected capacitor and resistor with impedance

1 −1 Z(ω) = + jωC , j = −i. (A29) R for the spectral density of the voltage noise we obtain FIG. 7. The correction ε as a function of the capacitance C. The R (U 2) = 4k T . (A30) behaviour of the correction confirms the analysis made in Sec. III. f B 1 + (ωRC)2 This spectral density is amplified by the 3 steps of the ampli- High school students are freely able to use this figure, fier Eq. (A27) while the university students can calculate it numerically. We do not recommend the analytical calculation of the integral 2 2 2 (U )f = |Υ(ω)| (U )f (A31) Eq. (A32). amp 2 After calculating Uamp(t) we can substitute it in Eq. (9) and for the mean square of the amplified voltage we have and obtain Z ∞ 2 2 2 4kB TR dω Uamp R1 + R2 U (t) = |Υ(ω)| . (A32) hU i = . amp 2 (A36) 0 1 + (ωRC) 2π Um R1 11

Using the described above correction factor, Eq. (11) now In our derivation of the Nyquist theorem we treat the av- reads erage energy of an inductance hELi participating in a high- quality 2 ZY R1 + R2 RV 2 V2 = U (t) , (A37) p Um R1 RV + Rav L/C Q ≡ 1 (B1) R with resonance circuit depicted in Fig. 8. We will compare the re- Z ∞ 4k TR dω k T U 2(t) = B = B . (A38) sult derived by analysis of the influence of the spectral density 2 2 0 1+(ωRC) 2π C of the random noise created by the resistor (E )f with the result of the Gibbs averaging of the energy Keeping the expression Eq. (14) for U ∗ unchanged, we finally arrive at the correct formula Eq. (17) used for the determina- 1 ~ω0 1 tion of the Boltzmann constant. ε¯ = ~ω0 coth , hEC i = hELi = ε¯ (B2) 2 2kB T 2 √ of an oscillator with resonance frequency ω0 = 1/ LC. The Appendix B: Johnson-Nyquist thermal noise and Callen-Welton impedance of the circuit Fluctuation-Dissipation Theorem 1 Z(ω) = jωL + R + . (B3) In order to understand the work of the described set-up, jωC only the equipartition theorem is necessary and the equipartition theorem is given in the all introductory courses in physics gives for the modulus of the conductivity σ = 1/Z in not only in the universities but also in many high-school 1 textbooks. Unfortunately speaking about the thermal fluc- |σ(ω)|2 = . (B4) 2 1 2 tuation of the voltage people erroneously suppose knowl- R + ωL − ωC edge of the frequency dependent spectral density described by Nyquist and Callen-Welton fluctuation-dissipation theorems. The spectral density of the current is given by the spectral These theorems already belong to the extended courses on density of the voltage statistical physics. Definitely those theorems cannot be re- 2 2 2 produced by at least of 51% of physics teachers. That is why (I )f = |σ(ω)| (U )f , (B5) in this appendix we will give a new derivation of all those theorems in the level corresponding to the introductory courses which for high-Q resonance circuit gives on physics read for the future physics teachers. This appendix 2 is oriented for readers willing to understand the work of the 2 π (E )f 2 (I )f = δ(ω − ω0), (E )f = 4Rε.¯ (B6) set-up, starting from the ideas and notions given in the univer- 2 RL sity courses on statistical physics and nothing is beyond this We will prove the Nyquist theorem frame. 2 (E )f = 4Rε.¯ (B7)

1. Nyquist theorem and its generalizations Expressing the dispersion of the current from its spectral density If the Nyquist theorem31 is derived from the statistical Z ∞ physics methods, it is viewed as an interesting applica- 2 2 dω I = (I )f , (B8) tion of the Callen-Welton Fluctuation and Dissipation The- 0 2π orem (FDT).32,33 In the Landau-Lifshtitz course in theoretical physics, both in the 5-th volume Statistical Physics33 and a simple integration of the δ-function in Eq. (B6) gives the first edition of the 8-th volume Electrodynamics of Con- 29 32 1 2 1 tinuous Media, a derivation following the original one is hELi = L I = ε.¯ (B9) given. In this approach a Gibbs averaging of the Dirac time- 2 2 dependent perturbation theory is performed. The practice in An analogous consideration of the energy of the capacitor teaching physics however, shows that this derivation meant for professionals is unreproducible by the students paying tuition Q2 Z ∞ dω (I2) fees and waiting to receive an educational service; the degra- hE i = , Q2 = (Q2) , (Q2) = f C 2C f 2π f ω2 dation of the physics education is global and the authors of the 0 present article are unable to point-out any contemporary text- gives books, where FDT is derived. That is why in this Appendix we give a new re-derivation of the FDT approbated by many 1 hE i = ε¯ (B10) recruits of students. C 2 12 in agreement with the virial theorem. In such a way we have has√ a sharp maximum at the resonance frequency ω0 = verified that the mean energy of an oscillator 1/ LC and is negligible far from the resonance and we can use the δ-function approximation Q2 LI2 hEC i + hELi = + =ε ¯ (B11) Z ∞ 2C 2 F (ω0) = F (ω)δ(ω − ω0)dω, ω0 > 0. (B17) 0 coincides with the Gibbs mean value and we have re-derived the Nyquist theorem Eq. (B7) as a consequence of detailed The coefficient in front of the delta function is given by the balance principle applied to an electric oscillator. In the clas- integral sical case of small frequencies ω k T , which covers ~ 0 B Z ∞ 2 almost all electronic applications, we obtain the equipartition ω dω π 1 2 2 2 2 = , (B18) theorem 0 R ω + (Lω − 1/C) 2 RL 1 which does not depend on the capacitance C and was solved hE i = hE i = k T (B12) C L 2 B by Schottky34 (Eq. (II)). Introducing the dimensionless variable x ≡ ωL/R, and dimensionless parameter Q the corre- and low frequency asymptotic of the Nyquist theorem for the sponding mathematical problem is spectral density of the noise of a resistor Z ∞ x2dx π 2 I(Q) = = , (B19) (E )f = 4kB TR. (B13) 2 2 2 2 0 x + (x − Q ) 2 One can say that the Nyquist theorem is a consequence of the which can be solved both analytically (given in Subsec. B 3) equipartition theorem and detailed balance principle. and numerically for different values of a. For Q = 0 we have R ∞ 2 a table integral −∞ dx/(1 + x ) = π. Using this δ-function approximation, we obtain for the 2. Nyquist theorem and its generalizations spectral density of the current

2 In Fig. 8 a resonance LC circuit with a resistor R creating 2 π (E )f 2 (I )f = δ(ω − ω0), (E )f = 4Rε.¯ (B20) random noise is depicted. The random thermal voltage Eω cre- 2 RL Experimentally such frequency dependence of the spectral density can be investigated using Fourier transformation of digital oscilloscopes and generators. Such an equipment is not typical for high-schools and that is why for physics teachers the present appendix in only an additional material complet- ing the theory from general introductory courses. Whence for the thermal averaged energies of the capacitor and inductance 1 Z ∞ dω hE i = C U 2 , U 2 = |Z |2 (I2) , C 2 c c C f 2π FIG. 8. Resonance circuit with inductance ZL = jωL, capaci- 0 Z ∞ tor ZC = 1/jωC and resistor ZR = R creating random electric 1 2 2 dω 2 hE i = L I , hIi = (I ) , voltage with spectral density given by Nyquist theorem (E )f = L f 2 0 2π 2~ωkB T coth(~ω/2kB T ). a trivial integration of the δ-functions ates a current amplitude Iω = σ(ω)E(ω), where the conduc- Z ∞ tivity σ(ω) = 1/Z(ω) is determined by the total impedance of F (ω0) = F (ω)δ(ω − ω0)dω (B21) the sequentially connected inductance, resistor and capacitor 0 1 gives Z(ω) = jωL + R + . (B14) jωC 1 1 ~ω0 hEC i = hELi = ε,¯ ε¯ = ~ω0 coth . (B22) As a gedanken experiment, let us analyse a high quality reso- 2 2 2kB T nance circuit, for which The Nyquist result for the spectral density of the noise can be r r 1 L pL/C presented in a more general form. R ωL = , Q ≡ 1. (B15) ωC C R In this case the square of the modulus of the conductivity 3. Solution of the Schottky integral

2 1 π 1 |σ(ω)| = ≈ δ(ω − ω0) (B16) Contemporary students use Mathematica or Maple but 100 2 1 2 2 RL R + ωL − ωC years ago physicists were able to perform elementary calculus 13 of integrals.34 In this subsection we give the calculation of the The solution of such a table integral is integral considered by Schottky when he analyzed the inﬂu- Z ∞ Z ∞ ence of noise the of an reconance circuit. Let us consider the dx 1 d(x/a) 1 π 2 2 = 2 = integrand from Eq. (B19) 0 x + a a 0 (x/a) + 1 a 2

x2 C D and hence for the solution of the integral we obtain = + , (B23) x2 + (x2 − Q2)2 x2 + a2 x2 + b2 1 π 1 π I = 2 2 (a − b) = . where ±a and ±b are the x values for which the denominator a − b 2 a + b 2 on the left is equal to 0 and C and D are coefficient we are Using the well known binomial theorem (a+b)2 = a2 +2ab+ going to find. Expanding the divisors both on the left and b2, we substitute the values for a2 +b2 and ab from Eqs. (B26) right sides into a polynomial, we obtain to obtain 2 2 2 2 4 2 2 2 2 2 (x + a )(x + b ) = x + (a + b )x + a b , (B24) (a + b)2 = 1 − 2Q2 + 2Q2 = 1, (B30) x2 + (x2 − Q2)2 = x4 + (1 − 2Q2)x2 + Q4. (B25) and hence for the solution we have Both polynomials have to be identical, meaning that the co- 1 π π efficients of the respective x degree terms should be equal. I = = , (B31) a + b 2 2 Clearly, the coefficients of the x4 terms are both equal to 1, and comparing the quadratic and zero degree terms, we ob- where we have taken only the positive value. tain a set of two equations (one per each degree term) √ a2 + b2 = 1 − 2Q2 (the terms with x2), b. Complex values of a and b, Q > 1/ 2 a2b2 = Q4 (the terms with x0). (B26) The solution of the integral in this case is analogous. The From the first√ equation it is evident√ that for real values of a parameters a and b are complex. However for the interesting and b, Q < 1/ 2, while for Q > 1/ 2 we have a2 + b2 < 0 for us case Q 1 we can recall the analytical continuation. and hence a and b are complex. If I(Q) = π/2 for 0 < Q < 1/2 this result has a unique analytical continuation and can be extended from a finite segment √ to the whole axis Q > 0. Now we address a new derivation of a. Real values of a and b, Q < 1/ 2 the FDT theorem.

First we are going for a solution in case of real values for a and b. Let us take a closer look at the ﬁrst of these two 4. Callen-Welton ﬂuctuation-dissipation theorem equations. Now going back to Eq. (B23), we combine both fractions on the right hand side into one and the nominator In the spectral density of the voltage noise becomes 2 (E )f = 2R(ω)~ω coth(~ω/2kB T ), (B32) Cx2 + Cb2 + Dx2 + Da2 = (C + D)x2 + (Cb2 + Da2). 1 R(ω) = <(Z(ω)),Z(ω) = , jωC Comparing it with the nominator on the left hand side x2, we obtain the following two equations the frequency dependent resistance R(ω) is the real part of the complex impedance which can be represented also by a 2 C + D = 1 (the terms with x ), (B27) frequency dependent capacitance Cb2 + Da2 = 0 (the terms with x0). (B28) i i C00 ωR(ω) = ω< = (C−1 − C−1∗) = , Solving these equations for C and D, we get ωC(ω) 2 |C|2

a2 b2 where C0 = <(C(ω)), C00 = =(C(ω)) and C = C0 + iC00. C = ,D = − . (B29) a2 − b2 a2 − b2 Here we suppose an arbitrary frequency of the impedance represented by a generalized capacitance. On the other hand, the Substituting these expressions for C and D into Eq. (B19), the spectral density of the charge is also given by the capacitance integral 2 2 2 (Q )f = |C(ω)| (E )f . (B33) Z ∞ x2dx Z ∞ C D I = 2 2 2 2 = 2 2 + 2 2 dx 0 x + (x − Q ) 0 x + a x + b In such a way, the Nyquist theorem for the thermal noise can a2 Z ∞ dx b2 Z ∞ dx be rewritten as = − . a2 − b2 x2 + a2 a2 − b2 x2 + b2 2 00 0 0 (Q )f = 2~C coth(~ω/2kB T ) (B34) 14 and for thermal averaged charge fluctuations we finally arrive In this case the Nyquist theorem can be rewritten as at Z ∞ ω x2 = ~ α00(ω) coth ~ dω, Vˆ = −f(t)ˆx. Z ∞ Z ∞ π 0 2kB T 2 2 dω ~ 00 ~ω Q = (Q )f = C (ω) coth dω. (B40) 0 2π π 0 2kB T cf. Eq. (V.124.10) and Eq. (V.123.1) of the Landau-Lifshitz (B35) 33 2 2 course of theoretical physics. The terminology also has to We use the spectral density (E )f = 2(E )ω which takes into be changed and the Nyquist theorem is now well-known as account the folding of positive and negative frequencies ω. fluctuation-dissipation theorem (FDT) for generalized suscep- Let us recall again the formula for the energy when the con- tibility α(ω) giving the general relations between the fluctua- sidered capacitor C is under the voltage E = Q0/C0 created tions of some quantum variable xˆ and dissipation (absorptive) 00 by a large capacitor (charge reservoir) C0 C part of the generalized susceptibility α = =(α(ω)). Callen and Welton rederived the Nyquist FDT applying Gibbs thermal averaging of the quantum mechanical second Q2 (Q − Q)2 Q2 E = + 0 = − EQ + const. (B36) order perturbation theory, as it is now described in every pro- 2C 2C0 2C fessionally written textbook on statistical physics. Even the history of physics is rewritten and in some textbooks one can read that the Nyquist theorem is an interesting application27 of One can consider that the FDT. The Sutherland-Einstein relations between the diffu- sion coefficient and mobility is one of the first examples of he V = −E(t)Q (B37) FDT.35–37 In short we represented all those theorems of the statistical physics as simple consequence of the principle of the de- is the interaction energy of the capacitor with the voltage tailed balance applied to the simplest physical system – the source. The formulae Eqs. (B35) and (B37) are so general harmonic oscillator. But each creation is a child of its own that it deserves to change the notations time, nowadays after more than a century of technology development, every teacher is able to illustrate the thermodynamic fluctuations of the voltage and the equipartition theorem af- ˆ αˆ(ω) = C(ω), xˆ = Q, f = E, V = V, hxîω = αωfω. ter only a single day work for the building of the experimental set-up. Tenths of students used the described experimen- Analogously in time representation tal set-up during the EPO5 and successfully determined the Boltzmann constant and the elementary formulae describing the electronic circuit operation as well. Z t 0 0 0 What has changed since the time of Habicht brothers till Q(t) = C(t − t )U(t )dt (B38) now? That is the invention/appearance of low-noise oper- −∞ ational amplifiers such as ADA4898 and multipliers of the AD633 type, which make possible RMS voltages measure- becomes ments of million times amplified thermal noise. Now could be motivated to follow our new derivation of the Nyquist and Z ∞ Callen-Welton theorems which is actually detailed balance x(t) = α(τ)f(t − τ)dτ. (B39) principle applied to the oscillator written in electric variables. 0