arxiv:physics/ v1 [physics.class-ph] 3 Oct 2005

Transcription

1 CESJ 1 (2003) 1 25 Central European Science Journals arxiv:physics/ v1 [physics.class-ph] 3 Oct 2005 The cubic period-distance relation for the Kater reversible pendulum M. Rossi 1, L. Zaninetti 2 1 Dipartimento di Matematica, Università degli Studi di Torino Via Carlo Alberto 10, Torino, Italy 2 Dipartimento di Fisica Generale, Università degli Studi di Torino via P.Giuria 1, Torino,Italy Received ; revised Abstract: We describe the correct cubic relation between the mass configuration of a Kater reversible pendulum and its period of oscillation. From an analysis of its solutions we conclude that there could be as many as three distinct mass configurations for which the periods of small oscillations about the two pivots of the pendulum have the same value. We also discuss a real compound Kater pendulum that realizes this property. c Central European Science Journals. All rights reserved. Keywords: reversible pendulum, physics of the pendulum, mathematics of the pendulum PACS (2000): Pa Laboratory experiments and apparatus; b General physics; Ud Linear algebra; 14Rxx Affine Geometry; b Computational methods in classical mechanics; 1 Introduction A well known consequence of the fundamental equation of rotational dynamics is that the period of small oscillations of a physical pendulum is given by T = 2π ω = 2π I mgh (1) where m is total mass of the pendulum, I its moment of inertia with respect to the center of oscillation O and h the distance of the center of mass from O. Then a physical michele.rossi@unito.it zaninetti@ph.unito.it

2 2 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 pendulum oscillates like a simple pendulum of length l = I mh = gt2 4π 2 (2) which is called the equivalent length of our physical pendulum. By the Huygens Steiner theorem (also known as the parallel axis theorem ) it is possible to write I = mh 2 +I 0 wherei 0 isthemoment ofinertiawithrespect tothecenter ofmass. Bysquaring equation (1) we get the following quadratic relation h 2 lh+ I 0 m = 0 (3) When l 2 4I 0 /m 0 that equation admits two real solutions h 1,h 2 such that h 1 +h 2 = l (4) In 1817 Captain H.Kater thought to use this last relation to empirically check the Huygens Steiner theorem. At this purpose he constructed his reversible pendulum consisting of a plated steel bar equipped with two weights, one of which can be moved along the bar. This pendulum is reversible because it can oscillate about two different suspension points realized by two knife edges symmetrically located on the bar. By adjusting the movable weight, it is possible to obtain a pendulum mass configuration such that the periods about the two pivots coincide, the equivalent length l is the distance between the two knife edges and condition (4) is satisfied. Themeasurement ofsuchacommonperiodt,ofthetotalmassmandofthedistancel between the two knife edges, gives then an easy way to perform an empirical measurement of the earth s (apparent) gravitational acceleration g by applying formula (2). This is why the Kater reversible pendulum is one of the favourite instrument to measure g in student labs. Anyway there is a subtle point in this procedure which is the determination of the right mass configuration of the pendulum. This problem gives rise to the following two questions: (1) how many possible positions of the movable weight determine a good mass configuration for which the periods of small oscillations about the two pivots coincide? (2) when a good mass configuration is realized, is the equivalent length l necessarily represented by the distance between the pivots? If the answer to the second question is assumed to be yes then the quadratic equation (3) gives precisely two possible good mass configurations since h depends linearly on the position x of the movable weight. These mass configurations can then be empirically obtained by the following standard procedure [1]: by varying the movable mass position x collect two series of data (x,t), one for each pivot,

3 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) make a parabolic fitting of the data by means of two parabolas of the following type: T = ax 2 +bx+c (5) these parabolas meet in at most two points (x 1,T), (x 2,T): positions x 1 and x 2 determine the two desired good mass configurations. Such a parabolic fitting is justified by two considerations. The first one is that we are looking for two good mass configurations, then the fitting curves have to admit at most two intersection points. The second one is the empiric observation of the data which apparently seem to be arranged just along two convex parabolas with vertical axis. This is what is usually done although the right answer to the second question should be no, aswas firstly pointed by Shedd and Birchby in 1907[2, 3, 4]. Their remark seems to have escaped general attention, perhaps due to the fact that, if the pendulum is well assembled, the previous parabolas meet at points whose abscissas give almost exactly the good mass configurations having the distance between pivots as equivalent length. The latter is much easier determined than any other equivalent length associated with further good mass configurations of the pendulum [5]! But what does mean well assembled? To fix ideas consider an ideal Kater pendulum consisting of an idealised massless rigid rod (x-axis) supporting two identical point masses, m f fixed at a and m m at a variable position x. The assembly has two distinct suspension points for the oscillations positionedat d/2andat+d/2, asketchofthis ideal pendulumisreportedinfigure1. Then d is the distance between the pivots, the center of mass is located at Fig. 1 Front view of idealized pendulum with two point masses and a massless rigid rod. b = am f +xm m m f +m m = x a 2

