Showing 1–50 of 329 results for author: Colombo, L

Symmetry reduction and reconstruction in contact geometry and Lagrange-Poincaré-Herglotz equations

Authors: Alexandre Anahory Simoes, Leonardo Colombo, Manuel de Leon, Modesto Salgado, Silvia Souto

Abstract: In this paper, we investigate the reduction process of a contact Lagrangian system whose Lagrangian is invariant under a group of symmetries. We give explicit coordinate expressions of the resulting reduced differential equations, the so-called Lagrange-Poincare-Herglotz equations. Our framework relied on the associated Herglotz vector field and its projected vector field, and the use of well-chos… ▽ More In this paper, we investigate the reduction process of a contact Lagrangian system whose Lagrangian is invariant under a group of symmetries. We give explicit coordinate expressions of the resulting reduced differential equations, the so-called Lagrange-Poincare-Herglotz equations. Our framework relied on the associated Herglotz vector field and its projected vector field, and the use of well-chosen quasi-velocities. Some examples are also discussed. △ Less

Submitted 13 August, 2024; originally announced August 2024.

Comments: 31 pages

Enhancing heat transport in MoS2 via defect-engineering

Authors: Riccardo Dettori, Francesco Siddi, Luciano Colombo, Claudio Melis

Abstract: MoS2 is one of the most investigated and promising transition-metal dichalcogenides. Its popularity stems from the interesting properties of the monolayer phase, which can serve as the fundamental block for numerous applications. In this paper, we propose an atomistic perspective on the modulation of thermal transport properties in monolayer MoS2 through strategic defect engineering, specifically… ▽ More MoS2 is one of the most investigated and promising transition-metal dichalcogenides. Its popularity stems from the interesting properties of the monolayer phase, which can serve as the fundamental block for numerous applications. In this paper, we propose an atomistic perspective on the modulation of thermal transport properties in monolayer MoS2 through strategic defect engineering, specifically the introduction of sulfur vacancies. Using a combination of molecular dynamics simulations and lattice dynamics calculations, we show how various distributions of sulfur vacancies -- ranging from random to periodically arranged configurations -- affect its thermal conductivity. Notably, we observe that certain periodic arrangements restore the thermal conductivity of the pristine system, due to a minimized interaction between acoustic and optical phonons facilitated by the imposed superperiodicity. This research deepens the understanding of phononic heat transport in two-dimensional materials and introduces a different point-of-view for phonon engineering in nanoscale devices, offering a pathway to enhance device performance and longevity through tailored thermal management strategies. △ Less

Submitted 29 July, 2024; originally announced July 2024.

NITRO-D: Native Integer-only Training of Deep Convolutional Neural Networks

Authors: Alberto Pirillo, Luca Colombo, Manuel Roveri

Abstract: Quantization has become increasingly pivotal in addressing the steadily increasing computational and memory requirements of Deep Neural Networks (DNNs). By reducing the number of bits used to represent weights and activations (typically from 32-bit floating-point to 16-bit or 8-bit integers), quantization reduces the memory footprint, energy consumption, and execution time of DNN models. However,… ▽ More Quantization has become increasingly pivotal in addressing the steadily increasing computational and memory requirements of Deep Neural Networks (DNNs). By reducing the number of bits used to represent weights and activations (typically from 32-bit floating-point to 16-bit or 8-bit integers), quantization reduces the memory footprint, energy consumption, and execution time of DNN models. However, traditional quantization methods typically focus on the inference of DNNs, while the training process still relies on floating-point operations. To date, only one work in the literature has addressed integer-only training for Multi-Layer Perceptron (MLP) architectures. This work introduces NITRO-D, a new framework for training arbitrarily deep integer-only Convolutional Neural Networks (CNNs) that operate entirely in the integer-only domain for both training and inference. NITRO-D is the first framework in the literature enabling the training of integer-only CNNs without the need to introduce a quantization scheme. Specifically, NITRO-D introduces a novel architecture integrating multiple integer local-loss blocks, which include the proposed NITRO Scaling Layer and the NITRO-ReLU activation function. Additionally, it introduces a novel integer-only learning algorithm derived from Local Error Signals (LES), utilizing IntegerSGD, an optimizer specifically designed to operate in an integer-only context. NITRO-D is implemented in an open-source Python library. Extensive experimental evaluations demonstrate its effectiveness across several state-of-the-art image recognition datasets. Results show significant performance improvements from 2.47% to 5.96% for integer-only MLP architectures over the state-of-the-art solution, and the capability of training integer-only CNN architectures with minimal accuracy degradation from -0.15% to -4.22% compared to floating-point LES. △ Less

Submitted 12 September, 2024; v1 submitted 16 July, 2024; originally announced July 2024.

Comments: 15 pages, 3 figures

ACM Class: I.2.6

Characterization of Acoustic Losses in Interdigitated VHF to mmWave Piezoelectric M/NEMS Resonators

Authors: Luca Colombo, Gabriel Giribaldi, Ryan Tetro, Jack Guida, Walter Gubinelli, Luca Spagnuolo, Nicol Maietta, Siddhartha Ghosh, Matteo Rinaldi

