BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Enrique Zuazua

Enrique Zuazua contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.OCmath.APmath.NANumerical AnalysisMachine LearningArtificial IntelligenceDistributed, Parallel, and Cluster Computingeess.SYnlin.AOquant-phSystems and Control

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

A Structure-Preserving Numerical Scheme for Optimal Control and Design of Mixing in Incompressible Flows

We develop a structure-preserving computational framework for optimal mixing control in incompressible flows. Our approach exactly conserves the continuous system's key invariants (mass and $L^2$-energy), while also maintaining discrete state-adjoint duality at every time step. These properties are achieved by integrating a centered finite-volume discretization in space with a time-symmetric Crank-Nicolson integrator for both the forward advection and its adjoint, all inside a gradient-based optimization loop. The result is a numerical solver that is faithful to the continuous optimality conditions and efficiently computes mixing-enhancing controls. In our numerical tests, the optimized time-dependent stirring produces a nearly exponential decay of a chosen mix-norm, achieving orders-of-magnitude faster mixing than any single steady flow. To our knowledge, this work provides the first evidence that enforcing physical structure at the discrete level can lead to both exact conservation and highly effective mixing outcomes in optimal flow design.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Moments, Time-Inversion and Source Identification for the Heat Equation

We address the initial source identification problem for the heat equation, a notably ill-posed inverse problem characterized by exponential instability. Departing from classical Tikhonov regularization, we propose a novel approach based on moment analysis of the heat flow, transforming the problem into a more stable inverse moment formulation. By evolving the measured terminal time moments backward through their governing ODE system, we recover the moments of the initial distribution. We then reconstruct the source by solving a convex optimization problem that minimizes the total variation of a measure subject to these moment constraints. This formulation naturally promotes sparsity, yielding atomic solutions that are sums of Dirac measures. Compared to existing methods, our moment-based approach reduces exponential error growth to polynomial growth with respect to the terminal time. We provide explicit error estimates on the recovered initial distributions in terms of moment order, terminal time, and measurement errors. In addition, we develop efficient numerical discretization schemes and demonstrate significant stability improvements of our approach through comprehensive numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Towards the Next Frontier of LLMs, Training on Private Data: A Cross-Domain Benchmark for Federated Fine-Tuning

The recent success of large language models (LLMs) has been largely driven by vast public datasets. However, the next frontier for LLM development lies beyond public data. Much of the world's most valuable information is private, especially in highly regulated sectors such as healthcare and finance, where data include patient histories or customer communications. Unlocking this data could represent a major leap forward, enabling LLMs with deeper domain expertise and stronger real-world utility. Yet, these data cannot be shared because they are distributed across institutions and constrained by privacy, regulatory, and organizational barriers. Moreover, institutional datasets are typically non-independent and identically distributed (non-IID), differing across sites in population characteristics, data modalities, documentation patterns, and task-specific label distributions. In this paper, we demonstrate a practical approach to unlocking private and distributed institutional data for LLM adaptation through federated collaboration across data silos. Built on the Sherpa.ai Federated Learning platform, our framework enables nodes to jointly fine-tune a shared LLM without exchanging private data. We evaluate this approach through a cross-domain benchmark in healthcare and finance, using four closed-ended question answering and classification datasets: MedQA, MedMCQA, FPB, and FiQA-SA. We compare three parameter-efficient fine-tuning (PEFT) strategies-LoRA, QLoRA, and IA3-across pretrained backbones under non-IID settings reflecting institutional data heterogeneity. Our results show that federated fine-tuning performs close to centralized training and outperforms isolated single-institution learning. From a Green AI perspective, QLoRA and IA3 improve efficiency with limited accuracy degradation, supporting federated PEFT as a viable approach for adapting LLMs where data cannot be shared.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Transmutation based Quantum Simulation for Non-unitary Dynamics