4 4 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 and the moment of inertia about the center of mass at b is given by I 0 = (b+a) 2 m f +(b x) 2 m m = m 4 (x+a)2 where m = 2m f = 2m m is the total mass of the pendulum. The moment of inertia I 1 and I 2 with respect to the two pivots are I 1 = ( a+x ) 2 m+( d 2 2 +b)2 m, (6) and I 2 = ( a+x ) 2 m+( d 2 2 b)2 m. (7) When x determines a good mass configuration the resulting periods T 1 and T 2 of small oscillations about the two pivots, respectively, have equal values. Equation (1) gives T 1 = 2π m(b+ d 2 )2 + m 4 (x+a)2 mg b+ d = 2π (x a+d)2 +(x+a) 2 2g x a+d 2 T 2 = 2π m(b d 2 )2 + m 4 (x+a)2 mg b d = 2π (x a d)2 +(x+a) 2 2g x a d 2 Then T 2 1 = T2 2 gives a cubic equation on the variable x. If it is assumed that (x a) 2 d 2 < 0 (8) which occurs, for instance, when suspension points are the end points of the pendulum bar, one finds that (x a) [ 2 ( x 2 +a 2) d 2] = 0 (9) Its solutions are then given by x = a, (10) d 2 x = ± 2 a2, (11) which represent all the possible positions of the movable weight giving a good mass configuration for the ideal Kater pendulum. The first solution, x = a, always exists. Furthermore, if d/ 2 a there are two additional positions which are symmetric with respect to the origin i.e. the middle point of the massless bar. Recall formula (2) to obtain the associated equivalent lengths. For the last two symmetric solutions it gives l = d But the equivalent length associated with the first solution is l = d 2 +2a2 d (12)

5 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) which in general does not coincide with the distance d between the two pivots. On the other hand if (8) is not assumed and we are in the more pathological case of a pendulum such that (x a) 2 d 2 > 0 then T1 2 = T2 2 reduces to a linear equation in the variable x whose solution is x = d2 4a and the associated equivalent length is l = a+ d2 4a which in general does not coincide with the distance d between the two pivots. Therefore for an ideal Kater pendulum the answers to the previous questions are: (1) there are at most three possible positions of the movable weight which determine a good mass configuration; (2) no; there always exists a good mass configuration whose associated equivalent length does not coincide with the distance between pivots. An immediate consequence is that a parabolic fitting of the empirical data (x,t) can t be the best fitting since two parabolas never meet at three points! Moreover in some particular cases a parabolic fitting may cause strong distortions in determining good mass configurations. For example: if either d/ 2 < a or (8) is not satisfied, the ideal Kater pendulum admits a unique good mass configuration; typically a parabolic fitting of data in this situation gives parabolas meeting only at imaginary points and the procedure stops; if a ±d/2 then the first solution of (9) is quite near to one of the two further symmetric solutions; a parabolic fitting of data gives only two intersection points but we do not know if one (and which one?) of them is nearer to the position associated with l than to the one associated with l ; in this situation also l l but they are not equal; then associating l with a so determined good mass configuration may cause a relevant error in the final value of g. One may object that we are discussing an empiric procedure by means of an ideal pendulum. In particular the position x = a for the movable mass gives the completely symmetric mass configuration with respect to the middle point of the ideal pendulum bar. When a physical pendulum with m f m m is considered, what is such a mass configuration? Does it occur again? The answer is yes. The key observation is that, for both pivots, the variable position x of the movable mass and the resulting period T of small oscillations are related by a cubic expressions of the following type (period-distance relations) ax 2 +bx+c = T 2 +dxt 2, (13) where the coefficients a,b,c,d depend on the pendulum parameters. This is precisely what Shedd and Birchby pointed out in their papers [2, 3, 4] giving theoretical and empirical evidences: they called the two (one for each pivot) equations (13) the equations

6 6 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 of the reversible pendulum (see equations (10) and (11) of their first paper). Here we will refer to (13) as the cubic period distance relation of the physical Kater pendulum considered. Note that only coefficients a, b, c, d depend on the pendulum parameters, whilethepolynomialtypeofequation(13)doesnotdependonthechoiceofthependulum. Thus, we can reduce the search for good mass configurations to a simple cubic equation similar to Eq. (9). A first point in the present paper is to give a mathematically rigorous proof of the following Theorem 1.1. Let p 1 (x,y),p 2 (x,y) be the following cubic polynomials p i (x,y) = A i x 2 +B i x+c i y 2 D i xy 2, i = 1,2. where A i,b i,c i,d i are real coefficients and D 1 D 2. Then they admit always two real common roots and two pairs of complex conjugated common roots which may be real under suitable conditions on coefficients A i,b i,c i,d i. Thinking them as points in the complex plane (x,y) they are symmetric three by three with respect to the x axis. Moreover these are all the common roots they can admit (that is: all the further common roots are at infinity ). This algebraic result leads to the following physical statement: Corollary 1.2. A physical Kater pendulum, with a sufficiently long bar, admits always a good mass configuration whose associated equivalent length does not in general coincide with the distance between the pivots. Under suitable conditions on the pendulum parameters, it may admit two further good mass configurations. They correspond to symmetric positions of the movable mass, with respect to the middle point of the bar (if also the pivots are symmetrically located). They admit a common associated equivalent length which is precisely the distance between the pivots. Moreover the pendulum can t admit any further good mass configuration. We will specify the meaning to the vague expressions sufficiently long and suitable conditions. Although Shedd and Birchby knew in practice the content of the previous statement (they actually wrote down all the three good mass configurations in period distance terms see formulas (27) of their first paper) they couldn t give a rigorous proof of it. They studied the geometry of the curves determined by the cubic period distance relations by means of an old and non standard Newton s classification. They then arrived to conclude that (see the bottom lines of p. 281 in their first paper): Of the nine possible intersections of two cubic curves, in the present case three are imaginary or at infinity, three belong to the condition that T is negative, and three belong to positive values of T, and can hence be experimentally realized.