Abstract: This work reports on a technology-agnostic and frequency-independent methodology combining a-priori modeling, Finite Element Analysis (FEA), and experimental results for the characterization of acoustic losses in interdigitated piezoelectric micro- and nano-electromechanical (M/NEMS) resonators. The proposed approach models the mechanical quality factor (Qm) and its dependency on piezoelectric (Qp… ▽ More This work reports on a technology-agnostic and frequency-independent methodology combining a-priori modeling, Finite Element Analysis (FEA), and experimental results for the characterization of acoustic losses in interdigitated piezoelectric micro- and nano-electromechanical (M/NEMS) resonators. The proposed approach models the mechanical quality factor (Qm) and its dependency on piezoelectric (Qpiezo) and metal (Qmetal) acoustic losses, as a function of the mode of vibration dispersion characteristics. The model is finally experimentally validated by exploiting the intrinsic on-chip multifrequency manufacturability of interdigitated devices. A broad range of available resonator technologies, frequencies, and piezoelectric materials are investigated for this purpose, including bulk X-cut Lithium Niobate (XLN) leaky surface acoustic wave resonators operating at Ultra High Frequency (UHF), thin film XLN Lamb Wave resonators spanning between Very High Frequency (VHF) and the Ku band, and Aluminum Nitride (AlN) and scandium-doped AlN (ScAlN) cross-sectional lame' mode resonators ranging from the Ku to Ka band (mmWave). △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: To be presented at IUS 2024

2-16 GHz Multifrequency X-Cut Lithium Niobate NEMS Resonators on a Single Chip

Authors: Ryan Tetro, Luca Colombo, Matteo Rinaldi

Abstract: This work presents the design, fabrication, and testing of X-Cut Lithium Niobate (LN) acoustic nanoelectromechanical (NEMS) Laterally Vibrating Resonators (LVRs) and Degenerate LVRs (d-LVRs) operating in the S0 (YZ30) and SH0 (YZ-10) modes between 2 to 16 GHz range, monolithically fabricated on a single chip. The NEMS topology is optimized to extend the aforementioned fundamental modes in the C-,… ▽ More This work presents the design, fabrication, and testing of X-Cut Lithium Niobate (LN) acoustic nanoelectromechanical (NEMS) Laterally Vibrating Resonators (LVRs) and Degenerate LVRs (d-LVRs) operating in the S0 (YZ30) and SH0 (YZ-10) modes between 2 to 16 GHz range, monolithically fabricated on a single chip. The NEMS topology is optimized to extend the aforementioned fundamental modes in the C-, X-, and Ku-bands while preserving performance and mass manufacturability. The devices present acoustic wavelengths (λ) varying between 1800 and 400 nm and are fabricated on a 100 nm ultra-thin LN film on high resistivity silicon with a 3-mask process. Experimental results highlighted quality factor at resonance (Qs) and mechanical quality factors (Qm) as high as 477 and 1750, respectively, and electromechanical coupling (kt2) as high as 32.7%. Large kt2 (>10%) are recorded over a broad range of frequencies (2 - 8 GHz), while Qm exceeding 100 are measured up to 15 GHz. Further enhancement to performance and range of operation on the same chip can be achieved by decreasing λ, refining the fabrication process, and optimizing device topology. These additional steps can help pave the way for manufacturing high-performance resonators on a single chip covering the entire 1 - 25 GHz spectrum. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 4 pages, 5 figures, 4 tables, to be presented at NEMS 2024 in Kyoto, Japan

Comparative analysis of spin wave imaging using nitrogen vacancy centers and time resolved magneto-optical measurements

Authors: Carolina Lüthi, Lukas Colombo, Franz Vilsmeier, Christian Back

Abstract: Spin waves, the fundamental excitations in magnetic materials, are promising candidates for realizing low-dissipation information processing in spintronics. The ability to visualize and manipulate coherent spin-wave transport is crucial for the development of spin wave-based devices. We use a recently discovered method utilizing nitrogen vacancy (NV) centers, point defects in the diamond lattice,… ▽ More Spin waves, the fundamental excitations in magnetic materials, are promising candidates for realizing low-dissipation information processing in spintronics. The ability to visualize and manipulate coherent spin-wave transport is crucial for the development of spin wave-based devices. We use a recently discovered method utilizing nitrogen vacancy (NV) centers, point defects in the diamond lattice, to measure spin waves in thin film magnetic insulators by detecting their magnetic stray field. We experimentally demonstrate enhanced contrast in the detected wavefront amplitudes by imaging spin waves underneath a reference stripline and phenomenologically model the results. By extracting the spin wave dispersion and comparing NV center based spin wave measurements to spin wave imaging conducted through the well-established time-resolved magneto-optical Kerr effect, we discuss the advantages and limitations of employing NV centers as spin wave sensors. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 7 pages, 6 figures

Comparison of two numerical methods for Riemannian cubic polynomials on Stiefel manifolds

Authors: Alexandre Anahory Simoes, Leonardo Colombo, Fátima Silva Leite

Abstract: In this paper we compare two numerical methods to integrate Riemannian cubic polynomials on the Stiefel manifold $\textbf{St}_{n,k}$. The first one is the adjusted de Casteljau algorithm, and the second one is a symplectic integrator constructed through discretization maps. In particular, we choose the cases of $n=3$ together with $k=1$ and $k=2$. The first case is diffeomorphic to the sphere and… ▽ More In this paper we compare two numerical methods to integrate Riemannian cubic polynomials on the Stiefel manifold $\textbf{St}_{n,k}$. The first one is the adjusted de Casteljau algorithm, and the second one is a symplectic integrator constructed through discretization maps. In particular, we choose the cases of $n=3$ together with $k=1$ and $k=2$. The first case is diffeomorphic to the sphere and the quasi-geodesics appearing in the adjusted de Casteljau algorithm are actually geodesics. The second case is an example where we have a pure quasi-geodesic different from a geodesic. We provide a numerical comparison of both methods and discuss the obtained results to highlight the benefits of each method. △ Less

Submitted 30 April, 2024; originally announced April 2024.

Comments: 12 pages, 3 figures

Scandium Aluminum Nitride Overmoded Bulk Acoustic Resonators for Future Wireless Communication

Authors: Walter Gubinelli, Pietro Simeoni, Ryan Tetro, Luca Colombo, Matteo Rinaldi

Abstract: This work reports on the modeling, fabrication, and experimental characterization of a 13 GHz 30% Scandium-doped Aluminum Nitride (ScAlN) Overmoded Bulk Acoustic Resonator (OBAR) for high-frequency Radio Frequency (RF) applications, notably in 5G technology and beyond. The Finite Element Analysis (FEA) optimization process targets the top and bottom metal electrode thicknesses, balancing the elect… ▽ More This work reports on the modeling, fabrication, and experimental characterization of a 13 GHz 30% Scandium-doped Aluminum Nitride (ScAlN) Overmoded Bulk Acoustic Resonator (OBAR) for high-frequency Radio Frequency (RF) applications, notably in 5G technology and beyond. The Finite Element Analysis (FEA) optimization process targets the top and bottom metal electrode thicknesses, balancing the electromechanical coupling coefficient and acoustic energy distribution to enhance device Figure of Merit (FOM). Experimental results on fabricated devices employing platinum and aluminum as bottom and top electrode, respectively, demonstrate a quality factor at resonance (Qs) of 210 and a coupling coefficient (kt2) of 5.2% at 13.3 GHz for the second bulk thickness overtone, effectively validating the simulation framework and hinting at the possible implementation of OBARs for advanced RF filters in 5G networks. △ Less

Submitted 23 April, 2024; originally announced April 2024.

Comments: 4 pages, 4 figures, submitted to 2024 Hilton Head Workshop Proceedings

The Interplay Between Symmetries and Impact Effects on Hybrid Mechanical Systems

Authors: William Clark, Leonardo Colombo, Anthony Bloch

Abstract: Hybrid systems are dynamical systems with continuous-time and discrete-time components in their dynamics. When hybrid systems are defined on a principal bundle we are able to define two classes of impacts for the discrete-time transition of the dynamics: interior impacts and exterior impacts. In this paper we define hybrid systems on principal bundles, study the underlying geometry on the switchin… ▽ More Hybrid systems are dynamical systems with continuous-time and discrete-time components in their dynamics. When hybrid systems are defined on a principal bundle we are able to define two classes of impacts for the discrete-time transition of the dynamics: interior impacts and exterior impacts. In this paper we define hybrid systems on principal bundles, study the underlying geometry on the switching surface where impacts occur and we find conditions for which both exterior and interior impacts are preserved by the mechanical connection induced in the principal bundle. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 6 pages. To be presented at a conference. Comments welcome

Exploring Dielectric Properties in Models of Amorphous Boron Nitride

Authors: Thomas Galvani, Ali K. Hamze, Laura Caputo, Onurcan Kaya, Simon Dubois, Luigi Colombo, Viet-Hung Nguyen, Yongwoo Shin, Hyeon-Jin Shin, Jean-Christophe Charlier, Stephan Roche

Abstract: We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about $100$ atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electro… ▽ More We report a theoretical study of dielectric properties of models of amorphous Boron Nitride, using interatomic potentials generated by machine learning. We first perform first-principles simulations on small (about $100$ atoms in the periodic cell) sample sizes to explore the emergence of mid-gap states and its correlation with structural features. Next, by using a simplified tight-binding electronic model, we analyse the dielectric functions for complex three dimensional models (containing about $10.000$ atoms) embedding varying concentrations of ${\rm sp^{1}, sp^{2}}$ and ${\rm sp^3}$ bonds between B and N atoms. Within the limits of these methodologies, the resulting value of the zero-frequency dielectric constant is shown to be influenced by the population density of such mid-gap states and their localization characteristics. We observe nontrivial correlations between the structure-induced electronic fluctuations and the resulting dielectric constant values. Our findings are however just a first step in the quest of accessing fully accurate dielectric properties of as-grown amorphous BN of relevance for interconnect technologies and beyond. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 27 pages, 10 figures

Amorphous Boron Nitride as a Diffusion Barrier to Cu Atoms

Authors: Onurcan Kaya, Hyeongjoon Kim, Byeongkyu Kim, Luigi Colombo, Hyeon-Jin Shin, Ivan Cole, Hyeon Suk Shin, Stephan Roche

Abstract: This study focuses on amorphous boron nitride ($\rm α$-BN) as a novel diffusion barrier for advanced semiconductor technology, particularly addressing the critical challenge of copper diffusion in back-end-of-logic (BEOL) interconnects. Owing to its ultralow dielectric constant and robust barrier properties, $\rm α$-BN is examined as an alternative to conventional low-k dielectrics. The investigat… ▽ More This study focuses on amorphous boron nitride ($\rm α$-BN) as a novel diffusion barrier for advanced semiconductor technology, particularly addressing the critical challenge of copper diffusion in back-end-of-logic (BEOL) interconnects. Owing to its ultralow dielectric constant and robust barrier properties, $\rm α$-BN is examined as an alternative to conventional low-k dielectrics. The investigation primarily employs theoretical modeling, using a Gaussian Approximation Potential, to simulate and understand the atomic-level interactions and barrier mechanisms of $\rm α$-BN. This machine learning-based approach allows for realistic simulations of its amorphous structure, enabling the exploration of the impact of different film morphologies on barrier efficacy. Complementing the theoretical study, experimental analyses are conducted on Plasma-Enhanced Chemical Vapor Deposition (PECVD) grown $\rm α$-BN films, evaluating their effectiveness in preventing copper diffusion in silicon-based substrates. The results from both the theoretical and experimental investigations highlight the potential of $\rm α$-BN as a highly effective diffusion barrier, suitable for integration in nanoelectronics. This research not only proposes $\rm α$-BN as a promising candidate for BEOL interconnects but also demonstrates the synergy of advanced computational models and experimental methods in material innovation for semiconductor applications. △ Less

Submitted 2 February, 2024; originally announced February 2024.

Comments: 16 pages, 9 figures

LSPE-STRIP on-sky calibration strategy using bright celestial sources

Authors: R. T. Génova-Santos, M. Bersanelli, C. Franceschet, M. Gervasi, C. López-Caraballo, L. Mandelli, M. Maris, A. Mennella, J. A. Rubiño-Martín, F. Villa, M. Zannoni, C. Baccigalupi, B. Caccianiga, L. Colombo, F. Cuttaia, F. Farsian, G. Morgante, S. Paradiso, G. Polenta, S. Ricciardi, M. Sandri, A. Taylor, L. Terenzi, M. Tomasi

Abstract: In this paper we describe the global on-sky calibration strategy of the LSPE-Strip instrument. Strip is a microwave telescope operating in the Q- and W-bands (central frequencies of 43 and 95 GHz respectively) from the Observatorio del Teide in Tenerife, with the goal to observe and characterise the polarised Galactic foreground emission, and complement the observations of the polarisation of the… ▽ More In this paper we describe the global on-sky calibration strategy of the LSPE-Strip instrument. Strip is a microwave telescope operating in the Q- and W-bands (central frequencies of 43 and 95 GHz respectively) from the Observatorio del Teide in Tenerife, with the goal to observe and characterise the polarised Galactic foreground emission, and complement the observations of the polarisation of the cosmic microwave background to be performed by the LSPE-SWIPE instrument and other similar experiments operating at higher frequencies to target the detection of the B-mode signal from the inflationary epoch of the Universe. Starting from basic assumptions on some of the instrument parameters (NET, 1/f noise knee frequency, beam properties, observing efficiency) we perform realistic simulations to study the level of accuracy that can be achieved through observations of bright celestial calibrators in the Strip footprint (sky fraction of 30 %) on the determination and characterisation of the main instrument parameters: global and relative gain factors (in intensity and in polarisation), polarisation direction, polarisation efficiency, leakage from intensity to polarisation, beams, window functions and pointing model. △ Less

Submitted 16 January, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: 46 pages, 20 figures. This paper is part of of the paper series "The LSPE/Strip instrument description and testing", to appear in Jinst. Comments are welcome

Reduction by Symmetry and Optimal Control with Broken Symmetries on Riemannian Manifolds

Authors: Jacob R. Goodman, Leonardo J. Colombo

Abstract: This paper studies the reduction by symmetry of variational problems on Lie groups and Riemannian homogeneous spaces. We derive the reduced equations of motion in the case of Lie groups endowed with a left-invariant metric, and on Lie groups that admits a bi-invariant metric. We repeated this analysis for Riemannian homogeneous spaces, where we derive the reduced equations by considering an altern… ▽ More This paper studies the reduction by symmetry of variational problems on Lie groups and Riemannian homogeneous spaces. We derive the reduced equations of motion in the case of Lie groups endowed with a left-invariant metric, and on Lie groups that admits a bi-invariant metric. We repeated this analysis for Riemannian homogeneous spaces, where we derive the reduced equations by considering an alternative variational problem written in terms of a connection on the horizontal bundle of the underlying Lie group. We study also the case that the underlying Lie group admits a bi-invariant metric, and consider the special case that the homogeneous space is in fact a Riemannian symmetric space. These ideas are applied to geodesics for a rigid body on $SO(3)$ to derive geodesic equations on the dual of its Lie algebra (a vector space), the heavy-top in $SE(3)$ to derive reduced equations of motion on the unit sphere $S^2$, geodesics on $S^2$ as a Riemannian symmetric space endowed with a bi-invariant metric and optimal control problems for applications to robotic manipulators. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 33 pages. arXiv admin note: text overlap with arXiv:2207.13574

Virtual Constraints on Lie groups

Authors: E. Stratoglou, A. Anahory Simoes, A. Bloch, L. Colombo

Abstract: This paper studies virtual constraints on Lie groups. Virtual constraints are invariant relations imposed on a control system via feedback. In this work, we introduce the notion of \textit{virtual constraints on Lie groups}, in particular, \textit{virtual nonholonomic constraints on Lie groups}, in a geometric framework. More precisely, this object is a controlled invariant subspace associated wit… ▽ More This paper studies virtual constraints on Lie groups. Virtual constraints are invariant relations imposed on a control system via feedback. In this work, we introduce the notion of \textit{virtual constraints on Lie groups}, in particular, \textit{virtual nonholonomic constraints on Lie groups}, in a geometric framework. More precisely, this object is a controlled invariant subspace associated with an affine connection mechanical control system on the Lie algebra associated with the Lie group which is the configuration space of the system. We demonstrate the existence and uniqueness of a control law defining a virtual nonholonomic constraint and we characterize the trajectories of the closed-loop system as solutions of a mechanical system associated with an induced constrained connection. Moreover, we characterize the reduced dynamics for nonholonomic systems in terms of virtual nonholonomic constraints, i.e., we characterize when can we obtain reduced nonholonomic dynamics from virtual nonholonomic constraints. We illustrate the theory with some examples and present simulation results for an application. △ Less

Submitted 29 December, 2023; originally announced December 2023.

Comments: 18 pages

On the integrability of hybrid Hamiltonian systems

Authors: Asier López-Gordón, Leonardo J. Colombo

Abstract: A hybrid system is a system whose dynamics are controlled by a mixture of both continuous and discrete transitions. The integrability of Hamiltonian systems is often identified with complete integrability or Liouville integrability, that is, the existence of as many independent integrals of motion in involution as the dimension of the phase space. Under certain regularity conditions, Liouville-Arn… ▽ More A hybrid system is a system whose dynamics are controlled by a mixture of both continuous and discrete transitions. The integrability of Hamiltonian systems is often identified with complete integrability or Liouville integrability, that is, the existence of as many independent integrals of motion in involution as the dimension of the phase space. Under certain regularity conditions, Liouville-Arnold theorem states that the invariant geometric structure associated with Liouville integrability is a fibration by Lagrangian tori, on which motions are linear. In this paper, we study an extension of the Liouville-Arnold theorem for hybrid systems whose continuous dynamics is given by a Hamiltonian vector field and we state conditions on the impact map and the switching surface under which a hybrid Hamiltonian system is, in a certain sense, completely integrable. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: 6 pages, preprint submitted to a conference, comments are welcome

MSC Class: 37J35; 37J39; 70F35; 70H06; 93C30

Parametric Frequency Divider based Ising Machines

Authors: Nicolas Casilli, Tahmid Kaisar, Luca Colombo, Siddhartha Ghosh, Philip X. -L. Feng, Cristian Cassella

Abstract: We report on a new class of Ising Machines (IMs) that rely on coupled parametric frequency dividers (PFDs) as macroscopic artificial spins. Unlike the IM counterparts based on subharmonic injection locking (SHIL), PFD IMs do not require strong injected continuous wave signals or applied DC voltages. Therefore, they show a significantly lower power consumption per spin compared to SHIL based IMs, m… ▽ More We report on a new class of Ising Machines (IMs) that rely on coupled parametric frequency dividers (PFDs) as macroscopic artificial spins. Unlike the IM counterparts based on subharmonic injection locking (SHIL), PFD IMs do not require strong injected continuous wave signals or applied DC voltages. Therefore, they show a significantly lower power consumption per spin compared to SHIL based IMs, making it feasible to accurately solve large scale combinatorial optimization problems (COPs) that are hard or even impossible to solve by using the current von Neumann computing architectures. Furthermore, using high quality (Q) factor resonators in the PFD design makes PFD IMs able to exhibit a nanoWatt level power per spin. Also, it remarkably allows a speed up of the phase synchronization among the PFDs, resulting in shorter time to solution and lower energy to solution despite the resonators' longer relaxation time. As a proof of concept, a 4 node PFD IM has been demonstrated. This IM correctly solves a set of MaxCut problems while consuming just 600 nanoWatts per spin. This power consumption is two orders of magnitude lower than the power per spin of state of the art SHIL based IMs operating at the same frequency. △ Less

Submitted 26 February, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

An energy-momentum method for ordinary differential equations with an underlying $k$-polysymplectic manifold

Authors: Leonardo Colombo, Javier de Lucas, Xavier Rivas, Bartosz M. Zawora

Abstract: This work presents a comprehensive review of the $k$-polysymplectic Marsden-Weinstein reduction theory, rectifying prior errors and inaccuracies in the literature while introducing novel findings. It also emphasises the genuine practical significance of seemingly minor technical details. On this basis, we introduce a novel $k$-polysymplectic energy-momentum method, new related stability analysis t… ▽ More This work presents a comprehensive review of the $k$-polysymplectic Marsden-Weinstein reduction theory, rectifying prior errors and inaccuracies in the literature while introducing novel findings. It also emphasises the genuine practical significance of seemingly minor technical details. On this basis, we introduce a novel $k$-polysymplectic energy-momentum method, new related stability analysis techniques, and apply them to Hamiltonian systems of ordinary differential equations relative to a $k$-polysymplectic manifold. We provide detailed examples of both physical and mathematical significance, including the study of complex Schwarz equations related to the Schwarz derivative, a series of isotropic oscillators, integrable Hamiltonian systems, quantum oscillators with dissipation, affine systems of differential equations, and polynomial dynamical systems. △ Less

Submitted 25 November, 2023; originally announced November 2023.

Comments: 40 pages

MSC Class: 34A26; 34D20; 37J39 (primary); 53B50; 53C15 (secondary)

Nonholonomic reduction for mechanical systems with collisions

Authors: Álvaro Rodríguez Abella, Leonardo J. Colombo

Abstract: This paper studies nonsmooth variational problems on principal bundles for nonholonomic systems with collisions taking place in the boundary of the manifold configuration space of the nonholonopmic system. In particular, we first extended to a nonsmooth context appropriate for collisions the variational principle for nonholonomic implicit Lagrangian systems, to obtain implicit Lagrange--d'Alembert… ▽ More This paper studies nonsmooth variational problems on principal bundles for nonholonomic systems with collisions taking place in the boundary of the manifold configuration space of the nonholonopmic system. In particular, we first extended to a nonsmooth context appropriate for collisions the variational principle for nonholonomic implicit Lagrangian systems, to obtain implicit Lagrange--d'Alembert--Pontryagin equations for nonholonomic systems with collisions, and after introducing the notion of connection on a principal bundle we consider Lagrange--Poincaré--Pointryagin reduction by symmetries for systems with collisions. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: arXiv admin note: text overlap with arXiv:2311.02614

MSC Class: 49J52; 49S05; 53B05; 53C05

Measuring the CMB primordial B-modes with Bolometric Interferometry

Authors: A. Mennella, P. Ade, A. Almela, G. Amico, L. H. Arnaldi, J. Aumont, S. Banfi, E. S. Battistelli, B. Bélier, L. Bergé, J. -Ph. Bernard, P. de Bernardis, M. Bersanelli, J. Bonaparte, J. D. Bonilla, E. Bunn, D. Buzi, F. Cacciotti, D. Camilieri, F. Cavaliere, P. Chanial, C. Chapron, L. Colombo, F. Columbro, A. Coppolecchia , et al. (89 additional authors not shown)

Abstract: The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest fo… ▽ More The Q&U Bolometric Interferometer for Cosmology (QUBIC) is the first bolometric interferometer designed to measure the primordial B-mode polarization of the Cosmic Microwave Background (CMB). Bolometric interferometry is a novel technique that combines the sensitivity of bolometric detectors with the control of systematic effects that is typical of interferometry, both key features in the quest for the faint signal of the primordial B-modes. A unique feature is the so-called "spectral imaging", i.e., the ability to recover the sky signal in several sub-bands within the physical band during data analysis. This feature provides an in-band spectral resolution of Δν/ν \sim 0.04 that is unattainable by a traditional imager. This is a key tool for controlling the Galactic foregrounds contamination. In this paper, we describe the principles of bolometric interferometry, the current status of the QUBIC experiment and future prospects. △ Less

Submitted 5 November, 2023; originally announced November 2023.

Comments: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Implicit Nonholonomic Mechanics with Collisions

Authors: Álvaro Rodríguez Abella, Leonardo Colombo

Abstract: In this paper, variational techniques are used to analyze the dynamics of nonholonomic mechanical systems with impacts. Implicit nonholonomic smooth Lagrangian and Hamiltonian systems are extended to a nonsmooth context appropriate for collisions. In particular, we provide a variational formulation for implicit nonholonomic mechanical systems with collisions, for those collisions that preserve ene… ▽ More In this paper, variational techniques are used to analyze the dynamics of nonholonomic mechanical systems with impacts. Implicit nonholonomic smooth Lagrangian and Hamiltonian systems are extended to a nonsmooth context appropriate for collisions. In particular, we provide a variational formulation for implicit nonholonomic mechanical systems with collisions, for those collisions that preserve energy and momentum at the impact. Lastly, the theoretical results are illustrated by examining the example of a rolling disk hitting a wall. △ Less

Submitted 5 November, 2023; originally announced November 2023.

Optimal Control with Obstacle Avoidance for Incompressible Ideal Flows of an Inviscid Fluid

Authors: Alexandre Anahory Simoes, Anthony Bloch, Leonardo Colombo

Abstract: It has been shown in previous works that an optimal control formulation for an incompressible ideal fluid flow yields Euler's fluid equations. In this paper we consider the modified Euler's equations by adding a potential function playing the role of a barrier function in the corresponding optimal control problem with the motivation of studying obstacle avoidance in the motion of fluid particles f… ▽ More It has been shown in previous works that an optimal control formulation for an incompressible ideal fluid flow yields Euler's fluid equations. In this paper we consider the modified Euler's equations by adding a potential function playing the role of a barrier function in the corresponding optimal control problem with the motivation of studying obstacle avoidance in the motion of fluid particles for incompressible ideal flows of an inviscid fluid From the physical point of view, imposing an artificial potential in the fluid context is equivalent to generating a desired pressure. Simulation results for the obstacle avoidance task are provided. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: 6 pages, conference

Reduction of Necessary Conditions for the Variational Collision Avoidance Problem

Authors: Jacob R. Goodman, Leonardo J. Colombo

Abstract: In this work, we study the reduction by a Lie group of symmetries of variational collision avoidance probelms of multiple agents evolving on a Riemannian manifold and derive necessary conditions for the reduced extremals. The problem consists of finding non-intersecting trajectories of a given number of agents, among a set of admissible curves, to reach a specified configuration, based on minimizi… ▽ More In this work, we study the reduction by a Lie group of symmetries of variational collision avoidance probelms of multiple agents evolving on a Riemannian manifold and derive necessary conditions for the reduced extremals. The problem consists of finding non-intersecting trajectories of a given number of agents, among a set of admissible curves, to reach a specified configuration, based on minimizing an energy functional that depends on the velocity, covariant acceleration and an artificial potential function used to prevent collision among the agents. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 7 pages, 8th IFAC Workshop on Lagrangian and Hamiltonian Methods for Nonlinear Control, LHMNC 2024. arXiv admin note: text overlap with arXiv:2310.18057, arXiv:2207.13574

Impact of Hydrogenation on the Stability and Mechanical Properties of Amorphous Boron Nitride

Authors: Onurcan Kaya, Luigi Colombo, Aleandro Antidormi, Marco A. Villena, Mario Lanza, Ivan Cole, Stephan Roche

Abstract: Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g/cm3), high thermal stability, and mechanical properties. The excellent properties o… ▽ More Interconnect materials with ultralow dielectric constant, and good thermal and mechanical properties are crucial for the further miniaturization of electronic devices. Recently, it has been demonstrated that ultrathin amorphous boron nitride (aBN) films have a very low dielectric constant, high density (above 2.1 g/cm3), high thermal stability, and mechanical properties. The excellent properties of aBN derive from the nature and degree of disorder, which can be controlled at fabrication, allowing tuning of the physical properties for desired applications. Here, we report an improvement in the stability and mechanical properties of amorphous boron nitride upon hydrogen doping. With the introduction of a Gaussian approximation potential (GAP) for atomistic simulations, we investigate the changing morphology of amorphous boron nitride with varying H doping concentrations. We found that for 8 at% of H doping, the concentration of $sp^3$-hybridized atoms reaches a maximum which leads to an improvement of thermal stability and mechanical properties by 20%. These results will be a guideline for experimentalists and process engineers to tune the growth conditions of amorphous boron nitride films for numerous applications. △ Less

