title: WISDOM Project -- X. The morphology of the molecular ISM in galaxy centres and its dependence on galaxy structure

abstract: We use high-resolution maps of the molecular interstellar medium (ISM) in the centres of eighty-six nearby galaxies from the millimetre-Wave Interferometric Survey of Dark Object Masses (WISDOM) and Physics at High Angular Resolution in Nearby GalaxieS (PHANGS) surveys to investigate the physical mechanisms setting the morphology of the ISM at molecular cloud scales. We show that early-type galaxies tend to have smooth, regular molecular gas morphologies, while the ISM in spiral galaxy bulges is much more asymmetric and clumpy when observed at the same spatial scales. We quantify these differences using non-parametric morphology measures (Asymmetry, Smoothness and Gini), and compare these measurements with those extracted from idealised galaxy simulations. We show that the morphology of the molecular ISM changes systematically as a function of various large-scale galaxy parameters, including galaxy morphological type, stellar mass, stellar velocity dispersion, effective stellar mass surface density, molecular gas surface density, star formation efficiency and the presence of a bar. We perform a statistical analysis to determine which of these correlated parameters best predicts the morphology of the ISM. We find the effective stellar mass surface (or volume) density to be the strongest predictor of the morphology of the molecular gas, while star formation and bars maybe be important secondary drivers. We find that gas self-gravity is not the dominant process shaping the morphology of the molecular gas in galaxy centres. Instead effects caused by the depth of the potential well such as shear, suppression of stellar spiral density waves and/or inflow affect the ability of the gas to fragment.


title: HR 6819 is a binary system with no black hole -- revisiting the source with infrared interferometry and optical integral field spectroscopy

abstract: Two scenarios have been proposed to match the existing observational constraints of the object HR 6819. The system could consist of a close inner B-type giant plus a black hole (BH) binary with an additional Be companion in a wide orbit. Alternatively, it could be a binary composed of a stripped B star and a Be star in a close orbit. Either scenario makes HR 6819 a cornerstone object as the stellar BH closest to Earth, or as an example of an important transitional, non-equilibrium phase for Be stars with solid evidence for its nature. We aimed to distinguish between the two scenarios for HR 6819. Both models predict two luminous stars but with very different angular separations and orbital motions. Therefore, the presence of bright sources in the 1-100 milliarcsec (mas) regime is a key diagnostic for determining the nature of the HR 6819 system. We obtained new high-angular resolution data with VLT/MUSE and VLTI/GRAVITY of HR 6819. The MUSE data are sensitive to bright companions at large scales, whilst the interferometric GRAVITY data are sensitive down to separations on mas scales and large magnitude differences. The MUSE observations reveal no bright companion at large separations and the GRAVITY observations indicate the presence of a stellar companion at an angular separation of ~1.2 mas that moves on the plane of the sky over a timescale compatible with the known spectroscopic 40-day period. We conclude that HR 6819 is a binary system and that no BH is present in the system. The unique nature of HR 6819, and its proximity to Earth make it an ideal system for quantitatively characterising the immediate outcome of binary interaction and probing how Be stars form.


title: Neural Galerkin Schemes with Active Learning for High-Dimensional Evolution Equations

abstract: Deep neural networks have been shown to provide accurate function approximations in high dimensions. However, fitting network parameters requires informative training data that are often challenging to collect in science and engineering applications. This work proposes Neural Galerkin schemes based on deep learning that generate training data with active learning for numerically solving high-dimensional partial differential equations. Neural Galerkin schemes build on the Dirac-Frenkel variational principle to train networks by minimizing the residual sequentially over time, which enables adaptively collecting new training data in a self-informed manner that is guided by the dynamics described by the partial differential equations. This is in contrast to other machine learning methods that aim to fit network parameters globally in time without taking into account training data acquisition. Our finding is that the active form of gathering training data of the proposed Neural Galerkin schemes is key for numerically realizing the expressive power of networks in high dimensions. Numerical experiments demonstrate that Neural Galerkin schemes have the potential to enable simulating phenomena and processes with many variables for which traditional and other deep-learning-based solvers fail, especially when features of the solutions evolve locally such as in high-dimensional wave propagation problems and interacting particle systems described by Fokker-Planck and kinetic equations.


title: Role of disorder in the electronic and magnetic properties of Ag$_3$LiIr$_2$O$_6$

abstract: The nature of magnetism in the intercalated honeycomb iridate Ag$_3$LiIr$_2$O$_6$ has been a subject of recent intensive debate, where the absence or presence of antiferromagnetic order has been reported to be related to possible structural disorder effects and, an enhanced Ir-O hybridization and itinerancy with respect to the parent Li$_2$IrO$_3$ has been suggested as the origin of distinct x-ray spectroscopy features. In the present work we investigate the microscopic nature of the electronic and magnetic properties of Ag$_3$LiIr$_2$O$_6$ via a combination of density functional theory combined with exact diagonalization of ab initio-derived models for various experimental and theoretical structures. We evaluate two possible scenarios, the itinerant quasimolecular framework (QMO) on the one hand, and the localized relativistic $j_{\rm eff} = 1/2$ and $j_{\rm eff} = 3/2$ picture on the other hand, and find that the latter description is still viable for this system. We further calculate resonant inelastic x-ray scattering spectra and show that agreement with experimental observations can be obtained if the presence of Ag vacancies leading to changes in Ir filling and structural disorder is assumed. Finally, we show that the experimentally observed antiferromagnetic spiral magnetic order is reproduced by our ab-initio derived magnetic models.


title: Simplified Stability Assessment of Power Systems with Variable-Delay Wide-Area Damping Control

abstract: Power electronic devices such as HVDC and FACTS can be used to improve the damping of poorly damped inter-area modes in large power systems. This involves the use of wide-area feedback signals, which are transmitted via communication networks. The performance of the closed-loop system is strongly influenced by the delay associated with wide-area signals. The random nature of this delay introduces a switched linear system model. The stability assessment of such a system requires linear matrix inequality based approaches. This makes the stability analysis more complicated as the system size increases. To address this challenge, this paper proposes a delay-processing strategy that simplifies the modelling and analysis in discrete-domain. In contrast to the existing stability assessment techniques, the proposed approach is advantageous because the stability, as well as damping performance, can be accurately predicted by a simplified analysis. The proposed methodology is verified with a case study on the 2-area 4-machine power system with a series compensated tie-line. The results are found to be in accordance with the predictions of the proposed simplified analysis.


title: Faking feature importance: A cautionary tale on the use of differentially-private synthetic data

abstract: Synthetic datasets are often presented as a silver-bullet solution to the problem of privacy-preserving data publishing. However, for many applications, synthetic data has been shown to have limited utility when used to train predictive models. One promising potential application of these data is in the exploratory phase of the machine learning workflow, which involves understanding, engineering and selecting features. This phase often involves considerable time, and depends on the availability of data. There would be substantial value in synthetic data that permitted these steps to be carried out while, for example, data access was being negotiated, or with fewer information governance restrictions. This paper presents an empirical analysis of the agreement between the feature importance obtained from raw and from synthetic data, on a range of artificially generated and real-world datasets (where feature importance represents how useful each feature is when predicting a the outcome). We employ two differentially-private methods to produce synthetic data, and apply various utility measures to quantify the agreement in feature importance as this varies with the level of privacy. Our results indicate that synthetic data can sometimes preserve several representations of the ranking of feature importance in simple settings but their performance is not consistent and depends upon a number of factors. Particular caution should be exercised in more nuanced real-world settings, where synthetic data can lead to differences in ranked feature importance that could alter key modelling decisions. This work has important implications for developing synthetic versions of highly sensitive data sets in fields such as finance and healthcare.


title: Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

abstract: Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations.


title: The Dark Matter Halo of M54

abstract: M54 is a prototype of a globular cluster embedded in a dark matter halo. Gaia EDR3 photometry and proper motions separate the old, metal-poor stars from the more metal rich and younger dwarf galaxy stars. The metal poor stars dominate the inner 50 pc, with a velocity dispersion profile that declines to a minimum around 30 pc then rises back to nearly the central velocity dispersion, as expected for a globular cluster at the center of a dark matter halo. The Jeans analysis of the three separate stellar populations give consistent masses that rise approximately linearly with radius to 1 kpc, implying a small core or cuspy halo. These data are compatible with an infalling CDM dark matter halo reduced to 3x10^8 M_sun at the 50 kpc apocenter 2.3 Gyr ago, with a central globular cluster surrounded by the remnant of a dwarf galaxy. Tides gradually remove material beyond 1 kpc but have little effect on the stars and dark matter within 300 pc of the center. M54 appears to be a transitional system between globular clusters with and without local dark halos, whose evolution within the galaxy depends sensitively on the time of accretion and orbital pericenter.


title: Estimating the coronal supra-arcade downflows radio emission: from centimetre through submillimetre wavelengths

abstract: Supra-arcade downflows (SADs) are infrequent, wiggly opaque structures observed to descend through the solar corona, mostly in EUV and soft X-ray frequencies. From their physical characteristics, SADs have been interpreted as voided (subdense) bubbles and are related to magnetic reconnection processes during long-term erupting flares. In this work we use numerical MHD simulations to compute flux density maps, which are convolved with telescope beams to synthesise images with the aim to assess the expected SADs emission at radio wavelengths and propose observing strategies, including the instruments that can be used. We assume that the emission is thermal bremsstrahlung from a fully ionised plasma without any appreciable gyroresonance contribution since magnetic fields are of the order of ~10 G. We find that SADs emission should be optically thin in the frequency [10-1000] GHz range, and the spatially integrated flux should be larger than 1 Jy. We conclude, therefore, that observing SADs in radio frequencies between [0.5-1000] GHz is feasible with present instrumentation. Moreover, since the emission is for the most part optically thin, the flux density is proportional to temperature, density and line-of-sight depth, and when combined with EUV and soft X-ray images, may allow a better density and temperature determination of SADs.