7 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) This conclusion does not exclude that, under some suitable conditions on the pendulum parameters, at least two of the three imaginary or at infinity intersections may become real and maybe physically realizable giving more than three good mass configurations. Actually we will see that these three intersections are not imaginary but definitely at infinity and they can never give physical results. A second aim of the present paper is to observe that the best fitting of empirical data (x,t) is then given by two cubic curves of type (13) instead of two parabolas of type (5). We will support this remark by experimental evidence for a real compound Kater pendulum. The paper is organised as follows. Section 2 is devoted to prove Theorem 1.1. and the physical statement of Corollary Here we set the main notation and describe the physics of a real Kater pendulum. The proof of Theorem 1.1. is based on elementary elements of complex algebraic and projective geometry. Anyway a non interested reader may skip it without losing any useful element to understand what follows. In Sec. 3 we describe an effective experiment. Section 4 is devoted to the analysis of experimental data by a linear fit of the period-distance cubics. An estimate of their intersection points gives the good mass configurations and the associated periods for our real Kater pendulum. Then the value of g. A comparison with a parabolic fitting of data is then given. Appendix A is devoted to discuss the suitable conditions on the pendulum parameters under which the pendulum admits all the possible good mass configurations (see Corollary 1.2.). In Appendix B we collect some further numerical methods to analyse our empirical data. 2 Physics of the Kater reversible pendulum Notation. Consider a physical Kater pendulum composed of a rigid bar equipped with two weights (see Fig. 2 and Fig. 3). The pendulum can be suspended by two knife edges, c 1 and c 2, symmetrically located on the bar. The weight m f is placed in a fixed position which is not between the knives. The other one, m m, can be moved along the bar. Small oscillations of the pendulum are parameterised by an angle ϕ such that ϕ sinϕ, that is, ϕ 3 0. The equation of motion of the pendulum is then given by ϕ+ mgh i I i ϕ = 0, (14) where g is the earth s apparent gravitational acceleration, m is the total mass of the pendulum, h i is the distance of the center of mass from the knife edge c i, and I i is the moment of inertia about c i. The Steiner s theorem [6] asserts that I i = I 0 +mh 2 i, (15) where I 0 is the moment of inertia with respect to the center of mass. The associated

8 8 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Fig. 2 Detailed side view of the Kater pendulum (not to scale). Fig. 3 Front view of the Kater pendulum (not to scale). The pendulum swings in the plane of the picture; its pivot can be inverted. period of small oscillations is T i = 2π Ii = 2π = 2π I 0 +mh 2 i. (16) ω i mgh i mgh i

9 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) Equation (16) implies that the Kater pendulum oscillates with the same period as a simple pendulum whose length is given by l i = I i mh i = I 0 +mh 2 i mh i. (17) Assume now that the movable mass m m is placed at a point x 0 on the bar such that T 1 = T 2 = T(x 0 ). (18) Such a point will be called a characteristic position of the pendulum. Equation (18) can be satisfied if and only if l 1 = l 2 = l(x 0 ). The length l = l(x 0 ) will be called the characteristic length of the pendulum associated with the characteristic position x 0. Analogously the associated periods T(x 0j ) will be the characteristic periods of the pendulum. The knowledge of l(x 0j ) and T(x 0j ) for each j = 1,2,3 yields the value of g from the relation and therefore l T = 2π g, (19) g = 4π2 T 2. (20) The variable position of m m is described by a linear coordinate x having origin at c 1. Then c 2 is the point x = d > 0 (see Fig. 2) while the fixed weight is placed at x f such that (d L)/2 < x f < 0. The movable and fixed weights are composed of disks whose radii are given respectively by r m and r f. L is the length of the pendulum bar. Therefore, the distance h between the pendulum center of mass and the origin c 1 depends on the position x of m m and is given by: h = where m b is the mass of the bar. Set d 2 m b +x f m f +xm m m b +m f +m m, (21) m = m b +m f +m m, (22) and Then h can be rewrite as follows K = d 2 m b +x f m f m. (23) The moment of inertia I 0 is then given by I 0 = (h x f ) 2 m f +(h x) 2 m m + h = K + m m x. (24) m ( h d 2 ) 2 m b +I 0, (25)

10 10 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 where Set and I 0 can be rewritten as follows: I 0 = r2 f 2 m f + r2 m 2 m m + L2 12 m b. (26) ( ) 2 d I 0 = I 0 +m f (x f K) 2 +m b 2 K +m m K 2, (27) I 0 = m m m m m m x2 2m m Kx+I 0. (28) From Eq. (17) the condition (18) is satisfied if and only if which is equivalent to requiring that From Eq. (24) we have I 0 +mh 2 1 mh 1 = I 0 +mh 2 2 mh 2, (29) (h 1 h 2 )(mh 1 h 2 I 0 ) = 0. (30) h 1 = h = K + m m m x (31) h 2 = d h = d K m m x, (32) m and we get the first characteristic position by imposing that h 1 = h 2, that is, x 01 = d 2 + m f 2m m (d 2x f ). (33) Two additional characteristic positions can be obtained by the second factor in Eq. (30). By letting mh 1 h 2 I 0 = 0 and expressing I 0 as in Eq. (28), we have x 2 dx mk2 mdk +I 0 m m = 0, (34) whose solutions are x 02 = d d mk2 mdk +I 0 (35) m m x 03 = d d 2 +4 mk2 mdk +I 0 m m. (36) To determine the associated characteristic lengths l(x 0j ) use Eqs. (17), (29), and (31). It follows that l(x 0j ) = I 0 +m ( K + mmx ) 2 m 0 j. (37) mk +m m x 0j