We present a quantum algorithm for simulating dissipative diffusion dynamics generated by positive semidefinite operators of the form $A=L^\dagger L$, a structure that arises naturally in standard discretizations of elliptic operators. Our main tool is the Kannai transform, which represents the diffusion semigroup $e^{-TA}$ as a Gaussian-weighted superposition of unitary wave propagators. This representation leads to a linear-combination-of-unitaries implementation with a Gaussian tail and yields query complexity $\tilde{\mathcal{O}}(\sqrt{\|A\| T \log(1/\varepsilon)})$, up to standard dependence on state-preparation and output norms, improving the scaling in $\|A\|, T$ and $\varepsilon$ compared with generic Hamiltonian-simulation-based methods. We instantiate the method for the heat equation and biharmonic diffusion under non-periodic physical boundary conditions, and we further use it as a subroutine for constant-coefficient linear parabolic surrogates arising in entropy-penalization schemes for viscous Hamilton--Jacobi equations. In the long-time regime, the same framework yields a structured quantum linear solver for $A\mathbf{x}=\mathbf{b}$ with $A=L^\dagger L$, achieving $\tilde{\mathcal{O}}(κ^{3/2}\log^2(1/\varepsilon))$ queries and improving the condition-number dependence over standard quantum linear-system algorithms in this factorized setting.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Universal Approximation of Dynamical Systems by Semi-Autonomous Neural ODEs and Applications

In this paper, we introduce semi-autonomous neural ordinary differential equations (SA-NODEs), a variation of the vanilla NODEs, employing fewer parameters. We investigate the universal approximation properties of SA-NODEs for dynamical systems from both a theoretical and a numerical perspective. Within the assumption of a finite-time horizon, under general hypotheses we establish an asymptotic approximation result, demonstrating that the error vanishes as the number of parameters goes to infinity. Under additional regularity assumptions, we further specify this convergence rate in relation to the number of parameters, utilizing quantitative approximation results in the Barron space. Based on the previous result, we prove an approximation rate for transport equations by their neural counterparts. Our numerical experiments validate the effectiveness of SA-NODEs in capturing the dynamics of various ODE systems and transport equations. Additionally, we compare SA-NODEs with vanilla NODEs, highlighting the superior performance and reduced complexity of our approach.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A framework for randomized time-splitting in linear-quadratic optimal control

Inspired by the successes of stochastic algorithms in the training of deep neural networks and the simulation of interacting particle systems, we propose and analyze a framework for randomized time-splitting in linear-quadratic optimal control. In our proposed framework, the linear dynamics of the original problem is replaced by a randomized dynamics. To obtain the randomized dynamics, the system matrix is split into simpler submatrices and the time interval of interest is split into subintervals. The randomized dynamics is then found by selecting randomly one or more submatrices in each subinterval. We show that the dynamics, the minimal values of the cost functional, and the optimal control obtained with the proposed randomized time-splitting method converge in expectation to their analogues in the original problem when the time grid is refined. The derived convergence rates are validated in several numerical experiments. Our numerical results also indicate that the proposed method can lead to a reduction in computational cost for the simulation and optimal control of large-scale linear dynamical systems.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Greedy search of optimal approximate solutions

In this paper we develop a procedure to deal with a family of parameter-dependent ill-posed problems, for which the exact solution in general does not exist. The original problems are relaxed by considering corresponding approximate ones, whose optimal solutions are well dfined, where the optimality is determined by the minimal norm requirement. The procedure is based upon greedy algorithms that preserve, at least asymptotically, Kolmogorov approximation rates. In order to provide a-priori estimates for the algorithm, a Tychonff-type regularization is applied, which adds an additional parameter to the model. The theory is developed in an abstract theoretical framework that allows its application to different kinds of problems. We present a specific example that considers a family of ill-posed elliptic problems. The required general assumptions in this case translate to rather natural uniform lower and upper bounds on coefficients of the considered operators.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Local Stability and Convergence of Unconstrained Model Predictive Control

The local stability and convergence for Model Predictive Control (MPC) of unconstrained nonlinear dynamics based on a linear time-invariant plant model is studied. Based on the long-time behavior of the solution of the Riccati Differential Equation (RDE), explicit error estimates are derived that clearly demonstrate the influence of the two critical parameters in MPC: the prediction horizon $T$ and the control horizon $τ$. In particular, if the MPC-controller has access to an exact (linear) plant model, the MPC-controls and the corresponding optimal state trajectories converge exponentially to the solution of an infinite-horizon optimal control problem when $T-τ\rightarrow \infty$. When the difference between the linear model and the nonlinear plant is sufficiently small in a neighborhood of the origin, the MPC strategy is locally stabilizing and the influence of modeling errors can be reduced by choosing the control horizon $τ$ smaller. The obtained convergence rates are validated in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the decay of one-dimensional locally and partially dissipated hyperbolic systems