title: Review on contraction analysis and computation of contraction metrics

abstract: Contraction analysis considers the distance between two adjacent trajectories. If this distance is contracting, then trajectories have the same long-term behavior. The main advantage of this analysis is that it is independent of the solutions under consideration. Using an appropriate metric, with respect to which the distance is contracting, one can show convergence to a unique equilibrium or, if attraction only occurs in certain directions, to a periodic orbit. Contraction analysis was originally considered for ordinary differential equations, but has been extended to discrete-time systems, control systems, delay equations and many other types of systems. Moreover, similar techniques can be applied for the estimation of the dimension of attractors and for the estimation of different notions of entropy (including topological entropy). This review attempts to link the references in both the mathematical and the engineering literature and, furthermore, point out the recent developments and algorithms in the computation of contraction metrics.


title: Conditional Reconstruction for Open-set Semantic Segmentation

abstract: Open set segmentation is a relatively new and unexploredtask, with just a handful of methods proposed to model suchtasks.We propose a novel method called CoReSeg thattackles the issue using class conditional reconstruction ofthe input images according to their pixelwise mask. Ourmethod conditions each input pixel to all known classes,expecting higher errors for pixels of unknown classes. Itwas observed that the proposed method produces better se-mantic consistency in its predictions, resulting in cleanersegmentation maps that better fit object boundaries. CoRe-Seg outperforms state-of-the-art methods on the Vaihin-gen and Potsdam ISPRS datasets, while also being com-petitive on the Houston 2018 IEEE GRSS Data Fusiondataset. Official implementation for CoReSeg is availableat:https://github.com/iannunes/CoReSeg.


title: ePA*SE: Edge-based Parallel A* for Slow Evaluations

abstract: Parallel search algorithms harness the multithreading capability of modern processors to achieve faster planning. One such algorithm is PA*SE (Parallel A* for Slow Expansions), which parallelizes state expansions to achieve faster planning in domains where state expansions are slow. In this work, we propose ePA*SE (Edge-based Parallel A* for Slow Evaluations) that improves on PA*SE by parallelizing edge evaluations instead of state expansions. This makes ePA*SE more efficient in domains where edge evaluations are expensive and need varying amounts of computational effort, which is often the case in robotics. On the theoretical front, we show that ePA*SE provides rigorous optimality guarantees. In addition, ePA*SE can be trivially extended to handle an inflation weight on the heuristic resulting in a bounded suboptimal algorithm w-ePA*SE (Weighted ePA*SE) that trades off optimality for faster planning. On the experimental front, we validate the proposed algorithm in two different planning domains: 1) motion planning for 3D humanoid navigation and 2) task and motion planning for a dual-arm robotic assembly task. We show that ePA*SE can be significantly more efficient than PA*SE and other alternatives. The open-source code for ePA*SE along with the baselines is available here: https://github.com/shohinm/parallel_search


title: Experiments in Adaptive Replanning for Fast Autonomous Flight in Forests

abstract: Fast, autonomous flight in unstructured, cluttered environments such as forests is challenging because it requires the robot to compute new plans in realtime on a computationally-constrained platform. In this paper, we enable this capability with a search-based planning framework that adapts sampling density in realtime to find dynamically-feasible plans while remaining computationally tractable. A paramount challenge in search-based planning is that dense obstacles both necessitate large graphs (to guarantee completeness) and reduce the efficiency of graph search (as heuristics become less accurate). To address this, we develop a planning framework with two parts: one that maximizes planner completeness for a given graph size, and a second that dynamically maximizes graph size subject to computational constraints. This framework is enabled by motion planning graphs that are defined by a single parameter, dispersion, which quantifies the maximum trajectory cost to reach an arbitrary state from the graph. We show through real and simulated experiments how the dispersion can be adapted to different environments in realtime, allowing operation in environments with varying density. The simulated experiment demonstrates improved performance over a baseline search-based planning algorithm. We also demonstrate flight speeds of up to 2.5m/s in real-world cluttered pine forests.


title: Micromechanics-based phase field fracture modelling of CNT composites

abstract: We present a novel micromechanics-based phase field approach to model crack initiation and propagation in carbon nanotube (CNT) based composites. The constitutive mechanical and fracture properties of the nanocomposites are first estimated by a mean-field homogenisation approach. Inhomogeneous dispersion of CNTs is accounted for by means of equivalent inclusions representing agglomerated CNTs. Detailed parametric analyses are presented to assess the effect of the main micromechanical properties upon the fracture behaviour of CNT-based composites. The second step of the proposed approach incorporates the previously estimated constitutive properties into a phase field fracture model to simulate crack initiation and growth in CNT-based composites. The modelling capabilities of the framework presented is demonstrated through three paradigmatic case studies involving mode I and mixed mode fracture conditions.


title: A persistent ultraviolet outflow from an accreting neutron star binary transient

abstract: All disc-accreting astrophysical objects produce powerful outflows. In binaries containing neutron stars (NS) or black holes, accretion often takes place during violent outbursts. The main disc wind signatures during these eruptions are blue-shifted X-ray absorption lines, which are preferentially seen in disc-dominated "soft states". By contrast,optical wind-formed lines have recently been detected in "hard states", when a hot corona dominates the luminosity. The relationship between these signatures is unknown, and no erupting system has revealed wind-formed lines between the X-ray and optical bands yet, despite the many strong resonance transitions in this ultraviolet (UV) region. Here, we show that the transient NS binary Swift J1858.6-0814 exhibits wind-formed, blue-shifted absorption associated with C IV, N V and He II in time-resolved UV spectroscopy during a luminous hard state. This represents the first evidence for a warm, moderately ionized outflow component in this state. Simultaneously observed optical lines also display transient blue-shifted absorption. Decomposing the UV data into constant and variable components, the blue-shifted absorption is associated with the former. This implies that the outflow is not connect to the luminous flares in the data. The joint presence of UV and optical wind features reveals a multi-phase and/or stratified outflow from the outer disc. This type of persistent mass loss across all accretion states has been predicted by radiation-hydrodynamic simulations and helps to explain the shorter-than-expected outbursts duration.


title: Privacy-Aware Crowd Labelling for Machine Learning Tasks

abstract: The extensive use of online social media has highlighted the importance of privacy in the digital space. As more scientists analyse the data created in these platforms, privacy concerns have extended to data usage within the academia. Although text analysis is a well documented topic in academic literature with a multitude of applications, ensuring privacy of user-generated content has been overlooked. Most sentiment analysis methods require emotion labels, which can be obtained through crowdsourcing, where non-expert individuals contribute to scientific tasks. The text itself has to be exposed to third parties in order to be labelled. In an effort to reduce the exposure of online users' information, we propose a privacy preserving text labelling method for varying applications, based in crowdsourcing. We transform text with different levels of privacy, and analyse the effectiveness of the transformation with regards to label correlation and consistency. Our results suggest that privacy can be implemented in labelling, retaining the annotational diversity and subjectivity of traditional labelling.


title: Personal Data Visualisation on Mobile Devices: A Systematic Literature Review

abstract: Personal data cover multiple aspects of our daily life and activities, including health, finance, social, Internet, Etc. Personal data visualisations aim to improve the user experience when exploring these large amounts of personal data and potentially provide insights to assist individuals in their decision making and achieving goals. People with different backgrounds, gender and ages usually need to access their data on their mobile devices. Although there are many personal tracking apps, the user experience when using these apps and visualisations is not evaluated yet. There are publications on personal data visualisation in the literature. Still, no systematic literature review investigated the gaps in this area to assist in developing new personal data visualisation techniques focusing on user experience. In this systematic literature review, we considered studies published between 2010 and 2020 in three online databases. We screened 195 studies and identified 29 papers that met our inclusion criteria. Our key findings are various types of personal data, and users have been addressed well in the found papers, including health, sport, diet, Driving habits, lifelogging, productivity, Etc. The user types range from naive users to expert and developers users based on the experiment's target. However, mobile device capabilities and limitations regarding data visualisation tasks have not been well addressed. There are no studies on the best practices of personal data visualisation on mobile devices, assessment frameworks for data visualisation, or design frameworks for personal data visualisations


title: From Digital Media to Empathic Reality: A Systematic Review of Empathy Research in Extended Reality Environments

abstract: Recent advances in extended reality (XR) technologies have enabled new and increasingly realistic empathy tools and experiences. In XR, all interactions take place in different spatial contexts, all with different features, affordances, and constraints. We present a systematic literature survey of recent work on empathy in XR. As a result, we contribute a research roadmap with three future opportunities in XR-enabled empathy research across both physical and virtual spaces.


title: Homogeneous ice nucleation in an ab initio machine learning model of water

abstract: Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine learning model trained on density-functional theory (DFT) energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation) which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.


title: Intrinsic spin-valley locking for conducting electrons in metal-semiconductor-metal lateral hetero-structures of $1H$-transition-metal dichalcogenides

abstract: Lateral-hetero structures of atomic layered materials alter the electronic properties of pristine crystals and provide a possibility to produce useful monolayer materials. We reveal that metal-semiconductor-metal lateral-hetero junctions of $1H$-transition-metal dichalcogenides intrinsically possess conducting channels of electrons with spin-valley locking, e.g., gate electrode. We theoretically investigate the electronic structure and transport properties of the lateral-hetero junctions and show that the hetero-structure produces conducting channels through the $K$ and $K'$ valleys in the semiconducting transition-metal dichalcogenide and restricts the spin of the conducting electrons in each valley due to the valley dependent charge transfer effect. Moreover, the theoretical investigation shows that the hetero-junction of WSe$_2$ realizes a high transmission probability for valley-spin locked electrons even in a long semiconducting region. The hetero-junction also provides a useful electronic transport property, a step-like I-V characteristic.


