BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Christian Lubich

Christian Lubich contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.OCphysics.plasm-phgr-qcmath-phmath.APmath.DSmath.MP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

On a large-stepsize integrator for charged-particle dynamics

Xiao and Qin [Computer Physics Comm., 265:107981, 2021] recently proposed a remarkably simple modification of the Boris algorithm to compute the guiding centre of the highly oscillatory motion of a charged particle with step sizes that are much larger than the period of gyrorotations. They gave strong numerical evidence but no error analysis. This paper provides an analysis of the large-stepsize modified Boris method in a setting that has a strong non-uniform magnetic field and moderately bounded velocities, considered over a fixed finite time interval. The error analysis is based on comparing the modulated Fourier expansions of the exact and numerical solutions, for which the differential equations of the dominant terms are derived explicitly. Numerical experiments illustrate and complement the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Polarized high-frequency wave propagation beyond the nonlinear Schrödinger approximation

This paper studies highly oscillatory solutions to a class of systems of semilinear hyperbolic equations with a small parameter, in a setting that includes Klein--Gordon equations and the Maxwell--Lorentz system. The interest here is in solutions that are polarized in the sense that up to a small error, the oscillations in the solution depend on only one of the frequencies that satisfy the dispersion relation with a given wave vector appearing in the initial wave packet. The construction and analysis of such polarized solutions is done using modulated Fourier expansions. This approach includes higher harmonics and yields approximations to polarized solutions that are of arbitrary order in the small parameter, going well beyond the known first-order approximation via a nonlinear Schrödinger equation. The given construction of polarized solutions is explicit, uses in addition a linear Schrödinger equation for each further order of approximation, and is accessible to direct numerical approximation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rank-$1$ matrix differential equations for structured eigenvalue optimization

A new approach to solving eigenvalue optimization problems for large structured matrices is proposed and studied. The class of optimization problems considered is related to computing structured pseudospectra and their extremal points, and to structured matrix nearness problems such as computing the structured distance to instability or to singularity. The structure can be a general linear structure and includes, for example, large matrices with a given sparsity pattern, matrices with given range and co-range, and Hamiltonian matrices. Remarkably, the eigenvalue optimization can be performed on the manifold of complex (or real) rank-1 matrices, which yields a significant reduction of storage and in some cases of the computational cost. The method relies on a constrained gradient system and the projection of the gradient onto the tangent space of the manifold of complex rank-$1$ matrices. It is shown that near a local minimizer this projection is very close to the identity map, and so the computationally favorable rank-1 projected system behaves locally like the %computationally expensive gradient system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rank-adaptive time integration of tree tensor networks

A rank-adaptive integrator for the approximate solution of high-order tensor differential equations by tree tensor networks is proposed and analyzed. In a recursion from the leaves to the root, the integrator updates bases and then evolves connection tensors by a Galerkin method in the augmented subspace spanned by the new and old bases. This is followed by rank truncation within a specified error tolerance. The memory requirements are linear in the order of the tensor and linear in the maximal mode dimension. The integrator is robust to small singular values of matricizations of the connection tensors. Up to the rank truncation error, which is controlled by the given error tolerance, the integrator preserves norm and energy for Schrodinger equations, and it dissipates the energy in gradient systems. Numerical experiments with a basic quantum spin system illustrate the behavior of the proposed algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Time-dependent electromagnetic scattering from thin layers

The scattering of electromagnetic waves from obstacles with wave-material interaction in thin layers on the surface is described by generalized impedance boundary conditions, which provide effective approximate models. In particular, this includes a thin coating around a perfect conductor and the skin effect of a highly conducting material. The approach taken in this work is to derive, analyse and discretize a system of time-dependent boundary integral equations that determines the tangential traces of the scattered electric and magnetic fields. In a familiar second step, the fields are evaluated in the exterior domain by a representation formula, which uses the time-dependent potential operators of Maxwell's equations. The time-dependent boundary integral equationis discretized with Runge--Kutta based convolution quadrature in time and Raviart--Thomas boundary elements in space. Using the frequency-explicit bounds from the well-posedness analysis given here together with known approximation properties of the numerical methods, the full discretization is proved to be stable and convergent, with explicitly given rates in the case of sufficient regularity. Taking the same Runge--Kutta based convolution quadrature for discretizing the time-dependent representation formulas, the optimal order of convergence is obtained away from the scattering boundary, whereas an order reduction occurs close to the boundary. The theoretical results are illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Finding the nearest passive or non-passive system via Hamiltonian eigenvalue optimization

