Dans cette thèse, on s'intéresse à deux aspects différents des systèmes dynamiques intégrables. La première partie est dévouée à l'étude de trois familles de systèmes hamiltoniens intégrables : les systèmes de pliage de Kapovich et Millson sur les espaces de modules de polygones 3D de longueurs de côtés fixées, les systèmes de Gelfand-Cetlin introduits par Guillemin et Sternberg sur les orbites coadjointes du groupe de Lie U(n), et une famille de systèmes définie par Nohara et Ueda sur la variété grassmannienne Gr(2,n). Dans chaque cas on montre que les fibres singulières de l'application moment sont des sous-variétés plongées et on en donne des modèles géométriques sous la forme de variétés quotients. La deuxième partie poursuit une étude initiée par Zung et Minh sur les actions totalement hyperboliques de Rn sur des variétés compactes de dimension n, qui apparaissent naturellement lors de l'étude des systèmes non-hamiltoniens intégrables dont toutes les singularités sont non-dégénérées. On s'intéresse au flot engendré par l'action d'un vecteur générique de Rn. On donne une définition d'indice pour ses singularités qu'on relie à la théorie de Morse classique, et on utilise ce flot pour obtenir des résultats sur le nombres d'orbites de dimension donnée. Une étude plus poussée est effectuée en dimension 2, et en particulier sur la sphère S2, où les orbites de l'action dessinent un graphe plongé dont on analyse la combinatoire. On termine en construisant explicitement des exemples d'actions hyperboliques en dimension 3 sur la sphère S3 et dans l'espace projectif RP3. / In this thesis, we are interested in two different aspects of integrable dynamical systems. The first part is devoted to the study of three families of integrable Hamiltonian systems: the systems of bending flows of Kapovich and Millson on the moduli spaces of 3D polygons with fixed side lengths, the Gelfand-Cetlin systems introduced by Guillemin and Sternberg on the coadjoint orbits of the Lie group U(n), and a family of integrable systems defined by Nohara and Ueda on the Grassmannian Gr(2,n). In each case we prove that the fibers of the momentum map are embedded submanifolds for which we give geometric models in terms of quotients manifolds. In the second part we carry on with a study initiated by Zung and Minh of the totally hyperbolic actions of R^n on compact n-dimensional manifolds that appear naturally when investigating integrable non-hamiltonian systems with nondegenerate singularities. We study the flow generated by the action of a generic vector in Rn. We define a notion of index for its singularities and we use this flow to obtain results on the number of orbits of given dimension. We investigate further the 2-dimensional case, and more particularly the case of the sphere S2, where the orbits of the action draw an embedded graph of which we analyse the combinatorics. Finally, we provide explicit examples of totally hyperbolic actions in dimension 3, on the sphere S3 and on the projective space RP3.

This paper describes LAPACK-based Fortran 77 subroutines for the reduction of a Hamiltonian matrix to square-reduced form and the approximation of all its eigenvalues using the implicit version of Van Loan's method. The transformation of the Hamilto- nian matrix to a square-reduced Hamiltonian uses only orthogonal symplectic similarity transformations. The eigenvalues can then be determined by applying the Hessenberg QR iteration to a matrix of half the order of the Hamiltonian matrix and taking the square roots of the computed values. Using scaling strategies similar to those suggested for algebraic Riccati equations can in some cases improve the accuracy of the computed eigenvalues. We demonstrate the performance of the subroutines for several examples and show how they can be used to solve some control-theoretic problems.

A new method is presented for the numerical computation of the generalized eigen- values of real Hamiltonian or symplectic pencils and matrices. The method is strongly backward stable, i.e., it is numerically backward stable and preserves the structure (i.e., Hamiltonian or symplectic). In the case of a Hamiltonian matrix the method is closely related to the square reduced method of Van Loan, but in contrast to that method which may suffer from a loss of accuracy of order sqrt(epsilon), where epsilon is the machine precision, the new method computes the eigenvalues to full possible accuracy.