title: The 2021 outburst of RS Oph: a pictorial atlas of the spectroscopic evolution. II. From day 19 to 102 (solar conjunction)

abstract: We present the second part (days 19 to 102) of our high-resolution Echelle atlas of the spectral evolution of RS Oph during the 2021 nova outburst. Some quick conclusions that can be drawn are: (1) the spectra evolution has progressed smoothly, suggesting that both the expansion of the ejecta through the pre-existing RG wind and the photo-ionization from the central source did not encountered sudden discontinuities prior to switch-off of the nuclear burning; (2) the 2021 spectral evolution is similar to that of the 2006 eruption, with a mirror-appearance of the triple-peaked line profiles (blue peak brighter than red in 2006, the opposite in 2021); (3) the evolution of the triple-peaked profiles progressed similarly through different excitation and ionization degrees, suggesting an origin in a bipolar outflow nested to an equatorial torus; (4) the time-behavior of coronal emission lines appears quite smooth, both in terms of radiated flux as well as evolution of the triple-peaked profiles, with a plateau at maximum lasting up to about day 87; (5) there seems to exist within RS Oph two distinct types of ejecta: one fast moving and producing Gaussian-like line profiles (primarely [NII] and [OIII]), that keep expanding at FWHM=1000 km/s with no sign of ongoing deceleration during the two months leading up to Solar conjunction, and slower moving ejecta giving rise to the triple-peaked line profiles (HeI, HeII, coronal lines, and also Balmer lines at later times), for which the separation in velocity of blue and red peaks is still shrinking at the time of Solar conjunction when it reaches an average of 330 km/s; (6) the "symbiotic band" at 6825 Ang, due to Raman scattering of OVI 1032 by neutral hydrogen, appears simultaneously with the coronal lines: its very presence indicates that a sizeable part of the wind of the red giant has remained neutral during the development of coronal lines.


title: Code Properties of the Holographic Sierpinski Triangle

abstract: We study the holographic quantum error correcting code properties of a Sierpinski Triangle-shaped boundary subregion in $AdS_4/CFT_3$. Due to existing no-go theorems in topological quantum error correction regarding fractal noise, this gives holographic codes a specific advantage over topological codes. We then further argue that a boundary subregion in the shape of the Sierpinski gasket in $AdS_5/CFT_4$ does not possess these holographic quantum error correction properties.


title: Over-Relaxation in Diffusive Integer Lattice Gas

abstract: One of the most striking draw-backs of standard lattice gas methods over lattice Boltzmann methods is a much more limited range of transport parameters that can be achieved. It is common for lattice Boltzmann methods to use over-relaxation to achieve arbitrarily small transport parameters in the hydrodynamic equations. Here, we show that it is possible to implement over-relaxation for integer lattice gases. For simplicity we focus here on lattice gases for the diffusion equation. We demonstrate that adding a flipping operation to lattice gases results in a multi-relaxation time lattice Boltzmann scheme with over-relaxation in the Boltzmann limit.


title: Stable and Semi-stable Sampling Approaches for Continuously Used Samples

abstract: Information retrieval systems are usually measured by labeling the relevance of results corresponding to a sample of user queries. In practical search engines, such measurement needs to be performed continuously, such as daily or weekly. This creates a trade-off between (a) representativeness of query sample to current query traffic of the product; (b) labeling cost: if we keep the same query sample, results would be similar allowing us to reuse their labels; and (c) overfitting caused by continuous usage of same query sample. In this paper we explicitly formulate this tradeoff, propose two new variants -- Stable and Semi-stable -- to simple and weighted random sampling and show that they outperform existing approaches for the continuous usage settings, including monitoring/debugging search engine or comparing ranker candidates.


title: Nemo: Guiding and Contextualizing Weak Supervision for Interactive Data Programming

abstract: Weak Supervision (WS) techniques allow users to efficiently create large training datasets by programmatically labeling data with heuristic sources of supervision. While the success of WS relies heavily on the provided labeling heuristics, the process of how these heuristics are created in practice has remained under-explored. In this work, we formalize the development process of labeling heuristics as an interactive procedure, built around the existing workflow where users draw ideas from a selected set of development data for designing the heuristic sources. With the formalism, we study two core problems of how to strategically select the development data to guide users in efficiently creating informative heuristics, and how to exploit the information within the development process to contextualize and better learn from the resultant heuristics. Building upon two novel methodologies that effectively tackle the respective problems considered, we present Nemo, an end-to-end interactive system that improves the overall productivity of WS learning pipeline by an average 20% (and up to 47% in one task) compared to the prevailing WS approach.


title: Disordered crystals reveal soft quasilocalized glassy excitations

abstract: Structural glasses formed by quenching a melt are known to host a population of low-energy quasilocalized (nonphononic) excitations whose frequencies $\omega$ follow a universal $\sim\!\omega^4$ distribution as $\omega\!\to\!0$, independently of the glass formation history, the interparticle interaction potential or spatial dimension. Here, we show that the universal quartic law of nonphononic excitations also holds in disordered crystals featuring finite long-range order, which is absent in their glassy counterparts. We thus establish that the degree of universality of the quartic law extends beyond structural glasses quenched from a melt. We further find that disordered crystals, whose level of disorder can be continuously controlled, host many more quasilocalized excitations than expected based on their degree of mechanical disorder -- quantified by the relative fluctuations of the shear modulus -- as compared to structural glasses featuring a similar degree of mechanical disorder. Finally, we show that the stability bound on nonlinear quasilocalized excitations -- previously established in structural glasses -- also holds in disordered crystals. Our results are related to glass-like anomalies experimentally observed in disordered crystals. More broadly, they constitute an important step towards tracing the essential ingredients necessary for the emergence of universal nonphononic excitations in disordered solids.

categories: cond-mat.dis-nn cond-mat.mtrl-sci cond-mat.soft cond-mat.stat-mech


title: Descending Price Auctions with Bounded Number of Price Levels and Batched Prophet Inequality

abstract: We consider descending price auctions for selling $m$ units of a good to unit demand i.i.d. buyers where there is an exogenous bound of $k$ on the number of price levels the auction clock can take. The auctioneer's problem is to choose price levels $p_1 > p_2 > \cdots > p_{k}$ for the auction clock such that auction expected revenue is maximized. The prices levels are announced prior to the auction. We reduce this problem to a new variant of prophet inequality, which we call \emph{batched prophet inequality}, where a decision-maker chooses $k$ (decreasing) thresholds and then sequentially collects rewards (up to $m$) that are above the thresholds with ties broken uniformly at random. For the special case of $m=1$ (i.e., selling a single item), we show that the resulting descending auction with $k$ price levels achieves $1- 1/e^k$ of the unrestricted (without the bound of $k$) optimal revenue. That means a descending auction with just 4 price levels can achieve more than 98\% of the optimal revenue. We then extend our results for $m>1$ and provide a closed-form bound on the competitive ratio of our auction as a function of the number of units $m$ and the number of price levels $k$.


title: Spectrally-modified frequency-swept pulses for optically-driven quantum light sources

abstract: We present a driving scheme for solid-state quantum emitters using frequency-swept pulses containing a spectral hole resonant with the optical transition in the emitter. Our scheme enables high-fidelity state inversion, exhibits robustness to variations in the laser pulse parameters and is immune to phonon-mediated excitation-induced dephasing, benefits that derive from the the insensitivity of the adiabaticity condition to variations in the experimental parameters. Our resonant driving approach could be combined with spectral filtering of the scattered pump light and photonic devices for enhanced collection efficiency to realize simultaneous high indistinguishability and brightness in single photon source applications.


title: Exploring Hierarchical Graph Representation for Large-Scale Zero-Shot Image Classification

abstract: The main question we address in this paper is how to scale up visual recognition of unseen classes, also known as zero-shot learning, to tens of thousands of categories as in the ImageNet-21K benchmark. At this scale, especially with many fine-grained categories included in ImageNet-21K, it is critical to learn quality visual semantic representations that are discriminative enough to recognize unseen classes and distinguish them from seen ones. We propose a \emph{H}ierarchical \emph{G}raphical knowledge \emph{R}epresentation framework for the confidence-based classification method, dubbed as HGR-Net. Our experimental results demonstrate that HGR-Net can grasp class inheritance relations by utilizing hierarchical conceptual knowledge. Our method significantly outperformed all existing techniques, boosting the performance by 7\% compared to the runner-up approach on the ImageNet-21K benchmark. We show that HGR-Net is learning-efficient in few-shot scenarios. We also analyzed our method on smaller datasets like ImageNet-21K-P, 2-hops and 3-hops, demonstrating its generalization ability. Our benchmark and code are available at https://kaiyi.me/p/hgrnet.html.


title: A Survey on Offline Reinforcement Learning: Taxonomy, Review, and Open Problems

abstract: With the widespread adoption of deep learning, reinforcement learning (RL) has experienced a dramatic increase in popularity, scaling to previously intractable problems, such as playing complex games from pixel observations, sustaining conversations with humans, and controlling robotic agents. However, there is still a wide range of domains inaccessible to RL due to the high cost and danger of interacting with the environment. Offline RL is a paradigm that learns exclusively from static datasets of previously collected interactions, making it feasible to extract policies from large and diverse training datasets. Effective offline RL algorithms have a much wider range of applications than online RL, being particularly appealing for real-world applications, such as education, healthcare, and robotics. In this work, we contribute with a unifying taxonomy to classify offline RL methods. Furthermore, we provide a comprehensive review of the latest algorithmic breakthroughs in the field using a unified notation as well as a review of existing benchmarks' properties and shortcomings. Additionally, we provide a figure that summarizes the performance of each method and class of methods on different dataset properties, equipping researchers with the tools to decide which type of algorithm is best suited for the problem at hand and identify which classes of algorithms look the most promising. Finally, we provide our perspective on open problems and propose future research directions for this rapidly growing field.