We study the time-asymptotic behavior of linear hyperbolic systems under partial dissipation which is localized in suitable subsets of the domain. More precisely, we recover the classical decay rates of partially dissipative systems satisfying the stability condition (SK) with a time-delay depending only on the velocity of each component and the size of the undamped region. To quantify this delay, we assume that the undamped region is a bounded space-interval and that the system without space-restriction on the dissipation satisfies the stability condition (SK). The former assumption ensures that the time spent by the characteristics of the system in the undamped region is finite and the latter that whenever the damping is active the solutions decay. Our approach consists in reformulating the system into n coupled transport equations and showing that the time-decay estimates are delayed by the sum of the times each characteristics spend in the undamped region.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reachable set for Hamilton-Jacobi equations with non-smooth Hamiltonian and scalar conservation laws

We give a full characterization of the range of the operator which associates, to any initial condition, the viscosity solution at time $T$ of a Hamilton-Jacobi equation with convex Hamiltonian. Our main motivation is to be able to treat the case of convex Hamiltonians with no further regularity assumptions. We give special attention to the case $H(p) = |p|$, for which we provide a rather geometrical description of the range of the viscosity operator by means of an interior ball condition on the sublevel sets. From our characterization of the reachable set, we are able to deduce further results concerning, for instance, sharp regularity estimates for the reachable functions, as well as structural properties of the reachable set. The results are finally adapted to the case of scalar conservation laws in dimension one.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turnpike in Lipschitz-nonlinear optimal control

We present a new proof of the turnpike property for nonlinear optimal control problems, when the running target is a steady control-state pair of the underlying system. Our strategy combines the construction of quasi-turnpike controls via controllability, and a bootstrap argument, and does not rely on analyzing the optimality system or linearization techniques. This in turn allows us to address several optimal control problems for finite-dimensional, control-affine systems with globally Lipschitz (possibly nonsmooth) nonlinearities, without any smallness conditions on the initial data or the running target. These results are motivated by applications in machine learning through deep residual neural networks, which may be fit within our setting. We show that our methodology is applicable to controlled PDEs as well, such as the semilinear wave and heat equation with a globally Lipschitz nonlinearity, once again without any smallness assumptions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turnpike in optimal control of PDEs, ResNets, and beyond

The \emph{turnpike property} in contemporary macroeconomics asserts that if an economic planner seeks to move an economy from one level of capital to another, then the most efficient path, as long as the planner has enough time, is to rapidly move stock to a level close to the optimal stationary or constant path, then allow for capital to develop along that path until the desired term is nearly reached, at which point the stock ought to be moved to the final target. Motivated in part by its nature as a resource allocation strategy, over the past decade, the turnpike property has also been shown to hold for several classes of partial differential equations arising in mechanics. When formalized mathematically, the turnpike theory corroborates the insights from economics: for an optimal control problem set in a finite-time horizon, optimal controls and corresponding states, are close (often exponentially), during most of the time, except near the initial and final time, to the optimal control and corresponding state for the associated stationary optimal control problem. In particular, the former are mostly constant over time. This fact provides a rigorous meaning to the asymptotic simplification that some optimal control problems appear to enjoy over long time intervals, allowing the consideration of the corresponding stationary problem for computing and applications. We review a slice of the theory developed over the past decade --the controllability of the underlying system is an important ingredient, and can even be used to devise simple turnpike-like strategies which are nearly optimal--, and present several novel applications, including, among many others, the characterization of Hamilton-Jacobi-Bellman asymptotics, and stability estimates in deep learning via residual neural networks.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Differentiability with respect to the initial condition for Hamilton-Jacobi equations

