Dipartimento di Matematica
E23
Via Ferrata 5
27100 Pavia, Italy
Contact
carlo.marcati (at) unipv.it

I am a researcher at the Dipartimento di Matematica of the University of Pavia.

## Research

My research activity focuses on

• the regularity of solutions to elliptic problems with point singularities,
• the numerical approximation of a priori electronic structure calculations in quantum chemistry,
• low-rank tensor representation,
• neural networks for PDEs.

### Papers and preprints

1. Carlo Marcati and Christoph Schwab. Exponential Convergence of Deep Operator Networks for Elliptic Partial Differential Equations. Submitted.bib | abstract | pdf ]
We construct deep operator networks (ONets) between infinite-​dimensional spaces that emulate with an exponential rate of convergence the coefficient-​to-solution map of elliptic second-​order PDEs. In particular, we consider problems set in d-​dimensional periodic domains, d=1,2,…, and with analytic right-​hand sides and coefficients. Our analysis covers diffusion-​reaction problems, parametric diffusion equations, and elliptic systems such as linear isotropic elastostatics in heterogeneous materials. We leverage the exponential convergence of spectral collocation methods for boundary value problems whose solutions are analytic. In the present periodic and analytic setting, this follows from classical elliptic regularity. Within the ONet branch and trunk construction of [Chen and Chen, 1993] and of [Lu et al., 2021], we show the existence of deep ONets which emulate the coefficient-​to-solution map to accuracy ε>0 in the H1 norm, uniformly over the coefficient set. We prove that the neural networks in the ONet have size O(|log(ε)|κ) for some κ>0 depending on the physical space dimension.

@Techreport{MS21_984,
author = {C. Marcati and Ch. Schwab},
title = {Exponential Convergence of Deep Operator Networks for Elliptic Partial Differential Equations},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2021-42},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2021/2021-42.pdf },
year = {2021}
}

2. Markus Faustmann, Carlo Marcati, Jens Markus Melenk, and Christoph Schwab. Weighted analytic regularity for the integral fractional Laplacian in polygons. Submitted.bib | abstract | pdf ]
We prove weighted analytic regularity of solutions to the Dirichlet problem for the integral fractional Laplacian in polygons with analytic right-​hand side. We localize the problem through the Caffarelli-​Silvestre extension and study the tangential differentiability of the extended solutions, followed by bootstrapping based on Caccioppoli inequalities on dyadic decompositions of vertex, edge, and edge-​vertex neighborhoods.

@Techreport{FMMS21_983,
author = {M. Faustmann and C. Marcati and J.M. Melenk and Ch. Schwab},
title = {Weighted analytic regularity for the integral fractional Laplacian in polygons},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2021-41},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2021/2021-41.pdf },
year = {2021}
}
3. Yanchen He, Carlo Marcati, and Christoph Schwab. Analytic regularity for the incompressible Navier-Stokes equations in polygons with mixed boundary conditions. Submitted.bib | abstract | pdf ]
We prove weighted analytic regularity of Leray-Hopf variational solutions for the stationary, incompressible Navier-Stokes Equations (NSE) in plane polygonal domains, subject to analytic body forces. We admit mixed boundary conditions which may change type at each vertex, under the assumption that homogeneous Dirichlet (''no-slip'') boundary conditions are prescribed on at least one side at each vertex of the domain. The weighted analytic regularity results are established in Hilbertian Kondrat'ev spaces with homogeneous corner weights. The proofs rely on a priori estimates for the corresponding linearized boundary value problem in sectors in corner-weighted Sobolev spaces and on an induction argument for the weighted norm estimates on the quadratic nonlinear term in the NSE, in a polar frame.
@Techreport{HMS21_971,
author = {Y. He and C. Marcati and Ch. Schwab},
title = {Analytic regularity for the Navier-Stokes equations in polygons with mixed boundary conditions},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2021-29},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2021/2021-29.pdf },
year = {2021}
}