title: Skew-Symmetric Adjacency Matrices for Clustering Directed Graphs

abstract: Cut-based directed graph (digraph) clustering often focuses on finding dense within-cluster or sparse between-cluster connections, similar to cut-based undirected graph clustering methods. In contrast, for flow-based clusterings the edges between clusters tend to be oriented in one direction and have been found in migration data, food webs, and trade data. In this paper we introduce a spectral algorithm for finding flow-based clusterings. The proposed algorithm is based on recent work which uses complex-valued Hermitian matrices to represent digraphs. By establishing an algebraic relationship between a complex-valued Hermitian representation and an associated real-valued, skew-symmetric matrix the proposed algorithm produces clusterings while remaining completely in the real field. Our algorithm uses less memory and asymptotically less computation while provably preserving solution quality. We also show the algorithm can be easily implemented using standard computational building blocks, possesses better numerical properties, and loans itself to a natural interpretation via an objective function relaxation argument.


title: Graph-based Multi-sensor Fusion for Consistent Localization of Autonomous Construction Robots

abstract: Enabling autonomous operation of large-scale construction machines, such as excavators, can bring key benefits for human safety and operational opportunities for applications in dangerous and hazardous environments. To facilitate robot autonomy, robust and accurate state-estimation remains a core component to enable these machines for operation in a diverse set of complex environments. In this work, a method for multi-modal sensor fusion for robot state-estimation and localization is presented, enabling operation of construction robots in real-world scenarios. The proposed approach presents a graph-based prediction-update loop that combines the benefits of filtering and smoothing in order to provide consistent state estimates at high update rate, while maintaining accurate global localization for large-scale earth-moving excavators. Furthermore, the proposed approach enables a flexible integration of asynchronous sensor measurements and provides consistent pose estimates even during phases of sensor dropout. For this purpose, a dual-graph design for switching between two distinct optimization problems is proposed, directly addressing temporary failure and the subsequent return of global position estimates. The proposed approach is implemented on-board two Menzi Muck walking excavators and validated during real-world tests conducted in representative operational environments.


title: Formalizing Elements of Probabilistic Mechanics

abstract: In this paper we create a model of particle motion on a three-dimensional lattice using discrete random walk with small steps. We rigorously construct a probability space of the particle trajectories. Unlike deterministic approach in classical mechanics, here we use probabilistic properties of particle movement to formally derive analogues of Newton's first and second laws of motion. Similar probabilistic models can potentially be applied to justify laws of thermodynamics in a consistent manner.


title: DDL-MVS: Depth Discontinuity Learning for MVS Networks

abstract: Traditional MVS methods have good accuracy but struggle with completeness, while recently developed learning-based multi-view stereo (MVS) techniques have improved completeness except accuracy being compromised. We propose depth discontinuity learning for MVS methods, which further improves accuracy while retaining the completeness of the reconstruction. Our idea is to jointly estimate the depth and boundary maps where the boundary maps are explicitly used for further refinement of the depth maps. We validate our idea and demonstrate that our strategies can be easily integrated into the existing learning-based MVS pipeline where the reconstruction depends on high-quality depth map estimation. Extensive experiments on various datasets show that our method improves reconstruction quality compared to baseline. Experiments also demonstrate that the presented model and strategies have good generalization capabilities. The source code will be available soon.


title: Algebraic hyperbolicity of very general hypersurfaces in products of projective spaces

abstract: We study the algebraic hyperbolicity of very general hypersurfaces in $\mathbb{P}^m \times \mathbb{P}^n$ by using three techniques that build on past work by Ein, Voisin, Pacienza, Coskun and Riedl, and others. As a result, we completely answer the question of whether or not a very general hypersurface of bidegree $(a,b)$ in $\mathbb{P}^m \times \mathbb{P}^n$ is algebraically hyperbolic, except in $\mathbb{P}^3 \times \mathbb{P}^1$ for the bidegrees $(a,b)= (7,3), (6,3)$ and $(5,b)$ with $b\geq 3.$ As another application of these techniques, we improve the known result that very general hypersurfaces in $\mathbb{P}^n$ of degree at least $2n-2$ are algebraically hyperbolic when $n\geq 6$ to $n \geq 5$, leaving $n=4$ as the only open case.


title: The role of band matching in the proximity effect between a superconductor and a normal metal

abstract: We combine density functional theory and scanning tunneling microscopy to study the proximity effects between a bulk Ru(0001) superconductor and an atomically thin overlayer of Co. We have identified that the Co monolayer can grow in two different stackings: the hcp and a reconstructed \textepsilon-like stacking. We analyze their electronic structure from both experiments and density functional theory. While the magnetic hcp stacking shows a weak proximity effect in combination with Shiba states and with almost no suppression of superconductivity of the substrate, the more complex \textepsilon-like stacking becomes almost fully superconducting and displays an edge state at the island rim. We identify this edge state as a trivial state caused by a local hcp rim around the \textepsilon-core. We explain the weak proximity effect between Ru and the magnetic hcp islands by a low transparency of the interface, while the large chemical unit cell of the non-magnetic \textepsilon-like stacking lifts the momentum conservation at the interface making it transparent and causing a clear proximity effect.


title: Detecting Chronic Kidney Disease(CKD) at the Initial Stage: A Novel Hybrid Feature-selection Method and Robust Data Preparation Pipeline for Different ML Techniques

abstract: Chronic Kidney Disease (CKD) has infected almost 800 million people around the world. Around 1.7 million people die each year because of it. Detecting CKD in the initial stage is essential for saving millions of lives. Many researchers have applied distinct Machine Learning (ML) methods to detect CKD at an early stage, but detailed studies are still missing. We present a structured and thorough method for dealing with the complexities of medical data with optimal performance. Besides, this study will assist researchers in producing clear ideas on the medical data preparation pipeline. In this paper, we applied KNN Imputation to impute missing values, Local Outlier Factor to remove outliers, SMOTE to handle data imbalance, K-stratified K-fold Cross-validation to validate the ML models, and a novel hybrid feature selection method to remove redundant features. Applied algorithms in this study are Support Vector Machine, Gaussian Naive Bayes, Decision Tree, Random Forest, Logistic Regression, K-Nearest Neighbor, Gradient Boosting, Adaptive Boosting, and Extreme Gradient Boosting. Finally, the Random Forest can detect CKD with 100% accuracy without any data leakage.


title: Neutrino mass measurements using cryogenic detectors

abstract: The determination of the absolute mass scale of neutrinos is one of the most important challenges in Particle Physics. The shape of the endpoint region of $\beta$-decay and electron capture (EC) spectra depends on the phase space factor, which, in turn, is function of the neutrino mass eigenstates. High energy resolution and high statistics measurements of $\beta$- and EC spectra are therefore considered a model-independent way for the determination of the neutrino mass scale. Since almost four decades, low temperature microcalorimeters are used for the measurement of low energy $\beta$- and EC spectra. The first efforts were focused on the development of large arrays for the measurement of the $^{187}$Re $\beta$-spectrum. In the last ten years, the attention moved to EC of $^{163}$Ho. This choice was mainly motivated by the very good performance which could be achieved with low temperature microcalorimeters enclosing $^{163}$Ho with respect to microcalorimeters with absorber containing $^{187}$Re. The development of low temperature microcalorimeters for the measurement of the finite neutrino mass is discussed and, in particular, the reasons for moving from $^{187}$Re to $^{163}$Ho. The possibility to reach sub-eV sensitivity on the effective electron neutrino mass with $^{163}$Ho, thanks to the multiplexing of large microcalorimeter arrays is demonstrated. In conclusion, an overview on other nuclides which have been proposed as good candidates, motivated by the excellent performance of low temperature microcalorimeters is presented.


title: Good things to do with extra Higgs doublets

abstract: In this contribution to the Snowmass 2021 process, we outline models with two or three Higgs doublets that address open questions of particle physics and cosmology. In particular, we show that with two additional Higgs doublets one can provide a mechanism for the generation of lepton asymmetry and hence baryon asymmetry, through CP violating Higgs decays, near weak scale temperatures. In another model with only one extra Higgs doublet, we illustrate that Yukawa couplings to quarks and neutrinos can lead to a viable mechanism for the generation of Dirac neutrino masses, sourced by the QCD chiral condensate of strange quarks. We adapt Spontaneous Flavor Violation -- a framework for coupling light fermions to new Higgs doublets while avoiding tree level flavor-changing neutral currents -- in constructing these models. In both cases, flavor data provide interesting constraints on the parameter space. Either scenario includes $\mathcal{O}{(1)}$ couplings of light quarks to the Higgs doublets which allow a future 100 TeV $pp$ collider to have reach for the new scalars up to $\mathcal{O}{(10~{\rm TeV})}$ masses, through resonant single production. In the neutrino mass model, collider data can shed light on the mass hierarchy of neutrinos. This article is based on work presented in Refs. [1,2].


title: The Asteroid-Comet Continuum