Submitted 22 December, 2023; v1 submitted 27 October, 2023; originally announced October 2023.

Comments: 15 pages, 7 figures

Reduction of Sufficient Conditions in Variational Obstacle Avoidance Problems

Authors: Jacob R. Goodman, Leonardo J. Colombo

Abstract: This paper studies sufficient conditions in a variational obstacle avoidance problem on complete Riemannian manifolds. That is, we minimize an action functional, among a set of admissible curves, which depends on an artificial potential function used to avoid obstacles. We provide necessary and sufficient conditions under which the resulting critical points, the so-called modified Riemannian cubic… ▽ More This paper studies sufficient conditions in a variational obstacle avoidance problem on complete Riemannian manifolds. That is, we minimize an action functional, among a set of admissible curves, which depends on an artificial potential function used to avoid obstacles. We provide necessary and sufficient conditions under which the resulting critical points, the so-called modified Riemannian cubics, are local minimizers. We then study the theory of reduction by symmetries of sufficient conditions for optimality in variational obstacle avoidance problems on Lie groups endowed with a left-invariant metric. This amounts to left-translating the Bi-Jacobi fields described to the Lie algebra, and studying the corresponding bi-conjugate points. New conditions are provided in terms of the invertibility of a certain matrix. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 9 pages, 8th IFAC Workshop on Lagrangian and Hamiltonian Methods for Nonlinear Control. arXiv admin note: text overlap with arXiv:2207.13574

