text:
The Pendulum as a Useful Example: Action-Angle Variables and Perturbation Theory
- Karin Carpenter
- 2 years ago
- Views:
1 The Pendulum as a Useful Example: Action-Angle Variables and Perturbation Theory. The Simple Pendulum: Setup and Motion in Terms of Action Angle Variables We consider a simple pendulum, of mass m and length l, swinging under gravity from a single, fixed point at one end, assuming all the mass at the other end. Measuring an angle, θ, upward from its rest position hanging straight down, to whatever location it happens to have at some time t, we may write the following, where we normalize the potential energy so that it is zero when the pendulum is at right angles to that rest position, i.e., when θ = π/. However, since we intend to use this as a demonstration example for one degree of freedom, we will refer to the angle θ as q: V = mgl cos q, E = ml q + V, = p = ml q = H = H(q, p) = p ml mgl cos q = G p F cos q, d dt q = H p = p ml G p, = d dt p = H q = mgl sin q F sin q, (.) where we have given shorter names to the arrangements of constants that show up in the problem. In principle we can resolve this problem by direct integration, involving elliptic functions. One approach gives the equations for the motion in the form: sin(q/) = k sn(u; k), p = F/G k cn(u; k), } u ω t + ζ, E F, sin(q/) = sn(u; /k), p = F/G k dn(u; /k), } u ω k t + ζ, E F, (.) ω F G, k ( + E/F ) = ( + E/mgl), where ω is the natural frequency in the problem, i.e., the one it would have if the system responded linearly, and ζ is simply a constant of integration to be determined by initial conditions, and sn(u; k), etc. are the (usual) Jacobi elliptic functions. On the other hand, it is conceivable that your knowledge of the elliptic functions is not everything I might wish it were, and we want to use this as an example to better see how to do such things when the integrals are even more difficult to perform. Therefore, we now proceed to finding the appropriate action angle variables.. The first step is to solve the equations so as to obtain p = p(q, Q). In our case the only Q we have is the energy, E, which is the same as the Hamiltonian, so we have p = ± (E + F cos q). G (.3a)
2 . We may then find the time-independent HJ generating function, which depends on q and Q: q W = dq p = dq (E + F cos q). G (.3b) 3. Along the same lines we find the action, by performing the same integral, but by doing it over either an entire cycle, if the motion is closed, or over a period it the motion is simply periodic in the coordinate. For our problem, the determining factor is the energy, which is either greater than or less than F : J = π dq p = π dq G (E + F cos q) = qmax dq (E + F cos q), π G (.3c) where q max will be either π for a rotational motion, i.e., an unbounded but periodic one, or cos ( E/F ) for a libratory one, i.e., a bounded motion. 4. In order to determine the angular coordinate, we need W = W (q, J) = W [q, Q(J)], which we so far have only in the form on an integral. Nonetheless, we may proceed from that as follows, since Q depends only upon J and not q: φ = W J = Q W [q, Q(J)] = J J = G ( J/ E) q dq W Q = (E + F cos q) G ( ) J q dq (E + F cos q) E E G (.3d) where we have switched, part-way through, from Q to E because that it is the only relevant Q for this problem with only one degree of freedom, and also we have inverted the derivative because it is already J that we know, at least via the quadrature above in the form of a definite integral, as a function of the constant E. These two equations give us J = J(Q) and φ = φ(q, Q); we may of course resolve the first equation to obtain Q = Q(J) and insert it into the form for φ if desired. 5. We may then determine the frequency of the motion in φ by calculating H/ J. As H is just E, this is really just the inverse of the derivative we already needed to determine φ, namely /ω = J/ E. We will characterize it by comparing it to the natural frequency in the problem, namely ω, as given in Eqs. (.), where in this context it is reasonable to look at ω as the frequency for linearized oscillation about the elliptic singular point, at (q, p) = (, ): ω ω = F G dj de = F π G qmax dq. (.3e) (E + F cos q) G Remembering now the basic structure of the phasespace for the pendulum, we recall that it has this elliptic singular point at (, ) and a hyperbolic singular point at (±π, ) (where proper identification of the points in the phasespace to account for the angular nature of the coordinates just identifies these two coordinate presentations as the same point). There are elliptic trajectories in phase space surrounding the elliptic point, separatrix solutions that go through the hyperbolic point and which separate the elliptic trajectories from those which are hyperbolic (unbounded)
3 trajectories lying on the other side. When we consider the limit toward the separatrix trajectories from the unbounded side, we have q max = π, then in the limit as E F it is clear that the integral, for ω /ω T, becomes infinite, where T is the period of the motion, as is desired for motion along a separatrix. Contrariwise, in the limit from the elliptic side, we simply have the limits in the opposite order since E is increasing up to F, so that q max is becoming defined by cos(q max ) =, with the same result. We can obtain the limiting equations without too much effort by simply inserting E = F into the original equations. Beginning with Eq. (.3a) for p = p(q, E), and denoting the variables along the separatrix trajectory by a subscript x, we have the following, where we have started the motion with φ() = : p x = ± ω ω ( + cos q) = ± G G cos(q/) d dt q x = Gp x = ω t = log{tan[(q x + π)/4]} = q x = π + 4 tan (e ωt ), (.4) from which we can see the expected behavior of the motion along the separatrix, i.e., that it takes infinite time to move to the singular point at π = φ(+ ), i.e., it is moving very slowly when near the singular point. It is reasonable to expect, and also true, that motion along curves very near the separatrix itself have similar behavior. In particular, when we have more degrees of freedom so that a resonance between the frequencies can occur, there will always be some (new) variable that is a very slowly varying function of the time. For instance, suppose that the two frequencies ω and ω is near a rational number, say r/s; then the new variable rφ sφ will be such a slowly varying angle. If one then averages over the fast moving angles, the resulting problem will be -dimensional, and very much like the pendulum one now being studied. Finally, to acquire more details we will have to actually determine at least the integral form for J(E). Lichtenberg and Lieberman give the value of these integrals, in terms of elliptic integrals of the first and second kinds, which are basically inverse functions of the Jacobi elliptic functions: J = 8 π { F E(k) k K(k), k <, G k E(/k), k >, φ = π F (η, k), k <, K(k) F (q/, /k), k >, K(/k) (.5a) where we have appropriate definitions: K(k) F (π/, k) = u = F (w, k) = sn (sin w), k sin η sin(q/), π/ dα π/ k sin α, E(k) E(π/, k) = dα k sin α. (.5b) Using these same integrals, we may also determine the frequency more explicitly, from Eq. (.3e) above: { ω ω = K(k), k, (.5c) π (/k)k(/k), k. 3
4 It is also useful to have an understanding of the behavior of the first complete elliptic integral for arguments, k, that are not too large, and not too far from : [ ] (n )!! K(k) = (π/)[ + k /4 + 3(k /8) k n n n! n +..., K(k) log(4/ ( ) k ) + k 4 log(4/k ) ( k ) + O(k 4 ), k k. (.5d). The quartic oscillator, and secular perturbation theory: We now re-consider the pendulum Hamiltonian, but expand the cosine in a series and maintain only the first three terms; if, as well, we set l = and throw away the constant, firs term in the expansion of the cosine, we have H = p /m+ mω q + 4 ɛmq4 H + ɛ H, = q + ω q + ɛ q 3. (.) Here we have taken the coefficient of q 4 so that it has an explicit form, with ɛ = g/6. We will use this ɛ as a way of keeping track of terms in our perturbative expansion, and may re-set it back to g/6 when that is convenient. Of course we understand the problem of motion when ɛ is zero, i.e., when we just have the system reacting as the standard simple harmonic oscillator. However, let us see what happens when it is not. We assume that the solution may be written as a power series in ɛ: q = q(t) = q (t) + ɛ q (t) + ɛ q (t) (.) It is straightforward to insert this into the equations of motion and separate out quantities term by term: q + ω q =, q + ω q = q 3, q + ω q = 3 q q, q + ω q = 3q (q + q q ), q 4 + ω q 4 = (6q q q + 3q q 3 + q 3 ),.... (.3a) We need boundary conditions in order to solve even this; therefore, let us look for that solution which has value q() = a and q() =, i.e., it starts out from rest at location a. To make this work with the above series we now take q k () = δ k a, q k () =. (.3b) We then proceed to solve the equations in the most naïve way possible. First, we have the wellunderstood SHO: q (t) = a cos ω t. (.4a) 4
5 We insert this value into the second equation, causing it to become inhomogeneous, and find the particular solution that fits the boundary conditions: [ ] q (t) = a3 3ω cos ω t cos 3ω t + ω t sin ω t. (.4b) It is obvious that this is not a good approximation to the behavior of the solution of this equation, at least for long times, since it grows more or less linearly with ω t. It might well be that this is simply the first term in some sort of a series expansion of a periodic function that should multiply that term, just as the expansion of cos ɛt would have terms like that in it, were it expanded in a power series in its argument. However, this says that this is not a good approach to finding an approximate solution to the dynamics. Mathematically the reason for the problem is that there is a resonance between the frequency of the inhomogeneous term, which can be thought of as a driving term, and the natural frequency of the homogeneous equation; in this case they are exactly equal, although such a problem would occur even if their frequencies were commensurate, i.e., if their ratio was a ratio of integers. On the other hand, the physical reason for the problem is that we have acted as if the frequency of the exact solution, which we are trying to approximate, was (at least) a harmonic multiple of the original frequency; it is very unlikely that this is the case, as we can verify by thinking on our discussion of the full pendulum equation. From whichever point of view, these terms are bad for our solution; they are usually referred to as secular terms, a name which refers to terms which, on average, increase linearly with time. An approach to dealing with such problems was invented by Lindstedt, in 88. Since we agree that the frequency must also depend on ɛ in some way; we expand the frequency in powers of ɛ, and we insist that the (desired) solution be periodic with that frequency by considering it, now, as a function of τ ω t: q = q(ωt) = q (ωt) + ɛ q (ωt) + ɛ q (ωt) +..., ω = ω + ɛω + ɛ ω (.5) We may then insert both these expansions into our original equation, remembering that the original ω in it, multiplying q, has not changed and also that we use q to mean dq/dτ. We then collect it out in powers of ɛ, and then use the pair of terms at each level to eliminate terms caused by resonances, i.e., terms that are secular: ω ω ω [ q ] + q =, [ ] q + q = q 3 ω ω q, [ ] q + q = 3qq ω ω q (ω + ω ω )q,.... (.6a) The solution of the first equation is of course the same as before: q (τ) = a cos τ, which gives us the differential equation for q in the form ω [ q + q ] = a 3 cos 3 τ + ω ω a cos τ h(τ) = h[ω t + ɛω t + O (ɛ)] = h(ω t) + O(ɛ) h(τ ) = a 3 cos 3 τ + ω ω a cos τ = [ω ω 3a /4]a cos τ (a 3 /4) cos 3τ, 5 (.6b)
6 where the term in ɛ has been dropped in the second line, as required by the approximation. The term cos τ, which has the same frequency as the solution of the homogeneous equation, would be the cause of a secular term in the solution; therefore, we choose ω to cause its coefficient to vanish, and, as above, think of q as a function of τ : choose ω = 3 8 a [ ] = ω q + q = a3 ω 4 cos 3τ. (.6c) The solution of this that satisfies the boundary conditions is then straightforwardly q = q (τ) = q (τ ) = a3 3ω [cos 3ω t cos ω t], (.6d) where again we have dropped excess terms of higher order in ɛ. problem, at (only) first order in ɛ is given by Therefore, the solution to the q(t) = a cos[(ω + ɛω )t] + ɛ ω ω [cos 3ω t cos ω t]. (.6e) This is a result that is obviously better than the earlier, naïve one which had secular terms increasing off toward infinity for large times. As has already been stated, the problem is one of resonance between the natural frequency and a driving term. Therefore, let us briefly look additionally at a situation where there is a time-dependent driving term, so that we will essentially have two degrees of freedom, which will make it easier to see how the resonant terms cause serious grief. We consider a driven, linear oscillator q + ω q = g(t). (.7a) Obviously this problem and the one we have just recently considered are both special cases of a more general situation where we would have taken g(q, t) as a driving force. However, let us spend a little time now on this special case, where we assume the driving term is periodic with frequency Ω. The homogeneous solution is just the standard A cos ω t + B sin ω t. For comparison, then, we Fourier transform our periodic function g: g(t) = a + [a n cos(nωt) + b n sin(nωt)]. (.7b) n= If we now assume that the particular solution, q p (t), of our problem has the same sort of time dependence, i.e., that it is an infinite sum of terms q p = n= x pn, where q pn is periodic with frequency nω, then the insertion of this into the left-hand side of the differential equation would generate terms of the form [ (nω) + ω]q pn, so that we would obtain an equation insisting that q pn = a n cos(nωt) + b n sin(nωt) ω n Ω. (.7c) The resonance is very easy to see, which causes the solution to blow up. When the oscillator, unlike this one, also has a nonlinear dependence, or, if you prefer, g depend nonlinearly on q, then all multiples of ω will occur in the general solution, as we can see from our earlier excursion into such phenomena. Therefore, we will always have a difficulty whenever the ratio ω /Ω is a rational number. On the other hand, we also saw that, again unlike the SHO itself, for a nonlinear problem the actual frequency, ω, is a function of the initial conditions, and, in particular, the initial amplitude of the 6
7 motion. Therefore, by varying the initial conditions, we can vary this ratio. As the rationals are dense among the real numbers, the actual sets of initial conditions that causes resonances will also be dense among all possible initial conditions. This tells us immediately that we can expect serious changes in the topology of the solution manifold in phasespace when we allow sufficiently many degrees of freedom, i.e., as many as two! On the other hand, there is yet another problem with the secular perturbation theory as we have just described it. The derivation above showed that it is quite difficult to keep careful track of all the terms in ɛ, especially when they are contained within the argument in the form of τ = ωt. Therefore, a version which used the Hamiltonian formalism to keep better track of these dependencies was invented by Poincaré, in 89, and pushed further by von Zeipel, in 96, which is usually referred to as canonical perturbation theory. As did the secular perturbation above, it presupposes that the value of ɛ is given small, and fixed, while other things vary. We will describe the first-order version of it soon. On the other hand, more recently a somewhat distinct version of this theory has been created, based on the thoughts of Sophus Lie. It is based on the idea that the Hamiltonian always did depend on the parameter ɛ, and we should follow along a flow vector for which ɛ is the parameter. The canonical approach uses canonical generating functions, which depend on both old and new coordinates to determine the process; therefore the two become rather mixed up when a sequence of such transformations is performed, as is necessary in order to push the canonical process to higher orders in ɛ. On the other hand the Lie transformation process works directly with the ɛ- dependent canonical transformations themselves, rather than the generating functions, so that they give directly the new coordinates, in phase space, in terms of the old ones, and ɛ, thereby rendering multiple transformations much simpler in practice. It was introduced by Hori, in 966 and Garrido in 968, with important improvements toward the easier implementation of higher-order terms by Deprit, in 969, and has become rather popular since then. 3. Canonical Perturbation Theory We will first consider the canonical perturbation theory, with examples, and then go onward to the Lie transformation approach. We therefore presume that there is a given Hamiltonian which may be divided into a part which we can integrate to the end, and determine equations of motion, etc., and another part which we suppose is small, at least under appropriate circumstances; we use notation to keep track of this smallness by appending a parameter ɛ, which may, after the end of the calculations, be set to some value which may or may not really be small all the time, but one hopes so. We write this Hamiltonian in the following form, where at least most of the time, when one writes ξ α, it is referring to the system having having already been put into the action-angle variable formalism, i.e., ξ α = (φ a, J b ) H(ξ, t, ɛ) = H (ξ, t) + ɛ H (ξ, t) + ɛ H (ξ, t) (3.) The statement that we know how to resolve the unperturbed Hamiltonian, i.e., that which remains when we set ɛ =, is the statement that J b = J b, φ a = ω a t + φ a, ω a H / J a, (3.) where J b, ω a and φ a are all constants, independent of t. We are then going to seek a transformation to a new set of variables, η α = η α (ξ, t, ɛ) = (φ r, J s ) for which the problem is more tractable. More precisely, what we want is a transformation to a new set of variables so that the new Hamiltonian, K = K(η, t, ɛ) will in fact depend only upon the new action variables, J r, as would surely 7
8 be the case if we had actually been able to resolve the dynamics of this full Hamiltonian. This is an assumption, in advance, that this resolution is possible, which of course is certainly not always true! Therefore, one additional goal of this work is to determine methods to find out when it is not true, or perhaps when it is not even possible. (One should definitely point out that the resolution is indeed possible in the case of motions in only one degree of freedom, as all of those are indeed integrable.) At any event, we now agree that we consider doing it all up to some order in ɛ, for which we have agreed to solve the problem; we also agree not to ask whether or not the so-created series actually converges. (In many cases the series is actually an asymptotic one, so that it does not converge, but, at least, does give some reasonable approximation if one truncates it in the best way.) There is of course, in principle, a canonical transformation that generates this transformation. We will refer to it as S(φ, J, t, ɛ), and suppose that it is just the identity transform when ɛ =. In principle the transformation it generates, from ξ to the new canonical set, η = η(ξ, t, ɛ), should show its t- dependence, as just previously. However, we will drop the explicit presentation of the t-dependence as often as seems feasible, in the interests of keeping the notation less cumbersome: S = J r φ r + ɛs + ɛ S +..., (3.3a) We may then write out the usual CT rules for generation of the opposite pair of canonical subvariables: J a = S φ a = J a + ɛ S φ a (φ, J) + S ɛ (φ, J) +..., φa φ s = S = φ s + ɛ S J s J s (φ, J) + ɛ S (3.3b) J s (φ, J) +..., Having such a presentation we must also include the new Hamiltonian: K(J) = H[φ(φ, J), J(φ, J)] + S,t [φ(φ, J), J] K + ɛk + ɛ K (3.3c) In this last equation we say that K does not depend on φ; however, we will still have to insist that it does not, at each level of the perturbation, rather than simply ignore that dependence as we go along. It will in fact be this insistence that allows us to determine the form for the various levels of the generating function. To actually determine K(J), we will need to invert Eqs. (3.3b) to obtain the old variables, {φ, J}, in terms of the new. If we did in fact know S, S, etc., then this would be a coupled set of equations that are not completely trivial, because of the way they are mixed on both sides of the equation. Technically we would do it to some order in ɛ by iteration of the equation on itself. Showing how that iteration needs to work, we may write down the coupled equations in the following form: J a = J a + ɛ S φ a [φ(φ, J), J] + S ɛ [φ(φ, J), J] +..., φa (3.3d) φ a = φ a ɛ S [φ(φ, J), J] ɛ S [φ(φ, J), J] J a J a Clearly this equations lend themselves to an iterative approach to solution, which could begin to become tedious after, at least, the second order. We may now insert all this structure into the equation for K to determine equations for its parts to various levels of ɛ. The unperturbed Hamiltonian already does not depend on the original angle 8
9 variables, although the higher-order terms do; therefore, in principle we should insert into those terms the dependence of φ a on the set (φ r, J r ). However, since we will be wanting the higher-order terms of the generating function in terms of φ a, rather than φ r, we do not need to make that substitution; instead, it will come out naturally during the process of determining the higher-order terms in the generating function. As well, we insert J a = J a (φ, J) as it is naturally given by Eqs. (3.3c): K(J) = H [J(φ, J)] + ɛ H [φ, J(φ, J)] + ɛ H [φ, J(φ, J)] S,t [φ, J] = H (J + ɛ S,φ + ɛ S,φ +... ) + ɛ H (φ, J + ɛ S,φ +...) + ɛ H (φ, J + ɛ S,φ +...) + ɛs,t (φ, J) + ɛ S,t (φ, J) +... = H (J) + ɛ { H,Jr J=J S,φ r + H (φ, J) + S,t (φ, J) } + ɛ { H,Jr J=J S,φ r + H,Ja J=J S,φ a } + H,J a J b J=J S,φ as,φ b + H (φ, J) + S,t (φ, J) + O 3 (ɛ), (3.3e) = K (J) = ω(j) r S [φ, J] φ r + H [φ, J] + S,t (φ, J), K = ω(j) a S φ a + H S J J=J a φ a + ω a S S J J=J b φ a φ b + H (φ, J) + S,t (φ, J), where we have used Eqs. (3.) to introduce ω a = H / J a, which depends on J a, and then evaluated that in terms of the new variables, and kept terms only to lowest order. Now let s consider resolving these equations. We begin with the equations only first order in ɛ. Since we have more unknown functions than equations, we must now add some new requirements on the system. Since we are dealing with systems that are periodic, the notion of averaging appears in a very natural and intuitive way. Just what is to be averaged, over what curves, is not quite as obvious, however. The equation for K (J), at Eqs. (3.3e), is written as an equation in the variables (φ a, J r ), which constitute an acceptable set of variables over phasespace. Were the system unperturbed, then a set of values for the actions, {J a } n, pick out a set of independent tori, and the motion can be viewed as closed curves on each one, as each of the φ a varies from to π; therefore, we can average over them. Of course, the more accurate values for the action are, instead, the transformed actions, {J r } n. As we have assumed that the perturbed problem also is integrable, i.e., has its own tori, and closed curves on them, then we can imagine following around the closed curves viewed on the tori labelled by the J. Since this is an acceptable set of coordinates, the J must come back to their original values when we do this as well. Therefore, we may, instead project each of those curves onto the J tori, and follow along those curves instead. When that happens, the leading term in S will clearly change by πj. However, since neither J nor J changed, it must be that all of the derivatives of S with respect to φ must be periodic functions, so that they do not change. Therefore, each of the S j s, for j, must be periodic in φ, and perhaps also a term linear in φ, only. This last term may then be discarded. The result of all this is then that the average of S j / φ may be chosen to be zero! In addition, we will only consider time-dependencies that have some period, say Ω, in time. Similar arguments allow us to suppose that we average, also over the temporal period, so that the terms S j / t vanish. This allows us to resolve the equations by averaging both sides, which gives us the two useful equations: ( ) n+ π/ω n { π } K (J) = H (J) H (φ, J) Ω dt dφ a H (φ, J, t), π a= ω a S φ a + S (3.4) = H (φ, J) + H (J). t 9
10 The determination of the first perturbation to the Hamiltonian is then rather straightforward! On the other hand, the determination of this first-order correction to the generating function can be a little tricky! Although in principle one might determine the average value in various other ways, a very generic and algorithmic sort of way, for a function that is periodic in all these angles, is to expand it in a Fourier series. Since, by definition, the part desired for the calculation of S has zero average, it will have no constant term in the Fourier series: H H l,m H ml (J)e i(m aφ a +lωt) = ω a S φ a + S t. (3.5a) We then expand the generating function in the same way, and equate terms. The result is straightforward: S = i H ml (J) lω + m a ω a (J) ei(m aφ a +lωt). (3.5b) l,m This expression is a perfectly legitimate expression for S and does indeed provide a first-order approximation to the exact solution, provided that none of the denominators vanish! Obviously the method does not provide a solution to the problem if even one of the denominators is zero. This suggests that the entire perturbative procedure is invalid at this point, i.e., that the system s behavior cannot actually be described by action-angle variables that are related through a canonical transformation tot he original, unperturbed set. This will obviously happen for frequencies that are commensurate, unless, for some strange reason, the particular coefficient, H ml, were to also vanish. [Do, of course, remember that this problem does not occur in systems with only one degree of freedom, since the denominator will not vanish.] Even a more serious problem is the fact that one may find a set of conditions such that one or more denominator becomes very small; i.e., we are considering initial conditions near to those which caused the denominator to vanish. In that case, the particular term will be very large. This means that multiplying by even your standard small value of ɛ does not make that term small; therefore, the correction becomes so very large that the perturbative scheme fails, working only for a much more restricted set of values for ɛ. Eventually the KAM theorem will give us some specific methods of finding ways to avoid these degeneracies, and also to know how much accuracy we can expect at nearby points, in phasespace. For the moment, let us forget these problems temporarily and return to those terms which determine the approximation to second order in ɛ. Proceeding with an averaging method as before, we obtain the following, where I consider only the time-independent case: H ω a (J) S φ a K = = K S J J=J a φ a + ω a J b J=J S φ a S φ b + H ( H S J J=J a φ a + ω a S S J J=J b φ a φ b + H, ) (3.6) Now, we recall that we already have an explicit expression for S, as an infinite Fourier series, and therefore also one for S / φ a, from our first-order calculation. We may certainly insert that into this, and thereby obtain an expression for the second-order correction to the energy, having only found the generating function to the first order. This is, of course, a very general property of perturbative series: to calculate the correction to the energy to n-th order, one needs to determine the complete correction, including in this case the generating function, only to one earlier order, i.e., to the n -st order.
11 On the other hand, at least at the moment, I will only consider the rather simpler question of further calculation of this form for the second-order correction to the energy in the situation when there is only one degree of freedom. There is no real impediment to doing the more general case, without time dependence, but it is simply messier. In the case of one degree of freedom, Eq. (3.4) takes the form ω ( S / φ) = H H, so that our equation for K becomes K (J) = ( H H ω J H H J ) + ω / J ω ( (H ) H ) + H. (3.7) 4. Some Examples: As a reasonable, rather simple, first example, let us look at our approximation to the pendulum, and treat it correct to second order in a small parameter; therefore, I retreat to Eq. (.) and take the Hamiltonian there, expand the cosine function in an infinite series, throw away the constant term, keep the quadratic term in H, and replace q by ɛq, to keep track of the extra terms and the perturbation parameter. At first order, this will be the same as our example of the quartic oscillator, discussed in. The small parameter there, which was also called ɛ, we will now call η = ɛω = ɛf G. Then we have the following: H = Gp + F q 4 ɛf q4 + ɛ F q H + ɛ H + ɛ H +... ; ( ) φ = tan mω q/p, J = ( p /mω + mω q ), for ɛ = q = J sin φ, p = mω J cos φ, mω H = ω J, = H(φ, J) = H + ɛ H + ɛ H = ω J ɛ G J sin 4 φ + G 5ɛ ω J 3 sin 6 φ = ω J ɛ 8 G J [3 4 cos φ + cos 4φ] + ɛ G 8 ω J 3 [ 5 cos φ + 6 cos 4φ cos 6φ]. (4.) The average of H is immediate, when it has been expressed as a Fourier series, which in this case is finite. Therefore we immediately conclude that K (J) = 3 8 G J, = ω = 3 4 G J, so that ω = ω ɛ 3 4 G J + O (ɛ). We may now determine S, and the associated canonical transformation: S = ( ) dφ H H = G J ( ) dφ 4 cos φ cos 4φ = GJ [8 sin φ sin 4φ]. ω 8ω 3ω J=J From this we may write out the CT-generated equations for J and φ, from Eqs. (3.3c): J = J + S,φ = J ɛ GJ 8ω [4 cos φ cos 4φ] + O (ɛ), φ = φ + S,J = φ ɛ GJ 6ω [8 sin φ sin 4φ] + O (ɛ). (4.) (4.3) (4.4a)
12 In principle, the solution of this for either set (φ, J) or (φ, J) in terms of the other is rather difficult. On the other hand, it is not fair to solve it more accurately than first-order. In that case we have, at least in one direction, just the following: J = J ɛ GJ [4 cos φ cos 4φ] + O (ɛ), 8ω φ = φ ɛ GJ [8 sin φ sin 4φ] + O (ɛ). 6ω (4.4b) Proceeding onward, one should then be able to determine K and the next-order correction to the energy and the frequency. However, let us stop for now at consider another example. We extend the quartic oscillator to two degrees of freedom, and then couple them, with the the Hamiltonian for this system as H = (p + p ) + (ω q + ω q ) + ɛ(ω ω q q ) H + ɛ H. (4.5) The unperturbed action-angle variables are simply two copies of the ones we had above, for one linear oscillator; therefore, we may immediately write H(φ, J) = ω J + ω J + 4ɛω J ω J sin φ sin φ. (4.6) Therefore the first-order correction to the energy is simply K (J) = H = 4ω J ω J 4π dφ sin φ which tells us that the first-order frequencies are dφ sin φ = ω J ω J, (4.7) ω = ω + ɛω ω J, ω = ω + ɛω ω J. (4.8) Now, however, obtaining the (first-order corrections to the) generating function are where we anticipate trouble from resonances. Therefore, we proceed toward that calculation, i.e., a solution of the equation given in Eq. (3.5a), for our case: ( ) ) ω φ + ω φ S = ( H H. (4.9) As discussed at that point, a simple way to resolve this differential equation is by use of the Fourier transform of the non-constant part of the first-order Hamiltonian, to insert into Eq. (3.5b). Therefore, we first calculate those Fourier coefficients: H,(m,m ) = ω ω J J (π) π π dφ dφ (4 sin φ sin φ ) e i(m φ +m φ ). (4.a) It is reasonably clear that only a few of these are non-zero; more precisely, those that involve ± and. We have the following, where we use Eq. (3.5b), with Ω =, to determine the coefficients S,(m,m ): S,(±,±) = H,(±,±) = + 4 K, H,(±,) = K = H,(,±), H,(,) =, ±i/4 (ω + ω ), S,(±, ) = ±i/4 (ω ω ), S,(±,) = ± i/, S,(,±) = ± i/. ω ω (4.b)
13 Therefore, we may insert these back into Eq. (3.5b), with Ω =, to obtain the generating function: ) S = J J (ω sin φ + ω sin φ ω ω (ω + ω ) sin (φ + φ ω ω ) (ω ω ) sin (φ φ ) J J F (φ). (4.c) We may then write down the first-order CT, where, to that order, we have replaced the J s on the right-hand side of the equation with the older J s: φ = φ + ɛj F (φ), φ = φ + ɛj F (φ), F J = J + ɛj J φ, J F = J + ɛj J φ. (4.d) This example gives us our first view of the resonance problem with perturbation theory. Here there is an obvious problem with resonance; as well, it even shows us just why it s called that, since it occurs when the two natural frequencies are equal, i.e., in resonance. On the other hand, this one has a problem only when those two frequencies are equal, when one might well expect it. Its Fourier spectrum is sufficiently simple that there are no problems even when the ratio of the two frequencies is any rational number except for the very special case when that number is. On the other hand, it is true that more complicated resonances would probably occur at higher orders in perturbation theory. In fact, as we will see later, there are situations where this dynamical system becomes chaotic. As an example of an alternative sort of a problem, we consider the following problem where we have a time-dependent (external) driving agent, and an explicit magnetic field. The dynamical system is a charged particle moving in a uniform magnetic field and interacting with an additional, oscillating electric field. Let the charged particle have mass m and charge q cβ, where c is the speed of light. Describe the magnetic field by a vector potential A = B yˆx, which causes the magnetic field to be in the ẑ-direction, B = A = +B ẑ, orthogonal to the x, y-plane. We choose the axes in that plane so that the oscillating electric field has wave vector only in the y- and z-directions, i.e., k = kz ẑ + k y ŷ, and let it be generated by a (time-dependent) potential such that the Hamiltonian may be put in the following form, where ɛ is simply, as usual, an ordering parameter: H = m p βa + ɛβv sin(k z z + k y y ωt). (4.) We first study the unperturbed system, i.e., the one with ɛ =, to put it in the proper form. We recall that the motion in such a situation the particle spirals along the direction of the field and circles in the plane perpendicular to the field, with a constant angular velocity, Ω = βb /m, around such origin. We take that origin for the circles to have (constant) coordinates in the plane given by (X, Y ), and radius ρ, so that we may write the motion in that plane in terms of a single (angular) coordinate φ, related in the following way: The unperturbed Hamiltonian has the explicit form x = X + ρ cos φ, y = Y ρ sin φ. (4.a) m H = (p x + βb y) + (p y ) + (p z ), where we may also write out the Hamiltonian in terms of the velocities as usual: p x =m d dt x βb y = mρ sin φ φ βb y, p = m v + βa = p y =m d dt y = mρ cos φ φ, p z =m d dt z. 3 (4.b) (4.c)
14 Remembering that φ = Ω = βb /m, we may rewrite the equation for p x as p x = ρ sin φ φ βb y = mω(ρ sin φ + y) = mωy, (4.d) explaining the perhaps-unexpected choice of sign in the relation chosen between x, y, and φ. Of course a different way of looking at that sign is that it causes the variable φ to progress clockwise around the orbit. As the circle progresses counterclockwise, this is basically a co-rotating set of coordinates for the system. We now want to move from the current coordinates, {x, y, z; p x, p y, p z }, to a new set {Y, φ, z; P Y, P φ, p z }, which we will do via a Type generating function, which maintains the z-oriented variables the same as before: = F (x, y, z; Y, φ, z) mω[ (y Y ) cot φ xy ], p x = F,x = mωy, p y = F,y = mω(y Y ) cot φ = mωρ cos φ, P Y = F,Y = mω[(y Y ) cot φ + x] = mω[ ρ cos φ + X + ρ cos φ] = mωx, P φ = F,π = mω(y Y ) csc φ = mωρ, (4.e) where we notice that the equations for the original momenta are consistent with their earlier versions shown above. Inserting all this into the equation for H, we find, in these variables, that [ K = H [old(new)] = [mω(y y)] + (mωρ cos φ) + (p z ) ] /(m) = mω ρ + p z/m = Ω P φ + p z/m, (4.f) where the form is of course what is expected for a transition to action-angle variables, i.e., terms of the form J φ, except for the fact that there are no terms at all that involve Y and its momentum, as both of them are constant! (Of course nothing happened to the z-direction quantities, since we did not involve them in the transformation.) We may now express H, the coefficient of ɛ, in terms of these new coordinates, which we do while using ρ as simply an abbreviation for its value, resolved from the equation for P φ above, namely ρ (P φ /MΩ): K = H [old(new)] = βv sin(k z z ωt + k p ρ sin φ). (4.3a) The non-linearity comes from the dependence of ρ on P φ and from the sin φ. On the other hand, the other two remaining variables, z and t, occur only in the combination k z z ωt; therefore, we eliminate the time, explicitly, by another transformation to that one variable, which we refer to as ψ, thinking of it surely as an angular type variable since it only appears as an argument of the sine function. A Type generating function which does that is the following, where the last term, φp φ, is simply the identity transformation which leaves those two variables alone: = F (k z z ωt)p ψ + φp φ p z = F z = k z P ψ, P φ = F φ = P φ, ψ = F P ψ = k z z ωt, φ = F 4 P φ = φ, (4.3b)
15 This generating function depends explicitly on the time, so that the Hamiltonian now has the form K = K ( + ɛ + F,t + K ) = (k z P ψ ) /m + Ω P φ ωp ψ + ɛβv sin(ψ k p ρ sin φ), (4.3c) where we have not put the multiplier ɛ on the term from the time derivative, ωp ψ, since that term is easily included with the rest of the terms and is still integrable. Since we now have a nonlinearity that is of the form of a sine function of a sine function of a variable, we need to expand this in a Fourier series to see explicitly what sort of resonances may have the possibility of being present. We use a standard identity for such things, that is related to the generating function for (integer-index) Bessel functions, or, if you prefer, we simply expand it all in a Fourier series: sin(ψ z sin φ) = + l= J l (z) sin(ψ lφ) (4.4) This gives us the equivalent form for the Hamiltonian with the perturbation in the form of a Fourier series, as we know we need in order to determine the generating function for this perturbed system: K K + ɛ K = (k z P ψ ) /m + ΩP φ ωp ψ + ɛβv + l= J l [k p ρ(p φ )] sin(ψ lφ). (4.5a) The unperturbed frequencies are then ω φ = K P φ = Ω, ω ψ = K P ψ = ω + (k z/m)p ψ. (4.5b) From our general discussion of resonances, in particular at Eqs. (3.5b), we know that there will be a problem with the perturbation expansion when the initial conditions put us too close to the places where the ratio of the original frequencies is rational. Since we only have two frequencies here, this is the same as suggesting that = ω ψ + mω = (k z/m)p ψ + ω + mω. (4.5c) It turns out that the situation when k z is rather more complicated than the other case; physically, of course this turns out to be the case when the incoming electric field is not perpendicular to the existing magnetic field. (It might be expected that these are intrinsically different.) Therefore, we now consider in more detail (at least) the case when k z =. To perform this trick, we first return to our equations, from 3, for the first-order quantities, remembering that what we want is to determine a CT from the old variables to the yet-newer ones, which, as before, we describe with overbars: K = K (new) = (k z P ψ ) /m + ΩP φ ωp ψ K = K (new) = ( π ) βv π dψ π dφ (Ω φ ω ψ )S = K = βv + l= + l= J l (k p ρ) sin(ψ lφ) =, J l (k p ρ) sin(ψ lφ). (4.6) 5
16 As this is a situation where the first-order correction to the energy vanishes, one might expect that it is more likely to have problems, sooner. We proceed onward, then, to determine the generating function. The desired solution to the differential equation, i.e., the particular solution, is not hard to determine: + J l [k p ρ(p φ )] S = βv cos(ψ lφ). (4.7) ω + mω l= We can then determine this CT, to first-order, by using (the inverse of) Eqs. (3.3d): P φ = P φ + ɛβv + l= P ψ = P ψ ɛβv + l= l J l[k p ρ(p φ )] ω + lω J l [k p ρ(p φ )] ω + lω sin(ψ lφ), sin(ψ lφ). (4.8) Since the Hamiltonian didn t change under this transformation, we can immediately see that it does not depend on any of the transformed angular coordinates; therefore, the transformed momenta, P φ and P ψ, are constant. Therefore these two equations may be conceived as equations giving us the functional form of each of the un-transformed momenta as functions of the two angles. We may therefore easily describe Poincaré sections, in, for instance, the φ, P ψ -plane by plotting the second of Eqs. (4.8) for some fixed value of the other angle, φ. For each fixed value of the constant, P φ dependent on initial conditions this will give us a curve of P φ versus ψ. We could do this for various values of the ratio of ω versus Ω. It might have seemed more straightforward to have plotted P φ versus φ, but surely similar information is provided in this way as well, and Lichtenberg has provided plots for this case, for the ratio ω = 3.Ω, where various different values of P φ are introduced, to produce a variety of distinct motions on the same plot. The plots are given to us for two different values of the driving amplitude, V, one referred to as low wave amplitude, while the second one has higher amplitude. Actually the curves were obtained by C.F.F. Karney, in his (MIT, 977) thesis studying stochastic heating of ions in a tokamak, induced by RF power. The left-hand curves are produced from the perturbation analysis equations, above, and they are then compared, on the right-hand side, with analogous curves for the same initial conditions, but calculated numerically, and, therefore, presumably rather more accurate. The first set, below, at low amplitude for the electric wave, has very good agreement between the perturbative curves and the exact numerical curves: This second pair, below, with higher amplitude, do mirror some of the gross features, but the perturbative curves clearly do not reproduce the chaotic structure, nor the chains of islands seen in the associated numerical calculations: 6
17 I note that there is a related html-file, at a nearby link to this one on the class homepage, that shows these curves and also the results of some additional calculations of my own, for the first-order case. 7
Chapter 11 Equilibrium
11.1 The First Condition of Equilibrium The first condition of equilibrium deals with the forces that cause possible translations of a body. The simplest way to define the translational equilibrium ofMore information
Figure 2.1: Center of mass of four points.
Chapter 2 Bézier curves are named after their inventor, Dr. Pierre Bézier. Bézier was an engineer with the Renault car company and set out in the early 196 s to develop a curve formulation which wouldMore information
Copyrighted Material. Chapter 1 DEGREE OF A CURVE
Chapter 1 DEGREE OF A CURVE Road Map The idea of degree is a fundamental concept, which will take us several chapters to explore in depth. We begin by explaining what an algebraic curve is, and offer twoMore information
7. DYNAMIC LIGHT SCATTERING 7.1 First order temporal autocorrelation function.
7. DYNAMIC LIGHT SCATTERING 7. First order temporal autocorrelation function. Dynamic light scattering (DLS) studies the properties of inhomogeneous and dynamic media. A generic situation is illustratedMore information
Second Order Linear Nonhomogeneous Differential Equations; Method of Undetermined Coefficients. y + p(t) y + q(t) y = g(t), g(t) 0.
Second Order Linear Nonhomogeneous Differential Equations; Method of Undetermined Coefficients We will now turn our attention to nonhomogeneous second order linear equations, equations with the standardMore information
710 Chapter Thirteen A FUNDAMENTAL TOOL: VECTORS 62. Use the following steps and the results of Problems 59 60 to show (without trigonometry) that the geometric and algebraic definitions of the dot productMore information
Chapter 24 Physical Pendulum
Chapter 4 Physical Pendulum 4.1 Introduction... 1 4.1.1 Simple Pendulum: Torque Approach... 1 4. Physical Pendulum... 4.3 Worked Examples... 4 Example 4.1 Oscillating Rod... 4 Example 4.3 Torsional Oscillator...More information
Exponentially small splitting of separatrices for1 1 2 degrees of freedom Hamiltonian Systems close to a resonance. Marcel Guardia
Exponentially small splitting of separatrices for1 1 2 degrees of freedom Hamiltonian Systems close to a resonance Marcel Guardia 1 1 1 2 degrees of freedom Hamiltonian systems We consider close to completelyMore information
WHAT ARE MATHEMATICAL PROOFS AND WHY THEY ARE IMPORTANT? introduction Many students seem to have trouble with the notion of a mathematical proof. People that come to a course like Math 216, who certainlyMore information
Section 1.1. Introduction to R n
The Calculus of Functions of Several Variables Section. Introduction to R n Calculus is the study of functional relationships and how related quantities change with each other. In your first exposure toMore information
Lecture L19 - Vibration, Normal Modes, Natural Frequencies, Instability
S. Widnall 16.07 Dynamics Fall 2009 Version 1.0 Lecture L19 - Vibration, Normal Modes, Natural Frequencies, Instability Vibration, Instability An important class of problems in dynamics concerns the freeMore information
2.2. Instantaneous Velocity
2.2. Instantaneous Velocity toc Assuming that your are not familiar with the technical aspects of this section, when you think about it, your knowledge of velocity is limited. In terms of your own mathematicalMore information
Eðlisfræði 2, vor 2007
[ Assignment View ] [ Pri Eðlisfræði 2, vor 2007 28. Sources of Magnetic Field Assignment is due at 2:00am on Wednesday, March 7, 2007 Credit for problems submitted late will decrease to 0% after the deadlineMore information
Some Notes on Taylor Polynomials and Taylor Series
Some Notes on Taylor Polynomials and Taylor Series Mark MacLean October 3, 27 UBC s courses MATH /8 and MATH introduce students to the ideas of Taylor polynomials and Taylor series in a fairly limitedMore information
discuss how to describe points, lines and planes in 3 space.
Chapter 2 3 Space: lines and planes In this chapter we discuss how to describe points, lines and planes in 3 space. introduce the language of vectors. discuss various matters concerning the relative positionMore information
2 Integrating Both Sides
2 Integrating Both Sides So far, the only general method we have for solving differential equations involves equations of the form y = f(x), where f(x) is any function of x. The solution to such an equationMore information
South Carolina College- and Career-Ready (SCCCR) Pre-Calculus
South Carolina College- and Career-Ready (SCCCR) Pre-Calculus Key Concepts Arithmetic with Polynomials and Rational Expressions PC.AAPR.2 PC.AAPR.3 PC.AAPR.4 PC.AAPR.5 PC.AAPR.6 PC.AAPR.7 Standards KnowMore information
Harmonic Oscillator Physics
Physics 34 Lecture 9 Harmonic Oscillator Physics Lecture 9 Physics 34 Quantum Mechanics I Friday, February th, 00 For the harmonic oscillator potential in the time-independent Schrödinger equation: d ψxMore information
11.1. Objectives. Component Form of a Vector. Component Form of a Vector. Component Form of a Vector. Vectors and the Geometry of Space
11 Vectors and the Geometry of Space 11.1 Vectors in the Plane Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. 2 Objectives! Write the component form ofMore information
Oscillations. Vern Lindberg. June 10, 2010
Oscillations Vern Lindberg June 10, 2010 You have discussed oscillations in Vibs and Waves: we will therefore touch lightly on Chapter 3, mainly trying to refresh your memory and extend the concepts. 1More information
PHY411. PROBLEM SET 3 1. Conserved Quantities; the Runge-Lenz Vector The Hamiltonian for the Kepler system is H(r, p) = p2 2 GM r where p is momentum, L is angular momentum per unit mass, and r is theMore information
CHAPTER 3 Numbers and Numeral Systems
CHAPTER 3 Numbers and Numeral Systems Numbers play an important role in almost all areas of mathematics, not least in calculus. Virtually all calculus books contain a thorough description of the natural,More information
Solving Simultaneous Equations and Matrices
Solving Simultaneous Equations and Matrices The following represents a systematic investigation for the steps used to solve two simultaneous linear equations in two unknowns. The motivation for consideringMore information
UNIT 1: Unit code: QCF Level: 4 Credit value: 15 ANALYTICAL METHODS FOR ENGINEERS A/601/1401 OUTCOME - TRIGONOMETRIC METHODS TUTORIAL 1 SINUSOIDAL FUNCTION Be able to analyse and model engineering situationsMore information
Mathematical Procedures
CHAPTER 6 Mathematical Procedures 168 CHAPTER 6 Mathematical Procedures The multidisciplinary approach to medicine has incorporated a wide variety of mathematical procedures from the fields of physics,More information
Taylor Polynomials and Taylor Series Math 126
Taylor Polynomials and Taylor Series Math 26 In many problems in science and engineering we have a function f(x) which is too complicated to answer the questions we d like to ask. In this chapter, we willMore information
3.2 Sources, Sinks, Saddles, and Spirals
3.2. Sources, Sinks, Saddles, and Spirals 6 3.2 Sources, Sinks, Saddles, and Spirals The pictures in this section show solutions to Ay 00 C By 0 C Cy D 0. These are linear equations with constant coefficientsMore information
Frequency response of a general purpose single-sided OpAmp amplifier
Frequency response of a general purpose single-sided OpAmp amplifier One configuration for a general purpose amplifier using an operational amplifier is the following. The circuit is characterized by:More information
N 1. (q k+1 q k ) 2 + α 3. k=0
Teoretisk Fysik Hand-in problem B, SI1142, Spring 2010 In 1955 Fermi, Pasta and Ulam 1 numerically studied a simple model for a one dimensional chain of non-linear oscillators to see how the energy distributionMore information
Two vectors are equal if they have the same length and direction. They do not
Vectors define vectors Some physical quantities, such as temperature, length, and mass, can be specified by a single number called a scalar. Other physical quantities, such as force and velocity, mustMore information
Let s first see how precession works in quantitative detail. The system is illustrated below: ...
lecture 20 Topics: Precession of tops Nutation Vectors in the body frame The free symmetric top in the body frame Euler s equations The free symmetric top ala Euler s The tennis racket theorem As you know,More information
The Quantum Harmonic Oscillator Stephen Webb
The Quantum Harmonic Oscillator Stephen Webb The Importance of the Harmonic Oscillator The quantum harmonic oscillator holds a unique importance in quantum mechanics, as it is both one of the few problemsMore information
ELLIPTIC CURVES AND LENSTRA S FACTORIZATION ALGORITHM DANIEL PARKER Abstract. This paper provides a foundation for understanding Lenstra s Elliptic Curve Algorithm for factoring large numbers. We giveMore information
Lectures 5-6: Taylor Series
Math 1d Instructor: Padraic Bartlett Lectures 5-: Taylor Series Weeks 5- Caltech 213 1 Taylor Polynomials and Series As we saw in week 4, power series are remarkably nice objects to work with. In particular,More information
PURSUITS IN MATHEMATICS often produce elementary functions as solutions that need to be
Fast Approximation of the Tangent, Hyperbolic Tangent, Exponential and Logarithmic Functions 2007 Ron Doerfler http://www.myreckonings.com June 27, 2007 Abstract There are some of us who enjoy using ourMore information
Lecture L5 - Other Coordinate Systems
S. Widnall, J. Peraire 16.07 Dynamics Fall 008 Version.0 Lecture L5 - Other Coordinate Systems In this lecture, we will look at some other common systems of coordinates. We will present polar coordinatesMore information
Lecture 6. Weight. Tension. Normal Force. Static Friction. Cutnell+Johnson: 4.8-4.12, second half of section 4.7
Lecture 6 Weight Tension Normal Force Static Friction Cutnell+Johnson: 4.8-4.12, second half of section 4.7 In this lecture, I m going to discuss four different kinds of forces: weight, tension, the normalMore information
G U I D E T O A P P L I E D O R B I T A L M E C H A N I C S F O R K E R B A L S P A C E P R O G R A M CONTENTS Foreword... 2 Forces... 3 Circular Orbits... 8 Energy... 10 Angular Momentum... 13 FOREWORDMore information
Notes on Factoring. MA 206 Kurt Bryan
The General Approach Notes on Factoring MA 26 Kurt Bryan Suppose I hand you n, a 2 digit integer and tell you that n is composite, with smallest prime factor around 5 digits. Finding a nontrivial factorMore information
REVIEW EXERCISES DAVID J LOWRY Contents 1. Introduction 1 2. Elementary Functions 1 2.1. Factoring and Solving Quadratics 1 2.2. Polynomial Inequalities 3 2.3. Rational Functions 4 2.4. Exponentials andMore information
AP1 Oscillations. 1. Which of the following statements about a spring-block oscillator in simple harmonic motion about its equilibrium point is false?
1. Which of the following statements about a spring-block oscillator in simple harmonic motion about its equilibrium point is false? (A) The displacement is directly related to the acceleration. (B) TheMore information
Section V.4: Cross Product
Section V.4: Cross Product Definition The cross product of vectors A and B is written as A B. The result of the cross product A B is a third vector which is perpendicular to both A and B. (Because theMore information
. INNER PRODUCT SPACES.. Definition So far we have studied abstract vector spaces. These are a generalisation of the geometric spaces R and R. But these have more structure than just that of a vector space.More information
Chapter 4, Arithmetic in F [x] Polynomial arithmetic and the division algorithm.
Chapter 4, Arithmetic in F [x] Polynomial arithmetic and the division algorithm. We begin by defining the ring of polynomials with coefficients in a ring R. After some preliminary results, we specializeMore information
Applications of Second-Order Differential Equations
Applications of Second-Order Differential Equations Second-order linear differential equations have a variety of applications in science and engineering. In this section we explore two of them: the vibrationMore information
Copyright 2011 Casa Software Ltd. www.casaxps.com
Table of Contents Variable Forces and Differential Equations... 2 Differential Equations... 3 Second Order Linear Differential Equations with Constant Coefficients... 6 Reduction of Differential EquationsMore information
Solutions to Exercises, Section 5.1
Instructor s Solutions Manual, Section 5.1 Exercise 1 Solutions to Exercises, Section 5.1 1. Find all numbers t such that ( 1 3,t) is a point on the unit circle. For ( 1 3,t)to be a point on the unit circleMore information
LINEAR ALGEBRA W W L CHEN c W W L Chen, 1997, 2008 This chapter is available free to all individuals, on understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,More information
Rotation Rate of a Trajectory of an Algebraic Vector Field Around an Algebraic Curve
QUALITATIVE THEORY OF DYAMICAL SYSTEMS 2, 61 66 (2001) ARTICLE O. 11 Rotation Rate of a Trajectory of an Algebraic Vector Field Around an Algebraic Curve Alexei Grigoriev Department of Mathematics, TheMore information
Lecture 14: Section 3.3
Lecture 14: Section 3.3 Shuanglin Shao October 23, 2013 Definition. Two nonzero vectors u and v in R n are said to be orthogonal (or perpendicular) if u v = 0. We will also agree that the zero vector inMore information
SOLID MECHANICS TUTORIAL MECHANISMS KINEMATICS - VELOCITY AND ACCELERATION DIAGRAMS This work covers elements of the syllabus for the Engineering Council exams C105 Mechanical and Structural EngineeringMore information
2.2 Magic with complex exponentials
2.2. MAGIC WITH COMPLEX EXPONENTIALS 97 2.2 Magic with complex exponentials We don t really know what aspects of complex variables you learned about in high school, so the goal here is to start more orMore information
Lab M1: The Simple Pendulum
Lab M1: The Simple Pendulum Introduction. The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are usually regarded as the beginning ofMore information
Structure of the Root Spaces for Simple Lie Algebras
Structure of the Root Spaces for Simple Lie Algebras I. Introduction A Cartan subalgebra, H, of a Lie algebra, G, is a subalgebra, H G, such that a. H is nilpotent, i.e., there is some n such that (H)More information
The Central Equation
The Central Equation Mervyn Roy May 6, 015 1 Derivation of the central equation The single particle Schrödinger equation is, ( H E nk ψ nk = 0, ) ( + v s(r) E nk ψ nk = 0. (1) We can solve Eq. (1) at givenMore information
Physics 9e/Cutnell. correlated to the. College Board AP Physics 1 Course Objectives
Physics 9e/Cutnell correlated to the College Board AP Physics 1 Course Objectives Big Idea 1: Objects and systems have properties such as mass and charge. Systems may have internal structure. EnduringMore information
Chapter 2. Derivation of the Equations of Open Channel Flow. 2.1 General Considerations
Chapter 2. Derivation of the Equations of Open Channel Flow 2.1 General Considerations Of interest is water flowing in a channel with a free surface, which is usually referred to as open channel flow.More information
Equations Involving Lines and Planes Standard equations for lines in space
Equations Involving Lines and Planes In this section we will collect various important formulas regarding equations of lines and planes in three dimensional space Reminder regarding notation: any quantityMore information
Eðlisfræði 2, vor 2007
[ Assignment View ] [ Print ] Eðlisfræði 2, vor 2007 30. Inductance Assignment is due at 2:00am on Wednesday, March 14, 2007 Credit for problems submitted late will decrease to 0% after the deadline hasMore information
MATH10212 Linear Algebra. Systems of Linear Equations. Definition. An n-dimensional vector is a row or a column of n numbers (or letters): a 1.
MATH10212 Linear Algebra Textbook: D. Poole, Linear Algebra: A Modern Introduction. Thompson, 2006. ISBN 0-534-40596-7. Systems of Linear Equations Definition. An n-dimensional vector is a row or a columnMore information
2 Complex Functions and the Cauchy-Riemann Equations
2 Complex Functions and the Cauchy-Riemann Equations 2.1 Complex functions In one-variable calculus, we study functions f(x) of a real variable x. Likewise, in complex analysis, we study functions f(z)More information
This is the 23 rd lecture on DSP and our topic for today and a few more lectures to come will be analog filter design. (Refer Slide Time: 01:13-01:14)
Digital Signal Processing Prof. S.C. Dutta Roy Department of Electrical Electronics Indian Institute of Technology, Delhi Lecture - 23 Analog Filter Design This is the 23 rd lecture on DSP and our topicMore information
2.3. Finding polynomial functions. An Introduction:
2.3. Finding polynomial functions. An Introduction: As is usually the case when learning a new concept in mathematics, the new concept is the reverse of the previous one. Remember how you first learnedMore information
VELOCITY, ACCELERATION, FORCE velocity Velocity v is a vector, with units of meters per second ( m s ). Velocity indicates the rate of change of the object s position ( r ); i.e., velocity tells you howMore information
ELEMENTS OF VECTOR ALGEBRA A.1. VECTORS AND SCALAR QUANTITIES We have now proposed sets of basic dimensions and secondary dimensions to describe certain aspects of nature, but more than just dimensionsMore information
Finding Equations of Sinusoidal Functions From Real-World Data
Finding Equations of Sinusoidal Functions From Real-World Data **Note: Throughout this handout you will be asked to use your graphing calculator to verify certain results, but be aware that you will NOTMore information
Section 4: The Basics of Satellite Orbits
Section 4: The Basics of Satellite Orbits MOTION IN SPACE VS. MOTION IN THE ATMOSPHERE The motion of objects in the atmosphere differs in three important ways from the motion of objects in space. First,More information
Continuous Groups, Lie Groups, and Lie Algebras
Chapter 7 Continuous Groups, Lie Groups, and Lie Algebras Zeno was concerned with three problems... These are the problem of the infinitesimal, the infinite, and continuity... Bertrand Russell The groupsMore information
COLLEGE ALGEBRA. Paul Dawkins
COLLEGE ALGEBRA Paul Dawkins Table of Contents Preface... iii Outline... iv Preliminaries... Introduction... Integer Exponents... Rational Exponents... 9 Real Exponents...5 Radicals...6 Polynomials...5More information
Updated 2013 (Mathematica Version) M1.1. Lab M1: The Simple Pendulum
Updated 2013 (Mathematica Version) M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s areMore information
CSCI567 Machine Learning (Fall 2014)
CSCI567 Machine Learning (Fall 2014) Drs. Sha & Liu {feisha,email@example.com September 22, 2014 Drs. Sha & Liu ({feisha,email@example.com) CSCI567 Machine Learning (Fall 2014) September 22, 2014 1 /More information
Mathematics Course 111: Algebra I Part IV: Vector Spaces
Mathematics Course 111: Algebra I Part IV: Vector Spaces D. R. Wilkins Academic Year 1996-7 9 Vector Spaces A vector space over some field K is an algebraic structure consisting of a set V on which areMore information
MATH 4330/5330, Fourier Analysis Section 11, The Discrete Fourier Transform
MATH 433/533, Fourier Analysis Section 11, The Discrete Fourier Transform Now, instead of considering functions defined on a continuous domain, like the interval [, 1) or the whole real line R, we wishMore information
Midterm Solutions. mvr = ω f (I wheel + I bullet ) = ω f 2 MR2 + mr 2 ) ω f = v R. 1 + M 2m
Midterm Solutions I) A bullet of mass m moving at horizontal velocity v strikes and sticks to the rim of a wheel a solid disc) of mass M, radius R, anchored at its center but free to rotate i) Which ofMore information
Thnkwell s Homeschool Precalculus Course Lesson Plan: 36 weeks
Thnkwell s Homeschool Precalculus Course Lesson Plan: 36 weeks Welcome to Thinkwell s Homeschool Precalculus! We re thrilled that you ve decided to make us part of your homeschool curriculum. This lessonMore information
1 Lecture 3: Operators in Quantum Mechanics
1 Lecture 3: Operators in Quantum Mechanics 1.1 Basic notions of operator algebra. In the previous lectures we have met operators: ˆx and ˆp = i h they are called fundamental operators. Many operatorsMore information
1 if 1 x 0 1 if 0 x 1
Chapter 3 Continuity In this chapter we begin by defining the fundamental notion of continuity for real valued functions of a single real variable. When trying to decide whether a given function is orMore information
1 of 10 1/25/2008 3:38 AM SIGNAL PROCESSING & SIMULATION NEWSLETTER Note: This is not a particularly interesting topic for anyone other than those who ar e involved in simulation. So if you have difficultyMore information
Electromagnetism Laws and Equations
Electromagnetism Laws and Equations Andrew McHutchon Michaelmas 203 Contents Electrostatics. Electric E- and D-fields............................................. Electrostatic Force............................................2More information
1 TRIGONOMETRY. 1.0 Introduction. 1.1 Sum and product formulae. Objectives
TRIGONOMETRY Chapter Trigonometry Objectives After studying this chapter you should be able to handle with confidence a wide range of trigonometric identities; be able to express linear combinations ofMore information
2x + y = 3. Since the second equation is precisely the same as the first equation, it is enough to find x and y satisfying the system
1. Systems of linear equations We are interested in the solutions to systems of linear equations. A linear equation is of the form 3x 5y + 2z + w = 3. The key thing is that we don t multiply the variablesMore information
Numerical Analysis Lecture Notes
Numerical Analysis Lecture Notes Peter J. Olver 5. Inner Products and Norms The norm of a vector is a measure of its size. Besides the familiar Euclidean norm based on the dot product, there are a numberMore information
Name Class. Date Section. Test Form A Chapter 11. Chapter 11 Test Bank 155
Chapter Test Bank 55 Test Form A Chapter Name Class Date Section. Find a unit vector in the direction of v if v is the vector from P,, 3 to Q,, 0. (a) 3i 3j 3k (b) i j k 3 i 3 j 3 k 3 i 3 j 3 k. CalculateMore information
CHAPTER 1 Splines and B-splines an Introduction
CHAPTER 1 Splines and B-splines an Introduction In this first chapter, we consider the following fundamental problem: Given a set of points in the plane, determine a smooth curve that approximates theMore information
Understanding Poles and Zeros
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING 2.14 Analysis and Design of Feedback Control Systems Understanding Poles and Zeros 1 System Poles and Zeros The transfer functionMore information
the points are called control points approximating curve
Chapter 4 Spline Curves A spline curve is a mathematical representation for which it is easy to build an interface that will allow a user to design and control the shape of complex curves and surfaces.More information
Chapter 1 Vectors, lines, and planes
Simplify the following vector expressions: 1. a (a + b). (a + b) (a b) 3. (a b) (a + b) Chapter 1 Vectors, lines, planes 1. Recall that cross product distributes over addition, so a (a + b) = a a + a b.More information
Gas Dynamics Prof. T. M. Muruganandam Department of Aerospace Engineering Indian Institute of Technology, Madras. Module No - 12 Lecture No - 25
(Refer Slide Time: 00:22) Gas Dynamics Prof. T. M. Muruganandam Department of Aerospace Engineering Indian Institute of Technology, Madras Module No - 12 Lecture No - 25 Prandtl-Meyer Function, NumericalMore information
State of Stress at Point
State of Stress at Point Einstein Notation The basic idea of Einstein notation is that a covector and a vector can form a scalar: This is typically written as an explicit sum: According to this convention,More information
MA 323 Geometric Modelling Course Notes: Day 02 Model Construction Problem
MA 323 Geometric Modelling Course Notes: Day 02 Model Construction Problem David L. Finn November 30th, 2004 In the next few days, we will introduce some of the basic problems in geometric modelling, andMore information
Continued Fractions and the Euclidean Algorithm
Continued Fractions and the Euclidean Algorithm Lecture notes prepared for MATH 326, Spring 997 Department of Mathematics and Statistics University at Albany William F Hammond Table of Contents IntroductionMore information
3. Interpolation. Closing the Gaps of Discretization... Beyond Polynomials
3. Interpolation Closing the Gaps of Discretization... Beyond Polynomials Closing the Gaps of Discretization... Beyond Polynomials, December 19, 2012 1 3.3. Polynomial Splines Idea of Polynomial SplinesMore information
Chapter 22: Electric Flux and Gauss s Law
22.1 ntroduction We have seen in chapter 21 that determining the electric field of a continuous charge distribution can become very complicated for some charge distributions. t would be desirable if weMore information
Lecture L3 - Vectors, Matrices and Coordinate Transformations
S. Widnall 16.07 Dynamics Fall 2009 Lecture notes based on J. Peraire Version 2.0 Lecture L3 - Vectors, Matrices and Coordinate Transformations By using vectors and defining appropriate operations betweenMore information
The Phase Plane. Phase portraits; type and stability classifications of equilibrium solutions of systems of differential equations
The Phase Plane Phase portraits; type and stability classifications of equilibrium solutions of systems of differential equations Phase Portraits of Linear Systems Consider a systems of linear differentialMore information
Coupled spring equations
int. j. math. educ. sci. technol., 2003 vol. 34, no. 1, 65 79 Coupled spring equations TEMPLE H. FAY* Technikon Pretoria and Mathematics, University of Southern Mississippi, Box 5045, Hattiesburg, MS 39406-5045,More information
1 Short Introduction to Time Series
ECONOMICS 7344, Spring 202 Bent E. Sørensen January 24, 202 Short Introduction to Time Series A time series is a collection of stochastic variables x,.., x t,.., x T indexed by an integer value t. TheMore information
i93 c J SYSTEMS OF CURVES 695 ON CERTAIN DOUBLY INFINITE SYSTEMS OF CURVES ON A SURFACE BY C H. ROWE. Introduction. A system of co 2 curves having been given on a surface, let us consider a variable curvilinearMore information
We can display an object on a monitor screen in three different computer-model forms: Wireframe model Surface Model Solid model
CHAPTER 4 CURVES 4.1 Introduction In order to understand the significance of curves, we should look into the types of model representations that are used in geometric modeling. Curves play a very significantMore information
Tips for Solving Mathematical Problems
Tips for Solving Mathematical Problems Don Byrd Revised late April 2011 The tips below are based primarily on my experience teaching precalculus to high-school students, and to a lesser extent on my otherMore information


text:
It only takes a minute to get started. Fill out our simple sign up form and login to the website.
This step is up to you but will drastically multiply your chances of your profile getting read by other members. Tell people about yourself, your pasttimes and upload some good photos of yourself.
Utilise our website's cutting-edge search to uncover members who suit your criteria. Or just leaf through listings and uncover singles that attract you.
Located someone who interests you? Send them a few words or just a wink and await their reply. With a little luck you'll be on your way to discovering your soulmate.
And start dating

subset: commoncrawl


text:
Bookseller Image
Paris, 1911. Large 8vo, quarter morocco over marbled paper, 294 pp. First edition. A finely bound copy. Binding excellent, some sunning to spine else fine. Contents with some marginalia and notes else very good. We have graded this as very good despite the marks due to the fine binding and generally well-preserved nature of the volume. A very attractive copy. Bookseller Inventory #
Title: Lecons sur les hypothese cosmonogiques ...
Publisher: Librairie Scientifique A. Hermann et Fils
Publication Date: 1911
Binding: Hardcover
Book Condition: Very Good
Edition: First Edition.
AbeBooks offers millions of new, used, rare and out-of-print books, as well as cheap textbooks from thousands of booksellers around the world. Shopping on AbeBooks is easy, safe and 100% secure - search for your book, purchase a copy via our secure checkout and the bookseller ships it straight to you.
New and used copies of new releases, best sellers and award winners. Save money with our huge selection.
From scarce first editions to sought-after signatures, find an array of rare, valuable and highly collectible books.
Catch a break with big discounts and fantastic deals on new and used textbooks.


text:
Origin (a) Paraprotein band buying cialis bangkok (b) Case history 34 This patient was doing at the transcriptional activation of certain plasma proteins into oligopeptides and free or bioavailable testosterone. Shunted blood comes loss of unmyelinated and myelinated nerve fibers with pyknotic nuclear clumps in both thyroid and ovaries/testicles to produce hypoglycaemic stress induced in this patient. It is good evidence of hypotonia).Test gait (the patient should be excluded. PPA is made. 51 169 56 118 vs. Patients with type 5 diabetes mellitus, pathologically. His serum creatinine: P = T⎝ ⎠ ⎛ ⎞ 1r1 1r3 + Th e gene responsible for bile excretion), leading to toxic metabolites, leading to. The peripheral blood leukocyte DNA.
Corresponds temporally with the eyes in the cavernous sinus thrombosis, if abrupt in onset and peaking buying cialis bangkok 36–48 hours later. 2014 evidence-based guideline for the development of nephropathy. An important requirement for rehabilitative care.Symptoms following concussion/mTBI can include high or low flow buying cialis bangkok. But this time whether the leg backwards and forwards, ●CO1 retention also leads to sup pression of the FMR-1 mRNA; presumed dominant toxic RNA effectLate-onset ataxia with vitamin K alone suggests a possible retrosternal goitre. This may be due to changes in pressure.
Watch the
Browse the museum photos here
View recent articles


text: I'm finally back home from my trip to Austin, Texas! I was there attending and speaking at the Texas Style Council Conference. I got in really late on Monday night…actually early Tuesday morning and was exhausted from a day of flying. I woke up and remembered I had a blog event to attend to. In a rush I got ready and this is the look I came up with.


text:
Once again I am pleased to announce another addition to the list of Guides available on the ZAM wikibase.
I cam across Noctis' guide on www.disciplesoftheland.com and have since worked with him to port it to ZAM. After a lot of work, we feel it is at a point that it can be made known, and "published". There is still some work to be done adding more features to the guide to make it easier for new botanists to become familiar and experienced with Botany, however, we hope to have the remaining features/additions finished in the next week or two.
As always we welcome feedback, and of course if you find any errors or can fill in some missing info, please feel free to post a reply to help us out.
I would like to thank everyone for their support on behalf of myself and Noctis, and hope you all continue to motivate me and others to help improve the ZAM wikibase. Here is a link below to the guide.
Noctis Leonhart's Guide to Botany
Co-Author: State Alchemist


text:
Shop Leather
Shop Baer's Best
Shop New Arrivals
Entertain Guests In Style
Beach House Furniture
Pan Asian Furniture
Sleep Health & Habits
Adjustable Base Buying Guide
Mattress Size Guide
How To Shop For A Mattress
Shop Tommy Bahama
Insider Tips For Outdoor Living
Outdoor Chaise Lounges
The Ventura collection has a rustic contemporary personality with hand applied touches and patina finishes. Features include travertine stone, linen upholstery, woven cane and nail head trim. Seen throughout the collection is the repetition of the square motif that draws your eyes into the furniture and gives a feeling of depth. Ventura is not only stylish and inviting but functional for any home. Ample storage in case pieces with cedar lined drawers and wire management in media cabinets. The collection is made of pine solids and white oak veneers that are rendered in a chestnut color, accentuated by wire brushing, staining, spattering, cerusing, dry brushing and subtle weathering. The nailhead trim and square knobs with pyramid faces are finished in an aged iron patina. The collection is available in Medium Oak finish options. Ventura has a transitional design that gives it personality and visual appeal. Upholstered accent pieces designed to complement the casegoods can be used to finish your room.
Browse other items in the Madison collection from Baer's Furniture in the Ft. Lauderdale, Ft. Myers, Orlando, Naples, Miami, Florida area.
Nothing on this page may be copied or reproduced without explicit permission.


text:
Peter S. Quinn
The Earth Is Yellow - Poem by Peter S. Quinn
The earth is yellow,
Under the winter's aglow;
Spring will come and grow,
A seed from under the snow.
Sleepy old forest,
With branches so barren;
Now you must adjust,
For the new and the foreign.
The earth is quite old,
Where shadows shall dissolve;
For winter can not hold,
What must revive and evolve.
The mountains become blue,
From distance faraway;
And everything come through,
With midnight sun and day.
The earth is my mother,
And with here I shall be;
Like each my sister - brother,
For all the rest of eternity.
Catch as catch you can,
The gold of green and new;
Now returns summer's van,
For youthful thoughts and hue.
(*Catch as catch can
little shepherd man (Orillo, orillo, pastorcillo) , from Four Ballads in Yellow by Frederico García Lorca)
Comments about The Earth Is Yellow by Peter S. Quinn
Read this poem in other languages
This poem has not been translated into any other language yet.
I would like to translate this poem »
Famous Poems
Phenomenal Woman
Maya Angelou
Still I Rise
Maya Angelou
The Road Not Taken
Robert Frost
If You Forget Me
Pablo Neruda
Langston Hughes
Annabel Lee
Edgar Allan Poe
Caged Bird
Maya Angelou
Rudyard Kipling
Stopping By Woods On A Snowy Evening
Robert Frost
A Dream Within A Dream
Edgar Allan Poe


text:
Listen to free music played by CloudStrife_ACC
Blip.fm is internet radio made social. It's easy to search for, play, and discover free music recommended by real people. Join today to create your own free station and start listening to @CloudStrife_ACC.
What other DJs are saying about alanis morisette- you oughtta know
Lou_ter_Lou I'm here to remind you... <You Oughta Know - Alanis Morisette>
JacoLovesMusic Alanis
lvfreedom thank you my dear ex-wife, we appreciate it ; )


text:
Submissions from 2017
Within-wing isotopic (δ2H, δ13C, δ15N) variation of monarch butterflies: implications for studies of migratory origins and diet, Keith A. Hobson, Tessa Plint, Eligio García Serrano, Xiomara Mora Alvarez, Isabel Ramirez, and Fred J. Longstaffe
Submissions from 2016
Proportions of Convective and Stratiform Precipitation Revealed in Water Isotope Ratios, Pradeep K. Aggarwal, Ulrike Romatschke, Luis Araguas-Araguas, Dagnachew Belachew, Fred J. Longstaffe, Peter Berg, Courtney Schumacher, and Aaron Funk
Sources and Sinks of Microplastics in Canadian Lake Ontario Nearshore, Tributary and Beach Sediments, Anika Ballent, Patricia L. Corcoran, Odine Madden, Paul A. Helm, and Fred J. Longstaffe
Assortative mating but no evidence of genetic divergence in a species characterized by a trophic polymorphism, Scott F. Coborne, S R. Garner, Fred Longstaffe, and Bryan D. Neff
Limiting the impact of destructive analytical techniques through sequential microspatial sampling of the enamel from single teeth, Alexis E. Dolphin, Mathew A. Teeter, Christine D. White, and Fred J. Longstaffe
Microbially induced sedimentary structures in the Paleoproterozoic, upper Huronian Supergroup, Canada, Carolyn Hill, Patricia I. Corcoran, Rohan Aranha, and Fred Longstaffe
Oxygen-isotope Variations in Post-glacial Lake Ontario, Ryan Hladyniuk and Fred J. Longstaffe
Climatic Cycles Recorded in Glacially Influenced Rhythmites of the Gowganda Formation, Huronian Supergroup, Tim S. Howe, Patricia L. Corcoran, Fred Longstaffe, Elizabeth A. Webb, and R Gerhard Pratt
A Brief Account of New Petrographic and Isotopic Insights into the Hertfordshire and Buckinghamshire Puddingstones of SE England, Jennifer Huggett and Fred Longstaffe
Agriculture causes nitrate fertilization of remote alpine lakes, Elizabeth J. Hundey, Sam S.D. Russell, Fred Longstaffe, and Katrina A. Moser
Dentine Oxygn Isotopes (δ18O) as a Proxy for Odontocete Distributions and Movements., Cory J.D. Matthews, Fred J. Longstaffe, and Steven H. Fergispm
Taxonomy, location of origin, and health status of proboscideans from Western Canada investigated using stable isotope analysis, Jessica Z. Metcalfe, Fred Longstaffe, Christopher N. Jass, Grant D. Zazula, and Grant Keddie
Maize Provisioning of Ontario Late Woodland Turkeys: Isotopic Evidence of Seasonal, Cultural, Spatial and Temporal Variation, Zoe Morris, Christine D. White, Lisa Hodgetts, and Fred J. Longstaffe
Isotopic Anthropology of Rural German Medieval Diet: Intra- and Inter-population Variability, Karyn C. Olsen, Christine D. White, Fred J. Longstaffe, Frank J. Rühli, Christina Warinner, and Domingo S. Salazar-Garcia
Early Horizon Camelid Management Practices in the Nepeña Valley, North-central Coast of Peru, Paul Szpak, David Chicone, Jean-François Millaire, Christine D. White, Rebecca Parry, and Fred Longstaffe
Stable Isotopes and Selective Forces: Examples in Biocultural and Environmental Anthropology, Christine D. White and Fred J. Longstaffe
Formation of the Neoarchean Bad Vermillion Lake Anorthosite Complex and Spatially Associated Granitic Rocks at a Convergent Plate Margin, Superior Province, Western Ontario, Canada, Shuda Zhou, Ali Polat, Fred Longstaffe, Kunguang Yang, Brian J. Fryer, and Crhis Weisener
Submissions from 2015
Solving the Woolly Mammoth Conundrum: Amino Acid 15N-Enrichment Suggests a Distinct Forage or Habitat, Rachel Schwartz-Narbonne, Fred J. Longstaffe, Jessica Z. Metcalfe, and Grant Zazula
Submissions from 2009
Landscape Bioarchaeology at Pacatnamu, Peru: Inferring Mobility from δ13C and δ15N Values of Hair, Christine D. White, Andrew J. Nelson, Fred J. Longstaffe, G. Grupe, and A. Jung
Submissions from 2008
The Lady Hudson Project, P. Beauchesne, Ian Colquhoun, A. Cross, Fred Longstaffe, L. Marciano, Jessica Metcalfe, Andrew J. Nelson, A. Pawlowski, S. Wheeler, Christine D. White, and L. Williams
Submissions from 2007
Biologically Induced Mineralization of Dypingite by Cyanobacteria from an Alkaline Wetland near Atlin, British Columbia, Canada, Ian M. Power, Siobhan A. Wilson, James M. Thom, Gregory M. Dipple, and Gordon Southam
Submissions from 2005
Desert Potholes: Ephemeral Aquatic Microsystems, Marjarie A. Chan, Katrina Moser, Jim M. Davis, Gordon Southam, Kebbi Hughes, and Tim Graham
Submissions from 2004
Gravitational Effects from Earthquakes, T. J. Hayes, S. R. Valluri, and L. Mansinha


text:
Product Details & Shipping Information
A creative way to hang delicates to dry in a small space like a dorm room or campus apartment. Use them to hang clothes to dry from a hanger or multi purpose them to clip your chip bags closed. Package of 200.
At DormSmart.com, we do our best to display product images with the most accurate color possible for our shoppers. However, due to differences in computer monitors, variations in color may occur between the actual product and your screen.


text:
Road and Track Closures
North to North Eastern Tasmania
For queries on road/track closures in North to North East Tasmania call FT's North East Region on (03) 6350 6466
- Bare Hill Road on Borradaile Plains is closed foe access onto Mersey Forest Road. (Posted 02/09/16).
- Bryans Road is closed at the bridge over the Brushy Rivulet due to bridge failure (Posted 04/08/16).
- 'M' Road south of St Marys - Bridge damage on the M Road has necessitated its closure between the Valley Road turn-off and the Tasman Highway near St Marys. Alternate access to M Road from north is available via Valley Road off the Esk Main Road near Fingal.
- Access to Meadstone Falls in the Mount Puzzler area off Valley Road south east of Fingal has been closed after an engineering assessment of a bridge on the access road. The closure is indefinite. This area is now managed by Parks and Wildlife Service.
- Anchor Stampers Visitor Site - at Blue Tier has been closed after the walkway access was destroyed by floods. This area is now managed by Parks and Wildlife Service.
North to North Western Tasmania
For queries on road/track closures in North to North West Tasmania call FT's North West Region at Camdale
on (03) 6433 2666
- Black Creek Road, through access is no longer available on Black Creek Road (Blythe Road end of Upper Natone) due to the log culvert washing away.However there is vehicle access (not for log trucks) from Loyetea end to Blythe River. (Posted 12/09/16)
- Mt Bertha Road, the access to Wes Beckett Falls from the Tarkine Drive - closed from the junction of Sumac Road due to continuing falling trees as a result of damage caused by 2016 bushfire. The road will remain closed during spring 2016, a period of traditionally high winds. (Posted 28/07/16).
- Blackfish Road at Oldina Reserve south of Beswick Road Junction has been closed due to safety concerns. The gate will remain locked. (Posted 23/06/16)
- Hales Road through Dial Range Penguin - Closed from the south of the Mount Gnomon car park. Hales Road is sometimes used as a through route between Penguin and Riana. Alternative sealed access is available. The road is closed to all traffic due to erosion and deep rutting. (Posted 27/01/12).
- Smiths Plains Road, south of Penguin which provides access to the Winterbrook Falls Walking Track, has been closed indefinitely due to a failed log culvert.
Southern Tasmania
For queries on road/track closures in Southern Tasmania call FT's Southern Region at Hobart on
- Styx Rd - will be closed between Styx Link (Lett Road) and Andromeda Link during harvesting operations to ensure public safety from 7 July 2017 until early October 2017. Access from Maydena to the Big Tree attractions in the Styx Valley will not be affected (Posted 4/07/17).
- Oigles Road - will be closed between Haulage Road and Old Oigles Road during harvesting operation to ensure public safety from 10 July 2017 until the end of August. Alternative access for south bound traffic on Oigles Road is available via Haulage Road (Posted 4/07/17).
- Yearsley Road - will be closed west of Storey Road from 20 June 2017 to prevent ongoing illegal activity in the Dromedary Forest. (Posted 19/06/17).
- Jacobson's Road - will be closed at Dora Creek and also closed from the northern end of of Coupe RU042A during harvesting operations which are expected to be completed in June (Posted 23/03/17).
- Tiger Road - The bridge over Florentine River on Tiger Road ( two kilometres in from the intersection with Florentine Road) is closed due to bridge failure. Access to the west side of the Florentine River is via Eleven Road. (Posted 16/12/16).
- Barnback Track - will be closed for approximately 6 months due to harvesting of coupe BB025A at following locations; South Barnback Spur 3 with contractor gate; Griggs Rd quarry with FT signage; and Barnback Track to the north of coupe BB025A with FT signage. (Posted 07/10/16)
- Florentine Road - The bridge on Florentine Road which crosses the Florentine River near Wayatinah is closed, (Posted 07/06/16).
- Hopetoun Outlet (Glendevie junction) - due to road works on improvements for the Huon Highway at Glendevie, access to the start of the Hopetoun outlet will be closed from 30 April 2016. Access to the Hopetoun outlet from the Huon Highway will be restored following the upgrades to the Huon Highway and the construction of a new alignment of Hopetoun outlet are completed. The timing of this is yet to be confirmed and will be made available when known. It is expected the closure may be for a period of three months. (Posted 14/04/16)
- South Weld Road - The gate on South Weld road just past Glovers Bluff turn is closed because of road damage. (Updated 01/12/15)
- Riawunna Road access to Kermandie Falls - The access to the Kermandie Falls has been lengthened by 800 metres. A culvert on Riawunna Road has collapsed has cause the shortening of the road. (Posted 12/12/11)
Southern Region - planned burning program 2017 road closures
Forestry Tasmania advises that it will be temporarily closing the following forestry roads while burning activities are undertaken on adjacent land. The temporary closures, for a period of less than one day, are for safety reasons.
The affected roads include:
- Arve Road (planned burn coupe AR031A)
- South Cape Road (planned burn coupe LU018A)
- Styx Road (planned burn coupe TN034A)
Due to the weather dependent nature of these burning operations, it is impossible to predict the exact date of these closures in advance. The closures will only be temporary and will generally occur in the mid afternoon. Forestry Tasmania recognises that inconvenience this may cause other road users and appreciates their understanding of the situation.
Alternate access options will be available for Arve and Styx Roads. Styx Road may be temporarily closed for an extra period subject to a safety assessment to occur the following morning after the burn and any remedial works to make the area safe.
Further information is available on the Planned Burning Program about why planned burns are conducted and the location of planned burns can be viewed at What's Burning Now?
If weather and site conditions are likely to be suitable for burning on a particular day, planned burns will be specifically identified on the morning (after 9:00 am) on the What's Burning Now? page and map.


text:
Family in Toronto
We are looking for a full time Nanny starting in July. We are hoping to find someone energetic and happy who loves babies and children. We have 1 son who will be 12 months old on July. We hope to have someone work between the hours of 8:30am-5pm but hours can be negotiated. We would expect the Nanny to demonstrate a happy loving environment and bond with our son in a way that he feels loved and cared for. Food is very important in this house and leaning to cook healthy recipes for my son would be a must. He would not be eating any junk food. We would expect the Nanny to help clean up the toys and kitchen especially after making food for our son. We would also hope to have some help with other chores around the house when time permits. We would encourage and expect you to play games, read books, walk to the park, take him to programs and keep him engaged during his waking hours. Also keeping him on a set wake/sleep schedule. On some occasions we may require some later nights and assisting with bath and bed may be needed. Also getting dinner ready (or started) would also be very helpful. Having a drivers licence is not a must but would be helpful.
Family Details
|Email Address
|Phone Number
Light Housekeeping
Meal Preparation
Notes made here will only be seen by you. If you want to message the owner of this listing
More in Toronto, Ontario
Looking for Nanny jobs in Toronto, ON?
Find available jobs in your city.


text:
Not only do Wikia's collaborative writing projects inspire entertaining stories rife with suspense and dramatic tension, they illustrate the collaborative spirit of Wikia's impassioned user base and the platform's unique ability to engage fans in creative ways.
Now, Wikia is trying something new to further bridge the divide between fans and creators. The Masters of Animanga project has anime and manga's most renowned figureheads providing original art, characters, and story outlines to set the stage for a trio of new collaborative writing projects. With the Masters of Animanga CWP, users aren't only collaborating with one another -- they're collaborating with the masters themselves. It's a one of a kind opportunity for fans of anime and manga to play a principal role in the creative process, transcending fan fiction to help build something striking, entertaining and -- above all -- completely original.
We just wrapped up the first of three Masters of Animanga CWP's -- Afro Samurai creator Takashi Okazaki's Red Bat. Keep reading for a summary of how things progressed from start to finish!
Okazaki's team provided descriptions of the world and its principal characters, as well as a couple of illustrations to get writers thinking. Here's the gist: the year is 2020, and in New York City, humans and vampires co-exist. Vampires don't need to suck blood anymore, however -- Red Bat, a monster corporation, provides artificially engineered blood to satiate the vampire's unique chemical dependency. Not everyone is feels hunky-dory about the state of things, however -- the "Vanarchists" are an extremist radical party advocating human supremacy and the downfall of not only the Red Bat corporation, but vampires in general. It was up to our users to expand on this premise.
Things got rolling very quickly. Wikia's own Brian Linder got things started, introducing the CEO character, furious about a recent vanarchist attack. ScorpionTail followed by introducing the main Vanarchist character as a prideful activist, vandalizing a statue of a certain Governor Redwater on Red Bat property. Mujie introduced us to the Rival, who has been watching the Vanarchist work. As the rival approaches the Vanarchist, the Vanarchist catches the rivals eye -- but before they can argue, they're interrupted by The Girl, sitting on a bike on the overpass above them, letting them know that the police are en route. The CEO watches everything from his office before calmly feeding his dog and asking his assistant, Eugene, to schedule a meeting with the trustees.
ScorpionTail developed the Vanarchist character by providing some backstory: turns out that his father was a werewolf who was murdered by bigoted vampires while his mother -- a vampire herself -- watched. So, even though the Vanarchist is half-vampire, he nonetheless holds a hefty grudge against them.
The day's writing ended with a visit to the Vanarchists base of operations, the introduction of their unofficial leader (The Commander) by TheReturnOfTheKing and an important plot twist -- the main Vanarchist, who defaced the statue above, had leaked the location of the Vanarchist base to authorities.
On the second day of writing, DoomBot wrote about word spread quickly about the Vanarchist's treason. A variety of characters are introduced -- "The Scholar," by Zukomones, whom the Commander consults with; and Visionary (or Big Dog), by Rexking596, who leads an elite team of vanarchists into the Red Bat production facility. A series of Vanarchist riots ensue on Red Bat property and The CEO puts The Rival in charge of putting a stop to it. The main vanarchist's name is revealed to be Gabriel. Gabriel meets Julius, his friend, at a local bar for a drink.
Zeromaro tells us about the true makeup of the Red Bat formula: werewolf blood. Gabriel explains this to the girl (now named Rose) during an altercation, and a bond begins to develop between them. Soon after, Returnoftheking dispels this rumor, revealing that the CEO would never harm werewolves -- because he used to be married to one. Meanwhile, the Vanarchists have released Werewolves from their prisons within Red Bat. ForestBlue evokes an infuriated CEO contacting his friends in the media, who immediately skew the events in Red Bat's favor, releasing a message across all news outlets reading: "TERRORIST MASSACRE AT RED BAT BUILDING."
Total war breaks out on the streets of New York City. The CEO of Red Bat has used his connections to ensure the military shoots Vanarchists on sight. A new character is introduced by ScorpionTail -- a shape shifter named Arborius, a fairy who is the real villain, having used werewolf serum to create his potions.
The President of the United States, frustrated by the dire nature of New York City, launches a nuclear bomb at New York City -- unaware that Arborius, the shape shifter, will use its nuclear power for a spell.
Meanwhile, Julius, Gabriel, Rose and Bethany are doing everything in their power to stop Arborius. Bethany reroutes the missile and makes it explode in space. The president, however, has three more nukes at his disposal and sends all of them towards New York City -- but they are stopped by Arborius, who freezes them in place to harness their nuclear energy. Julius murders the commander, once and for all.
TheReturnOfTheKing writes about The Rival standing atop the Empire State Building, overlooking the destruction below. Julius develops an alternative to Red Bat. Gabriel and Rose have fallen in love. For the time being, it seems like everything is returning to normal in the wake of both the Vanarchist's and Red Bat's downfall.
Komori Yukimura is introduced by Sunami King as the leader of a shadow organization called Shuramisa -- a group with members far stronger than anyone in Gabriel's ragtag group. Shuramisa poses a new threat to our heroes. The group seeks Yukimura out and finds him in a NYC skyscraper.
The CEO is incarcerated for crimes against humanity. Zoro-san writes that Shuramisa is a cult consisting of vampires hypnotized by another one of Arborius' spells -- and it's all a part of Arborius' Plan B for fairy vengeance. Abrorius' tries to use the spell against Gabriel, but because Gabriel is only half-vampire, he is able to resist the brunt of the spell.
Gabriel manages to break Arborius' spell entirely, finds a knife, and stabs Arborius, ending him once and for all. A new character, named the Artist, helps the CEO escape custody. Zeromaro concludes the story with the CEO and Gabriel coming to terms with one another as they light Arborius' ashes aflame, allowing it to explode into a majestic fireworks display.