abstract: The practical distinctions between asteroids and comets, viewed as products of accretion on either side of the snow line, are less clear-cut than previously understood. In this chapter, we discuss the numerous solar system populations which have physical and dynamical properties that conflict with any simple diagnosis of their nature and origin. Studies of these so-called "continuum" or "transition objects", which include many of the most intriguing bodies in the solar system, have implications for a broad range of scientific topics from the demise of comets and the activation of asteroids to the production of interplanetary debris and the origin of the terrestrial planet volatiles. We present an overview of the current state of knowledge concerning the asteroid-comet continuum and discuss the numerous physical processes behind the activity shown by small bodies in the solar system.


title: Structural Design and Impact Analysis of a 1.5U CubeSat on the Lunar Surface

abstract: Ahead of the United States' crewed return to the moon in 2024, Intuitive Machines, under a NASA Commercial Lunar Payload Services contract, will land their Nova-C lunar lander in October 2021. At 30 meters altitude during the terminal descent, EagleCam will be deployed, and will capture and transmit the first-ever third-person images of a spacecraft making an extraterrestrial landing. This paper will focus on the structural design, modeling, and impact analysis of a 1.5U CubeSat payload to withstand a ballistic, soft-touch landing on the lunar surface.


title: Optomechanical simulation of a parametric oscillator

abstract: We study an optomechhanical device supporting at least three optical modes in the infrared telecommunication band and three mechanical vibration modes. We model the coherent driving of each optical mode, independently of each other, to obtain an effective Hamiltonian showing the different types of parametric processes allowed in the device. We propose a bichromatic driving scheme, in the lossy optical cavity regime, under a mean field approximation, that provides the quantum simulation of a parametric oscillator with optical control of its parameters.


title: Adaptive Gradient Methods with Local Guarantees

abstract: Adaptive gradient methods are the method of choice for optimization in machine learning and used to train the largest deep models. In this paper we study the problem of learning a local preconditioner, that can change as the data is changing along the optimization trajectory. We propose an adaptive gradient method that has provable adaptive regret guarantees vs. the best local preconditioner. To derive this guarantee, we prove a new adaptive regret bound in online learning that improves upon previous adaptive online learning methods. We demonstrate the robustness of our method in automatically choosing the optimal learning rate schedule for popular benchmarking tasks in vision and language domains. Without the need to manually tune a learning rate schedule, our method can, in a single run, achieve comparable and stable task accuracy as a fine-tuned optimizer.


title: Optical precursors in waveguide quantum electrodynamics

abstract: When a broadband signal propagates through a dispersive medium, some frequency components move faster than the center of the pulse. This leads to the appearance of precursors, transient signals that emerge from the medium earlier than the main part of the pulse and seem to propagate superluminally. Here, we investigate the microscopic origin of precursors in a minimal setup: an array of qubits coupled to a waveguide. The linear transmission function only converges to that of a continuous medium for large qubit numbers. Nevertheless, the dispersion produced by only two qubits is enough to produce oscillatory transients. Precursors are best observed under conditions of electromagnetically-induced transparency, as the center of the pulse is significantly delayed. Under these conditions, just a single qutrit is enough to generate a precursor. Our results pave the way towards dispersion engineering of light with just a few qubits, and can be realized with superconducting qubits coupled to transmission lines or atoms coupled to optical waveguides.


title: Fixed-point-free pseudo-Anosov homeomorphisms, knot Floer homology and the cinquefoil

abstract: Given any genus-two, hyperbolic, fibered knot in $S^3$ with nonzero fractional Dehn twist coefficient, we show that its pseudo-Anosov representative has a fixed point. Combined with recent work of Baldwin--Hu--Sivek, this proves that knot Floer homology detects the cinquefoil knot $T(2,5)$, and that the cinquefoil is the only genus-two L-space knot in $S^3$. Our results have applications to Floer homology of cyclic branched covers over knots in $S^3$, to $\mathit{SU}(2)$-abelian Dehn surgeries, and to Khovanov and annular Khovanov homology. Along the way to proving our fixed point result, we describe a small list of train tracks carrying all pseudo-Anosov homeomorphisms in most strata on the punctured disk. As a consequence, we find a canonical track $\tau$ carrying all pseudo-Anosov homeomorphisms in a particular stratum $\mathcal{Q}_0$ on the genus-two surface, and describe every fixed-point-free pseudo-Anosov homeomorphism in $\mathcal{Q}_0$.


title: Modeling of an On-Orbit Maintenance Robotic Arm Test-Bed

abstract: This paper focuses on the development of a ground-based test-bed to analyze the complexities of contact dynamics between multibody systems in space. The test-bed consists of an air-bearing platform equipped with a 7 degrees-of-freedom (one degree per revolute joint) robotic arm which acts as the servicing satellite. The dynamics of the manipulator on the platform is modeled as an aid for the analysis and design of stabilizing control algorithms suited for autonomous on-orbit servicing missions. The dynamics are represented analytically using a recursive Newton-Euler multibody method with D-H parameters derived from the physical properties of the arm and platform. In addition, Product of Exponential (PoE) method is also employed to serve as a comparison with the D-H parameters approach. Finally, an independent numerical simulation created with the SimScape modeling environment is also presented as a means of verifying the accuracy of the recursive model and the PoE approach. The results from both models and SimScape are then validated through comparison with internal measurement data taken from the robotic arm itself.


title: Self-Supervised Online Learning for Safety-Critical Control using Stereo Vision

abstract: With the increasing prevalence of complex vision-based sensing methods for use in obstacle identification and state estimation, characterizing environment-dependent measurement errors has become a difficult and essential part of modern robotics. This paper presents a self-supervised learning approach to safety-critical control. In particular, the uncertainty associated with stereo vision is estimated, and adapted online to new visual environments, wherein this estimate is leveraged in a safety-critical controller in a robust fashion. To this end, we propose an algorithm that exploits the structure of stereo-vision to learn an uncertainty estimate without the need for ground-truth data. We then robustify existing Control Barrier Function-based controllers to provide safety in the presence of this uncertainty estimate. We demonstrate the efficacy of our method on a quadrupedal robot in a variety of environments. When not using our method safety is violated. With offline training alone we observe the robot is safe, but overly-conservative. With our online method the quadruped remains safe and conservatism is reduced.


title: Superconducting Phase of $\mathrm{Ti_xO_y}$ Thin Films Grown by Molecular Beam Epitaxy

abstract: We investigate the complex relationship between the growth conditions and the structural and transport properties of $\mathrm{Ti_xO_y}$ thin films grown by molecular beam epitaxy. Transport properties ranging from metallicity to superconductivity and insulating states are stabilized by effectively tuning the O/Ti ratio via the Ti flux rate and the O partial pressure, $P_{Ox}$, for films grown on (0001)-$\mathrm{Al_2 O_3}$ substrates at 850$^{\circ}$ C. A cubic $c-\mathrm{TiO_{1\pm\delta}}$ buffer layer is formed for low O/Ti ratios while a corundum cr-$\mathrm{Ti_2 O_3}$ layer is formed under higher oxidizing conditions. Metallicity is observed for c-$\mathrm{TiO_{1-\delta}}$ buffer layers. The superconducting $\mathrm{\gamma-Ti_3 O_5}$ Magn\'eli phase is found to nucleate on a c-$\mathrm{TiO_{1-\delta}}$ buffer for intermediate $P_{Ox}$ conditions and an insulator-superconducting transition is observed at 4.5 K (T$_C^{onset}=6 K$) for 85 nm thick films. Strain relaxation of the $\mathrm{\gamma-Ti_3 O_5}$ occurs with increasing film thickness and correlates with a thickness-dependent increase in T$_C$ observed for $\mathrm{Ti_xO_y}$ thin films.


title: The MKID Science Data Pipeline

abstract: We present The MKID Pipeline, a general use science data pipeline for the reduction and analysis of ultraviolet, optical and infrared (UVOIR) Microwave Kinetic Inductance Detector (MKID) data sets. This paper provides an introduction to the nature of MKID data sets, an overview of the calibration steps included in the pipeline, and an introduction to the implementation of the software.


title: The Mobile Production Vehicle Routing Problem: Using 3D Printing in Last Mile Distribution

abstract: We study a new variant of the vehicle routing problem, called the Mobile Production Vehicle Routing Problem (MoP-VRP). In this problem, vehicles are equipped with 3D printers, and production takes place on the way to the customer. The objective is to minimize the weighted cost incurred by travel and delay of service. We formulate a Mixed Integer Programming (MIP) model and develop an Adaptive Large Neighbourhood Search (ALNS) heuristic for this problem. To show the advantage of mobile production, we compare the problem with the Central Production Vehicle Routing Problem (CP-VRP), where production takes place in a central depot. We also propose an efficient ALNS for the CP-VRP. We generate benchmark instances based on Vehicle Routing Problem with Time Windows (VRPTW) benchmark instances, and realistic instances based on real-life data provided by the Danish Company 3D Printhuset. Overall, the proposed ALNS for both problems are efficient, and we solve instances up to 200 customers within a short computational time. We test different scenarios with varying numbers of machines in each vehicle, as well as different production time. The results show that these are the key factors that influence travel and delay costs. The key advantage of mobile production is flexibility: it can shorten the time span from the start of production to the delivery of products, and at the same time lower delivery costs. Moreover, long-term cost estimations show that this technology has low operation costs and thus is feasible in real life practice.


title: Further analysis on structure and spectral properties of symmetric graphs

abstract: Graph is an abstract representation commonly used to model networked systems and structure. In problems across various fields, including computer vision and pattern recognition, and neuroscience, graphs are often brought into comparison (a process is called \textit{graph matching}) or checked for symmetry. Friendliness property of the associated adjacency matrices, specified by their spectral properties, is important in deriving a convex relaxation of the (intractable) discrete graph matching problem. In this work, we study unfriendliness properties of symmetric graphs by studying its relation to the underlying graph structure. It is revealed that a symmetric graph has two or more subgraphs of the same topology, and are adjacent to the same set of vertices. We then show that if adjacency matrices of symmetric graphs have distinct eigenvalues then there exist eigenvectors orthogonal to the vector of all ones, making them unfriendly. Relation of graph symmetry to uncontrollability of multi-agent systems under agreement dynamics with one controlled node is revisited. Examples of both synthetic and real-world graphs are also given for illustrations.


