Publications
Monographs

M. Liero, M. Thomas, D. Peschka, eds., Special Issue on Energybased Mathematical Methods for Reactive Multiphase Flows, 103 of Zeitschrift für Angewandte Mathematik und Mechanik (ZAMM), WileyVCH Verlag, Weinheim, 2023, (Collection Published), DOI 10.1002/zamm.202302012 .
Articles in Refereed Journals

P. Bella, M. Kniely, Regularity of random elliptic operators with degenerate coefficients and applications to stochastic homogenization, Stochastic Partial Differential Equations. Analysis and Computations, published online on 27.02.2024, DOI https://doi.org/10.1007/s40072023003229 .
Abstract
We consider degenerate elliptic equations of second order in divergence form with a symmetric random coefficient field a. Extending the work of the first author, Fehrman, and Otto [Ann. Appl. Probab. 28 (2018), no. 3, 13791422], who established the largescale regularity of aharmonic functions in a degenerate situation, we provide stretched exponential moments for the minimal radius describing the minimal scale for this regularity. As an application to stochastic homogenization, we partially generalize results by Gloria, Neukamm, and Otto [Anal. PDE 14 (2021), no. 8, 24972537] on the growth of the corrector, the decay of its gradient, and a quantitative twoscale expansion to the degenerate setting. On a technical level, we demand the ensemble of coefficient fields to be stationary and subject to a spectral gap inequality, and we impose moment bounds on the coefficient field and its inverse. We also introduce the ellipticity radius, which encodes the minimal scale where these moments are close to their positive expectation value. 
M. Heida, Stochastic homogenization on perforated domains I: Extension operators, Networks and Heterogeneous Media, 18 (2023), pp. 1821/11821/78, DOI 10.3934/nhm.2023079 .
Abstract
This preprint is part of a major rewriting and substantial improvement of WIAS Preprint 2742. In this first part of a series of 3 papers, we set up a framework to study the existence of uniformly bounded extension and trace operators for W^{1,p}functions on randomly perforated domains, where the geometry is assumed to be stationary ergodic. We drop the classical assumption of minimaly smoothness and study stationary geometries which have no global John regularity. For such geometries, uniform extension operators can be defined only from W^{1,p} to W^{1,r} with the strict inequality r<p. In particular, we estimate the L^{r}norm of the extended gradient in terms of the L^{p}norm of the original gradient. Similar relations hold for the symmetric gradients (for ℝ^{d}valued functions) and for traces on the boundary. As a byproduct we obtain some Poincaré and Korn inequalities of the same spirit. Such extension and trace operators are important for compactness in stochastic homogenization. In contrast to former approaches and results, we use very weak assumptions: local (δ,M)regularity to quantify statistically the local Lipschitz regularity and isotropic cone mixing to quantify the density of the geometry and the mesoscopic properties. These two properties are sufficient to reduce the problem of extension operators to the connectivity of the geometry. In contrast to former approaches we do not require a minimal distance between the inclusions and we allow for globally unbounded Lipschitz constants and percolating holes. We will illustrate our method by applying it to the Boolean model based on a Poisson point process and to a Delaunay pipe process, for which we can explicitly estimate the connectivity terms. 
M. Heida, Stochastic homogenization on perforated domains III  General estimates for stationary ergodic random connected Lipschitz domains, Networks and Heterogeneous Media, 18 (2023), pp. 14101433, DOI 10.3934/nhm.2023062 .
Abstract
This is Part III of a series on the existence of uniformly bounded extension operators on randomly perforated domains in the context of homogenization theory. Recalling that randomly perforated domains are typically not John and hence extension is possible only from W ^{1,p} to W ^{1,r}, r < p, we will show that the existence of such extension operators can be guarantied if the weighted expectations of four geometric characterizing parameters are bounded: The local Lipschitz constant M, the local Lipschitz radius Δ , the mesoscopic Voronoi diameter ∂ and the local connectivity radius R. 
M. Mirahmadi, B. Friedrich, B. Schmidt, J. PérezRíos, Mapping atomic trapping in an optical superlattice onto the libration of a planar rotor in electric fields, New Journal of Physics, 25 (2023), pp. 023024/1023024/16, DOI 10.1088/13672630/acbab6 .
Abstract
We show that two seemingly unrelated problems  the trapping of an atom in a onedimensional optical superlattice (OSL) formed by the interference of optical lattices whose spatial periods differ by a factor of two, and the libration of a polar polarizable planar rotor (PR) in combined electric and optical fields  have isomorphic Hamiltonians. Since the OSL gives rise to a periodic potential that acts on atomic translation via the AC Stark effect, it is possible to establish a map between the translations of atoms in this system and the rotations of the PR due to the coupling of the rotor's permanent and induced electric dipole moments with the external fields. The latter system belongs to the class of conditionally quasiexactly solvable (CQES) problems in quantum mechanics and shows intriguing spectral properties, such as avoided and genuine crossings, studied in details in our previous works [our works]. We make use of both the spectral characteristics and the quasiexact solvability to treat ultracold atoms in an optical superlattice as a semifinitegap system. The band structure of this system follows from the eigenenergies and their genuine and avoided crossings obtained as solutions of the WhittakerHill equation. Furthermore, the mapping makes it possible to establish correspondence between concepts developed for the two eigenproblems individually, such as localization on the one hand and orientation/alignment on the other. This correspondence may pave the way to unraveling the dynamics of the OSL system in analytic form. 
P. Colli, G. Gilardi, A. Signori, J. Sprekels, CahnHilliardBrinkman model for tumor growth with possibly singular potentials, Nonlinearity, 36 (2023), pp. 44704500, DOI https://doi.org/10.1088/13616544/ace2a7 .
Abstract
We analyze a phase field model for tumor growth consisting of a CahnHilliardBrinkman system, ruling the evolution of the tumor mass, coupled with an advectionreactiondiffusion equation for a chemical species acting as a nutrient. The main novelty of the paper concerns the discussion of the existence of weak solutions to the system covering all the meaningful cases for the nonlinear potentials; in particular, the typical choices given by the regular, the logarithmic, and the double obstacle potentials are admitted in our treatise. Compared to previous results related to similar models, we suggest, instead of the classical noflux condition, a Dirichlet boundary condition for the chemical potential appearing in the CahnHilliardtype equation. Besides, abstract growth conditions for the source terms that may depend on the solution variables are postulated. 
P. Colli, G. Gilardi, A. Signori, J. Sprekels, Optimal control of a nonconserved phase field model of Caginalp type with thermal memory and double obstacle potential, Discrete and Continuous Dynamical Systems  Series S, 16 (2023), pp. 23052325, DOI 10.3934/dcdss.2022210 .
Abstract
In this paper, we investigate optimal control problems for a nonlinear state system which constitutes a version of the Caginalp phase field system modeling nonisothermal phase transitions with a nonconserved order parameter that takes thermal memory into account. The state system, which is a firstorder approximation of a thermodynamically consistent system, is inspired by the theories developed by Green and Naghdi. It consists of two nonlinearly coupled partial differential equations that govern the phase dynamics and the universal balance law for internal energy, written in terms of the phase variable and the socalled thermal displacement, i.e., a primitive with respect to time of temperature. We extend recent results obtained for optimal control problems in which the free energy governing the phase transition was differentiable (i.e., of regular or logarithmic type) to the nonsmooth case of a double obstacle potential. As is well known, in this nondifferentiable case standard methods to establish the existence of appropriate Lagrange multipliers fail. This difficulty is overcome utilizing of the socalled deep quench approach. Namely, the double obstacle potential is approximated by a family of (differentiable) logarithmic ones for which the existence of optimal controls and firstorder necessary conditions of optimality in terms of the adjoint state variables and a variational inequality are known. By proving appropriate bounds for the adjoint states of the approximating systems, we can pass to the limit in the corresponding firstorder necessary conditions, thereby establishing meaningful firstorder necessary optimality conditions also for the case of the double obstacle potential. 
P. Colli, G. Gilardi, A. Signori, J. Sprekels, On a CahnHilliard system with source term and thermal memory, Nonlinear Analysis. An International Mathematical Journal, 240 (2024), pp. 113461/1113461/16 (published online on 14.12.2023), DOI 10.1016/j.na.2023.113461 .
Abstract
A nonisothermal phase field system of CahnHilliard type is introduced and analyzed mathematically. The system constitutes an extension of the classical Caginalp model for nonisothermal phase transitions with a conserved order parameter. It couples a CahnHilliard type equation with source term for the order parameter with the universal balance law of internal energy. In place of the standard Fourier form, the constitutive law of the heat flux is assumed in the form given by the theory developed by Green and Naghdi, which accounts for a possible thermal memory of the evolution. This has the consequence that the balance law of internal energy becomes a secondorder in time equation for the thermal displacement or freezing index, that is, a primitive with respect to time of the temperature. Another particular feature of our system is the presence of the source term in the equation for the order parameter, which entails additional mathematical difficulties because the mass conservation of the order parameter is lost. We provide several mathematical results under general assumptions on the source term and the doublewell nonlinearity governing the evolution: existence and continuous dependence results are shown for weak and strong solutions to the corresponding initialboundary value problem. 
P. Colli, G. Gilardi, A. Signori, J. Sprekels, Optimal temperature distribution for a nonisothermal CahnHilliard system in two dimensions with source term and double obstacle potential, Annals of the Academy of Romanian Scientists. Mathematics and its Applications., 15 (2023), pp. 175204, DOI 10.56082/annalsarscimath.2023.12.175 .
Abstract
In this note, we study the optimal control of a nonisothermal phase field system of CahnHilliard type that constitutes an extension of the classical Caginalp model for nonisothermal phase transitions with a conserved order parameter. It couples a CahnHilliard type equation with source term for the order parameter with the universal balance law of internal energy. In place of the standard Fourier form, the constitutive law of the heat flux is assumed in the form given by the theory developed by Green and Naghdi, which accounts for a possible thermal memory of the evolution. This has the consequence that the balance law of internal energy becomes a secondorder in time equation for the thermal displacement or freezing index, that is, a primitive with respect to time of the temperature. Another particular feature of our system is the presence of the source term in the equation for the order parameter, which entails further mathematical difficulties because the mass conservation of the order parameter is no longer satisfied. In this paper, we study the case that the doublewell potential driving the evolution of the phase transition is given by the nondifferentiable double obstacle potential, thereby complementing recent results obtained for the differentiable cases of regular and logarithmic potentials. Besides existence results, we derive firstorder necessary optimality conditions for the control problem. The analysis is carried out by employing the socalled deep quench approximation in which the nondifferentiable double obstacle potential is approximated by a family of potentials of logarithmic structure for which meaningful firstorder necessary optimality conditions in terms of suitable adjoint systems and variational inequalities are available. Since the results for the logarithmic potentials crucially depend on the validity of the socalled strict separation property which is only available in the spatially twodimensional situation, our whole analysis is restricted to the twodimensional case. 
P. Colli, G. Gilardi, A. Signori, J. Sprekels, Optimal temperature distribution for a nonisothermal CahnHilliard system with source term, Applied Mathematics and Optimization. An International Journal with Applications to Stochastics, 88 (2023), pp. 68/168/31, DOI 10.1007/s00245023100399 .