The European Low Frequency Survey

Authors: Aniello Mennella, Kam Arnold, Susanna Azzoni, Carlo Baccigalupi, Anthony Banday, R. Belen Barreiro, Darcy Barron, Marco Bersanelli, Sean Casey, Loris Colombo, Elena de la Hoz, Cristian Franceschet, Michael E. Jones, Ricardo T. Genova-Santos, Roger J. Hoyland, Adrian T. Lee, Enrique Martinez-Gonzalez, Filippo Montonati, Jose-Alberto Rubino-Martin, Angela Taylor, Patricio Vielva

Abstract: In this paper we present the European Low Frequency Survey (ELFS), a project that will enable foregrounds-free measurements of primordial $B$-mode polarization to a level 10$^{-3}$ by measuring the Galactic and extra-Galactic emissions in the 5--120\,GHz frequency window. Indeed, the main difficulty in measuring the B-mode polarization comes not just from its sheer faintness, but from the fact tha… ▽ More In this paper we present the European Low Frequency Survey (ELFS), a project that will enable foregrounds-free measurements of primordial $B$-mode polarization to a level 10$^{-3}$ by measuring the Galactic and extra-Galactic emissions in the 5--120\,GHz frequency window. Indeed, the main difficulty in measuring the B-mode polarization comes not just from its sheer faintness, but from the fact that many other objects in the Universe also emit polarized microwaves, which mask the faint CMB signal. The first stage of this project will be carried out in synergy with the Simons Array (SA) collaboration, installing a 5.5--11 GHz coherent receiver at the focus of one of the three 3.5\,m SA telescopes in Atacama, Chile ("ELFS on SA"). The receiver will be equipped with a fully digital back-end based on the latest Xilinx RF System-on-Chip devices that will provide frequency resolution of 1\,MHz across the whole observing band, allowing us to clean the scientific signal from unwanted radio frequency interference, particularly from low-Earth orbit satellite mega-constellations. This paper reviews the scientific motivation for ELFS and its instrumental characteristics, and provides an update on the development of ELFS on SA. △ Less

Submitted 22 November, 2023; v1 submitted 25 October, 2023; originally announced October 2023.

Comments: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

On the Geometry of Virtual Nonlinear Nonholonomic Constraints

Authors: Efstratios Stratoglou, Alexandre Anahory Simoes, Anthony Bloch, Leonardo J. Colombo