title: Design and comparison of two linear controllers with precompensation gain for the Quadruple inverted pendulum

abstract: In this work we present a workflow for designing two linear control techniques applied to the dynamic system quadruple inverted pendulum mounted on a cart (QIP) where the steady state error on cart position is eliminated through a precompensation gain. The first control law designed was based on LQR, technique that stabilizes the system states based on a minimization of the total energy of the system states and the second control technique was pole placement by 2nd order system approximation, this method allowed to define a settling time and overshoot of the desired response of the plant, for this control the system was considered SISO, being the only output of the plant the horizontal position of the cart. To eliminate the steady-state error generated by both linear controllers, an analysis of the stationary response of the system under control in the Laplace Domain was performed. The simulation environment used to validate the results presented in this paper was Simscape of Matlab, it is a simulation environment that allowed to simulate the dynamic of nonlinear mechanical systems. The results obtained indicate that both controllers are able to stabilize the plant and move the cart position to the desired setpoint. Finally, we discuss about the performances of both linear control techniques.


title: Controllability Analysis of Motion of Artificial Satellite Under the Effect of Oblateness of the Earth

abstract: In this article we have studied the controllability of artificial satellite under the effect of zonal harmonic $J_2$ in cylindrical polar coordinates systems. Seven different cases of thrusters in various directions have been analyzed and it is found that the system is controllable if we apply thrusters in either $r$, $\theta$ and $z$ or $\theta$ and $z$ direction. The equations governing motion of satellite have been linearized and Kalman controllability test is applied to check the controllability of the system. We have also derived controller $u$ for the linearized system. The trajectory of the system have been plotted to show the controllability of the system.


title: Multi-system intervention optimization for interdependent infrastructure

abstract: The wellbeing of modern societies is dependent upon the functioning of their infrastructure networks. This paper introduces the 3C concept, an integrative multi-system and multi-stakeholder optimization approach for managing infrastructure interventions (e.g., maintenance, renovation, etc.). The proposed approach takes advantage of the benefits achieved by grouping (i.e., optimizing) intervention activities. Intervention optimization leads to substantial savings on both direct intervention costs (operator) and indirect unavailability costs (society) by reducing the number of system interruptions. The proposed optimization approach is formalized into a structured mathematical model that can account for the interactions between multiple infrastructure networks and the impact on multiple stakeholders (e.g., society and infrastructure operators), and it can accommodate different types of intervention, such as maintenance, removal, and upgrading. The different types of interdependencies, within and across infrastructures, are modeled using a proposed interaction matrix (IM). The IM allows integrating the interventions of different infrastructure networks whose interventions are normally planned independently. Moreover, the introduced 3C concept accounts for central interventions, which are those that must occur at a pre-established time moment, where neither delay nor advance is permitted. To demonstrate the applicability of the proposed approach, an illustrative example of a multi-system and multi-actor intervention planning is introduced. Results show a substantial reduction in the operator and societal costs. In addition, the optimal intervention program obtained in the analysis shows no predictable patterns, which indicates it is a useful managerial decision support tool.


title: Effect of Timing Error: A Case Study of Navigation Camera

abstract: We focus on the problem of timing errors in navigation camera as a case study in a broader problem of the effect of a timing error in cyber-physical systems. These systems rely on the requirement that certain things happen at the same time or certain things happen periodically at some period $T$. However, as these systems get more complex, timing errors can occur between the components thereby violating the assumption about events being simultaneous (or periodic). We consider the problem of a surgical navigation system where optical markers detected in the 2D pictures taken by two cameras are used to localize the markers in 3D space. A predefined array of such markers, known as a reference element, is used to navigate the corresponding CAD model of a surgical instrument on patient's images. The cameras rely on the assumption that the pictures from both cameras are taken exactly at the same time. If a timing error occurs then the instrument may have moved between the pictures. We find that, depending upon the location of the instrument, this can lead to a substantial error in the localization of the instrument. Specifically, we find that if the actual movement is $\delta$ then the observed movement may be as high as $5\delta$ in the operating range of the camera. Furthermore, we also identify potential issues that could affect the error in case there are changes to the camera system or to the operating range.


title: A 915-1220 TOPS/W Hybrid In-Memory Computing based Image Restoration and Region Proposal Integrated Circuit for Neuromorphic Vision Sensors in 65nm CMOS

abstract: In this work, we present a hybrid memory bit cell - collocated SRAM and DRAM (CRAM) consisting of 11 transistors for in-memory computing (IMC) based image restoration (IR) and region proposal (RP). A robust RP updated algorithm is proposed to improve the performance. This work demonstrates IMC based global parallel diffusion and column/row-wise projection to achieve a maximal energy efficiency of 1220 TOPS/W for image restoration and 915 TOPS/W when combined with region proposal.


title: ICARUS: A Specialized Architecture for Neural Radiance Fields Rendering

abstract: The practical deployment of Neural Radiance Fields (NeRF) in rendering applications faces several challenges, with the most critical one being low rendering speed on even high-end graphic processing units (GPUs). In this paper, we present ICARUS, a specialized accelerator architecture tailored for NeRF rendering. Unlike GPUs using general purpose computing and memory architectures for NeRF, ICARUS executes the complete NeRF pipeline using dedicated plenoptic cores (PLCore) consisting of a positional encoding unit (PEU), a multi-layer perceptron (MLP) engine, and a volume rendering unit (VRU). A PLCore takes in positions \& directions and renders the corresponding pixel colors without any intermediate data going off-chip for temporary storage and exchange, which can be time and power consuming. To implement the most expensive component of NeRF, i.e., the MLP, we transform the fully connected operations to approximated reconfigurable multiple constant multiplications (MCMs), where common subexpressions are shared across different multiplications to improve the computation efficiency. We build a prototype ICARUS using Synopsys HAPS-80 S104, a field programmable gate array (FPGA)-based prototyping system for large-scale integrated circuits and systems design. We evaluate the power-performance-area (PPA) of a PLCore using 40nm LP CMOS technology. Working at 400 MHz, a single PLCore occupies 16.5 $mm^2$ and consumes 282.8 mW, translating to 0.105 uJ/sample. The results are compared with those of GPU and tensor processing unit (TPU) implementations.


title: "Keep the fight unfair": Military rhetoric in quantum technology

abstract: Doing quantum ethics properly will require detailed socio-political analysis of the technologies and the organizations trying to build them. In this paper, I contribute to this task by analysing the public rhetoric of American military stakeholders in the quantum industry. I look at Air Force Research Laboratory involvement in the 2020 Quantum 2 Business conference, where they were the main sponsor. A critical thematic analysis shows a focus on enacting the violence of war, maintaining narratives that the Air Force provides a secure future for Americans, and marrying quantum technology with the aesthetics of war. I contextualize this with anti-imperialist theory, arguing that this rhetoric and the desire for quantum arms aligns with the reproduction of existing violent power structures. Insights about this example of military involvement in quantum spaces should help orient nascent critical quantum ethics interventions.


title: A Fully Memristive Spiking Neural Network with Unsupervised Learning

abstract: We present a fully memristive spiking neural network (MSNN) consisting of physically-realizable memristive neurons and memristive synapses to implement an unsupervised Spiking Time Dependent Plasticity (STDP) learning rule. The system is fully memristive in that both neuronal and synaptic dynamics can be realized by using memristors. The neuron is implemented using the SPICE-level memristive integrate-and-fire (MIF) model, which consists of a minimal number of circuit elements necessary to achieve distinct depolarization, hyperpolarization, and repolarization voltage waveforms. The proposed MSNN uniquely implements STDP learning by using cumulative weight changes in memristive synapses from the voltage waveform changes across the synapses, which arise from the presynaptic and postsynaptic spiking voltage signals during the training process. Two types of MSNN architectures are investigated: 1) a biologically plausible memory retrieval system, and 2) a multi-class classification system. Our circuit simulation results verify the MSNN's unsupervised learning efficacy by replicating biological memory retrieval mechanisms, and achieving 97.5% accuracy in a 4-pattern recognition problem in a large scale discriminative MSNN.


title: Two-fold theoretical and experimental misinterpretation of luminosity in the origin of Hubble tension

abstract: The comparison of the intrinsic luminosity of extragalactic objects detected at different redshifts requires the evaluation of the observed flux in a common rest-frame emission band. It is usually performed by adding a term to the observed magnitudes denominated K-correction. In this paper, the K-correction transformation has been carefully scrutinized. As a result of this inspection, it has been found that K-corrections systematically downgrades the measured energy flux by a factor of (1+z). C-correction, a more appropriate rest-frame magnitude transformation has been derived. C-correction has been applied to the re-analysis of JLA (Joint Light-curve analysis) supernovae sample. The best cosmological fit to JLA type Ia supernovae (SNIa) sample calibrated with CMB Planck $h=0.675$ for the common luminosity-angular distances relation $L2:D_L=D_A(1+z)^2$ corresponds to a closed Universe $L2(h_0=0.48,\Omega_M=2.25$). Nevertheless, considering the new luminosity-angular distances relation $L1:D_L=D_A(1+z)$ that accounts for expansion lensing, the new SNIa analysis converges to $L1(h_0=0.7558,\Omega_M=0.035,\Omega_\Lambda=0.965)$, providing a coherent explanation of the Hubble tension between Planck and SH0ES ($h=0.7422 \pm 0.0182$) results.