P.É. Druet, K. Hopf, A. Jüngel, Hyperbolicparabolic normal form and local classical solutions for crossdiffusion systems with incomplete diffusion, Communications in Partial Differential Equations, 48 (2023), pp. 863894, DOI 10.1080/03605302.2023.2212479 .
Abstract
We investigate degenerate crossdiffusion equations with a rankdeficient diffusion matrix that are considered to model populations which move as to avoid spatial crowding and have recently been found to arise in a meanfield limit of interacting stochastic particle systems. To date, their analysis in multiple space dimensions has been confined to the purely convective case with equal mobility coefficients. In this article, we introduce a normal form for an entropic class of such equations which reveals their structure of a symmetric hyperbolicparabolic system. Due to the statedependence of the range and kernel of the singular diffusive matrix, our way of rewriting the equations is different from that classically used for symmetric secondorder systems with a nullspace invariance property. By means of this change of variables, we solve the Cauchy problem for short times and positive initial data in H^s(mathbbT^d) for s>d/2+1. 
K. Fellner, J. Fischer, M. Kniely, B.Q. Tang, Global renormalised solutions and equilibration of reactiondiffusion systems with nonlinear diffusion, Journal of Nonlinear Science, 33 (2023), pp. 66/166/49, DOI 10.1007/s0033202309926w .
Abstract
The global existence of renormalised solutions and convergence to equilibrium for reactiondiffusion systems with nonlinear diffusion are investigated. The system is assumed to have quasipositive nonlinearities and to satisfy an entropy inequality. The difficulties in establishing global renormalised solutions caused by possibly degenerate diffusion are overcome by introducing a new class of weighted truncation functions. By means of the obtained global renormalised solutions, we study the largetime behaviour of complex balanced systems arising from chemical reaction network theory with nonlinear diffusion. When the reaction network does not admit boundary equilibria, the complex balanced equilibrium is shown, by using the entropy method, to exponentially attract all renormalised solutions in the same compatibility class. This convergence extends even to a range of nonlinear diffusion, where global existence is an open problem, yet we are able to show that solutions to approximate systems converge exponentially to equilibrium uniformly in the regularisation parameter. 
R. Finn, M. O'Donovan, P. Farrell, J. Moatti, T. Streckenbach, Th. Koprucki, S. Schulz, Theoretical study of the impact of alloy disorder on carrier transport and recombination processes in deep UV (Al,Ga)N light emitters, Applied Physics Letters, 122 (2023), pp. 241104/1241104/7, DOI 10.1063/5.0148168 .
Abstract
Aluminum gallium nitride [(Al,Ga)N] has gained significant attention in recent years due to its potential for highly efficient light emitters operating in the deep ultraviolet (UV) range (<280 nm). However, given that current devices exhibit extremely low efficiencies, understand ing the fundamental properties of (Al,Ga)Nbased systems is of key importance. Here, using a multiscale simulation framework, we study the impact of alloy disorder on carrier transport, radiative and nonradiative recombination processes in a cplane Al 0.7 Ga0.3 N/Al0.8 Ga0.2 N quantum well embedded in a pn junction. Our calculations reveal that alloy fluctuations can open "percolative" pathways that promote transport for the electrons and holes into the quantum well region. Such an effect is neglected in conventional and widely used transport sim ulations. Moreover, we find that the resulting increased carrier density and alloy induced carrier localization effects significantly increase nonradiative AugerMeitner recombination in comparison to the radiative process. Thus, to suppress such nonradiative process and poten tially related material degradation, a careful design (wider well, multiquantum wells) of the active region is required to improve the effi ciency of deep UV light emitters. 
G. Gilardi, A. Signori, J. Sprekels, Nutrient control for a viscous CahnHilliardKellerSegel model with logistic source describing tumor growth, Discrete and Continuous Dynamical Systems  Series S, 16 (2023), pp. 35523572, DOI 10.3934/dcdss.2023123 .
Abstract
In this paper, we address a distributed control problem for a system of partial differential equations describing the evolution of a tumor that takes the biological mechanism of chemotaxis into account. The system describing the evolution is obtained as a nontrivial combination of a CahnHilliard type system accounting for the segregation between tumor cells and healthy cells, with a KellerSegel type equation accounting for the evolution of a nutrient species and modeling the chemotaxis phenomenon. First, we develop a robust mathematical background that allows us to analyze an associated optimal control problem. This analysis forced us to select a source term of logistic type in the nutrient equation and to restrict the analysis to the case of two space dimensions. Then, the existence of an optimal control and firstorder necessary conditions for optimality are established. 
J. Riedel, P. Gelss, R. Klein, B. Schmidt, WaveTrain: A Python package for numerical quantum mechanics of chainlike systems based on tensor trains, Journal of Chemical Physics, 158 (2023), pp. 164801/1164801/15, DOI 10.1063/5.0147314 .
Abstract
WaveTrain is an opensource software for numerical simulations of chainlike quantum systems with nearestneighbor (NN) interactions only. The Python package is centered around tensor train (TT, or matrix product) format representations of Hamiltonian operators and (stationary or timeevolving) state vectors. It builds on the Python tensor train toolbox Scikit_tt, which provides efficient construction methods and storage schemes for the TT format. Its solvers for eigenvalue problems and linear differential equations are used in WaveTrain for the timeindependent and timedependent Schrödinger equations, respectively. Employing efficient decompositions to construct lowrank representations, the tensortrain ranks of state vectors are often found to depend only marginally on the chain length N. This results in the computational effort growing only slightly more than linearly with N, thus mitigating the curse of dimensionality. As a complement to the classes for full quantum mechanics, WaveTrain also contains classes for fully classical and mixed quantumclassical (Ehrenfest or mean field) dynamics of bipartite systems. The graphical capabilities allow visualization of quantum dynamics on the fly, with a choice of several different representations based on reduced density matrices. Even though developed for treating quasi onedimensional excitonic energy transport in molecular solids or conjugated organic polymers, including coupling to phonons, WaveTrain can be used for any kind of chainlike quantum systems, with or without periodic boundary conditions, and with NN interactions only. The present work describes version 1.0 of our WaveTrain software, based on version 1.2 of scikit_tt, both of which are freely available from the GitHub platform where they will also be further developed. Moreover, WaveTrain is mirrored at SourceForge, within the framework of the WavePacket project for numerical quantum dynamics. Workedout demonstration examples with complete input and output, including animated graphics, are available. 
J. Sprekels, F. Tröltzsch, Secondorder sufficient conditions for sparse optimal control of singular AllenCahn systems with dynamic boundary conditions, Discrete and Continuous Dynamical Systems  Series S, 16 (2023), pp. 37843812, DOI 10.3934/dcdss.2023163 .
Abstract
In this paper we study the optimal control of a parabolic initialboundary value problem of AllenCahn type with dynamic boundary conditions. Phase field systems of this type govern the evolution of coupled diffusive phase transition processes with nonconserved order parameters that occur in a container and on its surface, respectively. It is assumed that the nonlinear functions driving the physical processes within the bulk and on the surface are double well potentials of logarithmic type whose derivatives become singular at the boundary of their respective domains of definition. For such systems, optimal control problems have been studied in the past. We focus here on the situation when the cost functional of the optimal control problem contains a nondifferentiable term like the L^{1}norm leading to sparsity of optimal controls. For such cases, we derive secondorder sufficient conditions for locally optimal controls. 
A. Glitzky, M. Liero, A driftdiffusion based electrothermal model for organic thinfilm devices including electrical and thermal environment, ZAMM. Zeitschrift für Angewandte Mathematik und Mechanik, published online on 27.11.2023, DOI 10.1002/zamm.202300376 .
Abstract
We derive and investigate a stationary model for the electrothermal behavior of organic thinfilm devices including their electrical and thermal environment. Whereas the electrodes are modeled by Ohm's law, the electronics of the organic device itself is described by a generalized van Roosbroeck system with temperature dependent mobilities and using GaussFermi integrals for the statistical relation. The currents give rise to Joule heat which together with the heat generated by the generation/recombination of electrons and holes in the organic device occur as source terms in the heat flow equation that has to be considered on the whole domain. The crucial task is to establish that the quantities in the transfer conditions at the interfaces between electrodes and the organic semiconductor device have sufficient regularity. Therefore, we restrict the analytical treatment of the system to two spatial dimensions. We consider layered organic structures, where the physical parameters (total densities of transport states, LUMO and HOMO energies, disorder parameter, basic mobilities, activation energies, relative dielectric permittivity, heat conductivity) are piecewise constant. We prove the existence of weak solutions using Schauder's fixed point theorem and a regularity result for strongly coupled systems with nonsmooth data and mixed boundary conditions that is verified by Caccioppoli estimates and a Gehringtype lemma. 
TH. Eiter, M. Kyed, Y. Shibata, Periodic Lp estimates by Rboundedness: Applications to the NavierStokes equations, Acta Applicandae Mathematicae. An International Survey Journal on Applying Mathematics and Mathematical Applications, 188 (2023), pp. 1/11/43, DOI 10.1007/s10440023006123 .
Abstract
General evolution equations in Banach spaces are investigated. Based on an operatorvalued version of de Leeuw's transference principle, timeperiodic Lp estimates of maximal regularity type are established from Rbounds of the family of solution operators (Rsolvers) to the corresponding resolvent problems. With this method, existence of timeperiodic solutions to the NavierStokes equations is shown for two configurations: in a periodically moving bounded domain and in an exterior domain, subject to prescribed timeperiodic forcing and boundary data. 
TH. Eiter, M. Kyed, Y. Shibata, Falling drop in an unbounded liquid reservoir: Steadystate solutions, Journal of Mathematical Fluid Mechanics, 25 (2023), pp. 34/134/34, DOI 10.1007/s00021023007779 .
Abstract
The equations governing the motion of a threedimensional liquid drop moving freely in an unbounded liquid reservoir under the influence of a gravitational force are investigated. Provided the (constant) densities in the two liquids are sufficiently close, existence of a steadystate solution is shown. The proof is based on a suitable linearization of the equations. A setting of function spaces is introduced in which the corresponding linear operator acts as a homeomorphism. 
P. Farrell, J. Moatti, M. O'Donovan, S. Schulz, Th. Koprucki, Importance of satisfying thermodynamic consistency in optoelectronic device simulations for high carrier densities, Optical and Quantum Electronics, 55 (2023), pp. 978/1978/12, DOI 10.1007/s11082023052345 .
Abstract
We show the importance of using a thermodynamically consistent flux discretization when describing driftdiffusion processes within light emitting diode simulations. Using the classical ScharfetterGummel scheme with FermiDirac statistics is an example of such an inconsistent scheme. In this case, for an (In,Ga)N multi quantum well device, the Fermi levels show steep gradients on one side of the quantum wells which are not to be expected. This result originates from neglecting diffusion enhancement associated with FermiDirac statistics in the numerical flux approximation. For a thermodynamically consistent scheme, such as the SEDAN scheme, the spikes in the Fermi levels disappear. We will show that thermodynamic inconsistency has far reaching implications on the currentvoltage curves and recombination rates. 
M. Liero, A. Mielke, G. Savaré, Fine properties of geodesics and geodesic $lambda$convexity for the HellingerKantorovich distance, Archive for Rational Mechanics and Analysis, 247 (2023), pp. 112/1112/73, DOI 10.1007/s00205023019411 .
Abstract
We study the fine regularity properties of optimal potentials for the dual formulation of the HellingerKantorovich problem (HK), providing sufficient conditions for the solvability of the primal Monge formulation. We also establish new regularity properties for the solution of the HamiltonJacobi equation arising in the dual dynamic formulation of HK, which are sufficiently strong to construct a characteristic transportdilation flow driving the geodesic interpolation between two arbitrary positive measures. These results are applied to study relevant geometric properties of HK geodesics and to derive the convex behaviour of their Lebesgue density along the transport flow. Finally, exact conditions for functionals defined on the space of measures are derived that guarantee the geodesic lambdaconvexity with respect to the HellingerKantorovich distance. 
A. Mielke, R. Rossi, Balancedviscosity solutions to infinitedimensional multirate systems, Archive for Rational Mechanics and Analysis, 247 (2023), pp. 53/153/100, DOI 10.1007/s0020502301855y .
Abstract
We consider generalized gradient systems with rateindependent and ratedependent dissipation potentials. We provide a general framework for performing a vanishingviscosity limit leading to the notion of parametrized and true BalancedViscosity solutions that include a precise description of the jump behavior developing in this limit. Distinguishing an elastic variable $u$ having a viscous damping with relaxation time $eps^alpha$ and an internal variable $z$ with relaxation time $eps$ we obtain different limits for the three cases $alpha in (0,1)$, $alpha=1$ and $alpha>1$. An application to a delamination problem shows that the theory is general enough to treat nontrivial models in continuum mechanics. 
A. Mielke, Nonequilibrium steady states as saddle points and EDPconvergence for slowfast gradient systems, Journal of Mathematical Physics, 64 (2023), pp. 123502/1 123502/27, DOI 10.1063/5.0149910 .
Abstract
The theory of slowfast gradient systems leads in a natural way to nonequilibrium steady states, because on the slow time scale the fast subsystem stays in steady states that are driven by the interaction with the slow system. Using the theory of convergence of gradient systems in the sense of the energydissipation principle shows that there is a natural characterization of these nonequilibrium steady states as saddle points of a Lagrangian where the slow variables are fixed. We give applications to slowfast reactiondiffusion systems based on the socalled coshtype gradient structure for reactions. It is shown that two binary reaction give rise to a ternary reaction with a statedependent reaction coefficient. Moreover, we show that a reactiondiffusion equation with a thin membranelike layer convergences to a transmission condition, where the formerly quadratic dissipation potential for diffusion convergences to a coshtype dissipation potential for the transmission in the membrane limit. 
A. Mielke, On two coupled degenerate parabolic equations motivated by thermodynamics, Journal of Nonlinear Science, 33 (2023), pp. 42/142/55, DOI 10.1007/s00332023098923 .
Abstract
We discuss a system of two coupled parabolic equations that have degenerate diffusion constants depending on the energylike variable. The dissipation of the velocitylike variable is fed as a source term into the energy equation leading to conservation of the total energy. The motivation of studying this system comes from Prandtl's and Kolmogorov's one and twoequation models for turbulence, where the energylike variable is the mean turbulent kinetic energy. Because of the degeneracies there are solutions with timedependent support like in the porous medium equation, which is contained in our system as a special case. The motion of the free boundary may be driven by either selfdiffusion of the energylike variable or by dissipation of the velocitylike variable. The crossover of these two phenomena is exemplified for the associated planar traveling fronts. We provide existence of suitably defined weak and very weak solutions. After providing a thermodynamically motivated gradient structure we also establish convergence into steady state for bounded domains and provide a conjecture on the asymptotically selfsimilar behavior of the solutions in R^{d} for large times.
Contributions to Collected Editions