11 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) It is then easy to observe that l(x 02 ) and l(x 03 ) are equal and constant because x 02 and x 03 are symmetric. Precisely l(x 02 ) = l(x 03 ) = h 1 +h 2 = d, (38) and they do not depend on the other physical parameters of the pendulum. On the contrary, this is not true for l(x 01 ) because l(x 01 ) = d 2 I 0 +2 md + m f (m m +m f )(d 2x f ) 2. (39) 2m m md The reader may compare the characteristic positions (33), (35) and the associated characteristic lengths (39), (38), now obtained, with those given in Eq. (27) by Shedd and Birchby[2]. Moreover the period distance relations of the pendulum (what Shedd and Birchby called the equations of the Kater pendulum [2]) can be obtained by Eq. (29) when h 1 and h 2 are expressed as in Eqs. (31). When the pendulum oscillates about the pivot c i, the period T i and the distance x results to be related by the following cubic relations where A i x 2 +B i x+c i = T 2 i +D i xt 2 i, i = 1,2, (40) and A 1 = 4π2 m m gmk (41) B 1 = 0 C 1 = 4π2 ( I gmk 0 +mk 2) (42) D 1 = m m mk, (43) A 2 = 4π2 m m gm(d K) (44) B 2 = 8π2 m m d gm(d K) (45) 4π 2 ( C 2 = I gm(d K) 0 +m(d K) 2) (46) D 2 = m m m(d K). (47) All the possible characteristic positions are then given by all the common roots of Eqs. (40). Proof of Theorem For more details on the mathematics here involved see, for instance, Harris[7] or Shafarevich [8] among other introductory textbooks on algebraic geometry.

12 12 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Consider (x,y) as coordinates of points in the complex affine plane C 2. Then equations p 1 (x,y) = 0 and p 2 (x,y) = 0 give two cubic complex algebraic curves, C 1 and C 2, whose intersection points are precisely the common roots of p 1 and p 2. We can compactify C 2 by adding a line at infinity : this procedure produces the complex projective plane P 2 C. More precisely, we can consider our complex variables x and y to be a ratio of further variables, that is, x = X and y = Y Z Z. (48) The equations defining C 1 and C 2 multiplied by Z 3 become the following A 1 X 2 Z +B 1 XZ 2 +C 1 Z 3 = Y 2 Z +D 1 XY 2 (49) A 2 X 2 Z +B 2 XZ 2 +C 2 Z 3 = Y 2 Z +D 2 XY 2. (50) which are the defining equations of the projective completions C 1 and C 2, respectively. The main ingredient of the present proof is the following Theorem 2.1. (Bezout) Given two distinct irreducible complex algebraic plane curves of degree d 1 and d 2, their projective completions admits a finite number of intersection points. Precisely if every intersection point is counted with its algebraic multiplicity then this number is d 1 d 2. In particular the projective completions C 1 and C 2 meet in 9 points, counted with their algebraic multiplicities. The Bezout theorem is a consequence of the Fundamental Theorem of Algebra which asserts that on the field C of complex numbers every polynomial admits as many roots as its degree. The first step is to study the intersections at infinity, that is, which belong to the added line at infinity. The equation of this line is Z = 0 and by Eq. (49) it intersects both our cubics at y (that is, the point X = Z = 0 which is the infinity point of the affine y-axis x = 0) and in x (that is, the point Y = Z = 0 which is the infinity point of the affine x-axis y = 0). Both of these are inflection points for C 1 and C 2. In y the inflection tangent line of C 1 is given by t 1 : D 1 X +Z = 0, while the inflection tangent line of C 2 is t 2 : D 2 X +Z = 0. They cannot coincide since D 1 D 2. Therefore y is a simple intersection point of our cubics, that is, it admits intersection multiplicity 1. On the other hand, in x both C 1 and C 2 have the same inflection tangent line which is the infinity line Z = 0. Then x has intersection multiplicity 2. Consequently these infinity points count 3 of the 9 intersection points. The remaining 6 intersections must be affine, that is, they cannot belong to the compactifying line at infinity.

13 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) To find them note that for i = 1,2 t i C i = y, (51) with intersection multiplicity 3 because it is an inflection point for C i with tangent line t i. On the other hand t 1 C 2 = {y,p 1,P 2 } (52) t 2 C 1 = {y,q 1,Q 2 } (53) where P h =y, Q k =y, and P h =Q k, because t 1 and t 2 are always distinct. Therefore, the affine intersection points of C 1 and C 2 cannot belong to the lines t 1 and t 2, and they can be recovered by studying the common solutions to the following equations y 2 = A 1x 2 +B 1 x+c 1 1+D 1 x (54) y 2 = A 2x 2 +B 2 x+c 2 1+D 2 x, (55) because those points do not make the denominators vanish. So they are reduced to the roots of the following cubic equation (A 1 x 2 +B 1 x+c 1 )(1+D 2 x) = (A 2 x 2 +B 2 x+c 2 )(1+D 1 x). (56) It is a cubic equation with real coefficients. Therefore it admits 3 complex roots one of which is surely a real number. The remaining two roots are necessarily complex conjugated: their reality depends on the coefficients A i,b i,c i,d i. Proof of Corollary Recall the cubic period distance relations (40). Setting T 1 = T 2 = y they are represented by the two cubic curves C 1 and C whose coefficients are assigned by formulas (41) and (44), respectively. Note that they are real numbers and D 1 D 2 since d 0. The hypothesis of Theorem 1.1. are then satisfied and the characteristic positions of the pendulum must be represented by the real affine intersection points admitting y 0. To conclude the proof observe that Eq. (56) divided out by 4π 2 /gmk(d K) gives exactly the cubic equation (30). The real root is then given by (33) and it always occurs when m f m m d 2x f L The remaining two roots are then assigned by (35). A discussion of their reality is given in Appendix A. 3 The experiment The physical parameters characterising our pendulum are given in Table 1; the digits in parentheses indicate the uncertainties in the last digit.