Abstract: Virtual constraints are relations imposed on a control system that become invariant via feedback control, as opposed to physical constraints acting on the system. Nonholonomic systems are mechanical systems with non-integrable constraints on the velocities. In this work, we introduce the notion of virtual nonlinear nonholonomic constraints in a geometric framework which is a controlled invariant s… ▽ More Virtual constraints are relations imposed on a control system that become invariant via feedback control, as opposed to physical constraints acting on the system. Nonholonomic systems are mechanical systems with non-integrable constraints on the velocities. In this work, we introduce the notion of virtual nonlinear nonholonomic constraints in a geometric framework which is a controlled invariant submanifold and we show the existence and uniqueness of a control law preserving this submanifold. We illustrate the theory with various examples and present simulation results for an application. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 8 pages

Fundamental Antisymmetric Mode Acoustic Resonator in Periodically Poled Piezoelectric Film Lithium Niobate

Authors: Omar Barrera, Jack Kramer, Ryan Tetro, Sinwoo Cho, Vakhtang Chulukhadze, Luca Colombo, Ruochen Lu

Abstract: Radio frequency (RF) acoustic resonators have long been used for signal processing and sensing. Devices that integrate acoustic resonators benefit from their slow phase velocity (vp), in the order of 3 to 10 km/s, which allows miniaturization of the device. Regarding the subject of small form factor, acoustic resonators that operate at the so-called fundamental antisymmetric mode (A0), feature eve… ▽ More Radio frequency (RF) acoustic resonators have long been used for signal processing and sensing. Devices that integrate acoustic resonators benefit from their slow phase velocity (vp), in the order of 3 to 10 km/s, which allows miniaturization of the device. Regarding the subject of small form factor, acoustic resonators that operate at the so-called fundamental antisymmetric mode (A0), feature even slower vp (1 to 3 km/s), which allows for smaller devices. This work reports the design and fabrication of A0 mode resonators leveraging the advantages of periodically poled piezoelectricity (P3F) lithium niobate, which includes a pair of piezoelectric layers with opposite polarizations to mitigate the charge cancellation arising from opposite stress of A0 in the top and bottom piezoelectric layers. The fabricated device shows a quality factor (Q) of 800 and an electromechanical coupling (k2) of 3.29, resulting in a high figure of merit (FoM, Q times k2) of 26.3 at the resonant frequency of 294 MHz, demonstrating the first efficient A0 device in P3F platforms. The proposed A0 platform could enable miniature signal processing, sensing, and ultrasound transducer applications upon optimization. △ Less

Submitted 27 August, 2023; originally announced September 2023.

Comments: 4 pages, 6 figures, accepted by IEEE IUS 2023

Identifying frequency decorrelated dust residuals in B-mode maps by exploiting the spectral capability of bolometric interferometry

Authors: M. Regnier, E. Manzan, J. -Ch Hamilton, A. Mennella, J. Errard, L. Zapelli, S. A. Torchinsky, S. Paradiso, E. Battistelli, P. De Bernardis, L. Colombo, M. De Petris, G. D'Alessandro, B. Garcia, M. Gervasi, S. Masi, L. Mousset, N. Miron Granese, C. O'Sullivan, M. Piat, E. Rasztocky, G. E. Romero, C. G. Scoccola, M. Zannoni

Abstract: Astrophysical polarized foregrounds represent the most critical challenge in Cosmic Microwave Background (CMB) B-mode experiments. Multi-frequency observations can be used to constrain astrophysical foregrounds to isolate the CMB contribution. However, recent observations indicate that foreground emission may be more complex than anticipated. We investigate how the increased spectral resolution… ▽ More Astrophysical polarized foregrounds represent the most critical challenge in Cosmic Microwave Background (CMB) B-mode experiments. Multi-frequency observations can be used to constrain astrophysical foregrounds to isolate the CMB contribution. However, recent observations indicate that foreground emission may be more complex than anticipated. We investigate how the increased spectral resolution provided by band splitting in Bolometric Interferometry (BI) through a technique called spectral imaging can help control the foreground contamination in the case of unaccounted Galactic dust frequency decorrelation along the line-of-sight. We focus on the next generation ground-based CMB experiment CMB-S4, and compare its anticipated sensitivities, frequency and sky coverage with a hypothetical version of the same experiment based on BI. We perform a Monte-Carlo analysis based on parametric component separation methods (FGBuster and Commander) and compute the likelihood on the recovered tensor-to-scalar ratio. The main result of this analysis is that spectral imaging allows us to detect systematic uncertainties on r from frequency decorrelation when this effect is not accounted for in component separation. Conversely, an imager would detect a biased value of r and would be unable to spot the presence of a systematic effect. We find a similar result in the reconstruction of the dust spectral index, where we show that with BI we can measure more precisely the dust spectral index also when frequency decorrelation is present. The in-band frequency resolution provided by BI allows us to identify dust LOS frequency decorrelation residuals where an imager of similar performance would fail. This opens the prospect to exploit this potential in the context of future CMB polarization experiments that will be challenged by complex foregrounds in their quest for B-modes detection. △ Less

Submitted 21 February, 2024; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: 13 Pages, 15 figures, 4 tables. Accepted by A&A

Journal ref: A&A 686, A271 (2024)

Understanding coupled mass-heat transport in fluids by approach-to-equilibrium molecular dynamics

Authors: Antonio Cappai, Luciano Colombo, Claudio Melis

Abstract: We present a generalization of AEMD approach, routinely applied to estimate thermal conductivity, to the more general case in which Soret and Dufour effects determine a coupled heat-mass transfer. We show that, by starting from microscopical definitions of heat and mass currents, conservation laws dictates the form of the differential equations governing the time evolution. In particular, we focus… ▽ More We present a generalization of AEMD approach, routinely applied to estimate thermal conductivity, to the more general case in which Soret and Dufour effects determine a coupled heat-mass transfer. We show that, by starting from microscopical definitions of heat and mass currents, conservation laws dictates the form of the differential equations governing the time evolution. In particular, we focus to the well specific case in which a closed-form solution of the system is possible and derive the analytical form of time-evolution of temperature and concentration scalar fields in the case in which step-like initial profiles are imposed across a rectangular simulation cell. The validity of this new generalized expression is finally validated using as benchamrk system a two-component Lennard-Jones liquid system, for which generalized diffusivities are estimated in different reduced temperature and density region of phase diagram. △ Less

Submitted 30 August, 2023; originally announced August 2023.

Contact formalism for dissipative mechanical systems on Lie algebroids

Authors: Alexandre Anahory Simoes, Leonardo Colombo, Manuel de Leon, Modesto Salgado, Silvia Souto

Abstract: In this paper, we introduce a geometric description of contact Lagrangian and Hamiltonian systems on Lie algebroids in the framework of contact geometry, using the theory of prolongations. We discuss the relation between Lagrangian and Hamiltonian settings through a convenient notion of Legendre transformation. We also discuss the Hamilton-Jacobi problem in this framework and introduce the notion… ▽ More In this paper, we introduce a geometric description of contact Lagrangian and Hamiltonian systems on Lie algebroids in the framework of contact geometry, using the theory of prolongations. We discuss the relation between Lagrangian and Hamiltonian settings through a convenient notion of Legendre transformation. We also discuss the Hamilton-Jacobi problem in this framework and introduce the notion of a Legendrian Lie subalgebroid of a contact Lie algebroid. △ Less

Submitted 2 August, 2023; originally announced August 2023.

Comments: 36 pages. arXiv admin note: text overlap with arXiv:2204.11537 by other authors

MSC Class: 37J55; 53D10; 37C79; 37J37; 70H03; 70H05; 70H20

Cosmoglobe: Towards end-to-end CMB cosmological parameter estimation without likelihood approximations

Authors: J. R. Eskilt, K. Lee, D. J. Watts, V. Anshul, R. Aurlien, A. Basyrov, M. Bersanelli, L. P. L. Colombo, H. K. Eriksen, K. S. F. Fornazier, U. Fuskeland, M. Galloway, E. Gjerløw, L. T. Hergt, H. T. Ihle, J. G. S. Lunde, A. Marins, S. K. Nerval, S. Paradiso, F. Rahman, M. San, N. -O. Stutzer, I. K. Wehus

Abstract: We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. T… ▽ More We implement support for a cosmological parameter estimation algorithm as proposed by Racine et al. (2016) in Commander, and quantify its computational efficiency and cost. For a semi-realistic simulation similar to Planck LFI 70 GHz, we find that the computational cost of producing one single sample is about 20 CPU-hours and that the typical Markov chain correlation length is $\sim$100 samples. The net effective cost per independent sample is $\sim$2 000 CPU-hours, in comparison with all low-level processing costs of 812 CPU-hours for Planck LFI and WMAP in Cosmoglobe Data Release 1. Thus, although technically possible to run already in its current state, future work should aim to reduce the effective cost per independent sample by one order of magnitude to avoid excessive runtimes, for instance through multi-grid preconditioners and/or derivative-based Markov chain sampling schemes. This work demonstrates the computational feasibility of true Bayesian cosmological parameter estimation with end-to-end error propagation for high-precision CMB experiments without likelihood approximations, but it also highlights the need for additional optimizations before it is ready for full production-level analysis. △ Less

Submitted 31 October, 2023; v1 submitted 27 June, 2023; originally announced June 2023.

Comments: 10 pages, 8 figures. Accepted in A&A

Journal ref: A&A 678, A169 (2023)

Reduction by symmetries of contact mechanical systems on Lie groups

Authors: Alexandre Anahory Simoes, Leonardo Colombo, Manuel de León, Juan Carlos Marrero, David Martín de Diego, Edith Padrón

