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

Patrick Henning

Patrick Henning contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath-phmath.MPmath.SP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
6works
0followers
5topics
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

6 published item(s)

preprint2026arXiv

Nonlinear Inverse Iterations for Spin-Orbit Coupled Quantum Gases

This work concerns the computation of ground states of two-component spin-orbit coupled Bose-Einstein condensates (SO-coupled BECs), modelled by a coupled nonlinear eigenvalue problem of Gross-Pitaevskii type. Spin-orbit coupling gives rise to fascinating phenomena, including supersolid-like phases with spatially modulated densities. However, in such complex settings, conventional numerical approaches, such as generalized inverse iterations or gradient descent, often converge very slowly. To overcome this issue, we apply the concept of the J-method [E.~Jarlebring, S.~Kvaal, W.~Michiels. SIAM~J.~Sci.~Comput.~36-4,~2014] to construct a nonlinear inverse iteration scheme whose convergence can be accelerated through spectral shifting, analogous to techniques used for linear eigenproblems. For a fixed shift parameter, we establish local linear convergence rates determined by spectral gaps in the neighbourhood of each quasi-unique ground state. With adaptively chosen shifts, superlinear convergence is observed, which we verify through numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Shadow Lagrangian dynamics for superfluidity

Motivated by a similar approach for Born-Oppenheimer molecular dynamics, this paper proposes an extended "shadow" Lagrangian density for quantum states of superfluids. The extended Lagrangian contains an additional field variable that is forced to follow the wave function of the quantum state through a rapidly oscillating extended harmonic oscillator. By considering the adiabatic limit for large frequencies of the harmonic oscillator, we can derive the two equations of motions, a Schrödinger-type equation for the quantum state and a wave equation for the extended field variable. The equations are coupled in a nonlinear way, but each equation individually is linear with respect to the variable that it defines. The computational advantage of this new system is that it can be easily discretized using linear time stepping methods, where we propose to use a Crank-Nicolson-type approach for the Schrödinger equation and an extended leapfrog scheme for the wave equation. Furthermore, the difference between the quantum state and the extended field variable defines a consistency error that should go to zero if the frequency tends to infinity. By coupling the time-step size in our discretization to the frequency of the harmonic oscillator we can extract an easily computable consistency error indicator that can be used to estimate the numerical error without any additional costs. The findings are illustrated in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A note on optimal $H^1$-error estimates for Crank-Nicolson approximations to the nonlinear Schrödinger equation

In this paper we consider a mass- and energy--conserving Crank-Nicolson time discretization for a general class of nonlinear Schrödinger equations. This scheme, which enjoys popularity in the physics community due to its conservation properties, was already subject to several analytical and numerical studies. However, a proof of optimal $L^{\infty}(H^1)$-error estimates is still open, both in the semi-discrete Hilbert space setting, as well as in fully-discrete finite element settings. This paper aims at closing this gap in the literature. We also suggest a fixed point iteration to solve the arising nonlinear system of equations that makes the method easy to implement and efficient. This is illustrated by numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Computational homogenization of time-harmonic Maxwell's equations

In this paper we consider a numerical homogenization technique for curl-curl-problems that is based on the framework of the Localized Orthogonal Decomposition and which was proposed in [D. Gallistl, P. Henning, B. Verfürth. SIAM J. Numer. Anal. 56-3:1570-1596, 2018] for problems with essential boundary conditions. The findings of the aforementioned work establish quantitative homogenization results for the time-harmonic Maxwell's equations that hold beyond assumptions of periodicity, however, a practical realization of the approach was left open. In this paper, we transfer the findings from essential boundary conditions to natural boundary conditions and we demonstrate that the approach yields a computable numerical method. We also investigate how boundary values of the source term can effect the computational complexity and accuracy. Our findings will be supported by various numerical experiments, both in $2D$ and $3D$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Quantitative Anderson localization of Schrödinger eigenstates under disorder potentials

This paper concerns spectral properties of linear Schrödinger operators under oscillatory high-amplitude potentials on bounded domains. Depending on the degree of disorder, we prove the existence of spectral gaps amongst the lowermost eigenvalues and the emergence of exponentially localized states. We quantify the rate of decay in terms of geometric parameters that characterize the potential. The proofs are based on the convergence theory of iterative solvers for eigenvalue problems and their optimal local preconditioning by domain decomposition.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem: global convergence and computational efficiency

We propose a new normalized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem based on an energy inner product that depends on time through the density of the flow itself. The gradient flow is well-defined and converges to an eigenfunction. For ground states we can quantify the convergence speed as exponentially fast where the rate depends on spectral gaps of a linearized operator. The forward Euler time discretization of the flow yields a numerical method which generalizes the inverse iteration for the nonlinear eigenvalue problem. For sufficiently small time steps, the method reduces the energy in every step and converges globally in $H^1$ to an eigenfunction. In particular, for any nonnegative starting value, the ground state is obtained. A series of numerical experiments demonstrates the computational efficiency of the method and its competitiveness with established discretizations arising from other gradient flows for this problem.

0 citations0 reviewsNo code linkedOpen access