No abstract

We consider the time-dependent Landau-Lifshitz-Gilbert equation. We prove that each weak solution coincides with the (unique) strong solution, as long as the latter exists in time. Unlike available results in the literature, our analysis also includes the physically relevant lower-order terms like Zeeman contribution, anisotropy, stray field, and the Dzyaloshinskii-Moriya interaction (which accounts for the emergence of magnetic Skyrmions). Moreover, our proof gives a template on how to approach weakstrong uniqueness for even more complicated problems, where LLG is (nonlinearly) coupled to other (nonlinear) PDE systems.

We consider a linear elliptic PDE and a quadratic goal functional. The goal-oriented adaptive FEM algorithm (GOAFEM) solves the primal as well as a dual problem, where the goal functional is always linearized around the discrete primal solution at hand. We show that the marking strategy proposed in [M. Feischl, D. Praetorius and K. G. van der Zee, An abstract analysis of optimal goal-oriented adaptivity, SIAM J. Numer. Anal.54 (2016), 3, 1423–1448] for a linear goal functional is also optimal for quadratic goal functionals, i.e., GOAFEM leads to linear convergence with optimal convergence rates.

No abstract

No abstract

We consider an adaptive finite element method with arbitrary but fixed polynomial degree p ≥ 1, where adaptivity is driven by an edge-based residual error estimator. Based on the modified maximum criterion from [Diening et al, Found. Comput. Math. 16, 2016], we propose a goal-oriented adaptive algorithm and prove that it is instance optimal. Numerical experiments underline our theoretical findings.Date: July 31, 2019. 2010 Mathematics Subject Classification. 65N30, 41A25, 65N12, 65N50. Key words and phrases. Adaptive finite element method, goal-oriented algorithm, quantity of interest, maximum marking strategy, convergence, instance optimality.Acknowledgement: The authors acknowledge support through the Austrian Science Fund (FWF) through the doctoral school Dissipation and dispersion in nonlinear PDEs (grant W1245), the special research program Taming complexity in PDE systems (grant SFB F65), and the stand-alone project Optimal adaptivity for BEM and FEM-BEM coupling (grant P27005).

No abstract

We present a straightforward implementation scheme for solving the time-dependent Schrödinger equation for systems described by the Hubbard Hamiltonian with time-dependent hoppings. The computations can be performed for clusters of up to 14 sites with, in principle, general geometry. For the time evolution, we use the exponential midpoint rule, where the exponentials are computed via a Krylov subspace method, which only uses matrix-vector multiplication. The presented implementation uses standard libraries for constructing sparse matrices and for linear algebra. Therefore, the approach is easy to use on both desktop computers and computational clusters. We apply the method to calculate time evolution of double occupation and nonequilibrium spectral function of a photo-excited Mott-insulator. The results show that not only the double occupation increases due to creation of electron-hole pairs but also the Mott gap becomes partially filled.