Search

Turbulent transport driven by trapped electron modes (TEMs) is believed to drive significant heat and particle transport in quasihelically symmetric stellarators. Two three-dimensionally-shaped magnetic configurations with suppressed trapped-electron-mode (TEM)-driven turbulence were generated through optimization that targeted quasihelical symmetry and the available energy of trapped electrons. Initial equilibria have flux surface shapes with a helically rotating negative triangularity (NT) and positive triangularity (PT). In gyrokinetic simulations, TEMs are suppressed in the reduced-TEM NT and PT configurations, showing that negative triangularity does not have the same beneficial turbulence properties over positive triangularity as seen in tokamaks. Heat fluxes from TEMs are also suppressed. Without temperature gradients and with a strong density gradient, the most unstable modes at low $k_y$ were consistent with toroidal universal instabilities (UIs) in the NT case and slab UIs in the PT case. Nonlinear simulations show that UIs drive substantial heat flux in both the NT and PT configurations. A moderate increase in $\beta$ halves the heat flux in the NT configuration, while suppressing the heat flux in the PT geometry. Based on the present work, future optimizations aimed at reducing electrostatic drift wave-driven turbulent transport will need to consider UIs if $\beta$ is sufficiently small.

We present a method for obtaining a hierarchy of rigorous bounds on the time-evolution of a quantum mechanical system from an arbitrary initial state, systematically generalizing Mandelstam-Tamm-like relations. For any fixed level in the hierarchy, the bounds are tightest after short time-evolution and gradually loosen over time; we present evidence that for any fixed amount of time-evolution, the bounds can be made arbitrarily tight by moving up in the hierarchy. The computational effort to obtain the bounds scales polynomially with the number of degrees of freedom in the system being simulated. We demonstrate the method on both a single anharmonic oscillator and a system of two coupled anharmonic oscillators., Comment: 8 pages with appendices; 2 figures

We present a novel Automatic Target Recognition (ATR) system using open-vocabulary object detection and classification models. A primary advantage of this approach is that target classes can be defined just before runtime by a non-technical end user, using either a few natural language text descriptions of the target, or a few image exemplars, or both. Nuances in the desired targets can be expressed in natural language, which is useful for unique targets with little or no training data. We also implemented a novel combination of several techniques to improve performance, such as leveraging the additional information in the sequence of overlapping frames to perform tubelet identification (i.e., sequential bounding box matching), bounding box re-scoring, and tubelet linking. Additionally, we developed a technique to visualize the aggregate output of many overlapping frames as a mosaic of the area scanned during the aerial surveillance or reconnaissance, and a kernel density estimate (or heatmap) of the detected targets. We initially applied this ATR system to the use case of detecting and clearing unexploded ordinance on airfield runways and we are currently extending our research to other real-world applications., Comment: Accepted to the Thirty-Seventh Annual Conference on Innovative Applications of Artificial Intelligence (IAAI-25)

We describe the cryogenic system being developed for Taurus: a super-pressure balloon-borne microwave polarimeter scheduled to fly in 2027. The Taurus cryogenic system consists of a 660L liquid helium cryostat which achieves a base temperature of <100mK with the help of a capillary-fed superfluid tank and a closed cycle dilution refrigerator. The main tank is supported with fiberglass flexures and is encased in two layers of vapor-cooled shields which allow Taurus to make full use of the extended flight time offered by the super-pressure balloon platform. The Taurus cryostat is projected to hold for over 50 days while weighing under 1000lbs. We present the design, testing, and thermal analysis of the Taurus cryogenic systems.

The vanishing of the SO(N) beta-function of gauged N>4 supergravity has long seemed like an answer looking for a question. Here we conjecture that it is a consequence of a Montonon-Olive style electric/magnetic duality., Comment: Dedicated to the memory of Stanley Deser

