BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Philipp Bader appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We build an efficient and unitary (hence stable) method for the solution of the semi-classical Schrödinger equation subject with explicitly time-dependent potentials. The method is based on a combination of the Zassenhaus decomposition (Bader, Iserles, Kropielnicka & Singh 2014) with the Magnus expansion of the time-dependent Hamiltonian. We conclude with numerical experiments.
We consider the numerical integration of the matrix Hill's equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, the Hill's equations originate from a Hamiltonian function and the fundamental matrix solution is a symplectic matrix. This is a very important property to be preserved by the numerical integrators. In this work we present new sixth-and eighth-order symplectic exponential integrators that are tailored to the Hill's equation. The methods are based on an efficient symplectic approximation to the exponential of high dimensional coupled autonomous harmonic oscillators and yield accurate results for oscillatory problems at a low computational cost. Several numerical examples illustrate the performance of the new methods.
It is a classical theorem of Liouville that Hamiltonian systems preserve volume in phase space. Any symplectic Runge-Kutta method will respect this property for such systems, but it has been shown that no B-Series method can be volume preserving for all volume preserving vector fields (BIT 47 (2007) 351-378 and IMA J. Numer. Anal. 27 (2007) 381-405). In this paper we show that despite this result, symplectic Runge-Kutta methods can be volume preserving for a much larger class of vector fields than Hamiltonian systems, and discuss how some Runge-Kutta methods can preserve a modified measure exactly.
We propose splitting methods for the computation of the exponential of perturbed matrices which can be written as the sum $A=D+\varepsilon B$ of a sparse and efficiently exponentiable matrix $D$ with sparse exponential $e^D$ and a dense matrix $\varepsilon B$ which is of small norm in comparison with $D$. The predominant algorithm is based on scaling the large matrix $A$ by a small number $2^{-s}$, which is then exponentiated by efficient Padé or Taylor methods and finally squared in order to obtain an approximation for the full exponential. In this setting, the main portion of the computational cost arises from dense-matrix multiplications and we present a modified squaring which takes advantage of the smallness of the perturbed matrix $B$ in order to reduce the number of squarings necessary. Theoretical results on local error and error propagation for splitting methods are complemented with numerical experiments and show a clear improvement over existing methods when medium precision is sought.
The Schrödinger eigenvalue problem is solved with the imaginary time propagation technique. The separability of the Hamiltonian makes the problem suitable for the application of splitting methods. High order fractional time steps of order greater than two necessarily have negative steps and can not be used for this class of diffusive problems. However, there exist methods which use fractional complex time steps with positive real parts which can be used with only a moderate increase in the computational cost. We analyze the performance of this class of schemes and propose new methods which outperform the existing ones in most cases. On the other hand, if the gradient of the potential is available, methods up to fourth order with real and positive coefficients exist. We also explore this case and propose new methods as well as sixth-order methods with complex coefficients. In particular, highly optimized sixth-order schemes for near integrable systems using positive real part complex coefficients with and without modified potentials are presented. A time-stepping variable order algorithm is proposed and numerical results show the enhanced efficiency of the new methods.
We consider spherically trapped Bose gases in three dimensions with contact interactions, and investigate whether the Bose-Einstein condensate at zero temperature is stable against macroscopic fragmentation into a small number of mutually incoherent pieces. Our results are expressed in terms of a dimensionless interaction measure proportional to the Thomas-Fermi parameter. It is shown that while three-dimensional condensates are inherently much more stable against macroscopic fragmentation than their quasi-one- and quasi-two-dimensional counterparts, they fragment at a sufficiently large value of the dimensionless interaction measure, which we determine both fully numerically and semianalytically from a continuum limit of large particle numbers.
We present structure preserving integrators for solving linear quadratic optimal control problems. This problem requires the numerical integration of matrix Riccati differential equations whose exact solution is a symmetric positive definite time-dependent matrix which controls the stability of the equation for the state. This property is not preserved, in general, by the numerical methods. We propose second order exponential methods based on the Magnus series expansion which unconditionally preserve positivity for this problem and analyze higher order Magnus integrators. This method can also be used for the integration of nonlinear problems if they are previously linearized. The performance of the algorithms is illustrated with the stabilization of a quadrotor which is an unmanned aerial vehicle.
We consider the numerical integration of the Gross-Pitaevskii equation with a potential trap given by a time-dependent harmonic potential or a small perturbation thereof. Splitting methods are frequently used with Fourier techniques since the system can be split into the kinetic and remaining part, and each part can be solved efficiently using Fast Fourier Transforms. To split the system into the quantum harmonic oscillator problem and the remaining part allows to get higher accuracies in many cases, but it requires to change between Hermite basis functions and the coordinate space, and this is not efficient for time-dependent frequencies or strong nonlinearities. We show how to build new methods which combine the advantages of using Fourier methods while solving the timedependent harmonic oscillator exactly (or with a high accuracy by using a Magnus integrator and an appropriate decomposition).
We investigate the level population statistics and degree of coherence encoded in the single-particle density matrix of harmonically trapped low-dimensional [quasi-one-dimensional (quasi-1D) or quasi-two-dimensional (quasi-2D)] Bose gases with repulsive contact interactions. Using a variational analysis, we derive fragmentation of the condensate in the weakly confining directions into two (quasi-1D) respectively three (quasi-2D) mutually incoherent macroscopic pieces, upon increasing a dimensionless interaction measure beyond a critical value. Fragmented condensate many-body states in low-dimensional systems therefore occur well before the thermodynamic limit of infinite extension is reached, in which phase fluctuations of the matter wave field create an infinite number of nonmacroscopic fragments.