4. Carlo Marcati, Joost A. A. Opschoor, Philipp C. Petersen, and Christoph Schwab. Exponential ReLU neural network approximation rates for point and edge singularities. To appear in Found. Comp. Math..bib | abstract | pdf ]
We prove exponential expressivity with stable ReLU Neural Networks (ReLU NNs) in H1(Ω) for weighted analytic function classes in certain polytopal domains Ω, in space dimension d=2,3. Functions in these classes are locally analytic on open subdomains D⊂Ω, but may exhibit isolated point singularities in the interior of Ω or corner and edge singularities at the boundary ∂Ω. The exponential expression rate bounds proved here imply uniform exponential expressivity by ReLU NNs of solution families for several elliptic boundary and eigenvalue problems with analytic data. The exponential approximation rates are shown to hold in space dimension d=2 on Lipschitz polygons with straight sides, and in space dimension d=3 on Fichera-type polyhedral domains with plane faces. The constructive proofs indicate in particular that NN depth and size increase poly-logarithmically with respect to the target NN approximation accuracy ε>0 in H1(Ω). The results cover in particular solution sets of linear, second order elliptic PDEs with analytic data and certain nonlinear elliptic eigenvalue problems with analytic nonlinearities and singular, weighted analytic potentials as arise in electron structure models. In the latter case, the functions correspond to electron densities that exhibit isolated point singularities at the positions of the nuclei. Our findings provide in particular mathematical foundation of recently reported, successful uses of deep neural networks in variational electron structure algorithms.

@Techreport{MOPS20_938,
author = {C. Marcati and J. A. A. Opschoor and P. C. Petersen and Ch. Schwab},
title = {Exponential ReLU Neural Network Approximation Rates for Point and Edge Singularities},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2020-65},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2020/2020-65_rev1.pdf },
year = {2020}
}

5. Carlo Marcati, Maxim Rakhuba, and Johan E. M. Ulander. Low rank tensor approximation of singularly perturbed partial differential equations in one dimension. Calcolo 59 (2022).bib | abstract | pdf ]
We derive rank bounds on the quantized tensor train (QTT) compressed approximation of singularly perturbed reaction diffusion partial differential equations (PDEs) in one dimension. Specifically, we show that, independently of the scale of the singular perturbation parameter, a numerical solution with accuracy 0<ϵ<1 can be represented in QTT format with a number of parameters that depends only polylogarithmically on ϵ. In other words, QTT compressed solutions converge exponentially to the exact solution, with respect to a root of the number of parameters. We also verify the rank bound estimates numerically, and overcome known stability issues of the QTT based solution of PDEs by adapting a preconditioning strategy to obtain stable schemes at all scales. We find, therefore, that the QTT based strategy is a rapidly converging algorithm for the solution of singularly perturbed PDEs, which does not require prior knowledge on the scale of the singular perturbation and on the shape of the boundary layers.
@Techreport{MRU20_934,
author = {C. Marcati and M. Rakhuba and J. E. M. Ulander},
title = {Low rank tensor approximation of singularly perturbed partial differential equations in one dimension},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2020-61},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2020/2020-61.pdf },
year = {2020}
}  				
6. Yvon Maday, Carlo Marcati. Weighted analyticity of Hartree-Fock wavefunctions. Submitted.bib | abstract | pdf ]
We prove analytic-type estimates in weighted Sobolev spaces on the eigenfunctions of a class of elliptic and nonlinear eigenvalue problems with singular potentials, which includes the Hartree-Fock equations. Going beyond classical results on the analyticity of the wavefunctions away from the nuclei, we prove weighted estimates locally at each singular point, with precise control of the derivatives of all orders. Our estimates have far-reaching consequences for the approximation of the eigenfunctions of the problems considered, and they can be used to prove a priori estimates on the numerical solution of such eigenvalue problems.
@Techreport{MM20_932,
author = {Y. Maday and C. Marcati},
title = {Weighted analyticity of Hartree-Fock eigenfunctions},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2020-59},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2020/2020-59.pdf },
year = {2020}
}