Significant advancements have been made in understanding the physics of transition-edge sensors (TESs) over the past decade. However, key questions remain, particularly a detailed understanding of the current-dependent resistance of these detectors when biased within their superconducting transition. We use scanning superconducting quantum interference device (SQUID) microscopy (SSM) to image the local diamagnetic response of aluminum-manganese alloy (Al-Mn) transition-edge sensors (TESs) near their critical temperature of approximately 175 mK. By doing so, we gain insights into how the device dimensions influence TES transition width, which in turn affects device operation and informs optimal device design. Our images reveal that the Al-Mn thin film near the niobium (Nb) leads exhibits an excess diamagnetic response at temperatures several milli-Kelvin (mK) higher than the bulk of the film farther from the contacts. A possible origin of this behavior is a longitudinal proximity effect between the Nb and Al-Mn where the TES acts as a weak link between superconducting leads. We discuss how this effect shapes the temperature dependence of the resistance as the spacing between the leads decreases. This work demonstrates that magnetic imaging with SSM is a powerful tool for local characterization of superconducting detectors., Comment: 5 pages, 4 figures, Applied Superconductivity Conference (ASC 2024) proceedings, submitted to IEEE Transactions on Applied Superconductivity (IEEE TAS)

We derive the macroscopic charge and current densities of a Chern insulator initially occupying its electronic ground state as it responds to a finite-frequency electric field; we use a previously developed formalism based on microscopic polarization and magnetization fields in extended media. In a topologically trivial insulator, our result reduces to the familiar expression for the induced current density in linear response obtained from a Kubo analysis. But for a Chern insulator we find an extra "topological" term involving the (first) Chern number associated with the occupied bands, encoding the quantum anomalous Hall effect in the presence of a frequency-dependent electric field. While an analogous term has been introduced in the "modern theories of polarization and magnetization" for the linear response of finite-sized systems to static electric fields, our expression is valid for bulk Chern insulators in the presence of both static and finite-frequency electric fields, being derived analytically from a microscopic treatment of the electronic degrees of freedom, and can be generalized in a straightforward way to describe the response of a Chern insulator to electromagnetic fields that are not only frequency-dependent but also spatially inhomogeneous.

Model reduction is an active research field to construct low-dimensional surrogate models of high fidelity to accelerate engineering design cycles. In this work, we investigate model reduction for linear structured systems using dominant reachable and observable subspaces. When the training set $-$ containing all possible interpolation points $-$ is large, then these subspaces can be determined by solving many large-scale linear systems. However, for high-fidelity models, this easily becomes computationally intractable. To circumvent this issue, in this work, we propose an active sampling strategy to sample only a few points from the given training set, which can allow us to estimate those subspaces accurately. To this end, we formulate the identification of the subspaces as the solution of the generalized Sylvester equations, guiding us to select the most relevant samples from the training set to achieve our goals. Consequently, we construct solutions of the matrix equations in low-rank forms, which encode subspace information. We extensively discuss computational aspects and efficient usage of the low-rank factors in the process of obtaining reduced-order models. We illustrate the proposed active sampling scheme to obtain reduced-order models via dominant reachable and observable subspaces and present its comparison with the method where all the points from the training set are taken into account. It is shown that the active sample strategy can provide us $17$x speed-up without sacrificing any noticeable accuracy., Comment: 20 pages, 9 figures

Autonomous motorsports aim to replicate the human racecar driver with software and sensors. As in traditional motorsports, Autonomous Racing Vehicles (ARVs) are pushed to their handling limits in multi-agent scenarios at extremely high ($\geq 150mph$) speeds. This Operational Design Domain (ODD) presents unique challenges across the autonomy stack. The Indy Autonomous Challenge (IAC) is an international competition aiming to advance autonomous vehicle development through ARV competitions. While far from challenging what a human racecar driver can do, the IAC is pushing the state of the art by facilitating full-sized ARV competitions. This paper details the MIT-Pitt-RW Team's approach to autonomous racing in the IAC. In this work, we present our modular and fast approach to agent detection, motion planning and controls to create an autonomy stack. We also provide analysis of the performance of the software stack in single and multi-agent scenarios for rapid deployment in a fast-paced competition environment. We also cover what did and did not work when deployed on a physical system the Dallara AV-21 platform and potential improvements to address these shortcomings. Finally, we convey lessons learned and discuss limitations and future directions for improvement., Comment: Published in Journal of Field Robotics