We prove that the viscosity solution to a Hamilton-Jacobi equation with a smooth convex Hamiltonian of the form $H(x,p)$ is differentiable with respect to the initial condition. Moreover, the directional Gâteaux derivatives can be explicitly computed almost everywhere in $\mathbb{R}^N$ by means of the optimality system of the associated optimal control problem. We also prove that, in the one-dimensional case in space and in the quadratic case in any space dimension, these directional Gâteaux derivatives actually correspond to the unique duality solution to the linear transport equation with discontinuous coefficient, resulting from the linearization of the Hamilton-Jacobi equation. The motivation behind these differentiability results arises from the following optimal inverse-design problem: given a time horizon $T>0$ and a target function $u_T$, construct an initial condition such that the corresponding viscosity solution at time $T$ minimizes the $L^2$-distance to $u_T$. Our differentiability results allow us to derive a necessary first-order optimality condition for this optimization problem, and the implementation of gradient-based methods to numerically approximate the optimal inverse design.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The turnpike property in nonlinear optimal control -- A geometric approach

This paper presents, using dynamical system theory, a framework for investigating the turnpike property in nonlinear optimal control. First, it is shown that a turnpike-like property appears in general dynamical systems with hyperbolic equilibrium and then, apply it to optimal control problems to obtain sufficient conditions for the turnpike occurs. The approach taken is geometric and gives insights for the behaviors of controlled trajectories, allowing us to find simpler proofs for existing results on the turnpike properties. Attempts to remove smallness restrictions for initial and target states are also discussed based on the geometry of (un)stable manifold and exponential stabilizability of control systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A stochastic approach to the synchronization of coupled oscillators

This paper deals with an optimal control problem associated to the Kuramoto model describing the dynamical behavior of a network of coupled oscillators. Our aim is to design a suitable control function allowing us to steer the system to a synchronized configuration in which all the oscillators are aligned on the same phase. This control is computed via the minimization of a given cost functional associated with the dynamics considered. For this minimization, we propose a novel approach based on the combination of a standard Gradient Descent (GD) methodology with the recently-developed Random Batch Method (RBM) for the efficient numerical approximation of collective dynamics. Our simulations show that the employment of RBM improves the performances of the GD algorithm, reducing the computational complexity of the minimization process and allowing for a more efficient control calculation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Asymptotic behavior of scalar convection-diffusion equations

In these lecture notes, we address the problem of large-time asymptotic behaviour of the solutions to scalar convection-diffusion equations set in ${R}^N$. The large-time asymptotic behaviour of the solutions to many convection-diffusion equations is strongly linked with the behavior of the initial data at infinity. In fact, when the initial datum is integrable and of mass $M$, the solutions to the equations under consideration oftentimes behave like the associated self-similar profile of mass $M$, thus emphasising the role of scaling variables in these scenarios. However, these equations can also manifest other asymptotic behaviors, including weakly non-linear, linear or strongly non-linear behavior depending on the form of the convective term. We give an exhaustive presentation of several results and techniques, where we clearly distinguish the role of the spatial dimension and the form of the nonlinear convective term. Translation (English) by Borjan Geshovski

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Control under constraints for multi-dimensional reaction-diffusion monostable and bistable equations

Dynamic phenomena in social and biological sciences can often be modeled by employing reaction-diffusion equations. When addressing the control of these modes, from a mathematical viewpoint one of the main challenges is that, because of the intrinsic nature of the models under consideration, the solution, typically a proportion or a density function, needs to preserve given lower and upper bounds. Controlling the system to the desired final configuration then becomes complex, and sometimes even impossible. In the present work, we analyze the controllability to constant steady-states of spatially homogeneous semilinear heat equations, with constraints in the state, and using boundary controls, which is indeed a natural way of acting on the system in the present context. The nonlinearities considered are among the most frequent: monostable and bistable ones. We prove that controlling the system to a constant steady-state may become impossible when the diffusivity is too small (or when the domain is large), due to the existence of barrier functions. When such an obstruction does not arise, we build sophisticated control strategies combining the dissipativity of the system, the existence of traveling waves and some connectivity of the set of steady-states. This connectivity allows building paths that the controlled trajectories can follow, in a long time, with small oscillations, preserving the natural constraints of the system. This kind of strategy was successfully implemented in one-space dimension, where phase plane analysis techniques allowed to decode the nature of the set of steady-states. These techniques fail in the present multi-dimensional setting. We employ a fictitious domain technique, extending the system to a larger ball, and building paths of radially symmetric solution that can then be restricted to the original domain.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Model predictive control with random batch methods for a guiding problem