7. Carlo Marcati, Maxim Rakhuba, and Christoph Schwab. Tensor Rank bounds for Point Singularities in ℝ3. To appear in Adv. Comp. Math.. [ bib | abstract | pdf ]
We analyze rates of approximation by quantized, tensor-structured representations of functions with isolated point singularities in R3. We consider functions in countably normed Sobolev spaces with radial weights and analytic- or Gevrey-type control of weighted semi-norms. Several classes of boundary value and eigenvalue problems from science and engineering are discussed whose solutions belong to the countably normed spaces. It is shown that quantized, tensor-structured approximations of functions in these classes exhibit tensor ranks bounded polylogarithmically with respect to the accuracy ϵ∈(0,1) in the Sobolev space H1. We prove exponential convergence rates of three specific types of quantized tensor decompositions: quantized tensor train (QTT), transposed QTT and Tucker-QTT. In addition, the bounds for the patchwise decompositions are uniform with respect to the position of the point singularity. An auxiliary result of independent interest is the proof of exponential convergence of hp-finite element approximations for Gevrey-regular functions with point singularities in the unit cube Q=(0,1)3. Numerical examples of function approximations and of Schrödinger-type eigenvalue problems illustrate the theoretical results.
@Techreport{MRS19_872,
author = {C. Marcati and M. Rakhuba and Ch. Schwab},
title = {Tensor Rank bounds for Point Singularities in $\mathbb{R}^3$},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2019-68},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2019/2019-68.pdf },
year = {2019}
}