SPIDER is a balloon-borne instrument designed to map the cosmic microwave background at degree-angular scales in the presence of Galactic foregrounds. SPIDER has mapped a large sky area in the Southern Hemisphere using more than 2000 transition-edge sensors (TESs) during two NASA Long Duration Balloon flights above the Antarctic continent. During its first flight in January 2015, SPIDER observed in the 95 GHz and 150 GHz frequency bands, setting constraints on the B-mode signature of primordial gravitational waves. Its second flight in the 2022-23 season added new receivers at 280 GHz, each using an array of TESs coupled to the sky through feedhorns formed from stacks of silicon wafers. These receivers are optimized to produce deep maps of polarized Galactic dust emission over a large sky area, providing a unique data set with lasting value to the field. In this work, we describe the instrument's performance during SPIDER's second flight., Comment: Submitted to SPIE Astronomical Telescopes + Instrumentation 2024, JATIS

Traumatic brain injury (TBI) can cause a range of cognitive and communication challenges that negatively affect social participation in both face-to-face interactions and computer-mediated communication. In particular, individuals with TBI report barriers that limit access to participation on social media platforms. To improve access to and use of social media for users with TBI, we introduce the Social Media Accessibility and Rehabilitation Toolkit (\textbf{SMART-TBI}). The toolkit includes five aids (Writing Aid, Interpretation Aid, Filter Mode, Focus Mode, and Facebook Customization) designed to address the cognitive and communicative needs of individuals with TBI. We asked eight users with moderate-severe TBI and five TBI rehabilitation experts to evaluate each aid. Our findings revealed potential benefits of aids and areas for improvement, including the need for psychological safety, privacy control, and balancing business and accessibility needs; and overall mixed reactions among the participants to AI-based aids., Comment: In Proceedings The 26th International ACM SIGACCESS Conference on Computers and Accessibility (ASSETS'24)

The Simons Observatory (SO) is a cosmic microwave background experiment composed of three 0.42 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT) in the Atacama Desert of Chile. The Large Aperture Telescope Receiver (LATR) was integrated into the LAT in August 2023; however, because mirrors were not yet installed, the LATR optical chain was capped at the 4K stage. In this dark configuration we are able to characterize many elements of the instrument without contributions from atmospheric noise. Here we show this noise is below the required upper limit and its features are well described with a simple noise model. Maps produced using this noise model have properties that are in good agreement with the white noise levels of our dark data. Additionally, we show that our nominal scan strategy has a minimal effect on the noise when compared to the noise when the telescope is stationary

Taurus is a balloon-borne cosmic microwave background (CMB) experiment optimized to map the E-mode polarization and Galactic foregrounds at the largest angular scales ($\ell$ $\lt$ 30) and improve measurements of the optical depth to reionization ($\tau$). This will pave the way for improved measurements of the sum of neutrino masses in combination with high-resolution CMB data while also testing the $\Lambda CDM$ model on large angular scales and providing high-frequency maps of polarized dust foregrounds to the CMB community. These measurements take advantage of the low-loading environment found in the stratosphere and are enabled by NASA's super-pressure balloon platform, which provides access to 70% of the sky with a launch from Wanaka, New Zealand. Here we describe a general overview of Taurus, with an emphasis on the instrument design. Taurus will employ more than 10,000 100 mK transition edge sensor bolometers distributed across two low-frequency (150, 220 GHz) and one high-frequency (280, 350 GHz) dichroic receivers. The liquid helium cryostat housing the detectors and optics is supported by a lightweight gondola. The payload is designed to meet the challenges in mass, power, and thermal control posed by the super-pressure platform. The instrument and scan strategy are optimized for rigorous control of instrumental systematics, enabling high-fidelity linear polarization measurements on the largest angular scales.

