Dual S-matrix bootstrap. Part I. 2D theory

Using duality in optimization theory we formulate a dual approach to the S-matrix bootstrap that provides rigorous bounds to 2D QFT observables as a consequence of unitarity, crossing symmetry and analyticity of the scattering matrix. We then explain how to optimize such bounds numerically, and prove that they provide the same bounds obtained from the usual primal formulation of the S-matrix Bootstrap, at least once convergence is attained from both perspectives. These techniques are then applied to the study of a gapped system with two stable particles of different masses, which serves as a toy model for bootstrapping popular physical systems.

Iterative Solvers for EMI Models

This chapter deals with iterative solution algorithms for the four EMI formulations derived in (17, Chapter 10.1007/978-3-030-61157-6_5). Order optimal monolithic solvers robust with respect to material parameters, the number of degrees of freedom of discretization as well as the time-stepping parameter are presented and compared in terms of computational cost. Domain decomposition solver for the single-dimensional primal formulation is discussed.

Hybridization of Mixed High-Order Methods on General Meshes and Application to the Stokes Equations

This paper presents two novel contributions on the recently introduced Mixed High-Order (MHO) methods [`Arbitrary order mixed methods for heterogeneous anisotropic diffusion on general meshes', preprint (2013)]. We first address the hybridization of the MHO method for a scalar diffusion problem and obtain the corresponding primal formulation. Based on the hybridized MHO method, we then design a novel, arbitrary order method for the Stokes problem on general meshes. A full convergence analysis is carried out showing that, when independent polynomials of degree k are used as unknowns (at elements for the pressure and at faces for each velocity component), the energy-norm of the velocity and the L2-norm of the pressure converge with order (k + 1), while the L2-norm of the velocity (super-)converges with order (k + 2). The latter property is not shared by other methods based on a similar choice of unknowns. The theoretical results are numerically validated in two space dimensions on both standard and polygonal meshes.