8. Yvon Maday, Carlo Marcati. Analyticity and hp discontinuous Galerkin approximation of nonlinear Schrödinger eigenproblems. Submitted.bib | abstract | pdf ]
We study a class of nonlinear eigenvalue problems of Schrödinger type, where the potential is singular on a set of points. Such proble ms are widely present in physics and chemistry, and their analysis is of both theoretical and practical interest. In particular, we study the regularity of the eigenfunctions of the operators considered, and w e propose and analyze the approximation of the solution via an isotropically refined hp discontinuous Galerkin (dG) method. We show that, for weighted analytic potentials and for up-to-quartic nonlinear ities, the eigenfunctions belong to analytic-type non homogeneous weighted Sobolev spaces. We also prove quasi optimal a priori estimates on the error of the dG finite element method; when using an isotropical ly refined hp space the numerical solution is shown to converge with exponential rate towards the exact eigenfunction. In addition, we investigate the role of pointwise convergence in the doubling of the conve rgence rate for the eigenvalues with respect to the convergence rate of eigenfunctions. We conclude with a series of numerical tests to validate the theoretical results.
  @Techreport{MM19_873,
author = {Y. Maday and C. Marcati},
title = {Analyticity and hp discontinuous Galerkin approximation of nonlinear Schr\"odinger eigenproblems},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2019-69},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2019/2019-69.pdf },
year = {2019}
}

9. Alexey Chernov, Carlo Marcati, and Lorenzo Mascotto. p- and hp- virtual elements for the Stokes problem. Adv. Comput. Math., 47 (2021). [ bib | abstract | pdf ]
We analyse the p- and hp-versions of the virtual element method (VEM) for the the Stokes problem on a polygonal domain. The key tool in the analysis is the existence of a bijection between Poisson-like and Stokes-like VE spaces for the velocities. This allows us to re-interpret the standard VEM for Stokes as a VEM, where the test and trial discrete velocities are sought in Poisson-like VE spaces. The upside of this fact is that we inherit from [Beirão da Veiga, L., Chernov, A., Mascotto, L., Russo, A.: Exponential convergence of the hp virtual element method with corner singularity. Numer. Math. 138(3), 581–613 (2018)] an explicit analysis of best interpolation results in VE spaces, as well as stabilization estimates that are explicit in terms of the degree of accuracy of the method. We prove exponential convergence of the hp-VEM for Stokes problems with regular right-hand sides. We corroborate the theoretical estimates with numerical tests for both the p- and hp-versions of the method.
@ARTICLE{cmm2020,
author = {{Chernov}, Alexey and {Marcati}, Carlo and {Mascotto}, Lorenzo},
title = "{p- and hp- virtual elements for the Stokes problem}",
journal = {arXiv e-prints},
year = 2020,
month = jun,
eid = {arXiv:2006.10644},
pages = {arXiv:2006.10644},
archivePrefix = {arXiv},
eprint = {2006.10644},
primaryClass = {math.NA},
}

10. Dominik Schötzau, Carlo Marcati, and Christoph Schwab. Exponential convergence of mixed hp-DGFEM for the incompressible Navier-Stokes equations in ℝ2. IMA J. Numer. Anal.. [ bib | abstract | pdf ]
In a polygon Ω⊂ℝ2, we consider mixed hp-discontinuous Galerkin approximations of the stationary, incompressible Navier-Stokes equations, subject to no-slip boundary conditions. We use geometrically corner-refined meshes and hp spaces with linearly increasing polynomial degrees. Based on recent results on analytic regularity of velocity field and pressure of Leray solutions in Ω, we prove exponential rates of convergence of the mixed hp-discontinuous Galerkin finite element method (hp-DGFEM), with respect to the number of degrees of freedom, for small data which is piecewise analytic.
@Techreport{SMS20_888,
author = {D. Sch\"otzau and C. Marcati and Ch. Schwab},
title = {Exponential convergence of mixed hp-DGFEM for the incompressible Navier-Stokes equations in R²},
institution = {Seminar for Applied Mathematics, ETH Z{\"u}rich},
number = {2020-15},
url = {https://www.sam.math.ethz.ch/sam_reports/reports_final/reports2020/2020-15.pdf },
year = {2020}
}

11. Carlo Marcati, Christoph Schwab. Analytic Regularity for the incompressible Navier-Stokes Equations in polygons. SIAM J. Math. Anal. 52 (2020), no. 3, 2945–2968. bib | abstract | pdf ]
In a plane polygon P with straight sides, we prove analytic regularity of the Leray-Hopf solution of the stationary, viscous and incompressible Navier-Stokes equations. We assume small data, analytic volume force and no-slip boundary conditions. Analytic regularity is quantified in so-called countably normed, corner-weighted spaces with homogeneous norms. Implications of this analytic regularity include exponential convergence rates of so-called hp-Finite Element and Spectral Element discretizations.
@article{MarcatiSchwab2020,
AUTHOR = {Marcati, Carlo and Schwab, Christoph},
TITLE = {Analytic regularity for the incompressible {N}avier-{S}tokes
equations in polygons},
JOURNAL = {SIAM J. Math. Anal.},
FJOURNAL = {SIAM Journal on Mathematical Analysis},
VOLUME = {52},
YEAR = {2020},
NUMBER = {3},
PAGES = {2945--2968},
ISSN = {0036-1410},
MRCLASS = {35Q30 (35A20 76N10)},
MRNUMBER = {4113068},
DOI = {10.1137/19M1247334},
URL = {https://doi.org/10.1137/19M1247334},
}

12. Yvon Maday, Carlo Marcati. Regularity and hp discontinuous Galerkin finite element approximation of linear elliptic eigenvalue problems with singular potentials. Math. Models Methods Appl. Sci. 29 (2019), no. 8, 1585–1617. [ bib | abstract | pdf ]
We study the regularity in weighted Sobolev spaces of Schrödinger-type eigenvalue problems, and we analyse their approximation via a discontinuous Galerkin (dG) hp finite element method. In particular, we show that, for a class of singular potentials, the eigenfunctions of the operator belong to analytic-type non homogeneous weighted Sobolev spaces. Using this result, we prove that the an isotropically graded hp dG method is spectrally accurate, and that the numerical approximation converges with exponential rate to the exact solution. Numerical tests in two and three dimensions confirm the theoretical results and provide an insight into the the behaviour of the method for varying discretisation parameters.
@article {Maday2018,
AUTHOR = {Maday, Yvon and Marcati, Carlo},
TITLE = {Regularity and {$hp$} discontinuous {G}alerkin finite element
approximation of linear elliptic eigenvalue problems with
singular potentials},
JOURNAL = {Math. Models Methods Appl. Sci.},
FJOURNAL = {Mathematical Models and Methods in Applied Sciences},
VOLUME = {29},
YEAR = {2019},
NUMBER = {8},
PAGES = {1585--1617},
ISSN = {0218-2025},
MRCLASS = {65N25 (35J10 35P05 65N30)},
MRNUMBER = {3986800},
DOI = {10.1142/S0218202519500295},
URL = {https://doi.org/10.1142/S0218202519500295},
}

13. Paola F. Antonietti, Carlo Marcati, Ilario Mazzieri, and Alfio Quarteroni. High order discontinuous Galerkin methods on simplicial elements for the elastodynamics equation. Numerical Algorithms, 71(1):181–206, 2015. [ bib | abstract | http ]
In this work we apply the discontinuous Galerkin (dG) spectral element method on meshes made of simplicial elements for the approximation of the elastodynamics equation. Our approach combines the high accuracy of spectral methods, the geometrical flexibility of simplicial elements and the computational efficiency of dG methods. We analyze the dissipation, dispersion and stability properties of the resulting scheme, with a focus on the choice of different sets of basis functions. Finally, we apply the method on benchmark as well as realistic test cases.
@article{Antonietti2015,
author = {Antonietti, Paola F.
and Marcati, Carlo
and Mazzieri, Ilario
and Quarteroni, Alfio},
title = {High order discontinuous Galerkin methods on simplicial elements for the elastodynamics equation},
journal = {Numerical Algorithms},
year = {2015},
volume = {71},
number = {1},
pages = {181--206},
issn = {1572-9265},
doi = {10.1007/s11075-015-0021-7},
url = {http://dx.doi.org/10.1007/s11075-015-0021-7}
}


### PhD Thesis

Carlo Marcati. Discontinuous hp finite element methods for elliptic eigenvalue problems with singular potentials — with applications to quantum chemistry. Sorbonne Université, 2018. [ bib | abstract ]

In this thesis, we study elliptic eigenvalue problems with singular potentials, motivated by several models in physics and quantum chemistry, and we propose a discontinuous Galerkin hp finite element method for their solution. In these models, singular potentials occur naturally (associated with the interaction between nuclei and electrons). Our analysis starts from elliptic regularity in non homogeneous weighted Sobolev spaces. We show that elliptic operators with singular potential are isomorphisms in those spaces and that we can derive weighted analytic type estimates on the solutions to the linear eigenvalue problems. The isotropically graded hp method provides therefore approximations that converge with exponential rate to the solution of those eigenproblems. We then consider a wide class of nonlinear eigenvalue problems, and prove the convergence of numerical solutions obtained with the symmetric interior penalty discontinuous Galerkin method. Furthermore, when the non linearity is polynomial, we show that we can obtain the same analytic type estimates as in the linear case, thus the numerical approximation converges exponentially. We also analyze under what conditions the eigenvalue converges at an increased rate compared to the eigenfunctions. For both the linear and nonlinear case, we perform numerical tests whose objective is both to validate the theoretical results, but also evaluate the role of sources of errors not considered previously in the analysis, and to help in the design of hp/dG graded methods for more complex problems.
@phdthesis{marcati,
TITLE = {Discontinuous $hp$ finite element methods for elliptic eigenvalue problems with singular potentials},
AUTHOR = {Marcati, Carlo},
URL = {https://tel.archives-ouvertes.fr/tel-02072774},
SCHOOL = {{Sorbonne Universit{\'e}}},
YEAR = {2018},
MONTH = Oct,
KEYWORDS = {hp/dG graded finite element method ; discontinuous Galerkin ; nonlinear eigenvalue problem ; quantum chemistry ; weighted Sobolev spaces ; elliptic regularity},
TYPE = {Theses},
PDF = {https://tel.archives-ouvertes.fr/tel-02072774/file/these.pdf},
HAL_ID = {tel-02072774},
HAL_VERSION = {v1},
} 

### Posters

• YM60, Roscoff, France, May 2017
• Lions-Magenes Days, University of Pavia, Pavia (Italy), April 2015
• Electronic Structure Theory for Materials and (Bio)molecules, UCLA, Los Angeles (USA), July 2014

### Fellowships and grants

• Allocation doctorale DIM RDM-IdF, 2013
• Projet BOUM, SMAI (with Ludovick Gagnon)

## Education

2013-2018
PhD at LJLL under the supervision of prof. Yvon Maday.
2011-2013
Laurea Magistrale cum Laude in Ingegneria Matematica, Politecnico di Milano, Italy. My Master thesis.
2007-2011
Laurea cum Laude in Ingegneria Matematica, Politecnico di Milano, Italy.
My full CV.

### Summer schools

July-August 2015
Hands-on Summer School on Electronic Structure Theory for Materials and (Bio)molecules, UCLA, Los Angeles (USA)
June 2013
Summer School on Scientific Visualization, CINECA, Milano (Italy)