14 14 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Table 1 The physical parameters characterising the pendulum. m m (g) m f (g) m b (g) x f (cm) l(cm) d(cm) r f (cm) r m (cm) 1399(1) 1006(1) 1249(1) 26.73(1) 167.0(1) 99.3(1) 5.11(1) 5.12(1) The bar length is measured by means of a ruler whose accuracy is ±1mm. The radii r m and r f and the position x f are measured by a Vernier caliper accurate to ±0.1mm. The masses m b, m m, and m f are determined by means of a precision balance accurate to one gram. With reference to the structural conditions in Appendix A, we are in the case 3.b i.e. all the three characteristic positions occurs and precisely x 02,x 03 are placed between the knives while x 01 is on the opposite side of the bar with respect to m f. By recalling Eqs. (33) and (35), we expect that with associated characteristic lengths x 01 = (104.57±0.11)cm (57) x 02 = (61.74±0.40)cm (58) x 03 = (37.56±0.31)cm, (59) l(x 01 ) = (121.44±0.09)cm (60) l(x 02 ) = l(x 03 ) = d = (99.3±0.1)cm. (61) Throughout the experiment the movable mass m m will be placed in successive positions, generally 10 cm from each other, except near the theoretical characteristic positions (57) where the distances decrease (see the second column in Table 2). The period of small oscillation about the two pivots are measured for all those positions of m m. These periods are measured by recording the time of each of 9 consecutive oscillations when the pendulum starts from the angle ϕ 0 6 ±1. For this purpose we used a photogate timed by an electronic digital counter. We repeated the procedure for 18 positions of m m, at first with respect to c 1 and then c 2. The average of the 9 values is taken to be the period at the given position of m m whose error is given by half of its maximum excursion, that is, s. The initial angle ϕ 0 is sufficiently small that an equation similar to Eq. (14) is valid. By expanding an elliptic integral in a power series, it is possible to approximately express the period associated with the exact equation of pendulum motion ϕ+ mgh i I i sinϕ = 0 (62) by adding corrective terms,[9, 6] to the period expression given in Eq. (19). In the next section we will evaluate such a correction. The results are reported in Table 2. We did not choose positions too close to the estimated characteristic positions to prevent the casual occurrence of coincident period measures about the two pivots. In fact our distance measures are effected by an error of ±1mm. Such an error would cause a strong distortion in determining the empirical characteristic positions. One of them would be directly determined by direct measure and its error would not be lessened by the fitting procedure. The resolution of the LEYBOLD-LH model is 0.1ms.

15 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) x(cm) T 1 (s) T 2 (s) Table 2 The experimental data The linear fitting procedure We now describe a linear fitting procedure used to fit the experimental data listed in Table 2 and empirically determine the characteristic positions. The numerical computations were obtained using MAPLE and some FORTRAN code. From a numerical point of view we should fit the data by cubic polynomials like those in Eq. (40). Such a fitting can be treated linearly because the coefficients D 1 and D 2 may be determined a priori by Eqs. (41) and (44) which involve only the known physical parameters listed in Table 1. We obtain D 1 = (3.983±0.01) 10 2 cm 1 (63) D 2 = ( ±0.0047) 10 3 cm 1. (64) We used Maple V, Release 5.1 by Waterloo Maple Inc. The FORTRAN codes employed subroutines from Ref. [14] and the numerical package NAG-Mark 14. The plotting package is PGPLOT 5.2 developed by T. J. Pearson.

16 16 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Table 3 Coefficients of the cubic curve C 1 estimated by the linear method. A 1 ( ± )s 2 cm 2 B 1 0 s 2 cm 1 C 1 (7.641±0.011)s 2 Table 4 Coefficients of the cubic curve C 2 estimated by the linear method. A 2 ( ± )s 2 cm 2 B 2 ( ± )s 2 cm 1 C 2 (4.393±0.01) s 2 Table 5 χ 2 and critical p values for linear fitting by cubic curves C 1, degrees of freedom = 16 ;C 2 degrees of freedom = 15 χ χ 2 0 χ 2 (x,15)dx We can obtain the desired fitting of the data obtained in Sec. 3 by applying the least squares method to the following function Ξ i (A i,b i,c i ) = 18 h=1 T 2 h,i A ix 2 h +B ix h +C i 1+D i x h 2T h,i σ T 2, (65) where (x h,t h,i ) are the data of the ith set in Table 2. Two sources of error with period measurements were considered In any position and for both pivots, we considered the standard deviation of the 9 period electronic measurements, varying from s to s. The systematic error that formula 70 introduces on the data. For example when T=2.3s (the maximum period here analysed) and ϕ 0 12 the shift introduced on T is 0.006s. Afterthisanalysisweconsidered σ T = 0.006sasestimatederrorforperiodmeasurements. The obtained results are reported in Table 3, Table 4 and visualised in Fig.4. The merit function χ 2 and the associated p values are reported in Table 5 and each of them have to be understood as the maximum probability to obtain a better fitting. The estimated cubic coefficients of Table 3 and Table 4 allow us to evaluate the characteristic positions and the associated characteristic periods by intersecting their We also know that B 1 = 0 and the number of coefficients to be estimated by the fitting procedure can be reduced.