We consider a class of Heisenberg all-to-all coupled spin models inspired by neutrino interactions in a supernova with $N$ degrees of freedom. These models are characterized by a coupling matrix that is relatively simple in the sense that there are only a few, relative to $N$, non-trivial eigenvalues, in distinction to the classic Heisenberg spin-glass models, leading to distinct behavior in both the high-temperature and low-temperature regimes. When the momenta of the neutrinos are uniform and random in directions, we can calculate the large-$N$ partition function for the quantum Heisenberg model. In particular, the high-temperature partition function predicts a non-Gaussian density of states, providing interesting counter-examples showing the limits of general theorems on the density of states for quantum spin models. We can repeat the same argument for classical Heisenberg models, also known as rotor models, and we find the high-temperature expansion is completely controlled by the eigenvalues of the coupling matrix, and again predicts non-Gaussian behavior for the density of states as long as the number of eigenvalues does not scale linearly with $N$. Indeed, we derive the amusing fact that these \emph{thermodynamic} partition functions are essentially the generating function for counting permutations in the high-temperature regime. Finally, for the case relevant to neutrinos in a supernova, we identify the low-temperature phase as a unique state with the direction of the momenta of the neutrino dictating its coherent state in flavor-space, a state we dub the "flavor-momentum-locked" state., Comment: 26 pg., 3 app, multiple figs