R. Danabalan, M. Hintermüller, Th. Koprucki, K. Tabelow, MaRDI: Building research data infrastructures for mathematics and the mathematical sciences, in: 1st Conference on Research Data Infrastructure (CoRDI)  Connecting Communities, Y. SureVetter, C. Goble, eds., 1 of Proceedings of the Conference on Research Data Infrastructure, TIB Open Publishing, Hannover, 2023, pp. 69/169/4, DOI 10.52825/cordi.v1i.397 .
Abstract
MaRDI is building a research data infrastructure for mathematics and beyond based on semantic technologies (metadata, ontologies, knowledge graphs) and data repositories. Focusing on the algorithms, models and workflows, the MaRDI infrastructure will connect with other disciplines and NFDI consortia on data processing methods, solving real world problems and support mathematicians on research datamanagement 
A. Maltsi, A. Mielke, Th. Koprucki, Symmetries in Transmission Electron Microscopy images of semiconductor nanostructures with strain, in: 23nd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2023), P. Bardella, A. Tibaldi, eds., IEEE, 2023, pp. 111112, DOI 10.1109/NUSOD59562.2023.10273568 .
Abstract
Transmission electron microscopy is often used to image semiconductor nanostructures with strain. The resulting images exhibit symmetries, the source of which is not always known. We prove mathematically that the intensities are invariant under specific transformations, which allows us to distinguish between symmetries of the imaging process itself and symmetries of the inclusion. 
L. Ermoneit, B. Schmidt, J. Fuhrmann, Th. Koprucki, L.R. Schreiber, M. Kantner, Simulation of singleelectron shuttling for spinqubit transport in a SiGe quantum bus, in: Book of Abstracts of the International Workshop on Computational Nanotechnology 2023 (IWCN 2023), X. Orios Plaedvall, G. Abadal Berini, X. Cartoixà Soler, A. Cummings, C.F. Destefani, D. Jiménez Jiménez, J. Mart'in Mart'inez, R. Rodr'iguez Mart'inez, A. Benali, eds., pp. 8889.

