REVIEW 4 major objections 5 minor 30 references
Safe Quadrotor Navigation using Composite Control Barrier Functions
T0 review · 4 major / 5 minor · reviewed 2026-08-08 · deepseek-v4-flash
Pith's one-line read A single composite control barrier function, built by soft-minimizing all distance constraints, can guarantee quadrotor collision avoidance everywhere except on a zero-volume set of singular states.
desk verdict Strong empirical work and a promising idea, but the central safety-filter equation is internally inconsistent with the defined CBF, so the formal guarantee is unsupported as written. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the composite control barrier function $h_1$ in (15), a smooth soft-minimum of saturated third-order distance barriers $\nu_{i,2}$. It carries the argument by turning an arbitrary number of 'keep distance $\epsilon$ from point $x_i$' constraints into a single scalar constraint whose Lie-derivative regressor $L_g h_1$ has the interpretation of a virtual obstacle $\hat{x}(x) = \sum_i \lambda_i \cosh^{-2}(\nu_{i,2}(x)/\gamma) x_i$. The feasibility proof then reduces to analysing when this virtual obstacle is colinear with the thrust direction $Re_3$, which produces the null-set $X_{\mathrm{singular}}$.
What would settle it
Evaluate $L_f h_1 + \alpha_1 h_1$ and $L_g h_1$ over a grid of states near obstacles using the paper's parameters ($\kappa=20$, $\gamma=40$, $\alpha_1=1$, $p_0=-3$, $p_1=-2$) and search for a state where $L_g h_1 = 0$ while $L_f h_1 + \alpha_1 h_1 < 0$; any such state disproves CBF validity at the operating $\kappa$. In simulation with model (9), one can also drive the closed-loop filter toward the singular configuration $(x - \hat{x}) \times Re_3 = 0$ and check whether the QP (23) becomes infeasible.
Extended reading notes
Core claim
The paper proposes the safety filter in (23), which minimally corrects the geometric reference controller (22) by solving a quadratic program subject to the composite CBF constraint (16) and the thrust-rate CBF constraint (19). The composite function $h_1(x) = -\frac{\gamma}{\kappa} \log \sum_i e^{-\kappa \tanh(\nu_{i,2}(x)/\gamma)}$ aggregates all obstacle constraints, where $\nu_{i,2}$ is built from the squared distance to obstacle $i$ through the third-order chain (14). The proof that the infeasible set is negligible rests on the observation that conflicting constraints require the virtual obstacle $\hat{x}(x)$ -- a weighted average of obstacle positions -- to lie exactly on the thrust axis of the vehicle, and that this singular configuration has zero volume and can be left instantaneously by a suitable rotation command. The paper reports that evaluating $h_1$ and its derivatives for $10^4$ obstacles takes 3.4 ms on a laptop CPU, and the hardware experiments show the original distance constraints remaining satisfied over entire missions.
Load-bearing premise
The entire safety guarantee rests on assuming that the composite function $h_1$ in (15) qualifies as a control barrier function at the actual parameter value $\kappa=20$; the paper justifies this with a limit argument as $\kappa \to \infty$ and does not check the necessary and sufficient condition (8) at $\kappa=20$.
Editorial extensions
If this is right
- The safety filter can run at 100 Hz on an onboard computer with hundreds of active obstacles, since evaluating the composite function and its derivatives for $10^4$ obstacles takes 3.4 ms on a laptop CPU.
- Because the filter uses both attitude rates and thrust rate as inputs, it can avoid obstacles without a separate dynamic-inversion layer for the attitude dynamics.
- The only states where joint feasibility can fail form a zero-volume set, so the filter is safe almost everywhere in the formal sense; at those states one constraint must be slackened.
- Hardware flights, including references that actively steer into obstacles, show the original distance constraints $\nu_{i,0}$ remaining positive throughout the missions.
- The approach works with an onboard LiDAR and a 20 cm voxel map, using the 400 closest occupied voxels as obstacle constraints.
Reading between the lines
- As an extension of the paper's own logic, the $\kappa \to \infty$ feasibility argument should be backed by a numerical check of condition (8) at the operating $\kappa=20$; without it, the formal invariance guarantee does not strictly cover the hardware results.
- The zero-volume result suggests that, in continuous time with persistently exciting attitude commands, the probability of landing exactly in $X_{\mathrm{singular}}$ is zero; whether this survives sample-and-hold control at 100 Hz is a testable question the paper does not address.
- The virtual-obstacle interpretation points toward a map-free variant in which a signed distance field is fed directly into the composite barrier, potentially reducing the constraint count while keeping the same filter structure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a safety filter for rate-controlled quadrotors based on a single composite control barrier function h1 formed by a softmin (with tanh saturation) over many distance-based obstacle constraints, together with a thrust constraint h2. The filter solves a QP (23) that minimally corrects a geometric reference controller. The authors analyze recursive feasibility of the two stacked constraints, claim the infeasible configurations form a zero-volume set, report scalability to 10^4 obstacles, and validate the approach in indoor and outdoor hardware experiments against naive and adversarial references.
Significance. If the central derivation were correct, the paper would make a useful engineering contribution: a single composite CBF that scales to thousands of constraints and a feasibility analysis for a third-order quadrotor model, with hardware demonstration. The scalability numbers (Table I) and the adversarial experiments are strong practical evidence. However, as written the key inequality (16) does not match the derivative of the proposed h1, and the finite-kappa CBF validity is not established; the formal claims are therefore not supported by the presented derivation.
major comments (4)
- [Sec. IV-C, Eqs. (15)-(16)] Equation (16) is not the CBF condition for the h1 defined in Eq. (15). Differentiating (15) gives Lf h1 = sum_i [exp(-kappa tanh(nu_i,2/gamma)) / sum_j exp(-kappa tanh(nu_j,2/gamma))] sech^2(nu_i,2/gamma) Lf nu_i,2, and similarly for Lg h1, so the softmin weights are computed with the saturated values tanh(nu_i,2/gamma), not with nu_i,2, and the exponential weight is exp(-kappa(tanh(nu_i,2/gamma) - h1/gamma)), not e^{-kappa(nu_i,2 - h1)} as stated before Eq. (16). In addition, from Eq. (14b) the bracket in Lf h1 should be Lf^3 nu_0 - (p0+p1)Lf^2 nu_0 + p0 p1 Lf nu_0; the displayed expression Lf^3 nu_0 - p1 Lf^2 nu_0 + p1 p0 Lf nu_0 omits the -p0 Lf^2 nu_0 term. Since the proposed safety filter QP (23) uses Eq. (16), the implemented filter is not the one for which invariance is claimed. Please correct the transcription and supply a complete derivation, or state explicitly which h1 is actually used and prove the corresponding invariance condition.
- [Sec. IV-C, after Eq. (16); Table II] For the finite value kappa=20 used in the experiments, the validity of h1 as a CBF is not established. Theorem 2 requires verifying condition (8) at all states where Lg h1 = 0. The paper only invokes the kappa-to-infinity argument in Remark 1, which concerns the limiting softmin of nu_i,2 and does not imply condition (8) for finite kappa. The text even states that 'a large value of kappa is generally chosen to avoid evaluation of (8)', which is not a proof. Since Theorem 2 is stated as an if-and-only-if condition, the invariance guarantee of the set {h1 >= 0} is unsupported for the actual parameters. Please verify condition (8) for the specific h1 in Eq. (15), either symbolically or numerically over the operating region, or provide a theorem that covers the finite-kappa case with the tanh saturation.
- [Sec. IV-D, Proposition 1] The proof of Proposition 1 is a sketch rather than a rigorous zero-measure argument. It computes gradients of q = (x - xhat(x)) x R e3 with respect to r_i3 under kappa -> infinity and asserts that 'one can always choose Omega such that dq/dt != 0', but it does not show that the sublevel set has measure zero, and the phrase 'on either side in discontinuities of xhat(x)' is not a well-defined local argument. Please provide a complete proof, for example an explicit rank/submersion calculation or a measure-theoretic argument, and specify the exact definition of X_singular, including whether the case x = xhat(x) is included.
- [Sec. IV-E, QP (23)] The QP (23) has no slack variables, and the paper does not specify a fallback controller for the set X_singular where the constraints conflict. The statement that 'one of the constraints (16) or (19) must be slackened' is not part of the proposed algorithm. As written, the safety filter has no defined behavior on X_singular, so the claim 'guaranteed everywhere except in the zero volume set X_singular' leaves the algorithm incomplete at exactly the states where the QP is infeasible. Please either add a explicit fallback strategy with a formal guarantee (for example, a bounded-violation slack formulation) or clearly state and justify the operational assumption that X_singular is never reached.
minor comments (5)
- [Sec. V-C, Fig. 3] During experiment A the value of h1 crosses below zero; since a CBF safety filter is intended to enforce h1 >= 0 in the nominal model, please explain whether these crossings are due to the derivative mismatch in Eq. (16) or to disturbances and modeling errors, and quantify the maximum violation.
- [Sec. IV-C, Eq. (17)] The term 'weighted average' for the virtual obstacle is imprecise because the weights appearing in Eq. (16) are not normalized; please specify the normalized weights or adjust the wording.
- [Sec. IV-D, Proposition 1] The statement of Proposition 1 contains an incomplete clause ('…when.'); please complete the sentence and clarify the conditions under which X_singular is a null set.
- [Sec. V-A, Table I] The reported runtimes cover only the composition of h1 and its Lie derivatives, not the solution of the QP (23); it would be helpful to report the full filter-loop time, since the claimed 3.4 ms is not the end-to-end computation.
- [References] Reference [9] is missing author and title information; please complete the bibliographic entry.
Circularity Check
No significant circularity: the composite-CBF derivation rests on external theorems, the hardware self-citation is not load-bearing, and the unverified finite-kappa condition is a rigor gap rather than a circular step.
full rationale
The paper's derivation chain is not circular. The composite CBF construction (Eq. 15), the derivative identity (Eq. 7), and the validity condition (Theorem 2) are taken from the external reference [22] (Molnar and Ames), not from the authors' own prior work. The recursive-feasibility result (Proposition 1) is a new geometric sensitivity argument about the zero set of q = (x - xhat) x R e3; it does not assume the conclusion it is meant to prove, and it is explicitly derived for the 'idealized case (Remark 1)' from the same external theory. The main weakness flagged by a careful reading is that for the finite value kappa=20 used in experiments, the sufficient condition (8) is not verified; the text says 'we resort to the argument taken in Remark 1, i.e. that for one i, lim_{kappa -> infinity} h1(x) = gamma tanh(nu_{i,2}(x)/gamma), while nu_{i,2} is a CBF.' That is an appeal to a limiting argument to justify a finite-kappa design, and it is a correctness/rigor concern, not a circular reduction: the CBF candidate h1 is not defined in terms of the QP constraint (16), nor is any fitted parameter renamed as a prediction. The only self-citation is [16] for the physical quadrotor platform ('The experiments relied on a custom-built quadrotor from [16]'), which does not support any mathematical claim. The scalability numbers in Table I and the hardware experiments are external validation data rather than outputs forced by construction. The equation-level mismatch between the derivative of (15) and the regressor in (16) (e.g., the apparent dropped p0 Lf^2 nu_{i,0} term and the gamma/tanh exponent issue) would be a substantive bug if the implemented filter matches (16) literally, but a bug is not circularity: the central claim does not reduce to its inputs by definition. Score 0.
Assumptions & free parameters
free parameters (9)
- p0 =
-3
- p1 =
-2
- alpha1 =
1
- gamma =
40
- kappa =
20
- epsilon =
0.5
- alpha2 =
5
- epsilon_T =
7.5
- P (QP weight matrix) =
not specified
assumptions (5)
- standard math CBF theory: Theorem 1 (ECBF invariance) and Theorem 2 (softmin composite CBF condition) from [17,22] hold.
- domain assumption The quadrotor is exactly described by the rate-controlled third-order model (9), including thrust dynamics with no actuator delay or input magnitude bounds (except thrust lower bound).
- domain assumption Obstacles can be represented as static points xi with a uniform margin epsilon, using the closest 400 occupied voxels at 20cm resolution, updated at 10Hz.
- ad hoc to paper The condition (8) of Theorem 2 holds for the finite kappa used in practice, based on the kappa-to-infinity argument in Remark 1.
- ad hoc to paper The QP (23) is feasible at all times encountered in operation, i.e. the state never enters X_singular (a zero-volume set) and the solver always returns a solution.
Cite this review
Pith. "Pith review of Safe Quadrotor Navigation using Composite Control Barrier Functions." pith.science (2026). https://pith.science/paper/HYAWGMSQ
@misc{pith2026250204101,
author = {Pith},
title = {Pith review of: Safe Quadrotor Navigation using Composite Control Barrier Functions},
year = {2026},
howpublished = {\url{https://pith.science/paper/HYAWGMSQ}},
note = {Machine review of arXiv:2502.04101}
}
read the original abstract
This paper introduces a safety filter to ensure collision avoidance for multirotor aerial robots. The proposed formalism leverages a single Composite Control Barrier Function from all position constraints acting on a third-order nonlinear representation of the robot's dynamics. We analyze the recursive feasibility of the safety filter under the composite constraint and demonstrate that the infeasible set is negligible. The proposed method allows computational scalability against thousands of constraints and, thus, complex scenes with numerous obstacles. We experimentally demonstrate its ability to guarantee the safety of a quadrotor with an onboard LiDAR, operating in both indoor and outdoor cluttered environments against both naive and adversarial nominal policies.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
A reactive method for collision avoidance in industrial environments,
D. Alejo, J. A. Cobano, G. Heredia, and A. Ollero, “A reactive method for collision avoidance in industrial environments,” Journal of Intelligent & Robotic Systems , vol. 84, pp. 745–758, 2016
work page 2016
-
[2]
J. Dentler, S. Kannan, M. A. O. Mendez, and H. V oos, “A real- time model predictive position control with collision avoidance for commercial low-cost quadrotors,” in 2016 IEEE conference on control applications (CCA) . IEEE, 2016, pp. 519–525
work page 2016
-
[3]
Robust model predictive flight control of unmanned rotorcrafts,
K. Alexis, C. Papachristos, R. Siegwart, and A. Tzes, “Robust model predictive flight control of unmanned rotorcrafts,” Journal of Intelli- gent & Robotic Systems , vol. 81, pp. 443–469, 2016
work page 2016
-
[4]
Uncertainty- aware reinforcement learning for collision avoidance,
G. Kahn, A. Villaflor, V . Pong, P. Abbeel, and S. Levine, “Uncertainty- aware reinforcement learning for collision avoidance,” arXiv preprint arXiv:1702.01182, 2017
arXiv 2017
-
[5]
Reinforcement learning for collision- free flight exploiting deep collision encoding,
M. Kulkarni and K. Alexis, “Reinforcement learning for collision- free flight exploiting deep collision encoding,” arXiv preprint arXiv:2402.03947, 2024
arXiv 2024
-
[6]
Sampling-based algorithms for optimal motion planning,
S. Karaman and E. Frazzoli, “Sampling-based algorithms for optimal motion planning,” The international journal of robotics research , vol. 30, no. 7, pp. 846–894, 2011
2011
-
[7]
A real-time framework for kinodynamic planning with application to quadrotor obstacle avoidance,
R. Allen and M. Pavone, “A real-time framework for kinodynamic planning with application to quadrotor obstacle avoidance,” in AIAA Guidance, Navigation, and Control Conference , 2016, p. 1374
work page 2016
-
[8]
Barrier functions in cascaded controller: Safe quadrotor control,
M. Khan, M. Zafar, and A. Chatterjee, “Barrier functions in cascaded controller: Safe quadrotor control,” in 2020 American Control Confer- ence (ACC), 2020, pp. 1737–1742
work page 2020
Show all 30 references
-
[9]
Dynamic Systems and Control Conference, vol
Safety-Critical Control of a 3D Quadrotor With Range-Limited Sensing, ser. Dynamic Systems and Control Conference, vol. V olume 1, 10 2016. [Online]. Available: https://doi.org/10.1115/ DSCC2016-9913
2016
-
[10]
Comparative analysis of control barrier functions and artificial potential fields for obstacle avoidance,
A. Singletary, K. Klingebiel, J. Bourne, A. Browning, P. Tokumaru, and A. Ames, “Comparative analysis of control barrier functions and artificial potential fields for obstacle avoidance,” in 2021 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) , 2021...
2021
-
[11]
Safe control synthesis for multicopter via control barrier function backstepping,
J. Kim and Y . Kim, “Safe control synthesis for multicopter via control barrier function backstepping,” in 2023 62nd IEEE Conference on Decision and Control (CDC) , 2023, pp. 8720–8725
2023
-
[12]
Control barrier functions in dynamic uavs for kinematic obstacle avoidance: A col- lision cone approach,
M. Tayal, R. Singh, J. Keshavan, and S. Kolathaya, “Control barrier functions in dynamic uavs for kinematic obstacle avoidance: A col- lision cone approach,” in 2024 American Control Conference (ACC) , 2024, pp. 3722–3727
2024
-
[13]
Time-varying soft-maximum barrier functions for safety in unmapped and dynamic environments,
A. Safari and J. B. Hoagg, “Time-varying soft-maximum barrier functions for safety in unmapped and dynamic environments,” 2024. [Online]. Available: https://arxiv.org/abs/2409.01458
2024 arXiv
-
[14]
Control-barrier-aided teleoperation with visual-inertial slam for safe mav navigation in complex environments,
S. Zhou, S. Papatheodorou, S. Leutenegger, and A. P. Schoellig, “Control-barrier-aided teleoperation with visual-inertial slam for safe mav navigation in complex environments,” 2024. [Online]. Available: https://arxiv.org/abs/2403.04331
2024
-
[15]
Learning safe, generalizable perception-based hy- brid control with certificates,
C. Dawson et al. , “Learning safe, generalizable perception-based hy- brid control with certificates,” IEEE Robotics and Automation Letters , vol. 7, no. 2, pp. 1904–1911, 2022
1904
-
[16]
Neural control barrier functions for safe navigation,
M. Harms, M. Kulkarni, N. Khedekar, M. Jacquet, and K. Alexis, “Neural control barrier functions for safe navigation,” in 2024 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), 2024. [Online]. Available: https://arxiv.org/abs/2407.19907
2024 arXiv
-
[17]
Control barrier functions: Theory and applications,
A. D. Ames, S. Coogan, M. Egerstedt, G. Notomista, K. Sreenath, and P. Tabuada, “Control barrier functions: Theory and applications,” in 2019 18th European control conference (ECC) . IEEE, 2019, pp. 3420–3431
2019
-
[18]
Control barrier function based quadratic programs for safety critical systems,
A. D. Ames, X. Xu, J. W. Grizzle, and P. Tabuada, “Control barrier function based quadratic programs for safety critical systems,” IEEE Transactions on Automatic Control , vol. 62, no. 8, pp. 3861–3876, 2016
2016
-
[19]
Exponential control barrier functions for enforcing high relative-degree safety-critical constraints,
Q. Nguyen and K. Sreenath, “Exponential control barrier functions for enforcing high relative-degree safety-critical constraints,” in 2016 American Control Conference (ACC) , 2016, pp. 322–328
2016
-
[20]
High-order control barrier functions,
W. Xiao and C. Belta, “High-order control barrier functions,” IEEE Transactions on Automatic Control , vol. 67, no. 7, pp. 3655–3662, 2022
2022
-
[21]
Safe backstep- ping with control barrier functions,
A. J. Taylor, P. Ong, T. G. Molnar, and A. D. Ames, “Safe backstep- ping with control barrier functions,” in 2022 IEEE 61st Conference on Decision and Control (CDC) , 2022, pp. 5775–5782
2022
-
[22]
Composing control barrier functions for complex safety specifications,
T. G. Molnar and A. D. Ames, “Composing control barrier functions for complex safety specifications,” IEEE Control Systems Letters , 2023
2023
-
[23]
Composition of control barrier functions with differing relative degrees for safety under input constraints,
P. Rabiee and J. B. Hoagg, “Composition of control barrier functions with differing relative degrees for safety under input constraints,” in 2024 American Control Conference (ACC) , 2024, pp. 3692–3697
2024
-
[24]
Geometric tracking control of a quadrotor UA V on SE(3),
T. Lee, M. Leok, and N. H. McClamroch, “Geometric tracking control of a quadrotor UA V on SE(3),” in 49th IEEE Conference on Decision and Control (CDC) , 2010, pp. 5420–5425
2010
-
[25]
Control of complex maneuvers for a quadrotor uav using geometric methods on se(3),
——, “Control of complex maneuvers for a quadrotor uav using geometric methods on se(3),” 2011
2011
-
[26]
Complementary multi–modal sensor fusion for resilient robot pose estimation in subterranean environments,
S. Khattak, H. Nguyen, F. Mascarich, T. Dang, and K. Alexis, “Complementary multi–modal sensor fusion for resilient robot pose estimation in subterranean environments,” in 2020 International Con- ference on Unmanned Aircraft Systems (ICUAS) . IEEE, 2020, pp. 1024–1029
2020
-
[27]
V oxblox: Incremental 3d euclidean signed distance fields for on- board mav planning,
H. Oleynikova, Z. Taylor, M. Fehr, R. Siegwart, and J. Nieto, “V oxblox: Incremental 3d euclidean signed distance fields for on- board mav planning,” in 2017 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS) , 2017, pp. 1366–1373
2017
-
[28]
Accurate tracking of aggressive quadrotor trajectories using incremental nonlinear dynamic inversion and differ- ential flatness,
E. Tal and S. Karaman, “Accurate tracking of aggressive quadrotor trajectories using incremental nonlinear dynamic inversion and differ- ential flatness,” IEEE Transactions on Control Systems Technology , vol. 29, no. 3, pp. 1203–1218, 2021
2021
-
[29]
CasADi – A software framework for nonlinear optimization and optimal control,
J. A. E. Andersson, J. Gillis, G. Horn, J. B. Rawlings, and M. Diehl, “CasADi – A software framework for nonlinear optimization and optimal control,” Mathematical Programming Computation , vol. 11, no. 1, pp. 1–36, 2019
2019
-
[30]
qpOASES: A parametric active-set algorithm for quadratic program- ming,
H. Ferreau, C. Kirches, A. Potschka, H. Bock, and M. Diehl, “qpOASES: A parametric active-set algorithm for quadratic program- ming,” Mathematical Programming Computation , vol. 6, no. 4, pp. 327–363, 2014
2014
Reviewed August 8, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.