We propose and study an algorithm for computing a nearest passive system to a given non-passive linear time-invariant system (with much freedom in the choice of the metric defining `nearest', which may be restricted to structured perturbations), and also a closely related algorithm for computing the structured distance of a given passive system to non-passivity. Both problems are addressed by solving eigenvalue optimization problems for Hamiltonian matrices that are constructed from perturbed system matrices. The proposed algorithms are two-level methods that optimize the Hamiltonian eigenvalue of smallest positive real part over perturbations of a fixed size in the inner iteration, using a constrained gradient flow. They optimize over the perturbation size in the outer iteration, which is shown to converge quadratically in the typical case of a defective coalescence of simple eigenvalues approaching the imaginary axis. For large systems, we propose a variant of the algorithm that takes advantage of the inherent low-rank structure of the problem. Numerical experiments illustrate the behavior of the proposed algorithms.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large-stepsize integrators for charged-particle dynamics over multiple time scales

The Boris algorithm, a closely related variational integrator and a newly proposed filtered variational integrator are studied when they are used to numerically integrate the equations of motion of a charged particle in a non-uniform strong magnetic field, taking step sizes that are much larger than the period of the Larmor rotations. For the Boris algorithm and the standard (unfiltered) variational integrator, satisfactory behaviour is only obtained when the component of the initial velocity orthogonal to the magnetic field is filtered out. The particle motion shows varying behaviour over multiple time scales: fast Larmor rotation, guiding centre motion, slow perpendicular drift, near-conservation of the magnetic moment over very long times and conservation of energy for all times. Using modulated Fourier expansions of the exact and numerical solutions, it is analysed to which extent this behaviour is reproduced by the three numerical integrators used with large step sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A convergent algorithm for forced mean curvature flow driven by diffusion on the surface

The evolution of a closed two-dimensional surface driven by both mean curvature flow and a reaction--diffusion process on the surface is formulated into a system, which couples the velocity law not only to the surface partial differential equation but also to the evolution equations for the geometric quantities, namely the normal vector and the mean curvature on the surface. Two algorithms are considered for the obtained system. Both methods combine surface finite elements as a space discretisation and linearly implicit backward difference formulae for time integration. Based on our recent results for mean curvature flow, one of the algorithms directly admits a convergence proof for its full discretisation in the case of finite elements of polynomial degree at least two and backward difference formulae of orders two to five. Numerical examples are provided to support and complement the theoretical convergence results (demonstrating the convergence properties of the method without error estimate), and demonstrate the effectiveness of the methods in simulating a three-dimensional tumour growth model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A convergent evolving finite element algorithm for Willmore flow of closed surfaces

A proof of convergence is given for a novel evolving surface finite element semi-discretization of Willmore flow of closed two-dimensional surfaces, and also of surface diffusion flow. The numerical method proposed and studied here discretizes fourth-order evolution equations for the normal vector and mean curvature, reformulated as a system of second-order equations, and uses these evolving geometric quantities in the velocity law interpolated to the finite element space. This numerical method admits a convergence analysis in the case of continuous finite elements of polynomial degree at least two. The error analysis combines stability estimates and consistency estimates to yield optimal-order $H^1$-norm error bounds for the computed surface position, velocity, normal vector and mean curvature. The stability analysis is based on the matrix--vector formulation of the finite element method and does not use geometric arguments. The geometry enters only into the consistency estimates. Numerical experiments illustrate and complement the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Higher-order linearly implicit full discretization of the Landau--Lifshitz--Gilbert equation

For the Landau--Lifshitz--Gilbert (LLG) equation of micromagnetics we study linearly implicit backward difference formula (BDF) time discretizations up to order $5$ combined with higher-order non-conforming finite element space discretizations, which are based on the weak formulation due to Alouges but use approximate tangent spaces that are defined by $L^2$-averaged instead of nodal orthogonality constraints. We prove stability and optimal-order error bounds in the situation of a sufficiently regular solution. For the BDF methods of orders $3$ to~$5$, this requires %a mild time step restriction $τ\leqslant ch$ and that the damping parameter in the LLG equations be above a positive threshold; this condition is not needed for the A-stable methods of orders $1$ and $2$, for which furthermore a discrete energy inequality irrespective of solution regularity is proved.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measuring the stability of spectral clustering

As an indicator of the stability of spectral clustering of an undirected weighted graph into $k$ clusters, the $k$th spectral gap of the graph Laplacian is often considered. The $k$th spectral gap is characterized in this paper as an unstructured distance to ambiguity, namely as the minimal distance of the Laplacian to arbitrary symmetric matrices with vanishing $k$th spectral gap. As a conceptually more appropriate measure of stability, the structured distance to ambiguity of the $k$-clustering is introduced as the minimal distance of the Laplacian to Laplacians of graphs with the same vertices and edges but with weights that are perturbed such that the $k$th spectral gap vanishes. To compute a solution to this matrix nearness problem, a two-level iterative algorithm is proposed that uses a constrained gradient system of matrix differential equations in the inner iteration and a one-dimensional optimization of the perturbation size in the outer iteration. The structured and unstructured distances to ambiguity are compared on some example graphs. The numerical experiments show, in particular, that selecting the number $k$ of clusters according to the criterion of maximal stability can lead to different results for the structured and unstructured stability indicators.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time integration of tree tensor networks

Dynamical low-rank approximation by tree tensor networks is studied for the data-sparse approximation to large time-dependent data tensors and unknown solutions of tensor differential equations. A time integration method for tree tensor networks of prescribed tree rank is presented and analyzed. It extends the known projector-splitting integrators for dynamical low-rank approximation by matrices and Tucker tensors and is shown to inherit their favorable properties. The integrator is based on recursively applying the Tucker tensor integrator. In every time step, the integrator climbs up and down the tree: it uses a recursion that passes from the root to the leaves of the tree for the construction of initial value problems on subtree tensor networks using appropriate restrictions and prolongations, and another recursion that passes from the leaves to the root for the update of the factors in the tree tensor network. The integrator reproduces given time-dependent tree tensor networks of the specified tree rank exactly and is robust to the typical presence of small singular values in matricizations of the connection tensors, in contrast to standard integrators applied to the differential equations for the factors in the dynamical low-rank approximation by tree tensor networks.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Time-dependent acoustic scattering from generalized impedance boundary conditions via boundary elements and convolution quadrature

Generalized impedance boundary conditions are effective, approximate boundary conditions that describe scattering of waves in situations where the wave interaction with the material involves multiple scales. In particular, this includes materials with a thin coating (with the thickness of the coating as the small scale) and strongly absorbing materials. For the acoustic scattering from generalized impedance boundary conditions, the approach taken here first determines the Dirichlet and Neumann boundary data from a system of time-dependent boundary integral equations with the usual boundary integral operators, and then the scattered wave is obtained from the Kirchhoff representation. The system of time-dependent boundary integral equations is discretized by boundary elements in space and convolution quadrature in time. The well-posedness of the problem and the stability of the numerical discretization rely on the coercivity of the Calderón operator for the Helmholtz equation with frequencies in a complex half-plane. Convergence of optimal order in the natural norms is proved for the full discretization. Numerical experiments illustrate the behaviour of the proposed numerical method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Symplectic Integration of Post-Newtonian Equations of Motion with Spin

We present a non-canonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In the context of the splitting, it is possible to integrate the individual terms of the spin-orbit and spin-spin Hamiltonians analytically, exploiting the special structure of the underlying equations of motion. The outcome is a symplectic, time-reversible integrator, which can be raised to arbitrary order by composition. A fourth-order version is shown to give excellent behavior concerning error growth and conservation of energy and angular momentum in long-term simulations. Favorable properties of the integrator are retained in the presence of weak dissipative forces due to radiation damping in the full post-Newtonian equations.

0 citations0 reviewsNo code linkedOpen access