We model, simulate and control the guiding problem for a herd of evaders under the action of repulsive drivers. The problem is formulated in an optimal control framework, where the drivers (controls) aim to guide the evaders (states) to a desired region of the Euclidean space. The numerical simulation of such models quickly becomes unfeasible for a large number of interacting agents. To reduce the computational cost, we use the Random Batch Method (RBM), which provides a computationally feasible approximation of the dynamics. At each time step, the RBM randomly divides the set of particles into small subsets (batches), considering only the interactions inside each batch. Due to the averaging effect, the RBM approximation converges to the exact dynamics as the time discretization gets finer. We propose an algorithm that leads to the optimal control of a fixed RBM approximated trajectory using a classical gradient descent. The resulting control is not optimal for the original complete system, but rather for the reduced RBM model. We then adopt a Model Predictive Control (MPC) strategy to handle the error in the dynamics. While the system evolves in time, the MPC strategy consists in periodically updating the state and computing the optimal control over a long-time horizon, which is implemented recursively in a shorter time-horizon. This leads to a semi-feedback control strategy. Through numerical experiments we show that the combination of RBM and MPC leads to a significant reduction of the computational cost, preserving the capacity of controlling the overall dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Shape turnpike for linear parabolic PDE models

We introduce and study the turnpike property for time-varying shapes, within the viewpoint of optimal control. We focus here on second-order linear parabolic equations where the shape acts as a source term and we seek the optimal time-varying shape that minimizes a quadratic criterion. We first establish existence of optimal solutions under some appropriate sufficient conditions. We then provide necessary conditions for optimality in terms of adjoint equations and, using the concept of strict dissipativity, we prove that state and adjoint satisfy the measure-turnpike property, meaning that the extremal time-varying solution remains essentially close to the optimal solution of an associated static problem. We show that the optimal shape enjoys the exponential turnpike property in term of Hausdorff distance for a Mayer quadratic cost. We illustrate the turnpike phenomenon in optimal shape design with several numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Finite-Time Turnpike Phenomenon for Optimal Control Problems: Stabilization by Non-Smooth Tracking Terms

In this paper, problems of optimal control are considered where in the objective function, in addition to the control cost there is a tracking term that measures the distance to a desired stationary state. The tracking term is given by some norm and therefore it is in general not differentiable. In the optimal control problem, the initial state is prescribed. We assume that the system is either exactly controllable in the classical sense or nodal profile controllable. We show that both for systems that are governed by ordinary differential equations and for infinite-dimensional systems, for example for boundary control systems governed by the wave equation, under certain assumptions the optimal system state is steered exactly to the desired state after finite time.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The inverse problem for Hamilton-Jacobi equations and semiconcave envelopes

We study the inverse problem, or inverse design problem, for a time-evolution Hamilton-Jacobi equation. More precisely, given a target function $u_T$ and a time horizon $T>0$, we aim to construct all the initial conditions for which the viscosity solution coincides with $u_T$ at time $T$. As it is common in this kind of nonlinear equations, the target might not be reachable. We first study the existence of at least one initial condition leading the system to the given target. The natural candidate, which indeed allows determining the reachability of $u_T$, is the one obtained by reversing the direction of time in the equation, considering $u_T$ as terminal condition. In this case, we use the notion of backward viscosity solution, that provides existence and uniqueness for the terminal-value problem. We also give an equivalent reachability condition based on a differential inequality, that relates the reachability of the target with its semiconcavity properties. Then, for the case when $u_T$ is reachable, we construct the set of all initial conditions for which the solution coincides with $u_T$ at time $T$. Note that in general, such initial conditions are not unique. Finally, for the case when the target $u_T$ is not necessarily reachable, we study the projection of $u_T$ on the set of reachable targets, obtained by solving the problem backward and then forward in time. This projection is then identified with the solution of a fully nonlinear obstacle problem, and can be interpreted as the semiconcave envelope of $u_T$, i.e. the smallest reachable target bounded from below by $u_T$.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A parabolic approach to the control of opinion spreading

We analyze the problem of controlling to consensus a nonlinear system modeling opinion spreading. We derive explicit exponential estimates on the cost of approximately controlling these systems to consensus, as a function of the number of agents N and the control time-horizon T. Our strategy makes use of known results on the controllability of spatially discretized semilinear parabolic equations. Both systems can be linked through time-rescaling

0 citations0 reviewsNo code linkedOpen access