<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">Dear All,</span></font><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">&nbsp;&nbsp; &nbsp; Confusion reigns: &nbsp;</span></font><span class="Apple-style-span" style="font-size: 18px; ">As described in the papers</span></div><div><span class="Apple-style-span" style="font-size: 18px; "><br></span></div><div><span class="Apple-style-span" style="font-size: 18px; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">A. Dragt, F. Neri, et&nbsp;al., "LIE ALGEBRAIC TREATMENT OF LINEAR AND NONLINEAR BEAM DYNAMICS",</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">Annual Review of Nuclear and Particle Science 38, p. 455 (1988).</div></span></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px; "><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">A. Dragt, R. Gluckstern, et al., "THEORY OF EMITTANCE INVARIANTS", &nbsp;Lecture Notes in</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">Physics 343: Proceedings of the Joint US-CERN Capri School on Accelerator Physics,</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">Springer Verlag (1989).</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; "><br></div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; "><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">A. Dragt,&nbsp;F. Neri, et&nbsp;al., "GENERAL MOMENT INVARIANTS FOR LINEAR HAMILTONIAN SYSTEMS",&nbsp;</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; ">Physical Review A, 45, p. 2572 (1992).</div><div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font: normal normal normal 12px/normal Helvetica; "><br></div></div></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px; "></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">there are two kinds of invariants, which I call "kinematic" and "dynamic". &nbsp;It is important not to confound them.</span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">&nbsp;&nbsp; &nbsp; The Courant-Snyder invariants, generally treated in the linear approximation as described by Sergei below, &nbsp;are examples of dynamic invariants. &nbsp;Their use requires a knowledge of the equations of motion or, equivalently, a one-period or one-turn map. &nbsp;(They generally apply to periodic systems for which the linear part of the map has all its eigenvalues on the unit circle and distinct, but in principle could also be constructed when the eigenvalues are off the unit circle.) &nbsp;Contrary to the reservations that appear to be &nbsp;expressed by Sergei, they are extendable to the nonlinear case using normal form methods. &nbsp;See, for example, the first paper above and Sections 8.10 and 8.11 of the MaryLie manual available at the Web site</span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><a href="http://www.physics.umd.edu/dsat/">http://www.physics.umd.edu/dsat/</a></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">There are also papers in the literature that extend them, at least in the linear case, to non-periodic systems.</span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px; "><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">&nbsp;&nbsp; &nbsp;Examples of kinematic invariants include what I call "moment invariants" and "eigen emitances". &nbsp;They involve moments of particle distributions and, like Liouville's theorem, make no specific use of the equations of motion save for the symplectic condition. &nbsp;It is these invariants that are the current focus of interest with regard to questions of "emittance partitioning". &nbsp;See the draft paper</span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><span class="Apple-style-span" style="font-family: 'Lucida Grande', helvetica, arial, verdana, sans-serif; font-size: 13px; "><table summary="Additional metadata" style="position: static; z-index: auto; "><tbody><tr><td class="tablecell arxivid" style="padding-top: 0.1em; padding-right: 0.5em; padding-bottom: 0em; padding-left: 0em; vertical-align: top; font-weight: bold; "><a href="http://arxiv.org/abs/1010.1558v2" style="text-decoration: none; ">arXiv:1010.1558v2</a>&nbsp;[physics.acc-ph]</td></tr></tbody></table><br></span></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><span class="Apple-style-span" style="font-family: 'Lucida Grande', helvetica, arial, verdana, sans-serif; font-size: 13px; "><span class="Apple-style-span" style="font-family: Helvetica; font-size: 18px; ">Their computation and use are described in the papers above and Section 8.37 of the MaryLie manual.</span>&nbsp;</span></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">&nbsp;&nbsp; &nbsp; Moment invariants are currently known only for the case of linear maps, but there is some evidence that they should also exist in the &nbsp;nonlinear case. &nbsp;See Chapter 26 (only partially complete) and&nbsp;</span></font><span class="Apple-style-span" style="font-size: 18px; ">Section 6.8.2</span><span class="Apple-style-span" style="font-size: 18px; ">&nbsp;</span><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;">of "Lie Methods ...", also available at the Web site above. &nbsp;</span></font><span class="Apple-style-span" style="font-size: 18px; ">In the past I have tried to construct moment invariants for nonlinear maps by finding Casimir operators for the full Lie algebra of all symplectic maps. &nbsp;The problem is difficult because the usual method for constructing Casimir operators requires that the Killing metric be invertible. &nbsp;See Section 21.11 of "Lie Methods ...". &nbsp;This metric is not invertible in the nonlinear case, and therefore the standard machinery for constructing Casimirs fails. &nbsp;With considerable effort I found an alternate approach with some promising results, but was eventually overwhelmed by algebraic complexity. &nbsp;Moment invariants for the nonlinear case, if they exist at all, will require many pages even to write them down.</span></div><div><span class="Apple-style-span" style="font-size: 18px; "><br></span></div><div><span class="Apple-style-span" style="font-size: 18px; ">&nbsp;&nbsp; &nbsp; I hope these remarks provide some clarification.</span></div><div><span class="Apple-style-span" style="font-size: 18px; "><br></span></div><div><span class="Apple-style-span" style="font-size: 18px; ">Best,</span></div><div><span class="Apple-style-span" style="font-size: 18px; "><br></span></div><div><span class="Apple-style-span" style="font-size: 18px; ">Alex</span></div><div><font class="Apple-style-span" size="5"><span class="Apple-style-span" style="font-size: 18px;"><br></span></font></div><div><div><div>On Mar 10, 2011, at 7:32 PM, Sergei Nagaitsev wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Dear Kirk,<br><br>let me add my two cents to this discussion:<br><br>1. &nbsp;For the beam emittance to be a useful quantity, it needs to be conserved as the beam propagates along the beamline (consider non-interacting particles) through various external electromagnetic fields. &nbsp;To this end, the definition of emittance that reflects such a property is "an ensemble-averaged action". &nbsp;The particle motion is assumed to be integrable, i.e. there exists 3 functionally-independent constants of the motion in involution with the Hamiltonian. &nbsp;Particle actions, expressed through such constants of motion, are also constants of motion. &nbsp;Thus, the average actions (emittances) are conserved. &nbsp;In a linear-focusing transport channel with linear rf focusing, such constants of motion exist and are called Courant-Snyder invariants (2 transverse and 1 longitudinal). &nbsp;In a nonlinear (and generally time-dependent) focusing channel such constants of motion might not exist, therefore, the emittance (as a conserved quantity) is not defined. &nbsp;We may use an approximate expression for the emittance by using linear-only focusing to define actions and then by treating non-linearities as a perturbation. &nbsp;This leads to an effective emittance growth if a beamline has nonlinear elements even though the Liouville's theorem states that the phase-space density is conserved. &nbsp;In some cases this emittance growth is not "real" (irreversible) but just a reflection that we are using an incorrect action definition. &nbsp;Finally, my definition of the emittance (as the average action) is identical to the definition through eigen-values of a sigma matrix only in a case of a linear focusing channel. &nbsp;However, where possible (like in case of a bunch occupying a large portion of an rf bucket) we should use exact actions, not approximate. <br><br>2. When averaging particle actions over the distribution function at a given time t it is useful to remember that a time slice t=const in one frame is not t'=const in another frame because of Lorentz transformations.<br><br>Sergei &nbsp;<br><br>----- Original Message -----<br>From: Kirk T McDonald &lt;<a href="mailto:kirkmcd@Princeton.EDU">kirkmcd@Princeton.EDU</a>&gt;<br>Date: Thursday, March 10, 2011 8:49 pm<br>Subject: Re: [MAP] Liouville's theorem and electromagnetic fields<br>To: alex dragt &lt;<a href="mailto:dragtnb@comcast.net">dragtnb@comcast.net</a>&gt;<br>Cc: MAP List &lt;<a href="mailto:map-l@lists.bnl.gov">map-l@lists.bnl.gov</a>&gt;<br><br><br><blockquote type="cite">Folks,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I have added Alex' paper to DocDB 560. &nbsp;See Appendix A.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">It is gratifying to see that the fact that Liouville's theorem holds <br></blockquote><blockquote type="cite">for <br></blockquote><blockquote type="cite">both mechanical and canonical phase space is "well known to those who <br></blockquote><blockquote type="cite">know".<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The challenge now is to learn how best to use the "freedom" offered to <br></blockquote><blockquote type="cite">us by <br></blockquote><blockquote type="cite">this apparently nonintuitive result.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--Kirk<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--------------------------------------------------<br></blockquote><br></div></blockquote></div><br></div></body></html>