La première détection des ondes gravitationnelles en 2015 a ouvert un nouveau plan d'étude pour l'astrophysique des étoiles binaires compactes. En utilisant les données des détections faites par les détecteurs terrestres advanced LIGO et advanced Virgo, il est possible de contraindre les paramètres physiques de ces systèmes avec une analyse Bayésienne et ainsi approfondir notre connaissance physique des étoiles binaires compactes. Cependant, pour pouvoir être en mesure d'obtenir de tels résultats, il est essentiel d’avoir des algorithmes performants à la fois pour trouver les signaux de ces ondes gravitationnelles et pour l'estimation de paramètres. Le travail de cette thèse a ainsi été centré autour du développement d’algorithmes performants et adaptées au problème physique à la fois de la détection et de l'estimation des paramètres pour les ondes gravitationnelles. La plus grande partie de ce travail de thèse a ainsi été dédiée à l'implémentation d’un algorithme de type Hamiltonian Monte Carlo adapté à l'estimation de paramètres pour les signaux d’ondes gravitationnelles émises par des binaires compactes formées de deux étoiles à neutrons. L'algorithme développé a été testé sur une sélection de sources et a été capable de fournir de meilleures performances que d'autres algorithmes de type MCMC comme l'algorithme de Metropolis-Hasting et l'algorithme à évolution différentielle. L'implémentation d'un tel algorithme dans les pipelines d’analyse de données de la collaboration pourrait augmenter grandement l'efficacité de l'estimation de paramètres. De plus, il permettrait également de réduire drastiquement le temps de calcul nécessaire, ce qui est un facteur essentiel pour le futur où de nombreuses détections sont attendues. Un autre aspect de ce travail de thèse a été dédié à l'implémentation d'un algorithme de recherche de signaux gravitationnelles pour les binaires compactes monochromatiques qui seront observées par la future mission spatiale LISA. L'algorithme est une mixture de plusieurs algorithmes évolutionnistes, avec notamment l'inclusion d'un algorithme de Particle Swarm Optimisation. Cette algorithme a été testé dans plusieurs cas tests et a été capable de trouver toutes les sources gravitationnelles comprises dans un signal donné. De plus, l'algorithme a également été capable d'identifier des sources sur une bande de fréquence aussi grande que 1 mHz, ce qui n'avait pas été réalisé au moment de cette étude de thèse. / The first detection of gravitational waves in 2015 has opened a new window for the study of the astrophysics of compact binaries. Thanks to the data taken by the ground-based detectors advanced LIGO and advanced Virgo, it is now possible to constrain the physical parameters of compact binaries using a full Bayesian analysis in order to increase our physical knowledge on compact binaries. However, in order to be able to perform such analysis, it is essential to have efficient algorithms both to search for the signals and for parameter estimation. The main part of this thesis has been dedicated to the implementation of a Hamiltonian Monte Carlo algorithm suited for the parameter estimation of gravitational waves emitted by compact binaries composed of neutron stars. The algorithm has been tested on a selection of sources and has been able to produce better performances than other types of MCMC methods such as Metropolis-Hastings and Differential Evolution Monte Carlo. The implementation of the HMC algorithm in the data analysis pipelines of the Ligo/Virgo collaboration could greatly increase the efficiency of parameter estimation. In addition, it could also drastically reduce the computation time associated to the parameter estimation of such sources of gravitational waves, which will be of particular interest in the near future when there will many detections by the ground-based network of gravitational wave detectors. Another aspect of this work was dedicated to the implementation of a search algorithm for gravitational wave signals emitted by monochromatic compact binaries as observed by the space-based detector LISA. The developed algorithm is a mixture of several evolutionary algorithms, including Particle Swarm Optimisation. This algorithm has been tested on several test cases and has been able to find all the sources buried in a signal. Furthermore, the algorithm has been able to find the sources on a band of frequency as large as 1 mHz which wasn’t done at the time of this thesis study

Hamilton-Jacobi theory provides a powerful method for extracting the equations of motion out of some given systems in classical mechanics. On occasion it allows some systems to be solved by the method of separation of variables. If a system with n degrees of freedom has 2n - 1 constants of the motion that are polynomial in the momenta, then that system is called superintegrable. Such a system can usually be solved in multiple coordinate systems if the constants of the motion are quadratic in the momenta. All superintegrable two dimensional Hamiltonians of the form H = (p_x)sup2 + (p_y)sup2 + V(x,y), with constants that are quadratic in the momenta were classified by Kalnins et al [5], and the coordinate systems in which they separate were found. We discuss Hamilton-Jacobi theory and its development from a classical viewpoint, as well as superintegrability. We then proceed to use the theory to find equations of motion for some of the superintegrable Hamiltonians from Kalnins et al [5]. We also discuss some of the properties of the Poisson algebra of those systems, and examine the orbits.

In many applications from industry and technology computer simulations are performed using models which can be formulated by systems of differential equations. Often the equations underlie additional algebraic constraints. In this context we speak of descriptor systems. Very important characteristic values of such systems are the $L_\infty$-norms of the corresponding transfer functions. The main goal of this thesis is to extend a numerical method for the computation of the $L_\infty$-norm for standard state space systems to descriptor systems. For this purpose we develop a numerical method to check whether the transfer function of a given descriptor system is proper or improper and additionally use this method to reduce the order of the system to decrease the costs of the $L_\infty$-norm computation. When computing the $L_\infty$-norm it is necessary to compute the eigenvalues of certain skew-Hamiltonian/Hamiltonian matrix pencils composed by the system matrices. We show how we extend these matrix pencils to skew-Hamiltonian/Hamiltonian matrix pencils of larger dimension to get more reliable and accurate results. We also consider discrete-time systems, apply the extension strategy to the arising symplectic matrix pencils and transform these to more convenient structures in order to apply structure-exploiting eigenvalue solvers to them. We also investige a new structure-preserving method for the computation of the eigenvalues of skew-Hamiltonian/Hamiltonian matrix pencils and use this to increase the accuracy of the computed eigenvalues even more. In particular we ensure the reliability of the $L_\infty$-norm algorithm by this new eigenvalue solver. Finally we describe the implementation of the algorithms in Fortran and test them using two real-world examples.