We introduce a stochastic framework into the open--source Core Imaging Library (CIL) which enables easy development of stochastic algorithms. Five such algorithms from the literature are developed, Stochastic Gradient Descent, Stochastic Average Gradient (-Am\'elior\'e), (Loopless) Stochastic Variance Reduced Gradient. We showcase the functionality of the framework with a comparative study against a deterministic algorithm on a simulated 2D PET dataset, with the use of the open-source Synergistic Image Reconstruction Framework. We observe that stochastic optimisation methods can converge in fewer passes of the data than a standard deterministic algorithm.

We simulate a variety of optical systematics for Taurus, a balloon-borne cosmic microwave background (CMB) polarisation experiment, to assess their impact on large-scale E-mode polarisation measurements and constraints of the optical depth to reionisation {\tau}. We model a one-month flight of Taurus from Wanaka, New Zealand aboard a super-pressure balloon (SPB). We simulate night-time scans of both the CMB and dust foregrounds in the 150GHz band, one of Taurus's four observing bands. We consider a variety of possible systematics that may affect Taurus's observations, including non-gaussian beams, pointing reconstruction error, and half-wave plate (HWP) non-idealities. For each of these, we evaluate the residual power in the difference between maps simulated with and without the systematic, and compare this to the expected signal level corresponding to Taurus's science goals. Our results indicate that most of the HWP-related systematics can be mitigated to be smaller than sample variance by calibrating with Planck's TT spectrum and using an achromatic HWP model, with a preference for five layers of sapphire to ensure good systematic control. However, additional beam characterization will be required to mitigate far-sidelobe pickup from dust on larger scales., Comment: 21 pages, 6 tables, 7 figures, this version as published in JCAP after minor revisions to v1

The emergence of quantum computing technology over the last decade indicates the potential for a transformational impact in the study of quantum mechanical systems. It is natural to presume that such computing technologies would be valuable to large scientific institutions, such as United States national laboratories. However, detailed descriptions of what these institutions would like to use these computers for are limited. To help provide some initial insights into this topic, this report develops detailed use cases of how quantum computing technology could be utilized to enhance a variety of quantum physics research activities at Los Alamos National Laboratory, including quantum magnetic materials, high-temperature superconductivity and nuclear astrophysics simulations. The report discusses how current high-performance computers are used for scientific discovery today and develops detailed descriptions of the types of quantum physics simulations that Los Alamos National Laboratory scientists would like to conduct, if a sufficient computing technology became available. While the report strives to highlight the breadth of potential application areas for quantum computation, this investigation has also indicated that many more use cases exist at Los Alamos National Laboratory, which could be documented in similar detail with sufficient time and effort.

The Simons Observatory is a new ground-based cosmic microwave background experiment, which is currently being commissioned in Chile's Atacama Desert. During its survey, the observatory's small aperture telescopes will map 10% of the sky in bands centered at frequencies ranging from 27 to 280 GHz to constrain cosmic inflation models, and its large aperture telescope will map 40% of the sky in the same bands to constrain cosmological parameters and use weak lensing to study large-scale structure. To achieve these science goals, the Simons Observatory is deploying these telescopes' receivers with 60,000 state-of-the-art superconducting transition-edge sensor bolometers for its first five year survey. Reading out this unprecedented number of cryogenic sensors, however, required the development of a novel readout system. The SMuRF electronics were developed to enable high-density readout of superconducting sensors using cryogenic microwave SQUID multiplexing technology. The commissioning of the SMuRF systems at the Simons Observatory is the largest deployment to date of microwave multiplexing technology for transition-edge sensors. In this paper, we show that a significant fraction of the systems deployed so far to the Simons Observatory's large aperture telescope meet baseline specifications for detector yield and readout noise in this early phase of commissioning., Comment: 10 pages, 5 figures, 1 table. To be presented at SPIE Astronomical Telescopes + Instrumentation 2024

The Simons Observatory will map the temperature and polarization over half of the sky, at millimeter wavelengths in six spectral bands from the Atacama Desert in Chile. These data will provide new insights into the genesis, content, and history of our Universe; the astrophysics of galaxies and galaxy clusters; objects in our solar system; and time-varying astrophysical phenomena. This ambitious new instrument suite, initially comprising three 0.5 m small-aperture telescopes and one 6 m large aperture telescope, is designed using a common combination of new technologies and new implementations to realize an observatory significantly more capable than the previous generation. In this paper, we present the pre-deployment performance of the first mid-frequency "optics tube" which will be fielded on the large aperture telescope with sensitivity to the 90 and 150 GHz spectral bands. This optics tube contains lenses, filters, detectors, and readout components, all of which operate at cryogenic temperatures. It is one of seven that form the core of the large aperture telescope receiver in its initial deployment. We describe this optics tube, including details of comprehensive testing methods, new techniques for beam and passband characterization, and its measured performance. The performance metrics include beams, optical efficiency, passbands, and forecasts for the on-sky performance of the system. We forecast a sensitivity that exceeds the requirements of the large aperture telescope with greater than 30% margin in each spectral band, and predict that the instrument will realize diffraction-limited performance and the expected detector passbands.

We introduce a new family of minimal problems for reconstruction from multiple views. Our primary focus is a novel approach to autocalibration, a long-standing problem in computer vision. Traditional approaches to this problem, such as those based on Kruppa's equations or the modulus constraint, rely explicitly on the knowledge of multiple fundamental matrices or a projective reconstruction. In contrast, we consider a novel formulation involving constraints on image points, the unknown depths of 3D points, and a partially specified calibration matrix $K$. For $2$ and $3$ views, we present a comprehensive taxonomy of minimal autocalibration problems obtained by relaxing some of these constraints. These problems are organized into classes according to the number of views and any assumed prior knowledge of $K$. Within each class, we determine problems with the fewest -- or a relatively small number of -- solutions. From this zoo of problems, we devise three practical solvers. Experiments with synthetic and real data and interfacing our solvers with COLMAP demonstrate that we achieve superior accuracy compared to state-of-the-art calibration methods. The code is available at https://github.com/andreadalcin/MinimalPerspectiveAutocalibration, Comment: 8 pages main paper + 2 pages references + 8 pages supplementary; to be presented at CVPR 2024

The Simons Observatory (SO) is a cosmic microwave background (CMB) survey experiment that includes small-aperture telescopes (SATs) observing from an altitude of 5,200 m in the Atacama Desert in Chile. The SO SATs will cover six spectral bands between 27 and 280 GHz to search for primordial B-modes to a sensitivity of $\sigma(r)=0.002$, with quantified systematic errors well below this value. Each SAT is a self-contained cryogenic telescope with a 35$^\circ$ field of view, 42 cm diameter optical aperture, 40 K half-wave plate, 1 K refractive optics, and $<0.1$ K focal plane that holds $>12,000$ TES detectors. We describe the nominal design of the SATs and present details about the integration and testing for one operating at 93 and 145 GHz.

In this work, we consider the $\mathcal{H}_2$ optimal model reduction of dynamical systems that are linear in the state equation and up to quadratic nonlinearity in the output equation. As our primary theoretical contributions, we derive gradients of the squared $\mathcal{H}_2$ system error with respect to the reduced model quantities and, from the stationary points of these gradients, introduce Gramian-based first-order necessary conditions for the $\mathcal{H}_2$ optimal approximation of a linear quadratic output (LQO) system. The resulting $\mathcal{H}_2$ optimality framework neatly generalizes the analogous Gramian-based optimality framework for purely linear systems. Computationally, we show how to enforce the necessary optimality conditions using Petrov-Galerkin projection; the corresponding projection matrices are obtained from a pair of Sylvester equations. Based on this result, we propose an iteratively corrected algorithm for the $\mathcal{H}_2$ model reduction of LQO systems, which we refer to as LQO-TSIA (linear quadratic output two-sided iteration algorithm). Numerical examples are included to illustrate the effectiveness of the proposed computational method against other existing approaches., Comment: 18 pages, 4 figures

We revisit certain problems of pose estimation based on 3D--2D correspondences between features which may be points or lines. Specifically, we address the two previously-studied minimal problems of estimating camera extrinsics from $p \in \{ 1, 2 \}$ point--point correspondences and $l=3-p$ line--line correspondences. To the best of our knowledge, all of the previously-known practical solutions to these problems required computing the roots of degree $\ge 4$ (univariate) polynomials when $p=2$, or degree $\ge 8$ polynomials when $p=1.$ We describe and implement two elementary solutions which reduce the degrees of the needed polynomials from $4$ to $2$ and from $8$ to $4$, respectively. We show experimentally that the resulting solvers are numerically stable and fast: when compared to the previous state-of-the art, we may obtain nearly an order of magnitude speedup. The code is available at \url{https://github.com/petrhruby97/efficient\_absolute}, Comment: CVPR 2024, 11 pages, 8 figures, 5 tables

We introduce a detection algorithm for SAGBI basis in polynomial rings, analogous to a Gr\"obner basis detection algorithm previously proposed by Gritzmann and Sturmfels. We also present two accompanying software packages named SagbiGbDetection for Macaulay2 and Julia. Both packages allow the user to find one or more term orders for which a set of input polynomials form either Gr\"obner basis for the ideal they generate or a SAGBI basis for the subalgebra. Additionally, we investigate the computational complexity of homogeneous SAGBI detection and apply our implementation to several novel examples.

Using a novel optimization procedure it has been shown that the Helically Symmetric eXperiment (HSX) stellarator can be optimized for reduced trapped-electron-mode (TEM) instability [M.J.~Gerard et al., \textit{Nucl.~Fusion} \textbf{63} (2023) 056004]. Presently, with a set of 563 experimental candidate configurations, gyrokinetic simulations are performed to investigate the efficacy of available energy $E_\mathrm{A}$, quasi-helical symmetry, and flux-surface shaping parameters as metrics for TEM stabilization. It is found that lower values of $E_\mathrm{A}$ correlate with reduced growth rates, but only when separate flux-surface shaping regimes are considered. Moreover, configurations with improved quasi-helical symmetry demonstrate a similar reduction in growth rates and less scatter compared to $E_\mathrm{A}$. Regarding flux-surface shaping, a set of helical shaping parameters is introduced that show increased elongation is strongly correlated with reduced TEM growth rates, however, only when the quasi-helical symmetry is preserved. Using a newly derived velocity-space-averaged TEM resonance operator, these trends are analyzed to provide insights into the physical mechanism of the observed stabilization. For elongation, stabilization is attributed to geometric effects that reduce the destabilizing particle drifts across the magnetic field. Regarding quasi-helical symmetry, the TEM resonance in the maximally resonant trapping well is shown to increase as the quasi-helical symmetry is broken, and breaking quasi-helical symmetry increases the prevalence of highly resonant trapping wells. While these results demonstrate the limitations of using any single metric as a linear TEM proxy, it is shown that quasi-helical symmetry and plasma elongation are highly effective metrics for reducing TEM growth rates in helical equilibria.

The Simons Observatory (SO) is a cosmic microwave background instrumentation suite in the Atacama Desert of Chile. More than 65,000 polarization-sensitive transition-edge sensor (TES) bolometers will be fielded in the frequency range spanning 27 to 280 GHz, with three separate dichroic designs. The mid-frequency 90/150 GHz and ultra-high-frequency 220/280 GHz detector arrays, fabricated at NIST, account for 39 of 49 total detector modules and implement the feedhorn-fed orthomode transducer (OMT)-coupled TES bolometer architecture. A robust production-level fabrication framework for these detector arrays and the monolithic DC/RF routing wafers has been developed, which includes single device prototyping, process monitoring techniques, in-process metrology, and cryogenic measurements of critical film properties. Application of this framework has resulted in timely delivery of nearly 100 total superconducting focal plane components to SO with 88% of detector wafers meeting nominal criteria for integration into a detector module: a channel yield > 95% and Tc in the targeted range., Comment: 10 pages, 6 figures. Proceedings of the 20th International Conference on Low Temperature Detectors (LTD20). Submitted to JLTP

This work primarily focuses on an operator inference methodology aimed at constructing low-dimensional dynamical models based on a priori hypotheses about their structure, often informed by established physics or expert insights. Stability is a fundamental attribute of dynamical systems, yet it is not always assured in models derived through inference. Our main objective is to develop a method that facilitates the inference of quadratic control dynamical systems with inherent stability guarantees. To this aim, we investigate the stability characteristics of control systems with energy-preserving nonlinearities, thereby identifying conditions under which such systems are bounded-input bounded-state stable. These insights are subsequently applied to the learning process, yielding inferred models that are inherently stable by design. The efficacy of our proposed framework is demonstrated through a couple of numerical examples., Comment: 12 pages, 4 figures

This work discusses model reduction for differential-algebraic systems with quadratic output equations. Under mild conditions, these systems can be transformed into a Weierstra{\ss} canonical form and, thus, be decoupled into differential equations and algebraic equations. The corresponding decoupled states are referred to as proper and improper states. Due to the quadratic function of the state as an output, the proper and improper states are coupled in the output equation, which imposes a challenge from a model reduction viewpoint. Keeping the coupling in mind, our goal in this work is to find important subspaces of the proper and improper states and to reduce the system accordingly. To that end, we first propose the system's matrices, the so-called Gramians, to characterize the system's dominant subspaces. We pay particular attention to the computation of the observability Gramians that take into account the nonlinear coupling between the proper and the improper states. We furthermore show that the proposed Gramians are related to certain kernel functions, which are used to identify important subspaces. This allows us to propose a reduction algorithm to obtain reduced-order systems by removing the subspaces that are difficult to reach, as well as, difficult to observe. Moreover, we quantify the error between the full-order and reduced-order models and demonstrate the proposed methodology using three numerical experiments.

We revisit the problem of certifying the correctness of approximate solution paths computed by numerical homotopy continuation methods. We propose a conceptually simple approach based on a parametric variant of the Krawczyk method from interval arithmetic. Unlike most previous methods for certified path-tracking, our approach is applicable in the general setting of parameter homotopies commonly used to solve polynomial systems of equations. We also describe a novel preconditioning strategy and give theoretical correctness and termination results. Experiments using a preliminary implementation of the method indicate that our approach is competitive with specialized methods appearing previously in the literature, in spite of our more general setting., Comment: 19 pages, 4 figures, 4 tables. Accepted for the Proceedings of ISSAC 2024

The main results of this paper establish a partial correspondence between two previously-studied analogues of Groebner bases in the setting of algebras: namely, subalgebra (aka SAGBI) bases for quotients of polynomial rings and Khovanskii bases for valued algebras. We aim to bridge the gap between the concrete, computational aspects of the former and the more abstract theory of the latter. Our philosophy is that most interesting examples of Khovanskii bases can also be realized as subalgebra bases and vice-versa. We also discuss the computation of Newton-Okounkov bodies, illustrating how interpreting Khovanskii bases as subalgebra bases makes them more amenable to the existing computer algebra tools., Comment: 14 pages, 2 Figures

We systematically compile an exhaustive catalogue of multiview varieties and anchored multiview varieties arising from projections of points and lines in 1, 2, and 3-dimensional projective space. We say that two such varieties are ED-equivalent if there is a linear isomorphism between that that preserve ED-critical points. This gives rise to fourteen equivalence classes, and we determine various properties - dimension, set-theoretic equations, and multidegrees - for all varieties featured in our catalogue. In the case of points, we also present a complementary study of resectioning varieties and their singular loci. Finally, we propose conjectures for the Euclidean distance degrees of all varieties appearing in our comprehensive compilation.

Modern platforms leverage randomized experiments to make informed decisions from a given set of items (``treatments''). As a particularly challenging scenario, these items may (i) arrive in high volume, with thousands of new items being released per hour, and (ii) have short lifetime, say, due to the item's transient nature or underlying non-stationarity that impels the platform to perceive the same item as distinct copies over time. Thus motivated, we study a Bayesian multiple-play bandit problem that encapsulates the key features of the multivariate testing (or ``multi-A/B testing'') problem with a high volume of short-lived arms. In each round, a set of $k$ arms arrive, each available for $w$ rounds. Without knowing the mean reward for each arm, the learner selects a multiset of $n$ arms and immediately observes their realized rewards. We aim to minimize the loss due to not knowing the mean rewards, averaged over instances generated from a given prior distribution. We show that when $k = O(n^\rho)$ for some constant $\rho>0$, our proposed policy has $\tilde O(n^{-\min \{\rho, \frac 12 (1+\frac 1w)^{-1}\}})$ loss on a sufficiently large class of prior distributions. We complement this result by showing that every policy suffers $\Omega (n^{-\min \{\rho, \frac 12\}})$ loss on the same class of distributions. We further validate the effectiveness of our policy through a large-scale field experiment on {\em Glance}, a content-card-serving platform that faces exactly the above challenge. A simple variant of our policy outperforms the platform's current recommender by 4.32\% in total duration and 7.48\% in total number of click-throughs.

Galois/monodromy groups attached to parametric systems of polynomial equations provide a method for detecting the existence of symmetries in solution sets. Beyond the question of existence, one would like to compute formulas for these symmetries, towards the eventual goal of solving the systems more efficiently. We describe and implement one possible approach to this task using numerical homotopy continuation and multivariate rational function interpolation. We describe additional methods that detect and exploit a priori unknown quasi-homogeneous structure in symmetries. These methods extend the range of interpolation to larger examples, including applications with nonlinear symmetries drawn from vision and robotics., Comment: Extended journal version of conference paper published at ISSAC 2023

The Simons Observatory (SO) is a cosmic microwave background instrumentation suite being deployed in the Atacama Desert in northern Chile. The telescopes within SO use three types of dichroic transition-edge sensor (TES) detector arrays, with the 90 and 150 GHz Mid-Frequency (MF) arrays containing 65% of the approximately 68,000 detectors in the first phase of SO. All of the 26 required MF detector arrays have now been fabricated, packaged into detector modules, and tested in laboratory cryostats. Across all modules, we find an average operable detector yield of 84% and median saturation powers of (2.8, 8.0) pW with interquartile ranges of (1, 2) pW at (90, 150) GHz, respectively, falling within their targeted ranges. We measure TES normal resistances and superconducting transition temperatures on each detector wafer to be uniform within 3%, with overall central values of 7.5 mohm and 165 mK, respectively. Results on time constants, optical efficiency, and noise performance are also presented and are consistent with achieving instrument sensitivity forecasts., Comment: 8 pages, 3 figures. Proceedings of the 20th International Conference on Low Temperature Detectors (LTD20). Accepted to JLTP