M. Heida, Finite volumes for simulation of large molecules, in: Finite Volumes for Complex Applications X  Volume 1, Elliptic and Parabolic Problems: FVCA10, Strasbourg, France, October 30, 2023  November 03, 2023, Invited Contributions, E. Franck, J. Fuhrmann, V. MichelDansac, L. Navoret, eds., 432 of Springer Proceedings in Mathematics & Statistics, Springer International Publishing, Cham, 2023, pp. 305313, DOI 10.1007/9783031408649_25 .
Abstract
We study a finite volume scheme for simulating the evolution of large molecules within their reduced state space. The finite volume scheme under consideration is the SQRA scheme developed by Lie, Weber and Fackeldey. We study convergence of a more general family of FV schemes in up to 3 dimensions and provide a convergence result for the SQRAscheme in arbitrary space dimensions. 
M. O'Donovan, R. Finn, S. Schulz, Th. Koprucki, Atomistic study of Urbach tail energies in (Al,Ga)N quantum well systems, in: 23nd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2023), P. Bardella, A. Tibaldi, eds., IEEE, 2023, pp. 7980, DOI 10.1109/NUSOD59562.2023.10273479 .
Abstract
Aluminium gallium nitride is a system of interest for developing ultraviolet (UV) optoelectronic devices. Here Urbach tails induced by carrier localization effects play a key role in determining device behaviour. We study the electronic structure of Al x Ga 1x N/Al y Ga 1y N single quantum wells using an atomistic framework. Results show that the density of states exhibits a tail at low energies due to disorder in the alloy microstructure. Our analysis allows for insight into the orbital character of the states forming the Urbach tails, which can affect light polarization characteristics, and important quantity for deep UV light emitters. 
T. Boege, R. Fritze, Ch. Görgen, J. Hanselman, D. Iglezakis, L. Kastner, Th. Koprucki, T. Krause, Ch. Lehrenfeld, S. Polla, M. Reidelbach, Ch. Riedel, J. Saak, B. Schembera, K. Tabelow, M. Weber, Researchdata management planning in the German mathematical community, 130 of EMS Magazine, European Mathematical Society, 2023, pp. 4047, DOI 10.4171/MAG/152 .

R. Finn, M. O'Donovan, P. Farrell, T. Streckenbach, J. Moatti, Th. Koprucki, S. Schulz, Theoretical investigation of carrier transport and recombination processes for deep UV (Al,Ga)N light emitters, in: 23nd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2023), P. Bardella, A. Tibaldi, eds., IEEE, 2023, pp. 8384, DOI 10.1109/NUSOD59562.2023.10273485 .
Abstract
We present a theoretical study on the impact of alloy disorder on carrier transport and recombination rates in an (Al,Ga)N single quantum well based LED operating in the deep UV spectral range. Our calculations indicate that alloy fluctuations enable percolative pathways which can result in improved carrier injection into the well, but may also increase carrier leakage from the well. Additionally, we find that alloy disorder induces carrier localization effects, a feature particularly noticeable for the holes. These localization effects can lead to locally increased carrier densities when compared to a virtual crystal approximation which neglects alloy disorder. We observe that both radiative and nonradiative recombination rates are increased. Our calculations also indicate that AugerMeitner recombination increases faster than the radiative rate, based on a comparison with a virtual crystal approximation. 
B. Schembera, F. Wübbeling, Th. Koprucki, Ch. Biedinger, M. Reidelbach, B. Schmidt, D. Göddeke, J. Fiedler, Building ontologies and knowledge graphs for mathematics and its applications, in: 1st Conference on Research Data Infrastructure (CoRDI)  Connecting Communities, Y. SureVetter, C. Goble, eds., 1 of Proceedings of the Conference on Research Data Infrastructure, TIB Open Publishing, Hannover, 2023, pp. 29/129/5, DOI 10.52825/cordi.v1i.255 .
Abstract
Ontologies and knowledge graphs for mathematical algorithms and modelsare presented, that have been developed by the Mathematical Research Data Initiative.This enables FAIR data handling in mathematics and the applied disciplines. Moreover, challenges of harmonization during the ontology development are discussed. 
A. Glitzky, M. Liero, A bulksurface model for the electrothermal feedback in largearea organic lightemitting diodes, in: 93rd Annual Meeting 2023 of the International Association of Applied Mathematics and Mechanics (GAMM), 23 of Proc. Appl. Math. Mech. (Special Issue), WileyVCH Verlag, Weinheim, 2023, pp. e202300018/1e202300018/8, DOI 10.1002/pamm.202300018 .
Abstract
This work deals with an effective bulksurface thermistor model describing the electrothermal behavior of largearea thinfilm organic lightemitting diodes (OLEDs). This model was rigorously derived from a Laplace thermistor model by dimension reduction and consists of the heat equation in the threedimensional glass substrate and two semilinear equations describing the current flow through the electrodes coupled to algebraic equations that express the continuity of the electrical fluxes through the organic layers. The electrical problem lives on the surface of the glass substrate where the OLED is mounted. The source terms in the heat equation result from Joule heating and are concentrated on the part of the boundary where the currentflow problem is formulated. Schauder's fixedpoint theorem is used to establish the existence of weak solutions to this effective system. Since the heat source terms at the surface are a priori only in L1, the concept of entropy solutions for the heat equation is worked with. 
M. Landstorfer, M. Heida, Energie effizienter speichern, Spektrum der Wissenschaft, Spektrum der Wissenschaft Verlagsgesellschaft mbH, Heidelberg, 2023, pp. 7279.
Preprints, Reports, Technical Reports