17 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) Fig. 4 Theoretical cubics (dotted line), fitted cubics (full line). and experimental data (filled points). The experimental errors are much smaller than the filled points drawn, so they are not visible within this plot. Table 6 Estimated intersection points of fitting cubic curves. (x 01,T(x 01 )) ( cm, s) (x 02,T(x 02 )) (62.541cm, s) (x 03,T(x 03 )) (35.779cm, s) upper branches We obtain a cubic equation whose numerical solutions are reported in Table 6. Refer to Eqs. (20) and (60) to compute the associated values of g. We have g 1 = 4π 2 l(x 0 1 ) T(x 01 ) 2 (66) g 2 = 4π 2 l(x 0 2 ) T(x 02 ) 2 (67) These two cubic curves represent the period-distance relations in the plane (x,t) when oscillations are considered about c 1 or c 2 respectively. Then their common points coordinates give the characteristic positions of the pendulum and the associated periods. We have already observed in Sec. 2 that these cubics are symmetrical with respect to the x-axis. More precisely each of them is composed by two symmetrical branches. The branches lying under the x-axis are not physically interesting since their period coordinate T is negative. Therefore the only interesting common points of these two cubic curves are the intersection points of their upper branches.

18 18 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Table 7 Values of g obtained by formulas (66) and (69). g 1 (974.15±2.72) cms 2 g 2 g 3 (982.65±3.11) cms 2 (980.20±3.1) cms 2 g (979.00±1.72) cms 2 g 3 = 4π 2 l(x 0 3 ) T(x 03 ) 2, (68) and their numerical values are listed in Table 7. Their average gives g = (979.00±1.72)cms 2. (69) where the uncertainty is found implementing the error propagation equation (often called law of errors of Gauss) when the covariant terms are neglected(see equation(3.14) in[11]). We now consider the correction arising from the approximation of the exact equation of pendulum motion (62) already mentioned at the end of the previous section. This correction gives:[9, 6] l ( T = 2π ) 2 1+1/16ϕ0, (70) g and g = 4π 2 l ( ) /16ϕ0. (71) T 2 A small increase in the value of g is evident from Eq. (71), and we will refer to it as the finite amplitude correction (f.a.c.). With the data listed in (7) we obtain g + = (980.34±1.74)cms 2, (72) which is the gravity acceleration increased by the f.a.c.. An accurate measure of the value of g in Turin[10] gives g T = (4)cms 2. (73) This value will be considered as the true value of the acceleration due to the earth s apparentgravityfieldinturin. Bycomparingitwithg +, weseethatourmeasurement is 191ppm smaller than the true value. Note that the considered Kater pendulum admits characteristic positions sufficiently distant from each other (see formulas (57) and data collected in Table 6). Then it can be Further references for accurate measurements of g are the following. A world-wide survey of all the apparent gravity measurements (see < gives for Turin g = cms 2 ; this value differs from g T by 16ppm. An analytical formula provided by the U.S. Geological Survey [12] needs two input parameters, the local height above sea level and latitude, which in our case are 236m and respectively to give the local value of g. For Turin this formula gives g = cms 2, which differs from g T by 61ppm.

19 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) Table 8 Coefficients of the first fitting parabola P 1. A 1 ( ± )s cm 2 B 1 ( ± )s cm 1 C 1 (2.494±0.003)s Table 9 Coefficients of the second fitting parabola P 2. A 2 ( ± )s cm 2 B 2 ( ± )s cm 1 C 2 (2.113±0.004)s Table 10 χ 2 and critical p values for linear fitting by parabolas, first 13 data, degrees of freedom =10 P 1 P 2 χ χ 2 0 χ 2 (x,15)dx considered sufficiently well assembled, which means that a parabolic fitting (of type (5)), of the empirical data (x h,t h,i ) collected in Table 2, should give a sufficiently precise evaluation of characteristic positions x 02,x 03. As before we apply the least square method to the following function, but the sum is now extended to the first 13 entries of Table 2 in order to exclude the first intersection Θ i (A i,b i,c i ) = 13 h=1 ( Th,i (A i x 2 h +B ) 2 ix h +C i ). (74) The coefficients of the fitting parabolas are reported in Tables( 8) and ( 9) ; the χ 2 and the associated p values are reported in Table 10, in comparison the parabolic fit gives very bad results.. Their intersection points are given in Table ( 11). We get then the following two evaluations of g: σ T g = 4π 2 d T(x 02 ) = 2 (960.08±3.00)cms 2 (75)