Abstract: We study the dynamics of contact mechanical systems on Lie groups that are invariant under a Lie group action. Analogously to standard mechanical systems on Lie groups, existing symmetries allow for reducing the number of equations. Thus, we obtain Euler-Poincaré-Herglotz equations on the extended reduced phase space $\mathfrak{g}\times \R$ associated with the extended phase space $TG\times \R$, w… ▽ More We study the dynamics of contact mechanical systems on Lie groups that are invariant under a Lie group action. Analogously to standard mechanical systems on Lie groups, existing symmetries allow for reducing the number of equations. Thus, we obtain Euler-Poincaré-Herglotz equations on the extended reduced phase space $\mathfrak{g}\times \R$ associated with the extended phase space $TG\times \R$, where the configuration manifold $G$ is a Lie group and $\mathfrak{g}$ its Lie algebra. Furthermore, we obtain the Hamiltonian counterpart of these equations by studying the underlying Jacobi structure. Finally, we extend the reduction process to the case of symmetry-breaking systems which are invariant under a Lie subgroup of symmetries. △ Less

Submitted 13 June, 2023; v1 submitted 12 June, 2023; originally announced June 2023.

Comments: 38 pages

Jacobi structure for dissipative mechanical systems on Lie Algebroids

Authors: Alexandre Anahory Simoes, Leonardo Colombo, Manuel de Leon, Modesto Salgado, Silvia Souto

Abstract: We extend the Jacobi structure from $TQ\times \mathbb{R}$ and $T^{*}Q \times \mathbb{R}$ to $A\times \mathbb{R}$ and $A^{*}\times \mathbb{R}$, respectively, where $A$ is a Lie algebroid and $A^{*}$ carries the associated Poisson structure. We see that $A^*\times \mathbb{R}$ possesses a natural Jacobi structure from where we are able to model dissipative mechanical systems, generalizing previous mo… ▽ More We extend the Jacobi structure from $TQ\times \mathbb{R}$ and $T^{*}Q \times \mathbb{R}$ to $A\times \mathbb{R}$ and $A^{*}\times \mathbb{R}$, respectively, where $A$ is a Lie algebroid and $A^{*}$ carries the associated Poisson structure. We see that $A^*\times \mathbb{R}$ possesses a natural Jacobi structure from where we are able to model dissipative mechanical systems, generalizing previous models on $TQ\times \mathbb{R}$ and $\mathfrak{g}\times \mathbb{R}$. △ Less

Submitted 12 June, 2023; v1 submitted 30 May, 2023; originally announced May 2023.

Comments: 16 pages

Cosmoglobe DR1 results. II. Constraints on isotropic cosmic birefringence from reprocessed WMAP and Planck LFI data

Authors: J. R. Eskilt, D. J. Watts, R. Aurlien, A. Basyrov, M. Bersanelli, M. Brilenkov, L. P. L. Colombo, H. K. Eriksen, K. S. F. Fornazier, C. Franceschet, U. Fuskeland, M. Galloway, E. Gjerløw, B. Hensley, L. T. Hergt, D. Herman, H. T. Ihle, K. Lee, J. G. S. Lunde, S. K. Nerval, S. Paradiso, S. K. Patel, F. Rahman, M. Regnier, M. San , et al. (6 additional authors not shown)

Abstract: Cosmic birefringence is a parity-violating effect that might have rotated the plane of linearly polarized light of the cosmic microwave background (CMB) by an angle $β$ since its emission. This has recently been measured to be non-zero at a statistical significance of $3.6σ$ in the official Planck PR4 and 9-year WMAP data. In this work, we constrain $β$ using the reprocessed BeyondPlanck LFI and C… ▽ More Cosmic birefringence is a parity-violating effect that might have rotated the plane of linearly polarized light of the cosmic microwave background (CMB) by an angle $β$ since its emission. This has recently been measured to be non-zero at a statistical significance of $3.6σ$ in the official Planck PR4 and 9-year WMAP data. In this work, we constrain $β$ using the reprocessed BeyondPlanck LFI and Cosmoglobe DR1 WMAP polarization maps. These novel maps have both lower systematic residuals and a more complete error description than the corresponding official products. Foreground $EB$ correlations could bias measurements of $β$, and while thermal dust $EB$ emission has been argued to be statistically non-zero, no evidence for synchrotron $EB$ power has been reported. Unlike the dust-dominated Planck HFI maps, the majority of the LFI and WMAP polarization maps are instead dominated by synchrotron emission. Simultaneously constraining $β$ and the polarization miscalibration angle, $α$, of each channel, we find a best-fit value of $β=0.35^{\circ}\pm0.70^{\circ}$ with LFI and WMAP data only. When including the Planck HFI PR4 maps, but fitting $β$ separately for dust-dominated, $β_{>70\,\mathrm{GHz}}$, and synchrotron-dominated channels, $β_{\leq 70\,\mathrm{GHz}}$, we find $β_{\leq 70\,\mathrm{GHz}}=0.53^{\circ}\pm0.28^\circ$. This differs from zero with a statistical significance of $1.9σ$, and the main contribution to this value comes from the LFI 70 GHz channel. While the statistical significances of these results are low on their own, the measurement derived from the LFI and WMAP synchrotron-dominated maps agrees with the previously reported HFI-dominated constraints, despite the very different astrophysical and instrumental systematics involved in all these experiments. △ Less

Submitted 3 May, 2023; originally announced May 2023.

Comments: 10 pages, 7 figures, 2 tables. Submitted to A&A

Journal ref: A&A 679, A144 (2023)

Hamel equations and quasivelocities for nonholonomic systems with inequality constraints

Authors: Alexandre Anahory Simoes, Leonardo Colombo

Abstract: In this paper we derive Hamel equations for the motion of nonholonomic systems subject to inequality constraints in quasivelocities. As examples, the vertical rolling disk hitting a wall and the Chaplygin sleigh with a knife edge constraint hitting a circular table are shown to illustrate the theoretical results. In this paper we derive Hamel equations for the motion of nonholonomic systems subject to inequality constraints in quasivelocities. As examples, the vertical rolling disk hitting a wall and the Chaplygin sleigh with a knife edge constraint hitting a circular table are shown to illustrate the theoretical results. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Comments: 7 pages

Higher-order retraction maps and construction of numerical methods for optimal control of mechanical systems

Authors: Alexandre Anahory Simoes, Maria Barbero Liñán, Leonardo Colombo, David Martín de Diego

Abstract: Retractions maps are used to define a discretization of the tangent bundle of the configuration manifold as two copies of the configuration manifold where the dynamics take place. Such discretization maps can be conveniently lifted to a higher-order tangent bundle to construct geometric integrators for the higher-order Euler-Lagrange equations. Given a cost function, an optimal control problem for… ▽ More Retractions maps are used to define a discretization of the tangent bundle of the configuration manifold as two copies of the configuration manifold where the dynamics take place. Such discretization maps can be conveniently lifted to a higher-order tangent bundle to construct geometric integrators for the higher-order Euler-Lagrange equations. Given a cost function, an optimal control problem for fully actuated mechanical systems can be understood as a higher-order variational problem. In this paper we introduce the notion of a higher-order discretization map associated with a retraction map to construct geometric integrators for the optimal control of mechanical systems. In particular, we study applications to path planning for obstacle avoidance of a planar rigid body. △ Less

Submitted 31 March, 2023; originally announced March 2023.

Comments: 8 pages

Cosmoglobe DR1 results. I. Improved Wilkinson Microwave Anisotropy Probe maps through Bayesian end-to-end analysis

Authors: D. J. Watts, A. Basyrov, J. R. Eskilt, M. Galloway, L. T. Hergt, D. Herman, H. T. Ihle, S. Paradiso, F. Rahman, H. Thommesen, R. Aurlien, M. Bersanelli, L. A. Bianchi, M. Brilenkov, L. P. L. Colombo, H. K. Eriksen, C. Franceschet, U. Fuskeland, E. Gjerløw, B. Hensley, G. A. Hoerning, K. Lee, J. G. S. Lunde, A. Marins, S. K. Nerval , et al. (8 additional authors not shown)

Abstract: We present Cosmoglobe Data Release 1, which implements the first joint analysis of WMAP and Planck LFI time-ordered data, processed within a single Bayesian end-to-end framework. This framework builds directly on a similar analysis of the LFI measurements by the BeyondPlanck collaboration, and approaches the CMB analysis challenge through Gibbs sampling of a global posterior distribution, simultan… ▽ More We present Cosmoglobe Data Release 1, which implements the first joint analysis of WMAP and Planck LFI time-ordered data, processed within a single Bayesian end-to-end framework. This framework builds directly on a similar analysis of the LFI measurements by the BeyondPlanck collaboration, and approaches the CMB analysis challenge through Gibbs sampling of a global posterior distribution, simultaneously accounting for calibration, mapmaking, and component separation. The computational cost of producing one complete WMAP+LFI Gibbs sample is 812 CPU-hr, of which 603 CPU-hrs are spent on WMAP low-level processing; this demonstrates that end-to-end Bayesian analysis of the WMAP data is computationally feasible. We find that our WMAP posterior mean temperature sky maps and CMB temperature power spectrum are largely consistent with the official WMAP9 results. Perhaps the most notable difference is that our CMB dipole amplitude is $3366.2 \pm 1.4\ \mathrm{μK}$, which is $11\ \mathrm{μK}$ higher than the WMAP9 estimate and $2.5\ σ$ higher than BeyondPlanck; however, it is in perfect agreement with the HFI-dominated Planck PR4 result. In contrast, our WMAP polarization maps differ more notably from the WMAP9 results, and in general exhibit significantly lower large-scale residuals. We attribute this to a better constrained gain and transmission imbalance model. It is particularly noteworthy that the W-band polarization sky map, which was excluded from the official WMAP cosmological analysis, for the first time appears visually consistent with the V-band sky map. Similarly, the long standing discrepancy between the WMAP K-band and LFI 30 GHz maps is finally resolved, and the difference between the two maps appears consistent with instrumental noise at high Galactic latitudes. All maps and the associated code are made publicly available through the Cosmoglobe web page. △ Less