TH. Eiter, A.L. Silvestre, Representation formulas and farfield behavior of timeperiodic flow past a body, Preprint no. 3091, WIAS, Berlin, 2024, DOI 10.20347/WIAS.PREPRINT.3091 .
Abstract, PDF (325 kByte)
This paper is concerned with integral representations and asymptotic expansions of solutions to the timeperiodic incompressible NavierStokes equations for fluid flow in the exterior of a rigid body that moves with constant velocity. Using the timeperiodic Oseen fundamental solution, we derive representation formulas for solutions with suitable regularity. From these formulas, the decomposition of the velocity component of the fundamental solution into steadystate and purely periodic parts and their detailed decay rate in space, we deduce complete information on the asymptotic structure of the velocity and pressure fields. 
P. Colli, G. Gilardi, A. Signori, J. Sprekels, Curvature effects in pattern formation: Wellposedness and optimal control of a sixthorder CahnHilliard equation, Preprint no. 3085, WIAS, Berlin, 2024, DOI 10.20347/WIAS.PREPRINT.3085 .
Abstract, PDF (398 kByte)
This work investigates the wellposedness and optimal control of a sixthorder CahnHilliard equation, a higherorder variant of the celebrated and wellestablished CahnHilliard equation. The equation is endowed with a source term, where the control variable enters as a distributed mass regulator. The inclusion of additional spatial derivatives in the sixthorder formulation enables the model to capture curvature effects, leading to a more accurate depiction of isothermal phase separation dynamics in complex materials systems. We provide a wellposedness result for the aforementioned system when the corresponding nonlinearity of doublewell shape is regular and then analyze a corresponding optimal control problem. For the latter, existence of optimal controls is established, and the firstorder necessary optimality conditions are characterized via a suitable variational inequality. These results aim at contributing to improve the understanding of the mathematical properties and control aspects of the sixthorder CahnHilliard equation, offering potential applications in the design and optimization of materials with tailored microstructures and properties. 
J. Li, X. Liu, D. Peschka, Local wellposedness and global stability of onedimensional shallow water equations with surface tension and constant contact angle, Preprint no. 3084, WIAS, Berlin, 2024, DOI 10.20347/WIAS.PREPRINT.3084 .
Abstract, PDF (405 kByte)
We consider the onedimensional shallow water problem with capillary surfaces and moving contact lines. An energybased model is derived from the twodimensional water wave equations, where we explicitly discuss the case of a stationary force balance at a moving contact line and highlight necessary changes to consider dynamic contact angles. The moving contact line becomes our free boundary at the level of shallow water equations, and the depth of the shallow water degenerates near the free boundary, which causes singularities for the derivatives and degeneracy for the viscosity. This is similar to the physical vacuum of compressible flows in the literature. The equilibrium, the global stability of the equilibrium, and the local wellposedness theory are established in this paper. 
L. Ermoneit, B. Schmidt, Th. Koprucki, J. Fuhrmann, T. Breiten, A. Sala, N. Ciroth, R. Xue, L.R. Schreiber, M. Kantner, Optimal control of conveyormode spinqubit shuttling in a Si/SiGe quantum bus in the presence of charged defects, Preprint no. 3082, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3082 .
PDF (9473 kByte) 
A. Erhardt, D. Peschka, Ch. Dazzi, L. Schmeller, A. Petersen, S. Checa, A. Münch, B. Wagner, Modeling cellular selforganization in strainstiffening hydrogels, Preprint no. 3076, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3076 .
Abstract, PDF (14 MByte)
We develop a threedimensional mathematical model framework for the collective evolution of cell populations by an agentbased model (ABM) that mechanically interacts with the surrounding extracellular matrix (ECM) modeled as a hydrogel. We derive effective twodimensional models for the geometrical setup of a thin hydrogel sheet to study cellcell and cellhydrogel mechanical interactions for a range of external conditions and intrinsic material properties. We show that without any stretching of the hydrogel sheets, cells show the wellknown tendency to form long chains with varying orientations. Our results further show that external stretching of the sheet produces the expected nonlinear strainsoftening or stiffening response, with, however, little qualitative variation of the overall cell dynamics for all the materials considered. The behavior is remarkably different when solvent is entering or leaving from strain softening or stiffening hydrogels, respectively. 
D. Abdel, A. Glitzky, M. Liero, Analysis of a driftdiffusion model for perovskite solar cells, Preprint no. 3073, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3073 .
Abstract, PDF (541 kByte)
This paper deals with the analysis of an instationary driftdiffusion model for perovskite solar cells including FermiDirac statistics for electrons and holes and Blakemore statistics for the mobile ionic vacancies in the perovskite layer. The free energy functional is related to this choice of the statistical relations. Exemplary simulations varying the mobility of the ionic vacancy demonstrate the necessity to include the migration of ionic vacancies in the model frame. To prove the existence of weak solutions, first a problem with regularized state equations and reaction terms on any arbitrarily chosen finite time interval is considered. Its solvability follows from a time discretization argument and passage to the timecontinuous limit. Applying Moser iteration techniques, a priori estimates for densities, chemical potentials and the electrostatic potential of its solutions are derived that are independent of the regularization level, which in turn ensure the existence of solutions to the original problem. 
L. Araujo, C. Lasser, B. Schmidt, FSSH2: Fewest Switches Surface Hopping with robust switching probability, Preprint no. 3071, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3071 .
Abstract, PDF (199 kByte)
This study introduces the FSSH2 scheme, a redefined and numerically stable adiabatic Fewest Switches Surface Hopping (FSSH) method. It reformulates the standard FSSH hopping probability without nonadiabatic coupling vectors and allows for numerical time integration with larger step sizes. The advantages of FSSH2 are demonstrated by numerical experiments for five different model systems in one and two spatial dimensions with up to three electronic states. 
Y. Hadjimichael, Ch. Merdon, M. Liero, P. Farrell, An energybased finitestrain model for 3D heterostructured materials and its validation by curvature analysis, Preprint no. 3064, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3064 .
Abstract, PDF (6517 kByte)
This paper presents a comprehensive study of the intrinsic strain response of 3D het erostructures arising from lattice mismatch. Combining materials with different lattice constants induces strain, leading to the bending of these heterostructures. We propose a model for nonlinear elastic heterostructures such as bimetallic beams or nanowires that takes into account local prestrain within each distinct material region. The resulting system of partial differential equations (PDEs) in Lagrangian coordinates incorporates a nonlinear strain and a linear stressstrain relationship governed by Hooke?s law. To validate our model, we apply it to bimetallic beams and hexagonal heteronanowires and perform numerical simulations using finite element methods (FEM). Our simulations ex amine how these structures undergo bending under varying material compositions and crosssectional geometries. In order to assess the fidelity of the model and the accuracy of simulations, we compare the calculated curvature with analytically derived formula tions. We derive these analytical expressions through an energybased approach as well as a kinetic framework, adeptly accounting for the lattice constant mismatch present at each compound material of the heterostructures. The outcomes of our study yield valuable insights into the behavior of strained bent heterostructures. This is particularly significant as the strain has the potential to influence the electronic band structure, piezoelectricity, and the dynamics of charge carriers. 
M. Heida, On the computation of high dimensional Voronoi diagrams, Preprint no. 3041, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3041 .
Abstract, PDF (553 kByte)
We investigate a recently implemented new algorithm for the computation of a Voronoi diagram in high dimensions and generalize it to N nodes in general or nongeneral position using a geometric characterization of edges merging in a given vertex. We provide a mathematical proof that the algorithm is exact, convergent and has computational costs of O(E*nn(N)), where E is the number of edges and nn(N) is the computational cost to calculate the nearest neighbor among N points. We also provide data from performance tests in the recently developed Julia package ,,HighVoronoi.jl”. 
K.M. Gambaryan, O. Ernst, T. Boeck, O. Marquardt, Energy level alignment of confined hole states in $InAs_1xySb_xP_y$ double quantum dots, Preprint no. 3040, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3040 .
Abstract, PDF (310 kByte)
We present a combined experimental and theoretical study of uncapped In(As,Sb,P) double quantum dots (DQD), suited for application in novel resonant tunneling nanodiods or singlephoton nanooptical up and downconverters in the midinfrared spectral range. We provide details on the growth process using liquidphase epitaxy (LPE), as well as on the characterization using atomicforce microscopy (AFM) and scanning electron microscopy (SEM). We find that most DQDs exhibit an asymmetry such that the two QDs of each pair have different dimensions, giving rise to correspondingly different quantum confinement of hole states localized in each QD. Based on these data, we have performed systematic simulations based on an eightband $mathbfkcdotmathbfp$ model to identify the relationship between QD dimensions and the energy difference between corresponding confined hole states in the two QDs. Finally, we have determined the strength of an applied electric field required to energetically align the hole ground states of two QDs of different dimensions in order to facilitate hole tunneling. 
A. Mielke, R. Rossi, A. Stephan, On timesplitting methods for gradient flows with two dissipation mechanisms, Preprint no. 3033, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3033 .
Abstract, PDF (530 kByte)
We consider generalized gradient systems in Banach spaces whose evolutions are generated by the interplay between an energy functional and a dissipation potential. We focus on the case in which the dual dissipation potential is given by a sum of two functionals and show that solutions of the associated gradientflow evolution equation with combined dissipation can be constructed by a splitstep method, i.e. by solving alternately the gradient systems featuring only one of the dissipation potentials and concatenating the corresponding trajectories. Thereby the construction of solutions is provided either by semiflows, on the timecontinuous level, or by using Alternating Minimizing Movements in the timediscrete setting. In both cases the convergence analysis relies on the energydissipation principle for gradient systems. 
A. Stephan, Trottertype formula for operator semigroups on product spaces, Preprint no. 3030, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3030 .
Abstract, PDF (252 kByte)
We consider a Trottertypeproduct formula for approximating the solution of a linear abstract Cauchy problem (given by a strongly continuous semigroup), where the underlying Banach space is a product of two spaces. In contrast to the classical Trotterproduct formula, the approximation is given by freezing subsequently the components of each subspace. After deriving necessary stability estimates for the approximation, which immediately provide convergence in the natural strong topology, we investigate convergence in the operator norm. The main result shows that an almost optimal convergence rate can be established if the dominant operator generates a holomorphic semigroup and the offdiagonal coupling operators are bounded. 
W. van Oosterhout, M. Liero, Finitestrain poroviscoelasticity with degenerate mobility, Preprint no. 3027, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3027 .
Abstract, PDF (360 kByte)
A quasistatic nonlinear model for poroviscoelastic solids at finite strains is considered in the Lagrangian frame using the concept of secondorder nonsimple materials. The elastic stresses satisfy static frameindifference, while the viscous stresses satisfy dynamic frameindifference. The mechanical equation is coupled to a diffusion equation for a solvent or fluid content. The latter is pulledback to the reference configuration. To treat the nonlinear dependence of the mobility tensor on the deformation gradient, the result by Healey & Krömer is used to show that the determinant of the deformation gradient is bounded away from zero. Moreover, the focus is on the physically relevant case of degenerate mobilities. The existence of weak solutions is shown using a staggered timeincremental scheme and suitable energydissipation inequalities. 
A. Mielke, An introduction to the analysis of gradients systems, Preprint no. 3022, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3022 .
Abstract, PDF (862 kByte)
The present notes provide an extended version of a small lecture course (of 36 hours) given at the HumboldtUniversität zu Berlin in the Winter Term 2022/23. The material starting in Section 5.4 was added afterwards. The aim of these notes to give an introductory overview on the analytical approaches for gradientflow equations in Hilbert spaces, Banach spaces, and metric spaces and to show that on the first entry level these theories have a lot in common. The theories and their specific setups are illustrated by suitable examples and counterexamples. 
J. Sprekels, F. Tröltzsch, Secondorder sufficient conditions in the sparse optimal control of a phase field tumor growth model with logarithmic potential, Preprint no. 3020, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3020 .
Abstract, PDF (344 kByte)
his paper treats a distributed optimal control problem for a tumor growth model of viscous CahnHilliard type. The evolution of the tumor fraction is governed by a thermodynamic force induced by a doublewell potential of logarithmic type. The cost functional contains a nondifferentiable term in order to enhance the occurrence of sparsity effects in the optimal controls, i.e., of subdomains of the spacetime cylinder where the controls vanish. In the context of cancer therapies, sparsity is very important in order that the patient is not exposed to unnecessary intensive medical treatment. In this work, we focus on the derivation of secondorder sufficient optimality conditions for the optimal control problem. While in previous works on the system under investigation such conditions have been established for the case without sparsity, the case with sparsity has not been treated before. The results obtained in this paper also improve the known results on this phase field model for the case without sparsity. 
Y. Bokredenghel, M. Heida, Quenched homogenization of infinite range random conductance model on stationary point processes, Preprint no. 3017, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3017 .
Abstract, PDF (391 kByte)
We prove homogenization for elliptic longrange operators in the random conductance model on random stationary point processes in d dimensions with Dirichlet boundary conditions and with a jointly stationary coefficient field. Doing so, we identify 4 conditions on the point process and the coefficient field that have to be fulfilled at different stages of the proof in order to pass to the homogenization limit. The conditions can be clearly attributed to concentration of support, RellichPoincaré inequality, nondegeneracy of the homogenized matrix and ergodicity of the elliptic operator. 
A. Mielke, T. Roubíček, U. Stefanelli, A model of gravitational differentiation of compressible selfgravitating planets, Preprint no. 3015, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3015 .
Abstract, PDF (444 kByte)
We present a dynamic model for inhomogeneous viscoelastic media at finite strains. The model features a KelvinVoigt rheology, and includes a selfgenerated gravitational field in the actual evolving configuration. In particular, a fully Eulerian approach is adopted. We specialize the model to viscoelastic (barotropic) fluids and prove existence and a certain regularity of global weak solutions by a FaedoGalerkin semidiscretization technique. Then, an extension to multicomponent chemically reacting viscoelastic fluids based on a phenomenological approach by Eckart and Prigogine, is advanced and studied. The model is inspired by planetary geophysics. In particular, it describes gravitational differentiation of inhomogeneous planets and moons, possibly undergoing volumetric phase transitions. 
A. Mielke, S. Schindler, Convergence to selfsimilar profiles in reactiondiffusion systems, Preprint no. 3008, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3008 .
Abstract, PDF (380 kByte)
We study a reactiondiffusion system on the real line, where the reactions of the species are given by one reversible reaction pair satisfying the massaction law. We describe different positive limits at both sides of infinityand investigate the longtime behavior. Rescaling space and time according to the parabolic scaling, we show that solutions converge exponentially to a similarity profile when the scaled time goes to infinity. In the original variables, these profiles correspond to asymptotically selfsimilar behavior describing the phenomenon of diffusive mixing of the different states at infinity.Our method provides global exponential convergence for all initial states with finite relative entropy. For the case with equal stoichiometric coefficients, we can allow for selfsimilar profiles with arbitrary equilibrated states,while in the other case we need to assume that the two states atinfinity are sufficiently close such that the selfsimilar profile is relative flat. 
A. Mielke, S. Schindler, Existence of similarity profiles for diffusion equations and systems, Preprint no. 3007, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3007 .
Abstract, PDF (403 kByte)
We study the existence of selfsimilar profiles for diffusion equations and reactiondiffusion systems on the real line, where different nontrivial limits are imposed at both sides of infinity. The theses profiles solve a coupled system of nonlinear ODEs that can be treated by monotone operator theory. 
K. Hopf, A. Jüngel, Convergence of a finite volume scheme and dissipative measurevaluedstrong stability for a hyperbolicparabolic crossdiffusion system, Preprint no. 3006, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3006 .
Abstract, PDF (444 kByte)
This article is concerned with the approximation of hyperbolicparabolic crossdiffusion systems modeling segregation phenomena for populations by a fully discrete finitevolume scheme. It is proved that the numerical scheme converges to a dissipative measurevalued solution of the PDE system and that, whenever the latter possesses a strong solution, the convergence holds in the strong sense. Furthermore, the “parabolic density part” of the limiting measurevalued solution is atomic and converges to its constant state for long times. The results are based on Young measure theory and a weakstrong stability estimate combining Shannon and Rao entropies. The convergence of the numerical scheme is achieved by means of discrete entropy dissipation inequalities and an artificial diffusion, which vanishes in the continuum limit. 
TH. Eiter, Y. Shibata, Viscous flow past a translating body with oscillating boundary, Preprint no. 3000, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.3000 .
Abstract, PDF (326 kByte)
We study an incompressible viscous flow around an obstacle with an oscillating boundary that moves by a translational periodic motion, and we show existence of strong timeperiodic solutions for small data in different configurations: If the mean velocity of the body is zero, existence of timeperiodic solutions is provided within a framework of Sobolev functions with isotropic pointwise decay. If the mean velocity is nonzero, this framework can be adapted, but the spatial behavior of flow requires a setting of anisotropically weighted spaces. In the latter case, we also establish existence of solutions within an alternative framework of homogeneous Sobolev spaces. These results are based on the timeperiodic maximal regularity of the associated linearizations, which is derived from suitable Rbounds for the Stokes and Oseen resolvent problems. The pointwise estimates are deduced from the associated timeperiodic fundamental solutions. 
P. Gelss, S. Matera, R. Klein, B. Schmidt, Quantum dynamics of coupled excitons and phonons in chainlike systems: Tensor train approaches and higherorder propagators, Preprint no. 2995, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.2995 .
Abstract, PDF (546 kByte)
We investigate the use of tensortrain approaches to the solution of the timedependent Schrödinger equation for chainlike quantum systems with onsite and nearestneighbor interactions only. Using the efficient SLIM representation for lowrank tensor train representations of quantummechanical Hamiltonians, we aim at reducing the memory consumption as well as the computation costs, in order to mitigate the curse of dimensionality as much as possible. As an example, coupled excitons and phonons modeled in terms of FröhlichHolstein type Hamiltonians are studied here. By comparing our tensortrain based results with semianalytical results, we demonstrate the key role of the ranks of tensortrain representations for quantum state vectors. Both the computational effort of the propagations and the accuracy that can be reached crucially depend on the maximum number of ranks chosen. Typically, an excellent quality of the solutions is found only when the ranks exceeds a certain value. That threshold, however, is very different for excitons, phonons, and coupled systems. One class of propagation schemes used in the present work builds on splitting the Hamiltonian into two groups of interleaved nearestneighbor interactions which commutate within each of the groups. In addition to the first order LieTrotter and the second order StrangMarchuk splitting schemes, we have also implemented the 4th order YoshidaNeri and the 8th order KahanLi symplectic compositions. Especially the latter two are demonstrated to yield very accurate results, close to machine precision. However, due to the computational costs, currently their use is restricted to rather short chains. Another class of propagators involves explicit, timesymmetrized Euler integrators. Building on the traditional second order differencing method, we have also implemented higher order methods. Especially the 4th order variant is recommended for quantum simulations of longer chains, even though the high precision of the splitting schemes cannot be reached. Moreover, the scaling of the computational effort with the dimensions of the local Hilbert spaces is much more favorable for the differencing than for the splitting schemes. 
A. Mielke, S. Schindler, Selfsimilar pattern in coupled parabolic systems as nonequilibrium steady states, Preprint no. 2992, WIAS, Berlin, 2023, DOI 10.20347/WIAS.PREPRINT.2992 .
Abstract, PDF (370 kByte)
We consider reactiondiffusion systems and other related dissipative systems on unbounded domains which would have a Liapunov function (and gradient structure) when posed on a finite domain. In this situation, the system may reach local equilibrium on a rather fast time scale but the infinite amount of mass or energy leads to persistent mass or energy flow for all times. In suitably rescaled variables the system converges to a steady state that corresponds to asymptotically selfsimilar behavior in the original system.
Talks, Poster