20 20 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Table 11 Estimated intersection points of fitting parabolas. (x 02,T(x 02 )) (72.296cm : s) (x 03,T(x 03 )) (41.709cm : s) Their average gives g = 4π 2 d T(x 03 ) = 2 (989.11±3.14)cms 2. (76) g parabolic = (974.60±2.18)cms 2. (77) A comparison with g in (69) and g T in (73), gives a clear evidence of the better efficiency of a cubic fit with respect to a parabolic one. 5 Conclusions We summarise the main results of our theoretical and numerical analysis (1) The three solutions of the Kater pendulum concerning the distance period relationship discovered in 1907 by Shedd and Birchby in 1907 [2, 3, 4] are classified in a modern context (2) The first solution of the distance period relationship allows to deduce a new formula of g via the second equivalent length both in the idealized pendulum and in a commercial Kater pendulum, see respectively formula(12) and (39) (3) One of the main targets of our work the evaluation of g gives oscillating results : our best numerical fit to T 2 produces (the linear fit + non linear correction) a value of g that is 191 ppm smaller than the true vale : our worst fit to T 2 (the non linear fit + non linear correction) gives a value of g that is 1978 ppm smaller than the true vale (4) Concerning the fit to T through a parabola we obtain high values of χ 2 (χ 2 =893 for C 1 and χ 2 =12.51 for C 2 ) with respect to the linear fit to T 2 (χ 2 =2.69 for C 1 and χ 2 =1.57 for C 2 ). These high values of χ 2 allow to rule out the physical significance of this type of fit. Acknowledgements We would like to thank G.Maniscalco for his assistance in the preparation of experiment set up. We are also grateful to the anonymous referee for useful suggestions and improvements.

21 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) A Reality of characteristic positions and structural conditions The characteristic positions of our pendulum are given by Eqs. (33) and(35). The former, x 01, is always real. On the other hand x 02 and x 03 are real if and only if the square roots in Eq. (35) are real i.e.e if and only if m m d 2 +4mdK 4mK 2 4I 0 0 x 2 f dx f (m m +m b )d 2 4I 0 4m f 0. (A.1) The latter are the suitable conditions on the pendulum parameters of Corollary To avoid the overlapping of m f with c 1 we have to impose that x f r f. Then Eq. (A.1) is equivalent to requiring that md 2 4I 0 0 and d 2 1 md2 4I 0 x f r f. (A.2) 2 m f Note that the condition on the right in Eq. (A.2) is not empty if d 2 1 md2 4I 0 r f 2 m f d 2 m fr f ( ) m f r 2 f m f r + + f 2 +I 0. (A.3) m b +m m m b +m m m b +m m In particular, the latter ensures that the left condition in Eq.(A.2) is also satisfied because I md 2 4I 0 0 d 2 0 m. (A.4) Toavoidtheoverlapping ofm m withtheknife edges itfollowsthateither r m x d r m or d+r m x L+d. After some algebra we get the following results. 2 Assume that m m > m f and set M 1 = 2 m mr m +m f r f m m m f M 2 = 2 m fr f + m b +m m ( m f r f m b +m m ) 2 + m f r 2 f +I 0 m b +m m (A.5) (A.6) S 1 = m f m m d+ m m r m (A.7) 2m f m f S 2 = d 2 1 md2 4I 0 (A.8) 2 m f S 3 = d 2 1 m fd 2 +m b d 2 +4m m r m d 4m m rm 2 4I 0, (A.9) 2 m f Then we have the following possibilities:

22 22 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 (1) d < min(m 1,M 2 ): in this case the pendulum admits only one characteristic position given by x 01 because by Eqs. (A.2) and (A.3), x 02 and x 03 are not real; x 01 is not between the knives, but occurs on the opposite side of the bar with respect to m f ; the system is in an almost symmetrical mass configuration of the pendulum. (2) For min(m 1,M 2 ) d < max(m 1,M 2 ), we have the following possibilities: (2a) If M 1 < M 2 : we get only the characteristic position x 01 which is between the knives if and only if S 1 x f r f ; otherwise, x 01 is placed like in (1). (2b) If M 2 < M 1, we get all the characteristic positions; x 01 is like in (1) and x 02,x 03 occur between the knives if and only if S 2 x f min(s 3, r f ). (3)Formax(M 1,M 2 ) d,thependulumadmitsallthecharacteristicpositionsx 01,x 02,x 03 which are placed as follows: (3a) Only x 01 is placed between the knives when either S 1 < S 2 and S 1 x f S 2, (A.10) or S 3 < r f and max(s 1,S 3 ) < x f r f. (A.11) (In particular, if S 3 < S 1, we can also assume the position x f = S 1 for the fixed weight m f.) (3b) only x 02,x 03 are placed between the knives when S 2 < S 1 and S 2 x f < S 1. (3c) We obtain all the possible characteristic positions x 01,x 02,x 03 between the knives when max(s 1,S 2 ) x f min(s 3, r f ). (A.12) In the concrete case considered in the Section 3 we have M 1 = (62.61±0.25)cm M 2 = (70.87±0.05)cm S 1 = ( 12.28±0.08)cm S 2 = ( ±0.075)cm S 3 = ( 7.62±0.06)cm, (A.13) (A.14) (A.15) (A.16) (A.17) where the uncertainty is found by applying the law of errors of Gauss with the uncertainties listed in Table 1. Therefore we are in the case 3.b. B Further numerical methods We outline three additional numerical methods that may be applied to analyse experimental data. The final results obtained by means of each method are reported in Table B.1. B.1 The non-linear method In the fitting procedure of data reported in Table 2, all the coefficients A i,b i,c i and D i are considered as unknown parameters to be estimated. Therefore a fitting procedure