title: Third-order Analysis of Channel Coding in the Small-to-Moderate Deviations Regime

abstract: This paper studies the third-order characteristic of nonsingular discrete memoryless channels and the Gaussian channel with a maximal power constraint. The third-order term in our expansions employs a new quantity here called the \emph{channel skewness}, which affects the approximation accuracy more significantly as the error probability decreases. For the Gaussian channel, evaluating Shannon's (1959) random coding and sphere-packing bounds in the central limit theorem (CLT) regime enables exact computation of the channel skewness. For discrete memoryless channels, this work generalizes Moulin's (2017) bounds on the asymptotic expansion of the maximum achievable message set size for nonsingular channels from the CLT regime to include the moderate deviations (MD) regime, thereby refining Altu\u{g} and Wagner's (2014) MD result. For an example binary symmetric channel and most practically important $(n, \epsilon)$ pairs, including $n \in [100, 500]$ and $\epsilon \in [10^{-10}, 10^{-1}]$, an approximation up to the channel skewness is the most accurate among several expansions in the literature. A derivation of the third-order term in the type-II error exponent of binary hypothesis testing in the MD regime is also included; the resulting third-order term is similar to the channel skewness.


title: Electrostatic partners and zeros of orthogonal and multiple orthogonal polynomials

abstract: For a given polynomial $P$ with simple zeros, and a given semiclassical weight $w$, we present a construction that yields a linear second-order differential equation (ODE), and in consequence, an electrostatic model for zeros of $P$. The coefficients of this ODE are written in terms of a dual polynomial that we call the electrostatic partner of $P$. This construction is absolutely general and can be carried out for any polynomial with simple zeros and any semiclassical weight on the complex plane. An additional assumption of quasi-orthogonality of $P$ with respect to $w$ allows us to give more precise bounds on the degree of the electrostatic partner. In the case of orthogonal and quasi-orthogonal polynomials, we recover some of the known results and generalize others. Additionally, for the Hermite--Pad\'e or multiple orthogonal polynomials of type II, this approach yields a system of linear second-order differential equations, from which we derive an electrostatic interpretation of their zeros in terms of a vector equilibrium. More detailed results are obtained in the special cases of Angelesco, Nikishin, and generalized Nikishin systems. We also discuss the discrete-to-continuous transition of these models in the asymptotic regime, as the number of zeros tends to infinity, into the known vector equilibrium problems. Finally, we discuss how the system of obtained second-order ODEs yields a third-order differential equation for these polynomials, well described in the literature. We finish the paper by presenting several illustrative examples.


title: Minimax decision rules for planning under uncertainty

abstract: It is common to use minimax rules to make decisions for planning when there is great uncertainty on what will happen in the future. Minimax regret is one popular version of this. We give an analysis of the behaviour of minimax rules in the case with a finite set of possible future scenarios. The use of minimax rules avoids the need to determine probabilities for each scenario, which is an attractive feature in many public sector settings. However, minimax rules will have sensitivity to the choice of scenarios. In many cases using a minimax approach will mean the requirement for what may be regarded as arbitrary probabilities on scenarios is replaced by a similarly arbitrary choice of a very small number of specific scenarios. We investigate this phenomenon. When regret-based rules are used there are also problems arising since the independence of irrelevant alternatives property fails, which can lead to opportunities to game the process. Our analysis of these phenomena considers cases where the decision variables are chosen from a convex set in $R^n$, as well as cases with a finite set of decision choices.


title: Geometrically programmed self-limited assembly of tubules using DNA origami colloids

abstract: Self-assembly is one of the most promising strategies for making functional materials at the nanoscale, yet new design principles for making self-limiting architectures, rather than spatially unlimited periodic lattice structures, are needed. To address this challenge, we explore the trade-offs between addressable assembly and self-closing assembly of a specific class of self-limiting structures: cylindrical tubules. We make triangular subunits using DNA origami that have specific, valence-limited interactions and designed binding angles, and study their assembly into tubules that have a self-limited width that is much larger than the size of an individual subunit. In the simplest case, the tubules are assembled from a single component by geometrically programming the dihedral angles between neighboring subunits. We show that the tubules can reach many micrometers in length and that their average width can be prescribed through the dihedral angles. We find that there is a distribution in the width and the chirality of the tubules, which we rationalize by developing a model that considers the finite bending rigidity of the assembled structure as well as the mechanism of self-closure. Finally, we demonstrate that the distributions of tubules can be further sculpted by increasing the number of subunit species, thereby increasing the assembly complexity, and demonstrate that using two subunit species successfully reduces the number of available end states by half. These results help to shed light on the roles of assembly complexity and geometry in self-limited assembly and could be extended to other self-limiting architectures, such as shells, toroids, or triply-periodic frameworks.


title: Estimating Conditional Average Treatment Effects with Missing Treatment Information

abstract: Estimating conditional average treatment effects (CATE) is challenging, especially when treatment information is missing. Although this is a widespread problem in practice, CATE estimation with missing treatments has received little attention. In this paper, we analyze CATE estimation in the setting with missing treatments where unique challenges arise in the form of covariate shifts. We identify two covariate shifts in our setting: (i) a covariate shift between the treated and control population; and (ii) a covariate shift between the observed and missing treatment population. We first theoretically show the effect of these covariate shifts by deriving a generalization bound for estimating CATE in our setting with missing treatments. Then, motivated by our bound, we develop the missing treatment representation network (MTRNet), a novel CATE estimation algorithm that learns a balanced representation of covariates using domain adaptation. By using balanced representations, MTRNet provides more reliable CATE estimates in the covariate domains where the data are not fully observed. In various experiments with semi-synthetic and real-world data, we show that our algorithm improves over the state-of-the-art by a substantial margin.


title: Minimal sets on continua with a dense free interval

abstract: We study minimal sets on continua $X$ with a dense free interval $J$ and a locally connected remainder. This class of continua includes important spaces such as the topologist's sine curve or the Warsaw circle. In the case when minimal sets on the remainder are known and the remainder is connected, we obtain a full characterization of the topological structure of minimal sets. In particular, a full characterization of minimal sets on $X$ is given in the case when $X\setminus J$ is a local dendrite.


title: Ammonia, carbon dioxide and the non-detection of the 2152 cm$^{-1}$ CO band

abstract: CO is one of the most abundant ice components on interstellar dust grains. When it is mixed with amorphous solid water (ASW) or located on its surface, an absorption band of CO at 2152 cm$^{-1}$ is always present in laboratory measurements. This spectral feature is attributed to the interaction of CO with dangling-OH bonds (dOH) in ASW. However, this band is absent in observational spectra of interstellar ices. This raises the question whether CO forms a relatively pure layer on top of ASW or is in close contact with ASW, but not via dangling bonds. We aim to determine whether the incorporation of NH$_3$ or CO$_2$ into ASW blocks the dOH and therefore reduces the 2152 cm$^{-1}$ band. We performed laboratory experiments to simulate the layered structure of the ice mantle, that is, we grew CO ice on top of 1) pure ASW, 2) NH$_3$:H$_2$O=10:100 mixed ice, and 3) CO$_2$:H$_2$O=20:100 mixed ice. Infrared spectra were measured to quantify the strength of the 2152 cm$^{-1}$ band. In addition, a second set of experiments were performed to determine how the incorporation of NH$_3$ into ASW affects the dOH band. We found that annealing the ice reduces the 2152 cm$^{-1}$ band and that NH$_3$ blocks the dOH on ASW surface and therefore reduces the 2152 cm$^{-1}$ band more effectively than CO$_2$. We suggest that this difference between NH$_3$ and CO$_2$ can be ascribed to the polarity of the guest molecule (NH$_3$ is a polar species, whereas CO$_2$ is apolar). The polarity implies that the formation of an H-bond between the N atom of ammonia and the dOH is a barrier-less reaction. We also determined the pore surface area of the ice mixtures as a function of the annealing temperature, and found that the nondetection of 2152 cm$^{-1}$ band does not necessarily exclude the possibility of a porous ice mantle.


title: A Modern Gauss-Markov Theorem? Really?

abstract: We show that the theorems in Hansen (2021a) (the version accepted by Econometrica), except for one, are not new as they coincide with classical theorems like the good old Gauss-Markov or Aitken Theorem, respectively; the exceptional theorem is incorrect. Hansen (2021b) corrects this theorem. As a result, all theorems in the latter version coincide with the above mentioned classical theorems. Furthermore, we also show that the theorems in Hansen (2022) (the version published in Econometrica) either coincide with the classical theorems just mentioned, or contain extra assumptions that are alien to the Gauss-Markov or Aitken Theorem.


title: SPICEprop: Backpropagating Errors Through Memristive Spiking Neural Networks

abstract: We present a fully memristive spiking neural network (MSNN) consisting of novel memristive neurons trained using the backpropagation through time (BPTT) learning rule. Gradient descent is applied directly to the memristive integrated-and-fire (MIF) neuron designed using analog SPICE circuit models, which generates distinct depolarization, hyperpolarization, and repolarization voltage waveforms. Synaptic weights are trained by BPTT using the membrane potential of the MIF neuron model and can be processed on memristive crossbars. The natural spiking dynamics of the MIF neuron model are fully differentiable, eliminating the need for gradient approximations that are prevalent in the spiking neural network literature. Despite the added complexity of training directly on SPICE circuit models, we achieve 97.58% accuracy on the MNIST testing dataset and 75.26% on the Fashion-MNIST testing dataset, the highest accuracies among all fully MSNNs.


title: Chiral Phases on the Lattice

