0000000000417088
AUTHOR
Fernando Casas
High-order Runge–Kutta–Nyström geometric methods with processing
Abstract We present new families of sixth- and eighth-order Runge–Kutta–Nystrom geometric integrators with processing for ordinary differential equations. Both the processor and the kernel are composed of explicitly computable flows associated with non trivial elements belonging to the Lie algebra involved in the problem. Their efficiency is found to be superior to other previously known algorithms of equivalent order, in some case up to four orders of magnitude.
Efficient numerical integration of neutrino oscillations in matter
A special purpose solver, based on the Magnus expansion, well suited for the integration of the linear three neutrino oscillations equations in matter is proposed. The computations are speeded up to two orders of magnitude with respect to a general numerical integrator, a fact that could smooth the way for massive numerical integration concomitant with experimental data analyses. Detailed illustrations about numerical procedure and computer time costs are provided.
Unitary transformations depending on a small parameter
We formulate a unitary perturbation theory for quantum mechanics inspired by the LieDeprit formulation of canonical transformations. The original Hamiltonian is converted into a solvable one by a transformation obtained through a Magnus expansion. This ensures unitarity at every order in a small parameter. A comparison with the standard perturbation theory is provided. We work out the scheme up to order ten with some simple examples.
New Families of Symplectic Runge-Kutta-Nyström Integration Methods
We present new 6-th and 8-th order explicit symplectic Runge-Kutta-Nystrom methods for Hamiltonian systems which are more efficient than other previously known algorithms. The methods use the processing technique and non-trivial flows associated with different elements of the Lie algebra involved in the problem. Both the processor and the kernel are compositions of explicitly computable maps.
Geometric factors in the adiabatic evolution of classical systems
Abstract The adiabatic evolution of the classical time-dependent generalized harmonic oscillator in one dimension is analyzed in detail. In particular, we define the adiabatic approximation, obtain a new derivation of Hannay's angle requiring no averaging principle and point out the existence of a geometric factor accompanying changes in the adiabatic invariant.
The Magnus expansion and some of its applications
Approximate resolution of linear systems of differential equations with varying coefficients is a recurrent problem, shared by a number of scientific and engineering areas, ranging from Quantum Mechanics to Control Theory. When formulated in operator or matrix form, the Magnus expansion furnishes an elegant setting to build up approximate exponential representations of the solution of the system. It provides a power series expansion for the corresponding exponent and is sometimes referred to as Time-Dependent Exponential Perturbation Theory. Every Magnus approximant corresponds in Perturbation Theory to a partial re-summation of infinite terms with the important additional property of prese…
A pedagogical approach to the Magnus expansion
Time-dependent perturbation theory as a tool to compute approximate solutions of the Schrodinger equation does not preserve unitarity. Here we present, in a simple way, how the Magnus expansion (also known as exponential perturbation theory) provides such unitary approximate solutions. The purpose is to illustrate the importance and consequences of such a property. We suggest that the Magnus expansion may be introduced to students in advanced courses of quantum mechanics.
Floquet theory: exponential perturbative treatment
We develop a Magnus expansion well suited for Floquet theory of linear ordinary differential equations with periodic coefficients. We build up a recursive scheme to obtain the terms in the new expansion and give an explicit sufficient condition for its convergence. The method and formulae are applied to an illustrative example from quantum mechanics.
A Unifying Framework for Perturbative Exponential Factorizations
We propose a framework where Fer and Wilcox expansions for the solution of differential equations are derived from two particular choices for the initial transformation that seeds the product expansion. In this scheme, intermediate expansions can also be envisaged. Recurrence formulas are developed. A new lower bound for the convergence of theWilcox expansion is provided, as well as some applications of the results. In particular, two examples are worked out up to a high order of approximation to illustrate the behavior of the Wilcox expansion.
Non-Adiabatic Aspects of Time-Dependent Hamiltonian Systems
Extreme adiabatic behavior furnishes great simplification in the treatment of linear time-dependent Hamiltonian systems. But the actual time variation of the parameters is only finitely, rather than infinitely, slow. Then one is forced to consider corrections to the adiabatic limit.
Magnus and Fer expansions for matrix differential equations: the convergence problem
Approximate solutions of matrix linear differential equations by matrix exponentials are considered. In particular, the convergence issue of Magnus and Fer expansions is treated. Upper bounds for the convergence radius in terms of the norm of the defining matrix of the system are obtained. The very few previously published bounds are improved. Bounds to the error of approximate solutions are also reported. All results are based just on algebraic manipulations of the recursive relation of the expansion generators.
Analysis of multipactor RF breakdown in a waveguide containing a transversely magnetized ferrite
In this paper, the multipactor RF breakdown in a parallel-plate waveguide partially filled with a ferrite slab magnetized normal to the metallic plates is studied. An external magnetic field is applied along the vertical direction between the plates in order to magnetize the ferrite. Numerical simulations using an in-house 3-D code are carried out to obtain the multipactor RF voltage threshold in this kind of structures. The presented results show that the multipactor RF voltage threshold at certain frequencies becomes considerably lower than for the corresponding classical metallic parallel-plate waveguide with the same vacuum gap