Submitted 14 March, 2023; originally announced March 2023.

Comments: 65 pages, 61 figures. Data available at cosmoglobe.uio.no. Submitted to A&A

Journal ref: A&A 679, A143 (2023)

Liouville-Arnold theorem for contact Hamiltonian systems

Authors: Leonardo Colombo, Manuel de León, Manuel Lainz, Asier López-Gordón

Abstract: A Hamiltonian system is completely integrable (in the sense of Liouville) if there exist as many independent integrals of motion in involution as the dimension of the configuration space. Under certain regularity conditions, Liouville-Arnold theorem states that the invariant geometric structure associated with Liouville integrability is a fibration by Lagrangian tori (or, more generally, Abelian g… ▽ More A Hamiltonian system is completely integrable (in the sense of Liouville) if there exist as many independent integrals of motion in involution as the dimension of the configuration space. Under certain regularity conditions, Liouville-Arnold theorem states that the invariant geometric structure associated with Liouville integrability is a fibration by Lagrangian tori (or, more generally, Abelian groups), on which the motion is linear. In this paper, a Liouville-Arnold theorem for contact Hamiltonian systems is proven. In particular, it is shown that, given a $(2n+1)$-dimensional completely integrable contact system, one can construct a foliation by $(n+1)$-dimensional Abelian groups and induce action-angle coordinates in which the equations of motion are linearized. One important novelty with respect to previous attempts is that the foliation consists of $(n+1)$-dimensional coisotropic submanifolds given by the preimages of rays by the functions in involution. In order to prove the theorem, we first develop a version of Liouville-Arnold theorem for homogeneous functions in exact symplectic manifolds (which is of independent interest), and then apply symplectization to obtain the contact case. △ Less

Submitted 23 February, 2024; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 30 pages. Substantial changes have been made on the proof of Theorem 4 and an example has been added. Preprint submitted to a journal. Comments are welcome!

MSC Class: primary: 37J35; 37J55; 70H06; secondary: 53D10; 53E50; 53Z05

Nonholonomic systems with inequality constraints

Authors: Alexandre Anahory Simoes, Leonardo Colombo

Abstract: In this paper we derive the equations of motion for nonholonomic systems subject to inequality constraints, both, in continuous-time and discrete-time. The last is done by discretizing the continuous time-variational principle which defined the equations of motion for a nonholonomic system subject to inequality constraints. An example is show to illustrate the theoretical results. In this paper we derive the equations of motion for nonholonomic systems subject to inequality constraints, both, in continuous-time and discrete-time. The last is done by discretizing the continuous time-variational principle which defined the equations of motion for a nonholonomic system subject to inequality constraints. An example is show to illustrate the theoretical results. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Virtual Affine Nonholonomic Constraints

Authors: Efstratios Stratoglou, Alexandre Anahory Simoes, Anthony Bloch, Leonardo Colombo

Abstract: Virtual constraints are relations imposed in a control system that become invariant via feedback, instead of real physical constraints acting on the system. Nonholonomic systems are mechanical systems with non-integrable constraints on the velocities. In this work, we introduce the notion of virtual affine nonholonomic constraints in a geometric framework. More precisely, it is a controlled invari… ▽ More Virtual constraints are relations imposed in a control system that become invariant via feedback, instead of real physical constraints acting on the system. Nonholonomic systems are mechanical systems with non-integrable constraints on the velocities. In this work, we introduce the notion of virtual affine nonholonomic constraints in a geometric framework. More precisely, it is a controlled invariant affine distribution associated with an affine connection mechanical control system. We show the existence and uniqueness of a control law defining a virtual affine nonholonomic constraint. △ Less

Submitted 10 January, 2023; originally announced January 2023.

Comments: arXiv admin note: substantial text overlap with arXiv:2207.01299

Hamilton-Jacobi theory for nonholonomic and forced hybrid mechanical systems

Authors: Leonardo Colombo, Manuel de León, María Emma Eyrea Irazú, Asier López-Gordón

Abstract: A hybrid system is a system whose dynamics is given by a mixture of both continuous and discrete transitions. In particular, these systems can be utilised to describe the dynamics of a mechanical system with impacts. Based on the approach by Clark, we develop a geometric Hamilton-Jacobi theory for forced and nonholonomic hybrid dynamical systems. We state the corresponding Hamilton-Jacobi equation… ▽ More A hybrid system is a system whose dynamics is given by a mixture of both continuous and discrete transitions. In particular, these systems can be utilised to describe the dynamics of a mechanical system with impacts. Based on the approach by Clark, we develop a geometric Hamilton-Jacobi theory for forced and nonholonomic hybrid dynamical systems. We state the corresponding Hamilton-Jacobi equations for these classes of systems and apply our results to analyze some examples. △ Less

Submitted 13 May, 2024; v1 submitted 11 November, 2022; originally announced November 2022.

Comments: 31 pages, accepted for publication on Geometric Mechanics

MSC Class: 53Z05; 70F35; 70F40; 70H20; 93C30

Journal ref: Geom. Mech. 01 (02), July 2024

Multiscale modeling of resistive switching in gold nanogranular films

Authors: Miquel López-Suárez, Claudio Melis, Luciano Colombo, Walter Tarantino

Abstract: Metallic nanogranular films display a complex dynamical response to a constant bias, showing up as atypical resistive switching mechanism which could be used to create electrical components for neuromorphic applications. To model such a phenomenon we use a multiscale approach blending together an ab initio treatment of the electric current at the nanoscale, a molecular dynamical approach dictating… ▽ More Metallic nanogranular films display a complex dynamical response to a constant bias, showing up as atypical resistive switching mechanism which could be used to create electrical components for neuromorphic applications. To model such a phenomenon we use a multiscale approach blending together an ab initio treatment of the electric current at the nanoscale, a molecular dynamical approach dictating structural rearrangements, and a finite-element solution of the heat equation for heat propagation in the sample. We also consider structural changes due to electromigration which are modelled on the basis of experimental observations on similar systems. Within such an approach, we manage to describe some distinctive features of the resistive switching occurring in nanogranular film and provide a physical interpretation at the microscopic level. △ Less

Submitted 17 October, 2022; originally announced October 2022.

Discrete Mechanics and Optimal Control for Passive Walking with Foot Slippage

Authors: Alexandre Anahory Simoes, Asier López-Gordón, Anthony Bloch, Leonardo Colombo

Abstract: Forced variational integrators are given by the discretization of the Lagrange-d'Alembert principle for systems subject to external forces, and have proved useful for numerical simulation studies of complex dynamical systems. In this paper we model a passive walker with foot slip by using techniques of geometric mechanics, and we construct forced variational integrators for the system. Moreover, w… ▽ More Forced variational integrators are given by the discretization of the Lagrange-d'Alembert principle for systems subject to external forces, and have proved useful for numerical simulation studies of complex dynamical systems. In this paper we model a passive walker with foot slip by using techniques of geometric mechanics, and we construct forced variational integrators for the system. Moreover, we present a methodology for generating (locally) optimal control policies for simple hybrid holonomically constrained forced Lagrangian systems, based on discrete mechanics, applied to a controlled walker with foot slip in a trajectory tracking problem. △ Less

Submitted 28 September, 2022; originally announced September 2022.

Comments: 8 pages, preprint submitted to a conference. Comments are welcome!

Journal ref: 2023 American Control Conference (ACC), San Diego, CA, USA, 2023, pp. 4587-4592

Contracting Forced Lagrangian and Contact Lagrangian Systems: application to nonholonomic systems with symmetries

Authors: Alexandre Anahory Simoes, Leonardo Colombo

Abstract: In this paper we address the problem of identifying contracting systems among dynamical systems appearing in mechanics. First, we introduce a sufficient condition to identify contracting systems in a general Riemannian manifold. Then, we apply this technique to establish that a particular type of dissipative forced mechanical system is contracting, while stating immediate consequences of this fact… ▽ More In this paper we address the problem of identifying contracting systems among dynamical systems appearing in mechanics. First, we introduce a sufficient condition to identify contracting systems in a general Riemannian manifold. Then, we apply this technique to establish that a particular type of dissipative forced mechanical system is contracting, while stating immediate consequences of this fact for the stability of these systems. Finally, we use the previous results to study the stability of particular types of Contact and Nonholonomic systems. △ Less

Submitted 28 September, 2022; originally announced September 2022.

BeyondPlanck IV. On end-to-end simulations in CMB analysis -- Bayesian versus frequentist statistics

Authors: M. Brilenkov, K. S. F. Fornazier, L. T. Hergt, G. A. Hoerning, A. Marins, T. Murokoshi, F. Rahman, N. -O. Stutzer, Y. Zhou, F. B. Abdalla, K. J. Andersen, R. Aurlien, R. Banerji, A. Basyrov, A. Battista, M. Bersanelli, S. Bertocco, S. Bollanos, L. P. L. Colombo, H. K. Eriksen, J. R. Eskilt, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta , et al. (26 additional authors not shown)

Abstract: End-to-end simulations play a key role in the analysis of any high-sensitivity CMB experiment, providing high-fidelity systematic error propagation capabilities unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization… ▽ More End-to-end simulations play a key role in the analysis of any high-sensitivity CMB experiment, providing high-fidelity systematic error propagation capabilities unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data, or as a random realization independent from the data. We refer to these as Bayesian and frequentist simulations, respectively. We show that the two options lead to significantly different correlation structures, as frequentist simulations, contrary to Bayesian simulations, effectively include cosmic variance, but exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties, and we argue that they therefore also have different and complementary scientific uses, even if this dichotomy is not absolute. Before BeyondPlanck, most pipelines have used a mix of constrained and random inputs, and used the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BeyondPlanck represents the first end-to-end CMB simulation framework that is able to generate both types of simulations, and these new capabilities have brought this topic to the forefront. The Bayesian BeyondPlanck simulations and their uses are described extensively in a suite of companion papers. In this paper we consider one important applications of the corresponding frequentist simulations, namely code validation. That is, we generate a set of 1-year LFI 30 GHz frequentist simulations with known inputs, and use these to validate the core low-level BeyondPlanck algorithms; gain estimation, correlated noise estimation, and mapmaking. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Journal ref: A&A 675, A4 (2023)

BeyondPlanck X. Planck LFI frequency maps with sample-based error propagation

Authors: A. Basyrov, A. -S. Suur-Uski, L. P. L. Colombo, J. R. Eskilt, S. Paradiso, K. J. Andersen, R. Aurlien, R. Banerji, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, H. K. Eriksen, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, E. Gjerløw, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, H. T. Ihle , et al. (15 additional authors not shown)

Abstract: We present Planck LFI frequency sky maps derived within the BeyondPlanck framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncer… ▽ More We present Planck LFI frequency sky maps derived within the BeyondPlanck framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncertainties into higher-level science products. We show that the two dominant sources of LFI instrumental systematic uncertainties are correlated noise and gain fluctuations, and the products presented here support - for the first time - full Bayesian error propagation for these effects at full angular resolution. We compare our posterior mean maps with traditional frequency maps delivered by the Planck collaboration, and find generally good agreement. The most important quality improvement is due to significantly lower calibration uncertainties in the new processing, as we find a fractional absolute calibration uncertainty at 70 GHz of $δg_{0}/g_{0} =5 \cdot 10^{-5}$, which is nominally 40 times smaller than that reported by Planck 2018. However, the original Planck 2018 estimate has a non-trivial statistical interpretation, and this further illustrates the advantage of the new framework in terms of producing self-consistent and well-defined error estimates of all involved quantities without the need of ad hoc uncertainty contributions. We describe how low-resolution data products, including dense pixel-pixel covariance matrices, may be produced directly from the posterior samples without the need for computationally expensive analytic calculations or simulations. We conclude that posterior-based frequency map sampling provides unique capabilities in terms of low-level systematics modelling and error propagation, and may play an important role for future CMB B-mode experiments. (Abridged.) △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 32 pages, 23 figures, data available from https://www.cosmoglobe.uio.no/

Journal ref: A&A 675, A10 (2023)

BeyondPlanck XI. Bayesian CMB analysis with sample-based end-to-end error propagation

Authors: L. P. L. Colombo, J. R. Eskilt, S. Paradiso, H. Thommesen, K. J. Andersen, R. Aurlien, R. Banerji, M. Bersanelli, S. Bertocco, M. Brilenkov, M. Carbone, H. K. Eriksen, M. K. Foss, C. Franceschet, U. Fuskeland, S. Galeotta, M. Galloway, S. Gerakakis, E. Gjerløw, B. Hensley, D. Herman, M. Iacobellis, M. Ieronymaki, H. T. Ihle, J. B. Jewell , et al. (14 additional authors not shown)

Abstract: We present posterior sample-based cosmic microwave background (CMB) constraints from Planck LFI and WMAP observations derived through global end-to-end Bayesian processing. We use these samples to study correlations between CMB, foreground, and instrumental parameters, and we identify a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angul… ▽ More We present posterior sample-based cosmic microwave background (CMB) constraints from Planck LFI and WMAP observations derived through global end-to-end Bayesian processing. We use these samples to study correlations between CMB, foreground, and instrumental parameters, and we identify a particularly strong degeneracy between CMB temperature fluctuations and free-free emission on intermediate angular scales, which is mitigated through model reduction, masking, and resampling. We compare our posterior-based CMB results with previous Planck products, and find generally good agreement, but with higher noise due to exclusion of HFI data. We find a best-fit CMB dipole amplitude of $3362.7\pm1.4μK$, in excellent agreement with previous Planck results. The quoted uncertainty is derived directly from the sampled posterior distribution, and does not involve any ad hoc contribution for systematic effects. Similarly, we find a temperature quadrupole amplitude of $σ^{TT}_2=229\pm97μK^2$, in good agreement with previous results in terms of the amplitude, but the uncertainty is an order of magnitude larger than the diagonal Fisher uncertainty. Relatedly, we find lower evidence for a possible alignment between $\ell = 2$ and $\ell = 3$ than previously reported due to a much larger scatter in the individual quadrupole coefficients, caused both by marginalizing over a more complete set of systematic effects, and by our more conservative analysis mask. For higher multipoles, we find that the angular temperature power spectrum is generally in good agreement with both Planck and WMAP. This is the first time the sample-based asymptotically exact Blackwell-Rao estimator has been successfully established for multipoles up to $\ell\le600$, and it now accounts for the majority of the cosmologically important information. Cosmological parameter constraints are presented in a companion paper. (Abriged) △ Less

Submitted 30 August, 2022; originally announced August 2022.

Comments: 26 pages, 24 figures. Submitted to A&A. Part of the BeyondPlanck paper suite

Journal ref: A&A 675, A11 (2023)

Nonsmooth Herglotz variational principle

Authors: Asier López-Gordón, Leonardo Colombo, Manuel de León

Abstract: In this paper, the theory of smooth action-dependent Lagrangian mechanics (also known as contact Lagrangians) is extended to a non-smooth context appropriate for collision problems. In particular, we develop a Herglotz variational principle for non-smooth action-dependent Lagrangians which leads to the preservation of energy and momentum at impacts. By defining appropriately a Legendre transform,… ▽ More In this paper, the theory of smooth action-dependent Lagrangian mechanics (also known as contact Lagrangians) is extended to a non-smooth context appropriate for collision problems. In particular, we develop a Herglotz variational principle for non-smooth action-dependent Lagrangians which leads to the preservation of energy and momentum at impacts. By defining appropriately a Legendre transform, we can obtain the Hamilton equations of motion for the corresponding non-smooth Hamiltonian system. We apply the result to a billiard problem in the presence of dissipation. △ Less

Submitted 28 September, 2022; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 7 pages, 1 figure. Preprint submitted to a conference. Comments welcome!

MSC Class: 53Z05; 70F35; 70F40; 93C30

Journal ref: 2023 American Control Conference (ACC), San Diego, CA, USA, 2023, pp. 3376-3381

Reduction by Symmetry in Obstacle Avoidance Problems on Riemannian Manifolds

Authors: Jacob R. Goodman, Leonardo J. Colombo

Abstract: This paper studies the reduction by symmetry of a variational obstacle avoidance problem. We derive the reduced necessary conditions in the case of Lie groups endowed with a left-invariant metric, and for its corresponding Riemannian homogeneous spaces by considering an alternative variational problem written in terms of a connection on the horizontal bundle of the Lie group. A number of special c… ▽ More This paper studies the reduction by symmetry of a variational obstacle avoidance problem. We derive the reduced necessary conditions in the case of Lie groups endowed with a left-invariant metric, and for its corresponding Riemannian homogeneous spaces by considering an alternative variational problem written in terms of a connection on the horizontal bundle of the Lie group. A number of special cases where the obstacle avoidance potential can be computed explicitly are studied in detail, and these ideas are applied to the obstacle avoidance task for a rigid body evolving on SO$(3)$ and for the unit sphere $S^2$. △ Less

Submitted 27 July, 2022; originally announced July 2022.

Comments: 29 pages

Reduction in optimal control with broken symmetry for collision and obstacle avoidance of multi-agent system on Lie groups

Authors: Efstratios Stratoglou, Alexandre Anahory Simoes, Leonardo J. Colombo

Abstract: We study the reduction by symmetry for optimality conditions in optimal control problems of left-invariant affine multi-agent control systems, with partial symmetry breaking cost functions for continuous-time and discrete-time systems. We recast the optimal control problem as a constrained variational problem with a partial symmetry breaking Lagrangian and obtain the reduced optimality conditions… ▽ More We study the reduction by symmetry for optimality conditions in optimal control problems of left-invariant affine multi-agent control systems, with partial symmetry breaking cost functions for continuous-time and discrete-time systems. We recast the optimal control problem as a constrained variational problem with a partial symmetry breaking Lagrangian and obtain the reduced optimality conditions from a reduced variational principle via symmetry reduction techniques in both settings, continuous-time, and discrete-time. We apply the results to a collision and obstacle avoidance problem for multiple vehicles evolving on $SE(2)$ in the presence of a static obstacle. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 20 pages