Soit M une variété lisse fermée et considérons sont fibré cotangent T*M muni de la structure symplectique usuelle induite par la forme de Liouville. Une hypersurface S de T*M$ est dite étoilée fibre par fibre si pour tout point q de M, l'intersection Sq de S avec la fibre au dessus de q est le bord d'un domaine étoilé par rapport à l'origine 0q de la fibre T*qM. Un flot est naturellement associé à S, il s'agit de l'unique flot généré par le champ de Reeb le long de S, le flot de Reeb. <p><p>L'existence d'une orbite orbite fermée du flot de Reeb sur S fut annoncée par Weinstein dans sa conjecture en 1978. Indépendamment, Weinstein et Rabinowitz ont montré l'existence d'une orbite fermée sur les hypersurfaces de type étoilées dans l'espace réel de dimension 2n. Sous les hypothèses précédentes, l'existence d'une orbite fermée fut démontrée par Hofer et Viterbo. Dans le cas particulier du flot géodésique, l'existence de plusieurs orbites fermées fut notamment étudiée par Gromov, Paternain et Paternain-Petean. Dans cette thèse, ces résultats sont généralisés. <p><p>Les résultats principaux de cette thèse montrent que la structure topologique de la variété M implique, pour toute hypersurface étoilée fibre par fibre, l'existence de beaucoup d'orbites fermées du flot de Reeb. Plus précisément, une borne inférieure de la croissance du nombre d'orbites fermées du flot de Reeb en fonction de leur période est mise en évidence. /<p><p>Let M be a smooth closed manifold and denote by T*M the cotangent bundle over M endowed with its usual symplectic structure induced by the Liouville form. A hypersurface S of T*M is said to be fiberwise starshaped if for each point q in M the intersection Sq of S with the fiber at q bounds a domain starshaped with respect to the origin 0q in T*qM. There is a flow naturally associated to S, generated by the unique Reeb vector field R along S ,the Reeb flow. <p><p>The existence of one closed orbit was conjectured by Weinstein in 1978 in a more general setting. Independently, Weinstein and Rabinowitz established the existence of a closed orbit on star-like hypersurfaces in the 2n-dimensional real space. In our setting the Weinstein conjecture without the assumption was proved in 1988 by Hofer and Viterbo. The existence of many closed orbits has already been well studied in the special case of the geodesic flow, for example by Gromov, Paternain and Paternain-Petean. In this thesis we will generalize their results.<p><p>The main result of this thesis is to prove that the topological structure of $M$ forces, for all fiberwise starshaped hypersurfaces S, the existence of many closed orbits of the Reeb flow on S. More precisely, we shall give a lower bound of the growth rate of the number of closed Reeb-orbits in terms of their periods. / Doctorat en Sciences / info:eu-repo/semantics/nonPublished

Abstract In this thesis, the theory of static magnetic resonance spectral parameters of nuclear magnetic resonance (NMR) and electron spin resonance (ESR) spectroscopy is investigated in terms of the molecular Breit-Pauli Hamiltonian, which is obtained from the relativistic Dirac equation via the Foldy-Wouthuysen transformation. A leading-order perturbational relativistic theory of NMR nuclear shielding and spin-spin coupling tensors, and ESR electronic g-tensor, is presented. In addition, the possibility of external magnetic-field dependency of NMR parameters is discussed. Various first-principles methods of electronic structure theory and the role of one-electron basis sets and their performance in magnetic resonance properties in terms of their completeness profiles are discussed. The presented leading-order perturbational relativistic theories of NMR nuclear shielding tensors and ESR electronic g-tensors, as well as the theory of the magnetic-field dependent NMR shielding and quadrupole coupling are evaluated using first-principles wave function and density-functional theories.

This paper describes LAPACK-based Fortran 77 subroutines for the reduction of a Hamiltonian matrix to square-reduced form and the approximation of all its eigenvalues using the implicit version of Van Loan's method. The transformation of the Hamilto- nian matrix to a square-reduced Hamiltonian uses only orthogonal symplectic similarity transformations. The eigenvalues can then be determined by applying the Hessenberg QR iteration to a matrix of half the order of the Hamiltonian matrix and taking the square roots of the computed values. Using scaling strategies similar to those suggested for algebraic Riccati equations can in some cases improve the accuracy of the computed eigenvalues. We demonstrate the performance of the subroutines for several examples and show how they can be used to solve some control-theoretic problems.

The SHIRA method of Mehrmann and Watkins belongs among the structure preserving Krylov subspace methods for solving skew-Hamiltonian eigenvalue problems. It can also be applied to Hamiltonian eigenproblems by considering a suitable transformation. Structure induced shift-and-invert techniques are employed to steer the algorithm towards the interesting region of the spectrum. However, the shift cannot be altered in the middle of the computation without discarding the information that has been accumulated so far. This paper shows how SHIRA can be combined with ideas from Ruhe's Rational Krylov algorithm to yield a method that permits an adjustment of shift after every step of the computation, adding greatly to the flexibility of the algorithm. We call this new method rational SHIRA. A numerical example is presented to demonstrate its efficiency.