23 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) Table B.1 Average values g and corrected values g + (by f.a.c.). algorithm g g + linear fitting by parabolas (974.6±2.17)cms 2 (975.93±2.20)cms 2 linear fitting by cubics (979.00±1.72)cms 2 (980.34±1.74)cms 2 non-linear fit (977.25±1.71)cms 2 (978.25±1.74)cms 2 Cramer interpolation (980.06±4.88)cms 2 (981.40±4.89)cms 2 Spline interpolation (979.52±1.73)cms 2 (980.86±1.74)cms 2 performed by means of cubic polynomials like those in Eq. (40) is necessarily a non-linear one. We want to apply the least square method to minimise the following functions ( 18 X i (A i,b i,c i,d i ) = Th,i 2 A ix 2 ) h +B i x h +C i, (B.1) 1+D i x h h=1 which are non-linear in the unknown coefficients. The procedure is to apply the NAG- Mark14 subroutine E04FDF to find an unconstrained minimum of a sum of 18 nonlinear functions in 4 variables (see Ref. [13]). The final value of g is reported in the Table B.1. B.2 The Cramer interpolation method We present here a method that reduces our analysis in a local neighbourhood of the estimated characteristic positions where a cubic behaviour of the fitting curves is imposed. From Eqs. (41), (44), and (63) we know that D 1 and D 2 are completely determined by the pendulum parameters. Moreover, from Eqs. (41) we know that B 1 = 0. So to recover the remaining coefficients of C 1 and C 2, we need to interpolate two points of the first set of data in Table 2 and three points of the second one respectively. We have to solve a 2 2 and a 3 3 linear system by applying the Cramer theorem (which is the most practical method for solving a linear system of equations). If we choose data points that are close to a characteristic position, then the nearest point in C 1 C 2 to the chosen data will give an empirical estimation of such a characteristic position and its associated period. An iterated application of this procedure will produce a distribution of periods and we may obtain g from the mean value and its statistical error from the standard deviation (see the last line of Table B.2). The results obtained for every interpolation are reported in Table B.2; the chosen data points in the second and third columns are enumerated as they appear in Table 2. B.3 Cubic Spline Interpolation The last data analysis method to be proposed is the cubic spline interpolation (subroutine SPLINE and SPLINT from Numerical Recipes II). Once the three intersections are

24 24 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) 1 25 Table B.2 The Cramer interpolation method. Cramer Method Char. Chosen data Intersection Char. g position Series 1 Series 2 Position Period length x 01 15;18 15;17; cm t1,1= s cm cms 2 x 01 15;17 15;16; cm t1,2= s cm cms 2 x 01 16;18 16;17; cm t1,3= s cm cms 2 x 01 15;18 15;16; cm t1,4= s cm cms 2 x 02 7; 9 7; 8; cm t2,1= s 99.30cm cms 2 x 02 8;10 8; 9; cm t2,2= s 99.30cm cms 2 x 02 9;11 9;10; cm t2,3= s 99.30cm cms 2 x 02 10;12 10;11; cm t2,4= s 99.30cm cms 2 x 03 2; 4 2; 3; cm t3,1= s 99.30cm cms 2 x 03 3; 5 3; 4; cm t3,2= s 99.30cm cms 2 x 03 4; 6 4; 5; cm t3,3= s 99.30cm cms 2 x 03 5; 7 5; 6; cm t3,4= s 99.30cm cms 2 g=( ± 4.89)cms 2 obtained, the procedure was similar to the linear/nonlinear case and the final value of g is reported in Table B.1. References [1] D. Randall Peters, Student-friendly precision pendulum, Phys. Teach., Vol. 37, (1999),pp [2] J.C. Shedd, J.A. Birchby, A study of the reversible pendulum. Part I. Theoretical considerations, Phys. Rev. (Series I), Vol. 25, (1907),pp [3] J.C. Shedd, J.A. Birchby, A study of the reversible pendulum. Part II. Experimental verifications, Phys. Rev. (Series I), Vol 34, (1912), pp [4] J.C. Shedd, J.A. Birchby, A study of the reversible pendulum. Part III. A critique of captain Kater s paper of 1818, Phys. Rev. 1, Vol. 457, (1913), pp [5] D. Candela, K. M. Martini, R. V. Krotkov, K. H. Langley, Bessel s improved Kater pendulum in the teaching lab, Am. J. Phys., Vol. 69, (2001, pp [6] R. Resnick, D. Halliday, K. S. Krane, Physics,John Wiley & Sons, New York, [7] J. Harris, Algebraic Geometry, Springer-Verlag, New York, [8] I. R. Shafarevich, Basic Algebraic Geometry, Springer-Verlag, New York,1977.

25 M. Rossi and L. Zaninetti / Central European Science Journals 1 (2003) [9] R. A. Nelson, M. G. Olsson, The pendulum-rich physics from a simple system, Am. J. Phys.,Vol. 54,(1986), pp [10] G. Cerutti, P. DeMaria, Misure assolute dell accelerazione di gravità a Torino, In: Technical Report R432, Istituto di Metrologia G. Colonnetti, Turin,1996. [11] Philip R. Bevington, D. Keith Robinson, Data Reduction and Error Analysis for the Physical Sciences, McGraw-Hill, Inc., New York, [12] P. Moreland, Improving precision and accuracy in the g lab, Phys. Teach., Vol. 38, (2000). [13] NAG, [14] Numerical Recipes,