A. Mielke, On EVI flows for gradient systems on the (spherical) HellingerKantorovich space, Workshop ``Applications of Optimal Transportation'', February 5  9, 2024, Mathematisches Forschungsinstitut Oberwolfach, February 5, 2024.

A. Maltsi, Symmetries in Transmission Electron Microscopy images of semiconductor nanostructures with strain, 23rd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2023), September 18  21, 2023, Politecnico di Torino, Italy, September 21, 2023.

A. Maltsi, Symmetries in TEM imaging of semiconductor nanostructures with strain, 15th Annual Meeting Photonic Devices, March 29  31, 2023, ZuseInstitut Berlin, March 31, 2023.

A. Maltsi, Symmetries in TEM imaging of semiconductor nanostructures with strain, Leibniz MMS Days 2023, April 17  19, 2023, LeibnizInstitut für Agrartechnik und Bioökonomie (ATB), Potsdam, April 18, 2023.

A. Maltsi, Symmetries in transmission electron microscopy imaging of crystals with strain, 10th International Congress on Industrial and Applied Mathematics (ICIAM 2023), Tokyo, Japan, August 20  25, 2023.

L. Ermoneit, B. Schmidt, J. Fuhrmann, Th. Koprucki, M. Kantner, Coherent spinqubit shuttling in a SiGe quantum bus: Devicescale modeling, simulation and optimal control, Leibniz MMS Days 2023, Potsdam, April 17  19, 2023.

