2026-08-07 | | Total: 33
Increasing penetration of inverter-based resources in today's power systems requires substitution of the contribution from synchronous generators to dynamic voltage stability and inertial response. However, established approaches for power system planning are sequential, as stabilising measures are only considered at a later stage of the planning process. We investigate the advantages of an integrated approach for power system planning, where stabilising measures are considered simultaneously with the expansion of generation, transmission, and storage systems via simplified stability constraints on inertia and voltage stability. We find that system costs are reduced with the integrated approach and that the dual-use option of grid-forming battery energy storage systems is favoured over other stabilising measures like static synchronous compensators.
Fault diagnosis and prognosis in discrete event systems are studied in the scenario where the observations are possibly received with delay. To address this scenario, two conditions for diagnosis and prognosis with delayed observations are proposed, where we show that the state-of-the-art notion of prognosability must be revised to avoid conservativeness. Diagnosability and prognosability conditions are then verified by introducing a delay observer and a new verification function. Theoretical analysis indicates the effectiveness of the verification method for fault diagnosis and prognosis in the system.
We study stabilization of multidimensional continu- ity equations with source terms on bounded domains by means of localized interior flux feedback. The feedback is prescribed through the divergence of the flux and is chosen so that the error with respect to a reference profile satisfies a transport equation with localized damping. The main geometric condition is a finite-time characteristic damping inequality, requiring relevant characteristics to accumulate a uniform amount of damping over a time horizon. This condition is shown to yield exponential stability in L2 of the error, under a gain condition relating localized damping to compressive amplification of the transport field. Lyapunov-type entrance conditions ensure characteris- tic damping on support-restricted families of trajectories. A weighted Lyapunov functional provides a differential Lyapunov criterion and an input-to-state (ISS) estimate with respect to additive perturbations. We also discuss elliptic right-inverse realizations of the feedback flux and extend the characteristic damping argument to velocity fields depending nonlinearly on the state. A two-dimensional example finally illustrates the geometric, gain, and realization conditions.
We study Muon, a recently proposed matrix-aware optimization method, in the context of the Stiefel manifold. This manifold consists of matrices with orthonormal columns and is ubiquitous in machine learning and scientific computing. Existing extensions of Muon to this manifold rely on heuristic, approximate, or iterative updates with varying computational efficiency. We show that the corresponding Stiefel Muon update admits an exact closed-form solution and use this result to develop Skewon, a practical algorithm for orthogonality-constrained optimization with an efficient implementation. We further establish first-order convergence guarantees for Skewon in the smooth non-convex setting.
A framework is developed for the inference of dynamics described by a generalized system of ordinary differential equations. A stochastic gradient method is coined that infers dynamics from observed marginal probability density functions using the joint probability density function of the observable and latent variables. Diffusion and other irreversible processes observed in a low-dimensional state can be recast as deterministic, reversible flows in a sufficiently augmented state space, where the joint density satisfies the hyperbolic Liouville equation. The marginal distribution observed is the projection of these hyperbolic dynamics onto the observed coordinates, with the latent components carrying the randomness and memory. This reframing allows inference for irreversible or stochastic dynamics into the recovery of a deterministic Ordinary Differential Equation (ODE) from marginal observations. Instead of solving the high-dimensional Liouville equation for the joint density, the algorithm exploits its characteristic representation. Particles sampled from the initial distribution are transported along characteristic lines. The Eulerian sensitivity with respect to parameters is obtained by sensitivity propagation along the characteristic lines, with a crossed U-statistic producing an unbiased gradient estimator, which enables stochastic gradient descent. Four experiments validate the method: recovery of a three-mode linear system observed through the marginal of a single mode; a nonlinear Gompertz growth model with a hidden mode; a bistable system whose hidden mode turns a unimodal marginal bimodal; and Stokes--Oseen drag law recovery for particles in a cellular flow. Convergence behavior is analyzed across these settings.
We study stochastic extragradient (SEG) methods for solving monotone variational inequality problems (VIPs) over a feasible set. Although extragradient is a foundational algorithm for VIPs and its deterministic convergence theory is well developed, its stochastic counterpart remains less understood. Most existing analyses focus on independent-sample SEG (I-SEG) and assume either that the domain is compact or that the variance of the stochastic operator is uniformly bounded. The behavior of same-sample SEG (S-SEG), a natural variant with materially different properties, has received far less attention. In this work, we address these gaps in the literature. We first show that S-SEG is sensitive to samplewise Lipschitz parameters: mean Lipschitzness and bounded variance alone do not ensure convergence, even on a compact set. Then, for possibly unbounded domains, we establish a high-probability restricted-gap convergence for each SEG variant under a relaxed set of assumptions, and show that certain fundamental improvements to these results are impossible in general. Finally, we show that a known asymmetric double step-size selection that guarantees almost sure last-iterate convergence for I-SEG can fail for S-SEG: there exists a stochastic monotone VIP for which S-SEG diverges almost surely even under the modified step-sizes.
We determine the exact weak limited augmented Zarankiewicz numbers $z_{wL}(m,3)$ for $m=5,6,7,8,9$. We establish \[ z_{wL}(5,3)=10,\quad z_{wL}(6,3)=12,\quad z_{wL}(7,3)=14,\quad z_{wL}(8,3)=16,\quad z_{wL}(9,3)=18. \] Thus $z_{wL}(m,3)=2m$ for all $5\le m\le 9$. The lower bounds are given by explicit weak admissible constructions, while the matching upper bounds for $m=7,8,9$ are obtained by exact finite searches over canonical $C_4$-free base graphs. Together with the known small cases, this strongly supports the conjecture that \[ z_{wL}(m,3)=2m \quad \text{for all } m\ge 5. \] Since BSR$(m,n)\ge z_{wL}(m,n)$, these constructions imply \[ {\rm BSR}(m,3)\ge 2m \quad \text{for } m=5,6,7,8,9. \]
This note studies the distributed fault diagnosis problem in partially-observed discrete event systems, where the system is monitored by a group of agents to cooperatively diagnose faults within a finite number of steps. The novelty of this work is the creation of a methodology to verify when the faults can be diagnosed even in the presence of transmission delay impairments. To address this scenario, a new distributed diagnosability condition is proposed, which extends decentralized diagnosability conditions proposed in the literature. Such distributed diagnosability condition is then verified via a novel structure named delay recorder and a new diagnosis function. Theoretical analysis shows that the verification method can successfully determine whether the faults can be diagnosed.
Policy iteration (PI) is an important reinforcement learning tool for solving optimal control problems which includes an initialization stage, i.e., the search for an initial stabilizing controller. However, the initialization stage typically relies on complete model information, thereby imposing substantial constraints on the initialization of model-free PI. For stochastic systems with multiplicative noise dependent on state and control, the stability is not ensured by Hurwitz conditions as in the deterministic case, but rather by a Lyapunov-type inequality that incorporates both drift and diffusion terms. Therefore, the corresponding model-free PI initialization problem is more challenging. To this end, a novel spectrum assignment method is proposed to obtain an initial stabilizer for PI in continuous-time indefinite stochastic linear quadratic control. With the help of the Lyapunov-type operator's spectrum, the original system is gradually approximated from the stable auxiliary system by adjusting a cumulative factor, thereby obtaining a stabilizing control gain. Furthermore, by leveraging system data and adjusting the cumulative factor, we design a model-free algorithm that does not rely on an initial stabilizing policy and can achieve optimal control. Finally, simulation results are provided to validate the effectiveness of the proposed methods.
Bi-level interdiction models are frequently solved by metaheuristics whose solution quality cannot be assessed, because exact optima are unavailable at the scales tested. We supply them for the bi-level dual-defense attacker model (BDAM), which couples node interdiction, edge destruction and capacitated supply support, and which we previously solved by a hybrid metaheuristic. First, BDAM's dominant attacker-path term admits an exact single-level reformulation by lower-level dualization, a reduction available whenever arc lengths are linear in the defender's binary decisions; the resulting mixed-integer program certifies optimality on all eighteen three-row and five-row configurations of our earlier benchmark, seventeen in under ten seconds, with strong uncertified incumbents out to 15*30 grids. Our published averages sit 3.55% below that frontier, and the gap widens with scale. Second, the supply rule, like any rule priced on a single attacker shortest path, is ill-posed under ties, and our threat-corridor formulation is tie-invariant by construction. The tie-break moves the objective by under 10^(-3) but the realized supply cost by up to 0.83 units, so the defect is suppressed by the objective weight rather than absent. Third, MILP-DA pairs the exact anchor with the corridor decode and a feasibility repair; compared against our published figures with no re-implementation on either side, it wins on thirteen of eighteen certified configurations and on all eighteen larger ones, every loss falling on a three-row grid.
This article presents a priority-aware robust coverage optimization framework for an aerial sensor network under sensor location uncertainty. Each region is assigned a priority weight, and the objective is to maximize the weighted coverage while maintaining robustness against positional perturbations. A mathematical optimization model is developed by incorporating surveillance constraints and an RRF-based robustness formulation into the proposed framework. An efficient priority-aware robust orientation optimization (PAROO) algorithm is then proposed to determine the sensor orientations that maximize the weighted coverage objective. Experimental results on an airport-inspired surveillance scenario demonstrate that the proposed framework effectively directs sensing resources toward high-priority regions and achieves higher weighted coverage than representative baseline approaches, highlighting its practical applicability in security-sensitive environments.
The degree sequence optimization problem is to find a subgraph of a given graph which maximizes the sum over all vertices of a given function evaluated at the subgraph degree of that vertex. Here we study this problem and its complexity for quadratic functions. In particular, we introduce the critical roots of a graph, and show they define intervals over which the optimal value of the problem, as the quadratic root varies, is convex piecewise affine.
Constrained density functional theory (CDFT) provides a powerful framework for describing electronically excited and charge-localized states, which underlie a broad range of physical and chemical phenomena. However, the discretized optimization problems arising from CDFT calculations remain challenging, owing to the presence of both the Stiefel manifold constraint and additional nonconvex quadratic constraints. Existing algorithms either fail to enforce the quadratic constraints with high accuracy or face convergence issues due to double-loop iterative structures. In this paper, we first derive a subspace-splitting reformulation that decouples the two groups of constraints, by exploiting the inherent rotation invariance and introducing a nonlinear subspace alignment constraint. Based on this reformulation, we propose a single-loop damped alternating direction method of multipliers, called DASSP. To the best of our knowledge, DASSP is the first algorithm for CDFT calculations with rigorous convergence guarantees. Each iteration of DASSP comprises a spectral minimization step, a projected gradient step, and a damped dual ascent step, all of which admit efficient implementations. Numerical results on synthetic and realistic CDFT problems demonstrate that DASSP attains high feasibility accuracy and exhibits favorable efficiency without compromising robustness. We expect that this work will pave the way toward reliable and efficient large-scale CDFT applications.
We consider linear regression models fitted by minimizing Bregman losses of the form \[ \frac{1}{n}\sum_{i=1}^n \left[ φ(y_i)-φ(x_i^\topθ) -φ'(x_i^\topθ) \bigl(y_i-x_i^\topθ\bigr) \right], \] where \(y_i\) is the observed response and \(x_i^\topθ\) is the linear prediction. Even when the generating potential \(φ\) is strongly convex, the resulting regression objective may be nonconvex in \(θ\). This creates a gap between the convexity of the generating potential and the optimization geometry of the fitted model. The Hessian can be written as a weighted Gram matrix whose weights depend on the derivatives of the potential and the current residuals. This representation gives simple conditions for local strong convexity, smoothness, and conditional linear convergence of gradient descent. For the quadratic--quartic potential, we derive an exact scalar convexity condition, identify the interval of negative curvature, and obtain local and global sufficient conditions for positive curvature. Numerical experiments indicate that these scalar conditions may fail while the full Hessian remains positive definite at the evaluated points. They also indicate that the range of tested gradient-descent step sizes leading to convergence decreases as the quartic parameter grows. These results characterize how residual-dependent curvature interacts with the generating potential and the design matrix in reverse Bregman regression.
In this paper, we propose a fast gradient algorithm for the problem of minimizing a differentiable (possibly nonconvex) function in Hilbert spaces. We first extend the dry friction property for convex functions to what we call the dry-like friction property in a nonconvex setting, and then employ a line search technique to adaptively update parameters at each iteration. Depending on the choice of parameters, the proposed algorithm exhibits subsequential convergence to a critical point or full sequential convergence to an ``approximate'' critical point of the objective function. We also establish the full sequential convergence to a critical point under the Kurdyka--Lojasiewicz (KL) property of a merit function. Thanks to the parameters' flexibility, our algorithm can reduce to a number of existing inertial gradient algorithms with Hessian damping and dry friction. By exploiting variational properties of the Moreau envelope, the proposed algorithm is adapted to address weakly convex nonsmooth optimization problems. In particular, we extend the result on KL exponent for the Moreau envelope of a convex KL function to a broad class of KL functions that are not necessarily convex nor continuous. Simulation results illustrate the efficiency of our algorithm and demonstrate the potential advantages of combining dry-like friction with extrapolation and line search techniques.
Understanding how experts make decisions and being able to transfer that knowledge is important, especially in complex engineering applications. It is highly valuable for training novices, improving the performance of human-machine systems, and potentially enabling fully autonomous systems that perform as well as human experts. However, an expert's decision-making strategy, developed through years of experience, is often not directly accessible, since the implicit preferences and decision rules involved can be difficult to specify explicitly. This has motivated the use of observed decisions made by the expert to learn an interpretable model that captures the expert's decision-making process. In this work, we develop an inverse optimization approach to jointly learn the decision-maker's preferences (or perceived costs) and the decision rules governing their choices. We demonstrate the general applicability of our approach using three case studies that consider a shift assignment problem, a production planning problem, and a real-world routing problem, respectively. Across these case studies, modeling both perceived costs and decision rules leads to better predictions, highlighting the value of the proposed framework and its greater flexibility in capturing and replicating expert decision making.
We study continuous-time projection neural networks for variational inequalities on closed convex sets. A positive definite matrix that depends on the state preconditions the operator, and its inverse defines the projection metric. Existing convergence analyses of this flow cover Hessian-generated inverse metrics and state-dependent scalar metrics. In the first case, a Bregman distance eliminates metric-derivative terms. We treat a matrix metric whose inverse is of neither kind. We prove joint regularity of the projection in its argument and metric. Under common spectral bounds, the Euclidean Lipschitz estimate improves from the squared bound to the bound itself. For Lipschitz strongly monotone operators and Lipschitz metrics, we prove local exponential convergence with explicit rate and radius. On compact feasible sets, an explicit metric-variation bound yields global exponential convergence. A twice continuously differentiable, uniformly positive definite metric and a strongly monotone linear operator produce an annulus of periodic orbits. The inverse metric violates Hessian integrability throughout the annulus. Deterministic computations confirm the analytic formulas and quantify slack in both convergence certificates.
We present an open-source software package that implements a provably convergent Benders-type decomposition algorithm for multistage stochastic integer programs. In addition to standard cut families, such as Benders, strengthened Benders, and Lagrangian cuts, the algorithm incorporates rectified linear unit (ReLU) cuts, which provide convergence guarantees for general mixed-integer state variables. However, the dual problems used to generate these cuts often admit multiple optimal solutions. Although each solution yields a valid cut that separates the incumbent, the resulting cuts can differ in how well they approximate the subproblem cost. To strengthen these cuts, our package implements and evaluates two cut-selection strategies based on normalization and regularization of the dual problem. We also incorporate an alternating-cut criterion that uses cheaper Benders cuts when they are effective and invokes more expensive tight cuts only when necessary. Computational experiments on four classes of multistage stochastic integer programs benchmark these methods and provide insights on how problem structure affects their practical performance.
The Yau-Yau filter reformulates nonlinear state estimation as probability-density propagation governed by the Forward Kolmogorov equation (FKE). Applying it to target tracking, however, requires efficient FKE approximation on a finite computational domain. This paper proposes a Recentered-Domain Yau-Yau Filter (RD-YYF), which solves the FKE within a fixed-size local window centered at the latest state estimate. This design concentrates numerical resolution near the dominant posterior density. Offline, physics-informed neural networks (PINNs) generate FKE solution snapshots, while principal component analysis constructs a low-dimensional representation of density evolution. A lightweight residual surrogate maps the initial-condition coefficients and domain center to the terminal-solution coefficients. Online, the pretrained surrogate predicts density evolution within the recentered window, followed by observation update and state estimation. Experiments on two geometrically constrained target-tracking examples show that RD-YYF achieves lower tracking errors than the extended Kalman filter (EKF) and particle filter (PF), while retaining efficient per-timestep inference. Ablation results indicate that domain recentering improves density approximation in high-probability regions and accelerates offline PINN convergence. These results demonstrate the potential of RD-YYF for efficient nonlinear target tracking
Iterate convergence of Bregman proximal methods (BPMs) has long remained open, especially for nonconvex objectives. Recently, \citet{chen2026skl} made progress by establishing iterate convergence for a BPM via the so-called scaled Kurdyka-Łojasiewicz (SKŁ) property, but only for the Shannon entropy kernel and linearly constrained problems. In this paper, we develop a unified iterate convergence framework that applies to a broad group of kernels and composite objective functions. Our approach extends the analytical tools in \cite{chen2026skl}, in particular the SKŁ property, which plays a central role in ensuring convergence of the generated sequences. By introducing kernel-dependent parameterization functions, we show that the extended SKŁ property holds for all continuous subanalytic functions, particularly when the kernel has a closed domain. We then verify that the assumptions of the framework are satisfied by standard BPMs under mild regularity conditions, thereby establishing their iterate convergence for a wide range of objective functions. Furthermore, based on the parameterization functions, we show that the continuous-time BPM (mirror flow) converges to a stationary point for o-minimal definable objective functions, yielding the first trajectory convergence result for mirror flow without imposing convexity assumptions on the objective function or isolation assumptions on stationary points. Taken together, these discrete- and continuous-time convergence results provide a unified trajectory convergence theory for BPMs.
In this paper, we consider a class of multiblock nonconvex nonsmooth optimization problems, which covers many applications such as the analysis of pre-earthquake anomalies and machine learning. To solve this class of problems, we propose the inertial block proximal linearized method with two-phase adaptive momentum (IBPL$^+$-TP). Compared to the current methods, our method possesses three main advantages: (1) it introduces a two-phase adaptive momentum strategy to effectively update the extrapolation parameters, (2) it allows using two different extrapolation points to accelerate the convergence, (3) it allows the extrapolation parameters of these two extrapolation points to be independent of and unconstrained by all other parameters. While maintaining the above advantages, we prove that our method ensures the monotonic convergence of the objective function of this class of problems, and we also prove that the sequence generated by our method globally converges to a critical point, as well as establish the convergence rate of our method. To demonstrate the effectiveness of our method, we apply it to solve two nonconvex and nonsmooth machine learning problems, namely sparse nonnegative matrix factorization with $\ell_0$-constraints and sparse nonnegative CP decomposition with $\ell_0$-constraints. The numerical experimental results on solving these problems show that our method outperforms several state-of-the-art methods.
This work introduces a proximal stochastic subgradient method for minimizing the sum of an expected cost, whose integrand is potentially nonsmooth and nonconvex, and a lower semicontinuous, prox-bounded function. We target a broad class of integrands obeying a nonsmooth, localized variant of the descent lemma in the decision variable, a structural assumption that simultaneously covers smooth losses with Lipschitz gradient and differences of such losses with convex functions. At each iteration the expected cost is replaced by a sample average that is progressively refined, and the proximal-subgradient stepsize is selected by an Armijo-type line search enforcing a sufficient-decrease property up to stochastic errors induced by the sample-based approximation. This framework accommodates substantially more general problem formulations than existing methods, in particular, it requires neither (weak) convexity of the regularizer nor a uniform bound on the variance of the stochastic oracle, and our analysis yields convergence guarantees that are new even in the smooth setting. Specifically, we establish almost sure convergence of the sequence of function values and stationarity of every accumulation point of the trajectories under the relaxed requirement that the sample-size sequence be merely nondecreasing and unbounded, with no prescribed growth rate. Leveraging the Kurdyka-Lojasiewicz (KL) property, we further upgrade this subsequential guarantee to convergence of the whole trajectory to a single stationary point. Finally, for exponential-type KL desingularizing functions and polynomially growing sample sizes, we derive explicit polynomial convergence rates, up to a logarithmic factor, for both the function values and the iterates.
This study addresses the critical yet under researched area of supported evacuation for vulnerable populations during wildfires, such as hospital patients and long term care residents, by developing a two stage stochastic optimization model that optimizes facility location, fleet sizing, and vehicle routing under strict time windows. To overcome the problem NP hard complexity, the authors propose an innovative solution methodology leveraging Logic Based Benders Decomposition, featuring Combinatorial Benders Cuts and logic based inequalities. Extensive numerical experiments and real world data from a community wildfire drill in Roxborough Park, Colorado, demonstrate that the proposed approach yields high quality solutions, significantly improving shelter placement, vehicle utilization, and overall evacuation efficiency compared to alternative policies.
A Hider distributes two identical balls among three boxes whose search costs satisfy $a\ge b\ge c>0$. A Searcher opens boxes adaptively until both balls are found; every opening incurs the corresponding box cost and recovers at most one ball. We give an exact solution of this heterogeneous multi-look search-cost game. The value is the maximum of three explicit rational functions. In the three parameter regimes, an optimal Hider strategy is the product-form distribution restricted respectively to three, five, or all six placements. Every non-wasteful deterministic Searcher policy is payoff-equivalent to one of $72$ elementary decision trees, which induce only $42$ distinct opening-count profiles. An analytic argument settles the first regime. Two exact Searcher mixtures settle the second, and four settle the third. Feasibility of the parameterized mixtures over the full cost region is established by exact rational Bernstein-basis certificates. Independent implementations reproduce the policy set, the zero-sum linear-program values, and all $154$ dyadic nodes examined by the Bernstein verifier.
This paper introduces a contour integral method (CIM) for efficiently computing outputs of parametric linear systems in control form over specified time intervals and to a user-prescribed accuracy. The CIM approximates the inverse Laplace transform via a quadrature rule applied along a modified integration contour. For parametric systems, we show how CIM integrates effectively with projection-based model order reduction (MOR) where a greedy algorithm builds the projection spaces following an error estimate we derive for this setting. We additionally demonstrate that the developed projection framework naturally enforces Hermite interpolation conditions. This combination substantially lowers the cost of evaluating the input-output relations across the parameter domain, for a wide range of input functions, and for initial conditions well captured by a low-dimensional subspace.. We demonstrate the accuracy and efficiency of the approach on benchmark non-parametric and parametric control systems, comparing against state-of-the-art projection-based MOR methods.