abstract: While chiral quantum field theories (QFTs) describe a wide range of physical systems, from the standard model to topological quantum matter, the realization of chiral QFTs on a lattice has proved to be difficult due to the Nielsen-Ninomiya theorem and the possible presence of quantum anomalies. In this thesis, we use the connection between chiral phases of matter and chiral quantum field theories (QFTs) to define chiral QFTs on a lattice and allow a huge class of exotic field theories to be simulated numerically. Our work builds on the 'mirror fermion' approach to the problem of defining chiral theories on a lattice, which defines chiral field theories as the edge modes of chiral phases. We begin by reviewing the deep connections between chiral phases of matter, chiral field theories, and anomalies. We then develop numerical treatments of an $SU(2)$ chiral field theory, and provide a semiclassically solvable definition of Abelian $2+1$ chiral topological orders. This leads to an exactly solvable definition of chiral $U(1)$ SPT phases with zero correlation length, which we use to extract the edge chiral field theories exactly. These zero-correlation length models are vastly more simple than previous approaches to defining chiral field theories on the lattice.


title: The case of equality in geometric instances of Barthe's reverse Brascamp-Lieb inequalit

abstract: The works of Bennett, Carbery, Christ, Tao and of Valdimarsson have clarified when equality holds in the Brascamp-Lieb inequality. Here we characterize the case of equality in the Geometric case of Barthe's reverse Brascamp-Lieb inequality.


title: SMTNet: Hierarchical cavitation intensity recognition based on sub-main transfer network

abstract: With the rapid development of smart manufacturing, data-driven machinery health management has been of growing attention. In situations where some classes are more difficult to be distinguished compared to others and where classes might be organised in a hierarchy of categories, current DL methods can not work well. In this study, a novel hierarchical cavitation intensity recognition framework using Sub-Main Transfer Network, termed SMTNet, is proposed to classify acoustic signals of valve cavitation. SMTNet model outputs multiple predictions ordered from coarse to fine along a network corresponding to a hierarchy of target cavitation states. Firstly, a data augmentation method based on Sliding Window with Fast Fourier Transform (Swin-FFT) is developed to solve few-shot problem. Secondly, a 1-D double hierarchical residual block (1-D DHRB) is presented to capture sensitive features of the frequency domain valve acoustic signals. Thirdly, hierarchical multi-label tree is proposed to assist the embedding of the semantic structure of target cavitation states into SMTNet. Fourthly, experience filtering mechanism is proposed to fully learn a prior knowledge of cavitation detection model. Finally, SMTNet has been evaluated on two cavitation datasets without noise (Dataset 1 and Dataset 2), and one cavitation dataset with real noise (Dataset 3) provided by SAMSON AG (Frankfurt). The prediction accurcies of SMTNet for cavitation intensity recognition are as high as 95.32%, 97.16% and 100%, respectively. At the same time, the testing accuracies of SMTNet for cavitation detection are as high as 97.02%, 97.64% and 100%. In addition, SMTNet has also been tested for different frequencies of samples and has achieved excellent results of the highest frequency of samples of mobile phones.


title: On-Demand Directional Microwave Photon Emission Using Waveguide Quantum Electrodynamics

abstract: Routing quantum information between non-local computational nodes is a foundation for extensible networks of quantum processors. Quantum information transfer between arbitrary nodes is generally mediated either by photons that propagate between them, or by resonantly coupling nearby nodes. The utility is determined by the type of emitter, propagation channel, and receiver. Conventional approaches involving propagating microwave photons have limited fidelity due to photon loss and are often unidirectional, whereas architectures that use direct resonant coupling are bidirectional in principle, but can generally accommodate only a few local nodes. Here we demonstrate high-fidelity, on-demand, directional, microwave photon emission. We do this using an artificial molecule comprising two superconducting qubits strongly coupled to a bidirectional waveguide, effectively creating a chiral microwave waveguide. Quantum interference between the photon emission pathways from the molecule generates single photons that selectively propagate in a chosen direction. This circuit will also be capable of photon absorption, making it suitable for building interconnects within extensible quantum networks.


title: Electro-optic modulator based on Vanadium dioxide epsilon-near-zero vertical nanowire triple-cavity in Silicon Waveguide

abstract: We present an electro-optic modulator exploiting an array of Vanadium dioxide nanorods operating in the epsilon near regime as the active switching material integrated in a silicon waveguide. The modulator takes advantage of the insulator-to-metal transition of Vanadium dioxide to achieve a robust modulation depth of 19.7 dB/micron with a broad wavelength of operation. Using simulations, we demonstrate how the effective permittivity of the metamaterial can be tuned to a near-zero value by varying the nanorod geometry. The paper also proposes a novel hexagonal array design which achieves low insertion losses while retaining a strong modulation depth. The results provide insight into the design of ultra-compact modulators with high operation frequencies and low insertion losses.


title: Extinction of multiple populations and a team of Die-out Lyapunov functions

abstract: The extinction of species is a major problem of concern with a large literature. Our investigation gives insight into when species extinctions must occur, with an emphasis on determining which species might possibly die out and on how fast they die out. We investigate a differential equations model for population interactions with the goal of determining when several species (\ie, coordinates of a bounded solution) must die out or ``go extinct'' and must do so exponentially fast. Typically each coordinate represents the population density of a different species. For our main tool, we create what we call ``die-out'' Lyapunov functions. A given system may have several or many such functions, each of which is a function of a different set of coordinates. That die-out function implies that one of the species in its subset must die out exponentially fast -- for almost every choice of coefficients of the system. We create a ``team'' of die-out functions that work together to show that $k$ species must die, where $k$ is determined separately. Secondly, we present a ``trophic'' condition for generalized Lotka-Volterra systems that guarantees that there is a trapping region that is globally attracting. That implies that all solutions are bounded.


title: Queueing Systems with Preferred Service Delivery Times and Multiple Customer Classes

abstract: Motivated by the operational problems in click and collect systems, such as curbside pickup programs, we study a joint admission control and capacity allocation problem. We consider a system where arriving customers have preferred service delivery times and gauge the service quality based on the service provider's ability to complete the service as close as possible to the preferred time. Customers can be of different priority classes, and their priority may increase as they wait longer in the queue. The service provider can reject customers upon their arrival if the system is overloaded or outsource the service (alternatively work overtime) when the capacity is not enough. The service provider's goal is to find the minimum-cost admission and capacity allocation policy to dynamically decide when to serve and whom to serve. We model this problem as a Markov Decision Process. Our structural results partially characterize a set of suboptimal solutions, and we develop solution methods using these results. We also develop a problem-specific approximation method that is based on state aggregation to overcome the computational challenges. We present extensive computational results and discuss the impact of problem parameters on the optimal policy.


title: A mathematical study of the role of tBregs in breast cancer

abstract: A model for the mathematical study of immune response to breast cancer is proposed and studied, both analytically and numerically. It is a simplification of a complex one, recently introduced by two of the present authors. It serves for a compact study of the dynamical role in cancer promotion of a relatively recently described subgroup of regulatory B cells, which are evoked by the tumour.


title: A general approximation to nested Bayes factors with informed priors

abstract: A staple of Bayesian model comparison and hypothesis testing, Bayes factors are often used to quantify the relative predictive performance of two rival hypotheses. The computation of Bayes factors can be challenging, however, and this has contributed to the popularity of convenient approximations such as the BIC. Unfortunately, these approximations can fail in the case of informed prior distributions. Here we address this problem by outlining an approximation to informed Bayes factors for a focal parameter $\theta$. The approximation is computationally simple and requires only the maximum likelihood estimate $\hat\theta$ and its standard error. The approximation uses an estimated likelihood of $\theta$ and assumes that the posterior distribution for $\theta$ is unaffected by the choice of prior distribution for the nuisance parameters. The resulting Bayes factor for the null hypothesis $\mathcal{H}_0: \theta = \theta_0$ versus the alternative hypothesis $\mathcal{H}_1: \theta \sim g(\theta)$ is then easily obtained using the Savage--Dickey density ratio. Three real-data examples highlight the speed and closeness of the approximation compared to bridge sampling and Laplace's method. The proposed approximation facilitates Bayesian reanalyses of standard frequentist results, encourages application of Bayesian tests with informed priors, and alleviates the computational challenges that often frustrate both Bayesian sensitivity analyses and Bayes factor design analyses. The approximation is shown to suffer under small sample sizes and when the posterior distribution of the focal parameter is substantially influenced by the prior distributions on the nuisance parameters. The proposed methodology may also be used to approximate the posterior distribution for $\theta$ under $\mathcal{H}_1$.


title: CAMERA: A Method for Cost-aware, Adaptive, Multifidelity, Efficient Reliability Analysis

abstract: Estimating probability of failure in aerospace systems is a critical requirement for flight certification and qualification. Failure probability estimation involves resolving tails of probability distribution, and Monte Carlo sampling methods are intractable when expensive high-fidelity simulations have to be queried. We propose a method to use models of multiple fidelities that trade accuracy for computational efficiency. Specifically, we propose the use of multifidelity Gaussian process models to efficiently fuse models at multiple fidelity, thereby offering a cheap surrogate model that emulates the original model at all fidelities. Furthermore, we propose a novel sequential \emph{acquisition function}-based experiment design framework that can automatically select samples from appropriate fidelity models to make predictions about quantities of interest in the highest fidelity. We use our proposed approach in an importance sampling setting and demonstrate our method on the failure level set estimation and probability estimation on synthetic test functions as well as two real-world applications, namely, the reliability analysis of a gas turbine engine blade using a finite element method and a transonic aerodynamic wing test case using Reynolds-averaged Navier--Stokes equations. We demonstrate that our method predicts the failure boundary and probability more accurately and computationally efficiently while using varying fidelity models compared with using just a single expensive high-fidelity model.