L. Ermoneit, M. Kantner, Th. Koprucki, B. Schmidt, Coherent spinqubit shuttling for scalable quantum processors: Modeling, simulation and optimal control, MATH+ Day, HumboldtUniversität zu Berlin, October 20, 2023.

M. Heida, Finite volumes for simulation of large molecules, Finite Volumes for Complex Applications (FVCA10), Strasbourg, France, October 30  November 3, 2023.

M. Heida, Diskrete Operatoren in Modellbildung und Numerik, Universität der Bundeswehr München, Institut für Mathematik und Informatik, July 13, 2023.

M. Heida, Perspectives for homogenization on randomly perforated domains, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), SPP 2256 ``Variational Methods for Predicting Complex Phenomena in Engineering Structures and Materials'', May 30  June 2, 2023, Technische Universität Dresden, June 1, 2023.

M. O'Donovan, Atomistic study of Urbach tail energies in (Al,Ga)N quantum well systems, 23rd International Conference on Numerical Simulation of Optoelectronic Devices (NUSOD 2023), September 18  21, 2023, Politecnico di Torino, Italy, September 21, 2023.

M. O'Donovan, Impact of alloy disorder on carrier transport and recombination in (Al,Ga)Nbased UVC emitters, The 6th International Workshop on Ultraviolet Materials and Devices (IWUMD 2023), June 5  8, 2023, Metz Congrès Robert Schuman, France, June 7, 2023.

M. O'Donovan, Modeling random alloy fluctuations in carrier transport simulations of IIIN based light emitting diodes  Connecting atomistic tightbinding to driftdiffusion, 15th Annual Meeting Photonic Devices, March 29  31, 2023, ZuseInstitut Berlin, March 31, 2023.

M. O'Donovan, Tight binding simulations of localization in alloy fluctuations in nitride based LEDs, Seminar zu Physik der Gruppe IIINitridHalbleiter und nanophotonischer Bauelemente und Advanced IIINitride Materials and Photonic Devices (IIINMPD), Technische Universität Berlin, AG Experimentelle Nanophysik und Photonik, May 17, 2023.

L. Schütz, M. Heida, M. Thomas, Materials with discontinuities on many scales, SCCS Days 2023 of the Collaborative Research Center  CRC 1114 ``Scaling Cascades in Complex Systems'', November 13  15, 2023.

L. Schütz, An existence theory for solitary waves on a ferrofluid jet, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), Session 14 ``Applied Analysis'', May 30  June 2, 2023, Technische Universität Dresden, May 30, 2023.

L. Schütz, Towards stochastic homogenization of a rateindependent delamination model, Hausdorff School ``Analysis of PDEs: Variational and Geometric Perspectives'', Bonn, July 10  14, 2023.

M. Kniely, A thermodynamically correct framework for electroenergyreactiondiffusion systems, 22nd ECMI Conference on Industrial and Applied Mathematics, June 26  30, 2023, Wrocław University of Science and Technology, Poland, June 30, 2023.

M. Kniely, On a thermodynamically consistent electroenergyreactiondiffusion system, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), Session 14 ``Applied Analysis'', May 30  June 2, 2023, Technische Universität Dresden, June 1, 2023.

J. Rehberg, Maximal parabolic regularity for the treatment of real world problems, Oberseminar für Optimale Steuerung und Inverse Probleme, Universität DuisburgEssen, Fakultät für Mathematik, May 4, 2023.

J. Rehberg, A view on the KohnSham system from the perspective of functional analysists, Technische Universität Braunschweig Institut für Analysis und Algebra, November 13, 2023.

J. Rehberg, Nonsmooth elliptic and parabolic regularity, Technische Universität Darmstadt, Fachbereich Mathematik, November 27, 2023.

J. Sprekels, Sparse optimal control of singular AllenCahn systems with dynamic boundary conditions, Kolloquium, Università di Pavia, Dipartimento di Matematica ``F. Casorati'', Italy, April 18, 2023.

A. Glitzky, An effective bulksurface thermistor model for largearea organic lightemitting diodes, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), Session 14 ``Applied Analysis'', May 30  June 2, 2023, Technische Universität Dresden, May 30, 2023.

K. Hopf, Normal form and the Cauchy problem for crossdiffusive mixtures, Workshop ``Variational Methods for Evolution'', December 3  8, 2023, Mathematisches Forschungsinstitut Oberwolfach, December 4, 2023.

K. Hopf, On crossdiffusive coupling of hyperbolicparabolic type, Variational and Geometric Structures for Evolution, October 9  13, 2023, Centro Internazionale per la Ricerca Matematica (CIRM), Levico Terme, Italy, October 13, 2023.

K. Hopf, Structure and approximation of crossdiffusive mixtures with incomplete diffusion, Universität Kassel, Fachbereich Mathematik und Naturwissenschaften, September 21, 2023.

K. Hopf, Structure, dynamics, and approximation of crossdiffusive mixtures with incomplete diffusion, Universität Hamburg, Fachbereich Mathematik, May 10, 2023.

K. Hopf, The Cauchy problem for multicomponent systems with strong crossdiffusion, Johannes GutenbergUniversität Mainz, Fachbereich Physik, Mathematik und Informatik, January 11, 2023.

M. Thomas, Damage in viscoelastic materials at finite strains, Workshop ``Variational Methods for Evolution'', December 3  8, 2023, Mathematisches Forschungsinstitut Oberwolfach, December 7, 2023.

M. Eigel, M. Heida, M. Landstorfer, A. Selahi, Recovery of battery ageing dynamics with multiple timescales, MATH+ Day, HumboldtUniversität zu Berlin, October 20, 2023.

TH. Eiter, R. Lasarzik, Analysis of energyvariational solutions for hyperbolic conservation laws, Presentation of project proposals in SPP 2410 ``Hyperbolic Balance Laws in Fluid Mechanics: Complexity, Scales, Randomness'', Bad Honnef, April 28, 2023.

TH. Eiter, Artificial boundary conditions for timeperiodic flow past a body, 10th International Congress on Industrial and Applied Mathematics (ICIAM 2023), Minisymposium 00558 ``Bifurcations, Periodicity and Stability in Fluidstructure Interactions'', August 20  25, 2023, Waseda University, Tokyo, Japan, August 21, 2023.

TH. Eiter, Energyvariational solutions for a class of hyperbolic conservation laws, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), Session 14 ``Applied Analysis'', May 30  June 2, 2023, Technische Universität Dresden, June 2, 2023.

TH. Eiter, The concept of energyvariational solutions for hyperbolic conservation laws, Seminar on Partial Differential Equations, Czech Academy of Sciences, Institute of Mathematics, Prague, Czech Republic, March 28, 2023.

M. Kantner, L. Ermoneit, B. Schmidt, J. Fuhrmann, A. Sala, L.R. Schreiber, Th. Koprucki, Optimal control of a SiGequantum bus for coherent electron shuttling in scalable quantum computing architectures, Silicon Quantum Electronics Workshop 2023, Kyoto, Japan, October 31  November 2, 2023.

TH. Koprucki, Building research data services for the community, Leibniz MMS Days 2023, April 17  19, 2023, LeibnizInstitut für Agrartechnik und Bioökonomie (ATB), Potsdam, April 17, 2023.

TH. Koprucki, MaRDI  The Mathematical Research Data Initiative within the German National Research Data Infrastructure (NFDI), SIAM Conference on Computational Science and Engineering (CSE23), Minisymposium MS301 ``Interfaces, Workflows, and Knowledge Graphs for FAIR CSE'', February 27  March 3, 2023, Amsterdam, Netherlands, March 2, 2023.

TH. Koprucki, MaRDI  The Mathematical Research Data Initiative within the German National Research Data Infrastructure (NFDI), Kolloquium der AG ``Modellierung, Numerik, Differentialgleichungen'', Technische Universität Berlin, January 31, 2023.

M. Liero, Variational modeling of biomechanical systems, 10th International Conference on Computational Bioengineering (ICCB 2023), Minisymposium 223 ``Continuum Biomechanics of Active Biological Systems'', September 20  22, 2023, Technische Universität Wien, Austria, September 22, 2023.

M. Liero, Analysis for thermomechanical models with internal variables, Presentation of project proposals in SPP 2256 ``Variational Methods for Predicting Complex Phenomena in Engineering Structures and Materials'', Bad Honnef, March 27, 2023.

M. Liero, Balancedviscosity solutions for a PenroseFife phasefield model with friction, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), SPP 2256 ``Variational Methods for Predicting Complex Phenomena in Engineering Structures and Materials'', May 30  June 2, 2023, Technische Universität Dresden, June 1, 2023.

M. Liero, EDPconvergence for evolutionary systems with gradient flow structure, 29th Nordic Congress of Mathematicians with EMS, July 3  7, 2023, Aalborg University, Department of Mathematical Sciences, Denmark, July 4, 2023.

M. Liero, Luminance inhomogeneities in largearea OLEDs due to electrothermal feedback, Hybride Optoelektronische Materialsysteme (HYD Seminar), Integrative Research Institute for the Sciences (IRIS Adlershof), Hybrid Devices Group, Berlin, April 20, 2023.

M. Liero, On the geometry of the HellingerKantorovich space (hybrid talk), Seminar ``Modern Methods in Applied Stochastics and Nonparametric Statistics'', WIAS Berlin, January 31, 2023.

A. Mielke, Variational and geometric structures for thermomechanical systems, Variational and Geometric Structures for Evolution, October 9  13, 2023, Centro Internazionale per la Ricerca Matematica (CIRM), Levico Terme, Italy, October 11, 2023.

A. Mielke, Asymptotic selfsimilar behavior in reactiondiffusion systems on the real line, Minisymposium ``Interacting Particle Systems and Variational Methods'', Einhoven University of Technology, Department of Mathematics and Computer Science, Netherlands, February 3, 2023.

A. Mielke, Balancedviscosity solutions as limits in generalized gradient systems under slow loading, Hausdorff School ``Analysis of PDEs: Variational and Geometric Perspectives'', July 10  14, 2023, Universität Bonn, Hausdorff School for Advanced Studies in Mathematics.

A. Mielke, EDPconvergence for gradient systems and nonequilibrium steady states, Nonlinear Diffusion and Nonlocal Interaction Models  Entropies, Complexity, and MultiScale Structures, May 28  June 2, 2023, Universidad de Granada, Spain, May 30, 2023.

A. Mielke, Nonequilibrium steady states and EDPconvergence for slowfast gradient systems, In Search of Model Structures for Nonequilibrium Systems, April 24  28, 2023, Westfälische WilhelmsUniversität Münster, April 25, 2023.

A. Mielke, Viscoelastic fluid models for geodynamic processes in the lithosphere, ``SPP Meets TP'' Workshop: Variational Methods for Complex Phenomena in Solids, February 21  24, 2023, Universität Bonn, Hausdorff Institute for Mathematics, February 24, 2023.

A. Stephan, Gradient systems and timesplitting methods (online talk), PDE & Applied Mathematics Seminar, University of California, Riverside, Department of Mathematics, USA, November 8, 2023.

A. Stephan, On timesplitting methods for gradient flows with two dissipation mechanisms, Gradient Flows facetoface 3, September 11  14, 2023, Université Claude Bernard Lyon 1, France, September 11, 2023.

A. Stephan, On timesplitting methods for gradient flows with two dissipation mechanisms, PDE Afternoon, Technische Universität Wien, Austria, December 13, 2023.

A. Stephan, Positivity and polynomial decay of energies for squarefield operators, Variational and Geometric Structures for Evolution, October 9  13, 2023, Centro Internazionale per la Ricerca Matematica (CIRM), Levico Terme, Italy, October 13, 2023.

A. Stephan, Fastslow chemical reaction systems: Gradient systems and EDPconvergence, Oberseminar Dynamics, Technische Universität München, Department of Mathematics, April 17, 2023.

A. Stephan, On timesplitting methods for gradient flows with two dissipation mechanisms, In Search of Model Structures for Nonequilibrium Systems, April 24  28, 2023, Westfälische WilhelmsUniversität Münster, April 28, 2023.

A. Stephan, On timesplitting methods for gradient flows with two dissipation mechanisms, 10th International Congress on Industrial and Applied Mathematics (ICIAM 2023), Minisymposium 01181 ``Variational Methods for Multiscale Dynamics'', August 20  25, 2023, Waseda University, Tokyo, Japan, August 24, 2023.

W. van Oosterhout, Analysis of poroviscoelastic solids at finite strains, 93rd Annual Meeting of the International Association of Applied Mathematics and Mechanics (GAMM 2023), Session 14 ``Applied Analysis'', May 30  June 2, 2023, Technische Universität Dresden, June 2, 2023.

W. van Oosterhout, Poroviscoelastic solids at finite strains with degenerate mobilities, Nonlinear PDEs: Recent Trends in the Analysis of Continuum Mechanics, July 17  21, 2023, Universität Bonn, Hausdorff School for Advanced Studies in Mathematics, July 19, 2023.
External Preprints

R. Finn, M. O'Donovan, P. Farrell, J. Moatti, T. Streckenbach, Th. Koprucki, S. Schulz, Theoretical study of the impact of alloy disorder on carrier transport and recombination processes in deep UV (Al, Ga)N light emitters, Preprint no. hal04037215, Hyper Articles en Ligne (HAL), 2023.
Abstract
Aluminium gallium nitride ((Al,Ga)N) has gained significant attention in recent years due to its potential for highly efficient light emitters operating in the deep ultraviolet (UV) range (< 280 nm). However, given that current devices exhibit extremely low efficiencies, understanding the fundamental properties of (Al,Ga)Nbased systems is of key importance. Here, using a multiscale simulation framework, we study the impact of alloy disorder on carrier transport, radiative and nonradiative recombination processes in a cplane Al0.7Ga0.3N/Al0.8Ga0.2N quantum well embedded in a pin junction. Our calculations reveal that alloy fluctuations can open "percolative" pathways that promote transport for the electrons and holes into the quantum well region. Such an effect is neglected in conventional, and widely used transport simulations. Moreover, we find also that the resulting increased carrier density and alloy induced carrier localization effects significantly increase nonradiative AugerMeitner recombination in comparison to the radiative process. Thus, to avoid such nonradiative process and potentially related material degradation, a careful design (wider well, multi quantum wells) of the active region is required to improve the efficiency of deep UV light emitters. 
B. Schembera, F. Wübbeling, H. Kleikamp, Ch. Bledinger, J. Fiedler, M. Reidelbach, A. Shehu, B. Schmidt, Th. Koprucki, D. Iglezakis, D. Göddeke, Ontologies for models and algorithms in applied mathematics and related disciplines, Preprint no. arXiv:2310.20443, Cornell University, 2023, DOI 10.48550/arXiv.2310.20443 .
Abstract
In applied mathematics and related disciplines, the modelingsimulationoptimization workflow is a prominent scheme, with mathematical models and numerical algorithms playing a crucial role. For these types of mathematical research data, the Mathematical Research Data Initiative has developed, merged and implemented ontologies and knowledge graphs. This contributes to making mathematical research data FAIR by introducing semantic technology and documenting the mathematical foundations accordingly. Using the concrete example of microfracture analysis of porous media, it is shown how the knowledge of the underlying mathematical model and the corresponding numerical algorithms for its solution can be represented by the ontologies. 
P. Farrell, J. Moatti, M. O'Donovan, S. Schulz, Th. Koprucki, Importance of satisfying thermodynamic consistency in light emitting diode simulations, Preprint no. hal04012467, Hyper Articles en Ligne (HAL), 2023.
Abstract
We show the importance of using a thermodynamically consistent flux discretization when describing driftdiffusion processes within light emitting diode simulations. Using the classical ScharfetterGummel scheme with FermiDirac statistics is an example of such an inconsistent scheme. In this case, for an (In,Ga)N multi quantum well device, the Fermi levels show steep gradients on one side of the quantum wells which are not to be expected. This result originates from neglecting diffusion enhancement associated with FermiDirac statistics in the numerical flux approximation. For a thermodynamically consistent scheme, such as the SEDAN scheme, the spikes in the Fermi levels disappear. We will show that thermodynamic inconsistency has far reaching implications on the currentvoltage curves and recombination rates.
Research Groups
 Partial Differential Equations
 Laser Dynamics
 Numerical Mathematics and Scientific Computing
 Nonlinear Optimization and Inverse Problems
 Interacting Random Systems
 Stochastic Algorithms and Nonparametric Statistics
 Thermodynamic Modeling and Analysis of Phase Transitions
 Nonsmooth Variational Problems and Operator Equations