title: PRIMA: General and Precise Neural Network Certification via Scalable Convex Hull Approximations

abstract: Formal verification of neural networks is critical for their safe adoption in real-world applications. However, designing a precise and scalable verifier which can handle different activation functions, realistic network architectures and relevant specifications remains an open and difficult challenge. In this paper, we take a major step forward in addressing this challenge and present a new verification framework, called PRIMA. PRIMA is both (i) general: it handles any non-linear activation function, and (ii) precise: it computes precise convex abstractions involving multiple neurons via novel convex hull approximation algorithms that leverage concepts from computational geometry. The algorithms have polynomial complexity, yield fewer constraints, and minimize precision loss. We evaluate the effectiveness of PRIMA on a variety of challenging tasks from prior work. Our results show that PRIMA is significantly more precise than the state-of-the-art, verifying robustness to input perturbations for up to 20%, 30%, and 34% more images than existing work on ReLU-, Sigmoid-, and Tanh-based networks, respectively. Further, PRIMA enables, for the first time, the precise verification of a realistic neural network for autonomous driving within a few minutes.


title: Symmetric decompositions, triangulations and real-rootedness

abstract: Polynomials which afford nonnegative, real-rooted symmetric decompositions have been investigated recently in algebraic, enumerative and geometric combinatorics. Br\"and\'en and Solus have given sufficient conditions under which the image of a polynomial under a certain operator associated to barycentric subdivision has such a decomposition. This paper gives a new proof of their result which generalizes to subdivision operators in the setting of uniform triangulations of simplicial complexes, introduced by the first named author. Sufficient conditions under which these decompositions are also interlacing are described. Applications yield new classes of polynomials in geometric combinatorics which afford nonnegative, real-rooted symmetric decompositions. Some interesting questions in $f$-vector theory arise from this work.


title: Flag Hilbert-Poincar\'e series of hyperplane arrangements and their Igusa zeta functions

abstract: We introduce and study a class of multivariate rational functions associated with hyperplane arrangements, called flag Hilbert-Poincar\'e series. These series are intimately connected with Igusa local zeta functions of products of linear polynomials, and their motivic and topological relatives. Our main results include a self-reciprocity result for central arrangements defined over fields of characteristic zero. We also prove combinatorial formulae for a specialization of the flag Hilbert-Poincar\'e series for irreducible Coxeter arrangements of types $\mathsf{A}$, $\mathsf{B}$, and $\mathsf{D}$ in terms of total partitions of the respective types. We show that a different specialization of the flag Hilbert-Poincar\'e series, which we call the coarse flag Hilbert-Poincar\'e series, exhibits intriguing nonnegativity features and - in the case of Coxeter arrangements - connections with Eulerian polynomials. For numerous classes and examples of hyperplane arrangements, we determine their (coarse) flag Hilbert-Poincar\'e series. Some computations were aided by a SageMath package we developed.


title: Experimental Validation of a Dynamic Equivalent Model for Microgrids

abstract: The goal of this paper is the experimental validation of a gray-box equivalent modeling approach applied to microgrids. The main objective of the equivalent modeling is to represent the dynamic response of a microgrid with a simplified model. The main contribution of this work is the experimental validation of a two-step process, composed by the definition of a nonlinear equivalent model with operational constraints, adapted to the microgrid environment, and the identification procedure used to define the model parameters. Once the parameters are identified, the simplified model is ready to reproduce the microgrid behavior to voltage and frequency variations, in terms of active and reactive power exchanges at the point of common coupling. To validate the proposed approach, a set of experimental tests have been carried out on a real LV microgrid considering different configurations, including both grid-connected and islanded operating conditions. Results show the effectiveness of the proposed technique and the applicability of the model to perform dynamic simulations.


title: Topology-Aware Correlations Between Relations for Inductive Link Prediction in Knowledge Graphs

abstract: Inductive link prediction -- where entities during training and inference stages can be different -- has been shown to be promising for completing continuously evolving knowledge graphs. Existing models of inductive reasoning mainly focus on predicting missing links by learning logical rules. However, many existing approaches do not take into account semantic correlations between relations, which are commonly seen in real-world knowledge graphs. To address this challenge, we propose a novel inductive reasoning approach, namely TACT, which can effectively exploit Topology-Aware CorrelaTions between relations in an entity-independent manner. TACT is inspired by the observation that the semantic correlation between two relations is highly correlated to their topological structure in knowledge graphs. Specifically, we categorize all relation pairs into several topological patterns, and then propose a Relational Correlation Network (RCN) to learn the importance of the different patterns for inductive link prediction. Experiments demonstrate that TACT can effectively model semantic correlations between relations, and significantly outperforms existing state-of-the-art methods on benchmark datasets for the inductive link prediction task.


title: Relative field line helicity of a large eruptive solar active region

abstract: Context. Magnetic helicity is a physical quantity of great importance in the study of astrophysical and natural plasmas. Although a density for helicity cannot be defined, a good proxy for it is field line helicity. The appropriate quantity for use in solar conditions is relative field line helicity (RFLH). Aims. This work aims to study in detail the behaviour of RFLH, for the first time, in a solar active region (AR). Methods. The target active region is the large, eruptive AR 11158. In order to compute RFLH and all other quantities of interest we use a non-linear force-free reconstruction of the AR coronal magnetic field of excelent quality. Results. We find that the photospheric morphology of RFLH is quite different than that of the magnetic field or of the electrical current, and this is not sensitive to the chosen gauge in the computation of RFLH. The value of helicity experiences a large decrease, 25% of its pre-flare value, during an X-class flare of the AR, a change that is also depicted in the photospheric morphology of RFLH. Moreover, the area of this change coincides with the area that encompasses the flux rope, the magnetic structure that later erupted. Conclusions. The use of RFLH can provide important information about the value and location of the magnetic helicity expelled from the solar atmosphere during eruptive events.


title: Rediscussion of eclipsing binaries. Paper IV. The evolved G-type system AN Camelopardalis

abstract: AN Cam is a little-studied eclipsing binary containing somewhat evolved components in an orbit with a period of 21.0 d and an eccentricity of 0.47. A spectroscopic orbit based on photoelectric radial velocities was published in 1977. AN Cam has been observed using the TESS satellite in three sectors: the data were obtained in long-cadence mode and cover nine eclipses. By modelling these data and published radial velocities we obtain masses of 1.380 +/- 0.021 Msun and 1.402 +/- 0.025 Msun, and radii of 2.159 +/- 0.012 Rsun and 2.646 +/- 0.014 Rsun. We also derive a precise orbital ephemeris from these data and recent times of minimum light, but find that the older times of minimum light cannot be fitted assuming a constant orbital period. This could be caused by astrophysical or instrumental effects; forthcoming TESS observations will help the investigation of this issue. We use the Gaia EDR3 parallax and optical/infrared apparent magnitudes to measure effective temperatures of 6050 +/- 150 K and 5750 +/- 150 K: the primary star is hotter but smaller and less massive than its companion. A comparison with theoretical models indicates that the system has an approximately solar chemical composition and an age of 3.3 Gyr. Despite the similarity of their masses the two stars are in different evolutionary states: the primary is near the end of its main-sequence lifetime and the secondary is now a subgiant. AN Cam is a promising candidate for constraining the strength of convective core overshooting in 1.4 Msun stars.


title: The OTELO survey: Faint end of the luminosity function of [OII] emitters at <z>= 1.43

abstract: In this paper, we aim to study the main properties and luminosity function (LF) of the [OII] emitters detected in the OTELO survey in order to characterise the star formation processes in low-mass galaxies at $z\sim1.43$ and to constrain the faint-end of the LF. Here, we describe the selection method and analysis of the emitters obtained from narrow-band scanning techniques. In addition, we present several relevant properties of the emitters and discuss the selection biases and uncertainties in the determination of the LF and the star formation rate density (SFRD). We confirmed a total of 60 sources from a preliminary list of 332 candidates as [OII] emitters. Approximately 93% of the emitters have masses in the range of $10^{8}<M_{*}/{\rm M_{\odot}}<10^{9}$. All of our emitters are classified as late-type galaxies, with a lower value of $(u-v)$\, when compared with the rest of the emitters of the OTELO survey. We find that the cosmic variance strongly affects the normalisation ($\phi^*$) of the LF and explains the discrepancy of our results when compared with those obtained from surveys of much larger volumes. However, we are able to determine the faint-end slope of the LF, namely, $\alpha=-1.42\pm0.06$, by sampling the LF down to $\sim1$\,dex lower than in previous works. We present our calculation of the SFRD of our sample and compare it to the value obtained in previous studies from the literature.


title: Puiseux Series and Algebraic Solutions of First Order Autonomous AODEs -- A MAPLE Package

abstract: There exist several methods for computing exact solutions of algebraic differential equations. Most of the methods, however, do not ensure existence and uniqueness of the solutions and might fail after several steps, or are restricted to linear equations. The authors have presented in previous works a method to overcome this problem for autonomous first order algebraic ordinary differential equations and formal Puiseux series solutions and algebraic solutions. In the first case, all solutions can uniquely be represented by a sufficiently large truncation and in the latter case by its minimal polynomial. The main contribution of this paper is the implementation, in a MAPLE-package named FirstOrderSolve, of the algorithmic ideas presented therein. More precisely, all formal Puiseux series and algebraic solutions, including the generic and singular solutions, are computed and described uniquely. The computation strategy is to reduce the given differential equation to a simpler one by using local parametrizations and the already known degree bounds.


title: sparta: Sparse Tables and their Algebra with a View Towards High Dimensional Graphical Models

abstract: A graphical model is a multivariate (potentially very high dimensional) probabilistic model, which is formed by combining lower dimensional components. Inference (computation of conditional probabilities) is based on message passing algorithms that utilize conditional independence structures. In graphical models for discrete variables with finite state spaces, there is a fundamental problem in high dimensions: A discrete distribution is represented by a table of values, and in high dimensions such tables can become prohibitively large. In inference, such tables must be multiplied which can lead to even larger tables. The sparta package meets this challenge by implementing methods that efficiently handles multiplication and marginalization of sparse tables. The package was written in the R programming language and is freely available from the Comprehensive R Archive Network (CRAN). The companion package jti, also on CRAN, was developed to showcase the potential of sparta in connection to the Junction Tree Algorithm. We show, that jti is able to handle highly complex graphical models which are otherwise infeasible due to lack of computer memory, using sparta as a backend for table operations.


title: Software Startup Practices -- Software Development in Startups through the Lens of the Essence Theory of Software Engineering

abstract: Software startups continue to be important drivers of economy globally. As the initial investment required to found a new software company becomes smaller and smaller resulting from technological advances such as cloud technology, increasing numbers of new software startups are born. Typically, the main argument for studying software startups is that they differ from mature software organizations in various ways, thus making the findings of many existing studies not directly applicable to them. How, exactly, software startups really differ from other types of software organizations as an on-going debate. In this paper, we seek to better understand how software startups differ from mature software organizations in terms of development practices. Past studies have primarily studied method use, and in comparison, we take on a more atomic approach by focusing on practices. Utilizing the Essence Theory of Software Engineering as a framework, we split these practices into categories for analysis while simultaneously evaluating the suitability of the theory for the context of software startups. Based on the results, we propose changes to the Essence Theory of Software Engineering for it to better fit the startup context.


title: Fast Outage Analysis of Large-scale Production Clouds with Service Correlation Mining

abstract: Cloud-based services are surging into popularity in recent years. However, outages, i.e., severe incidents that always impact multiple services, can dramatically affect user experience and incur severe economic losses. Locating the root-cause service, i.e., the service that contains the root cause of the outage, is a crucial step to mitigate the impact of the outage. In current industrial practice, this is generally performed in a bootstrap manner and largely depends on human efforts: the service that directly causes the outage is identified first, and the suspected root cause is traced back manually from service to service during diagnosis until the actual root cause is found. Unfortunately, production cloud systems typically contain a large number of interdependent services. Such a manual root cause analysis is often time-consuming and labor-intensive. In this work, we propose COT, the first outage triage approach that considers the global view of service correlations. COT mines the correlations among services from outage diagnosis data. After learning from historical outages, COT can infer the root cause of emerging ones accurately. We implement COT and evaluate it on a real-world dataset containing one year of data collected from Microsoft Azure, one of the representative cloud computing platforms in the world. Our experimental results show that COT can reach a triage accuracy of 82.1%~83.5%, which outperforms the state-of-the-art triage approach by 28.0%~29.7%.


title: Revisiting the role of CP-conserving processes in cosmological particle-antiparticle asymmetries

abstract: We point out qualitatively different possibilities on the role of CP-conserving processes in generating cosmological particle-antiparticle asymmetries, with illustrative examples from models in leptogenesis and asymmetric dark matter production. In particular, we consider scenarios in which the CP-violating and CP-conserving processes are either both decays or both scatterings, thereby being naturally of comparable rates. This is in contrast to the previously considered CP-conserving processes in models of leptogenesis in different see-saw mechanisms, in which the CP-conserving scatterings typically have lower rates compared to the CP-violating decays, due to a Boltzmann suppression. We further point out that the CP-conserving processes can play a dual role if the asymmetry is generated in the mother sector itself, in contrast to the conventional scenarios in which it is generated in the daughter sector. This is because, the CP-conserving processes initially suppress the asymmetry generation by controlling the out-of-equilibrium number densities of the bath particles, but subsequently modify the ratio of particle anti-particle yields at the present epoch by eliminating the symmetric component of the bath particles through pair-annihilations, leading to a competing effect stemming from the same process at different epochs. We find that the asymmetric yields for relevant particle-antiparticle systems can vary by orders of magnitude depending upon the relative size of the CP-conserving and violating reaction rates.


title: Fine-Grained Off-Road Semantic Segmentation and Mapping via Contrastive Learning

abstract: Road detection or traversability analysis has been a key technique for a mobile robot to traverse complex off-road scenes. The problem has been mainly formulated in early works as a binary classification one, e.g. associating pixels with road or non-road labels. Whereas understanding scenes with fine-grained labels are needed for off-road robots, as scenes are very diverse, and the various mechanical performance of off-road robots may lead to different definitions of safe regions to traverse. How to define and annotate fine-grained labels to achieve meaningful scene understanding for a robot to traverse off-road is still an open question. This research proposes a contrastive learning based method. With a set of human-annotated anchor patches, a feature representation is learned to discriminate regions with different traversability, a method of fine-grained semantic segmentation and mapping is subsequently developed for off-road scene understanding. Experiments are conducted on a dataset of three driving segments that represent very diverse off-road scenes. An anchor accuracy of 89.8% is achieved by evaluating the matching with human-annotated image patches in cross-scene validation. Examined by associated 3D LiDAR data, the fine-grained segments of visual images are demonstrated to have different levels of toughness and terrain elevation, which represents their semantical meaningfulness. The resultant maps contain both fine-grained labels and confidence values, providing rich information to support a robot traversing complex off-road scenes.


title: Democratic Forking: Choosing Sides with Social Choice

abstract: Any community in which membership is optional may eventually break apart, or fork. For example, forks may occur in political parties, business partnerships, social groups, cryptocurrencies, and federated governing bodies. Forking is typically the product of informal social processes or the organized action of an aggrieved minority, and it is not always amicable. Forks usually come at a cost, and can be seen as consequences of collective decisions that destabilize the community. Here, we provide a social choice setting in which agents can report preferences not only over a set of alternatives, but also over the possible forks that may occur in the face of disagreement. We study this social choice setting, concentrating on stability issues and concerns of strategic agent behavior.


title: GraphMineSuite: Enabling High-Performance and Programmable Graph Mining Algorithms with Set Algebra

abstract: We propose GraphMineSuite (GMS): the first benchmarking suite for graph mining that facilitates evaluating and constructing high-performance graph mining algorithms. First, GMS comes with a benchmark specification based on extensive literature review, prescribing representative problems, algorithms, and datasets. Second, GMS offers a carefully designed software platform for seamless testing of different fine-grained elements of graph mining algorithms, such as graph representations or algorithm subroutines. The platform includes parallel implementations of more than 40 considered baselines, and it facilitates developing complex and fast mining algorithms. High modularity is possible by harnessing set algebra operations such as set intersection and difference, which enables breaking complex graph mining algorithms into simple building blocks that can be separately experimented with. GMS is supported with a broad concurrency analysis for portability in performance insights, and a novel performance metric to assess the throughput of graph mining algorithms, enabling more insightful evaluation. As use cases, we harness GMS to rapidly redesign and accelerate state-of-the-art baselines of core graph mining problems: degeneracy reordering (by up to >2x), maximal clique listing (by up to >9x), k-clique listing (by 1.1x), and subgraph isomorphism (by up to 2.5x), also obtaining better theoretical performance bounds.


title: Effects of Image Compression on Face Image Manipulation Detection: A Case Study on Facial Retouching

abstract: In the past years, numerous methods have been introduced to reliably detect digital face image manipulations. Lately, the generalizability of these schemes has been questioned in particular with respect to image post-processing. Image compression represents a post-processing which is frequently applied in diverse biometric application scenarios. Severe compression might erase digital traces of face image manipulation and hence hamper a reliable detection thereof. In this work, the effects of image compression on face image manipulation detection are analyzed. In particular, a case study on facial retouching detection under the influence of image compression is presented. To this end, ICAO-compliant subsets of two public face databases are used to automatically create a database containing more than 9,000 retouched reference images together with unconstrained probe images. Subsequently, reference images are compressed applying JPEG and JPEG 2000 at compression levels recommended for face image storage in electronic travel documents. Novel detection algorithms utilizing texture descriptors and deep face representations are proposed and evaluated in a single image and differential scenario. Results obtained from challenging cross-database experiments in which the analyzed retouching technique is unknown during training yield interesting findings: (1) most competitive detection performance is achieved for differential scenarios employing deep face representations; (2) image compression severely impacts the performance of face image manipulation detection schemes based on texture descriptors while methods utilizing deep face representations are found to be highly robust; (3) in some cases, the application of image compression might as well improve detection performance.


title: Foundations of Population-Based SHM, Part IV: The Geometry of Spaces of Structures and their Feature Spaces

abstract: One of the requirements of the population-based approach to Structural Health Monitoring (SHM) proposed in the earlier papers in this sequence, is that structures be represented by points in an abstract space. Furthermore, these spaces should be metric spaces in a loose sense; i.e. there should be some measure of distance applicable to pairs of points; similar structures should then be close in the metric. However, this geometrical construction is not enough for the framing of problems in data-based SHM, as it leaves undefined the notion of feature spaces. Interpreting the feature values on a structure-by-structure basis as a type of field over the space of structures, it seems sensible to borrow an idea from modern theoretical physics, and define feature assignments as sections in a vector bundle over the structure space. With this idea in place, one can interpret the effect of environmental and operational variations as gauge degrees of freedom, as in modern gauge field theories. This paper will discuss the various geometrical structures required for an abstract theory of feature spaces in SHM, and will draw analogies with how these structures have shown their power in modern physics. In the second part of the paper, the problem of determining the normal condition cross section of a feature bundle is addressed. The solution is provided by the application of Graph Neural Networks (GNN), a versatile non-Euclidean machine learning algorithm which is not restricted to inputs and outputs from vector spaces. In particular, the algorithm is well suited to operating directly on the sort of graph structures which are an important part of the proposed framework for PBSHM. The solution of the normal section problem is demonstrated for a heterogeneous population of truss structures for which the feature of interest is the first natural frequency.


title: Automatic Exploration Process Adjustment for Safe Reinforcement Learning with Joint Chance Constraint Satisfaction

abstract: In reinforcement learning (RL) algorithms, exploratory control inputs are used during learning to acquire knowledge for decision making and control, while the true dynamics of a controlled object is unknown. However, this exploring property sometimes causes undesired situations by violating constraints regarding the state of the controlled object. In this paper, we propose an automatic exploration process adjustment method for safe RL in continuous state and action spaces utilizing a linear nominal model of the controlled object. Specifically, our proposed method automatically selects whether the exploratory input is used or not at each time depending on the state and its predicted value as well as adjusts the variance-covariance matrix used in the Gaussian policy for exploration. We also show that our exploration process adjustment method theoretically guarantees the satisfaction of the constraints with the pre-specified probability, that is, the satisfaction of a joint chance constraint at every time. Finally, we illustrate the validity and the effectiveness of our method through numerical simulation.


title: XMMU J050722.1-684758: Discovery of a new Be X-ray binary pulsar likely associated with the supernova remnant MCSNR J0507-6847

abstract: We report the discovery of a new high mass X-ray binary pulsar, XMMU J050722.1-684758, possibly associated with the supernova remnant MCSNR J0507-6847 in the Large Magellanic Cloud, using XMM-Newton X-ray observations. Pulsations with a periodicity of 570 s are discovered from the Be X-ray binary XMMU J050722.1-684758 confirming its nature as a HMXB pulsar. The HMXB is located near the geometric centre of the supernova remnant MCSNR J0507-6847 (0.9 arcmin from the centre) which supports the XRB-SNR association. The estimated age of the supernova remnant is 43-63 kyr which points to a middle aged to old supernova remnant. The large diameter of the supernova remnant combined with the lack of distinctive shell counterparts in optical and radio indicates that the SNR is expanding into the tenous environment of the superbubble N103. The estimated magnetic field strength of the neutron star is $B\gtrsim10^{14}$ G assuming a spin equilibrium condition which is expected from the estimated age of the parent remnant and assuming that the measured mass-accretion rate remained constant throughout.


title: Highly accurate operator factorization methods for the integral fractional Laplacian and its generalization

abstract: In this paper, we propose a new class of operator factorization methods to discretize the integral fractional Laplacian $(-\Delta)^\frac{\alpha}{2}$ for $\alpha \in (0, 2)$. The main advantage of our method is to easily increase numerical accuracy by using high-degree Lagrange basis functions, but remain the scheme structure and computer implementation unchanged. Moreover, our discretization of the fractional Laplacian results in a symmetric (multilevel) Toeplitz differentiation matrix, which not only saves memory cost in simulations but enables efficient computations via the fast Fourier transforms. The performance of our method in both approximating the fractional Laplacian and solving the fractional Poisson problems was detailedly examined. It shows that our method has an optimal accuracy of ${\mathcal O}(h^2)$ for constant or linear basis functions, while ${\mathcal O}(h^4)$ if quadratic basis functions are used, with $h$ a small mesh size. Note that this accuracy holds for any $\alpha \in (0, 2)$ and can be further increased if higher-degree basis functions are used. If the solution of fractional Poisson problem satisfies $u \in C^{m, l}(\bar{\Omega})$ for $m \in {\mathbb N}$ and $0 < l < 1$, then our method has an accuracy of ${\mathcal O}\big(h^{\min\{m+l,\, 2\}}\big)$ for constant and linear basis functions, while ${\mathcal O}\big(h^{\min\{m+l,\, 4\}}\big)$ for quadratic basis functions. Additionally, our method can be readily applied to study generalized fractional Laplacians with a symmetric kernel function, and numerical study on the tempered fractional Poisson problem demonstrates its efficiency.


title: Rapid microwave-only characterization and readout of quantum dots using multiplexed gigahertz-frequency resonators

abstract: Superconducting resonators enable fast characterization and readout of mesoscopic quantum devices. Finding ways to perform measurements of interest on such devices using resonators only is therefore of great practical relevance. We report the experimental investigation of an InAs nanowire multi-quantum dot device by probing GHz resonators connected to the device. First, we demonstrate accurate extraction of the DC conductance from measurements of the high-frequency admittance. Because our technique does not rely on DC calibration, it could potentially obviate the need for DC measurements in semiconductor qubit devices. Second, we demonstrate multiplexed gate sensing and the detection of charge tunneling on microsecond time scales. The GHz detection of dispersive resonator shifts allows rapid acquisition of charge-stability diagrams, as well as resolving charge tunneling in the device with a signal-to-noise ratio of up to 15 in one microsecond. Our measurements show that GHz-frequency resonators may serve as a universal tool for fast tune-up and high-fidelity readout of semiconductor qubits.


title: Deformations of K\"ahler manifolds to normal bundles and restricted volumes of big classes

abstract: The deformation of a variety $X$ to the normal cone of a subvariety $Y$ is a classical construction in algebraic geometry. In this paper we study the case when $(X,\omega)$ is a compact K\"ahler manifold and $Y$ is a submanifold. The deformation space $\mathcal{X}$ is fibered over $\mathbb{P}^1$ and all the fibers $X_{\tau}$ are isomorphic to $X$, except the zero-fiber, which has the projective completion of the normal bundle $N_{Y|X}$ as one of its components. The first main result of this paper is that one can find K\"ahler forms on modifications of $\mathcal{X}$ which restricts to $\omega$ on $X_1$ and which makes the volume of the normal bundle in the zero-fiber come arbitrarily close to the volume of $X$. Phrased differently, we find K\"ahler deformations of $(X,\omega)$ such that almost all of the mass ends up in the normal bundle. The proof relies on a general result on the volume of big cohomology classes, which is the other main result of the paper. A $(1,1)$ cohomology class on a compact K\"ahler manifold $X$ is said to be big if it contains the sum of a K\"ahler form and a closed positive current. A quantative measure of bigness is provided by the volume function, and there is also a related notion of restricted volume along a submanifold. We prove that if $Y$ is a smooth hypersurface which intersects the K\"ahler locus of a big class $\alpha$ then up to a dimensional constant, the restricted volume of $\alpha$ along $Y$ is equal to the derivative of the volume at $\alpha$ in the direction of the cohomology class of $Y$. This generalizes the corresponding result on the volume of line bundles due to Boucksom-Favre-Jonsson and independently Lazarsfeld-Musta\c{t}\u{a}.


title: Nonlinear versions of Korovkin's abstract theorems

abstract: In this paper we prove Korovkin type theorems for sequences of sublinear, monotone and weak additive operators acting on function spaces C(X); where X is a compact or a locally compact metric space. Our results are illustrated by a series of examples


title: MAMBPO: Sample-efficient multi-robot reinforcement learning using learned world models

abstract: Multi-robot systems can benefit from reinforcement learning (RL) algorithms that learn behaviours in a small number of trials, a property known as sample efficiency. This research thus investigates the use of learned world models to improve sample efficiency. We present a novel multi-agent model-based RL algorithm: Multi-Agent Model-Based Policy Optimization (MAMBPO), utilizing the Centralized Learning for Decentralized Execution (CLDE) framework. CLDE algorithms allow a group of agents to act in a fully decentralized manner after training. This is a desirable property for many systems comprising of multiple robots. MAMBPO uses a learned world model to improve sample efficiency compared to model-free Multi-Agent Soft Actor-Critic (MASAC). We demonstrate this on two simulated multi-robot tasks, where MAMBPO achieves a similar performance to MASAC, but requires far fewer samples to do so. Through this, we take an important step towards making real-life learning for multi-robot systems possible.


title: Real-time RGBD-based Extended Body Pose Estimation

abstract: We present a system for real-time RGBD-based estimation of 3D human pose. We use parametric 3D deformable human mesh model (SMPL-X) as a representation and focus on the real-time estimation of parameters for the body pose, hands pose and facial expression from Kinect Azure RGB-D camera. We train estimators of body pose and facial expression parameters. Both estimators use previously published landmark extractors as input and custom annotated datasets for supervision, while hand pose is estimated directly by a previously published method. We combine the predictions of those estimators into a temporally-smooth human pose. We train the facial expression extractor on a large talking face dataset, which we annotate with facial expression parameters. For the body pose we collect and annotate a dataset of 56 people captured from a rig of 5 Kinect Azure RGB-D cameras and use it together with a large motion capture AMASS dataset. Our RGB-D body pose model outperforms the state-of-the-art RGB-only methods and works on the same level of accuracy compared to a slower RGB-D optimization-based solution. The combined system runs at 30 FPS on a server with a single GPU. The code will be available at https://saic-violet.github.io/rgbd-kinect-pose


title: ASC-Net : Adversarial-based Selective Network for Unsupervised Anomaly Segmentation

abstract: We introduce a neural network framework, utilizing adversarial learning to partition an image into two cuts, with one cut falling into a reference distribution provided by the user. This concept tackles the task of unsupervised anomaly segmentation, which has attracted increasing attention in recent years due to their broad applications in tasks with unlabelled data. This Adversarial-based Selective Cutting network (ASC-Net) bridges the two domains of cluster-based deep learning methods and adversarial-based anomaly/novelty detection algorithms. We evaluate this unsupervised learning model on BraTS brain tumor segmentation, LiTS liver lesion segmentation, and MS-SEG2015 segmentation tasks. Compared to existing methods like the AnoGAN family, our model demonstrates tremendous performance gains in unsupervised anomaly segmentation tasks. Although there is still room to further improve performance compared to supervised learning algorithms, the promising experimental results shed light on building an unsupervised learning algorithm using user-defined knowledge.


title: A Machine Learning Approach for Predicting Human Preference for Graph Layouts

abstract: Understanding what graph layout human prefer and why they prefer is significant and challenging due to the highly complex visual perception and cognition system in human brain. In this paper, we present the first machine learning approach for predicting human preference for graph layouts. In general, the data sets with human preference labels are limited and insufficient for training deep networks. To address this, we train our deep learning model by employing the transfer learning method, e.g., exploiting the quality metrics, such as shape-based metrics, edge crossing and stress, which are shown to be correlated to human preference on graph layouts. Experimental results using the ground truth human preference data sets show that our model can successfully predict human preference for graph layouts. To our best knowledge, this is the first approach for predicting qualitative evaluation of graph layouts using human preference experiment data.


title: Tree of Knowledge: an Online Platform for Learning the Behaviour of Complex Systems

abstract: Many social sciences such as psychology and economics try to learn the behaviour of complex agents such as humans, organisations and countries. The current statistical methods used for learning this behaviour try to infer generally valid behaviour, but can only learn from one type of study at a time. Furthermore, only data from carefully designed studies can be used, as the phenomenon of interest has to be isolated and confounding factors accounted for. These restrictions limit the robustness and accuracy of insights that can be gained from social/economic systems. Here we present the online platform TreeOfKnowledge which implements a new methodology specifically designed for learning complex behaviours from complex systems: agent-based behaviour learning. With agent-based behaviour learning it is possible to gain more accurate and robust insights as it does not have the restriction of conventional statistics. It learns agent behaviour from many heterogenous datasets and can learn from these datasets even if the phenomenon of interest is not directly observed, but appears deep within complex systems. This new methodology shows how the internet and advances in computational power allow for more accurate and powerful mathematical models.


title: SasCsvToolkit -- A versatile parallel 'bag-of-tasks' job submission application on heterogeneous and homogeneous platforms for Big Data Analytics such as for Biomedical Informatics

abstract: Background: The need for big data analysis requires being able to process large data which are being held fine-tuned for usage by corporate. It is only very recently that the need for big data has caught attention for low budget corporate groups and academia who typically do not have money and resources to buy expensive licenses of big data analysis platforms such as SAS. The corporate continue to work on SAS data format largely because of systemic organizational history and that the prior codes have been built on them. The data-providers continue to thus provide data in SAS formats. Acute sudden need has arisen because of this gap of data being in SAS format and the coders not having a SAS expertise or training background as the economic and inertial forces acting of having shaped these two class of people have been different. Method: We analyze the differences and thus the need for SasCsvToolkit which helps to generate a CSV file for a SAS format data so that the data scientist can then make use of his skills in other tools that can process CSVs such as R, SPSS, or even Microsoft Excel. At the same time, it also provides conversion of CSV files to SAS format. Apart from this, a SAS database programmer always struggles in finding the right method to do a database search, exact match, substring match, except condition, filters, unique values, table joins and data mining for which the toolbox also provides template scripts to modify and use from command line. Results: The toolkit has been implemented on SLURM scheduler platform as a `bag-of-tasks` algorithm for parallel and distributed workflow though serial version has also been incorporated.


title: Towards a Broad-Band Coherent Perfect Absorption in systems without Scale-Invariance

abstract: We experimentally and theoretically challenge the concept of coherent perfect absorption (CPA) as a narrow frequency resonant mechanism associated with scattering processes that respect scale-invariance. Using a microwave platform, consisting of a lossy nonlinear resonator coupled to two interrogating antennas, we show that a coherent incident excitation can trigger a self-induced perfect absorption once its intensity exceeds a critical value. Importantly, a (near) perfect absorption persists for a broad band frequency range around the nonlinear CPA condition. Its origin is traced to a quartic behavior that the absorbance spectrum acquires in the proximity of a CPA associated with a new kind of exceptional point degeneracy related to the zeros of the nonlinear scattering operator.


title: Enumerating all bilocal Clifford distillation protocols through symmetry reduction

abstract: Entanglement distillation is an essential building block in quantum communication protocols. Here, we study the class of near-term implementable distillation protocols that use bilocal Clifford operations followed by a single round of communication. We introduce tools to enumerate and optimise over all protocols for up to $n=5$ (not necessarily equal) Bell-diagonal states using a commodity desktop computer. Furthermore, by exploiting the symmetries of the input states, we find all protocols for up to $n=8$ copies of a Werner state. For the latter case, we present circuits that achieve the highest fidelity with perfect operations and no decoherence. These circuits have modest depth and number of two-qubit gates. Our results are based on a correspondence between distillation protocols and double cosets of the symplectic group, and improve on previously known protocols.


title: Cosmological Constraints on the Coupling Model from Observational Hubble Parameter and Baryon Acoustic Oscillation Measurements

abstract: In the paper, we consider two models in which dark energy is coupled with either dust matter or dark matter, and discuss the conditions that allow more time for structure formation to take place at high redshifts. These models are expected to have a larger age of the universe than that of $\Lambda$CDM [universe consists of cold dark matter (CDM) and dark energy (a cosmological constant, $\Lambda$)], so it can explain the formation of high redshift gravitationally bound systems which the $\Lambda$CDM model cannot interpret. We use the observational Hubble parameter data (OHD) and Hubble parameter obtained from cosmic chronometers method ($H(z)$) in combination with baryon acoustic oscillation (BAO) data to constrain these models. With the best-fitting parameters, we discuss how the age, the deceleration parameter, and the energy density parameters evolve in the new universes, and compare them with that of $\Lambda$CDM.


title: Trotter-Kato approximations of semilinear stochastic evolution equations in Hilbert spaces

abstract: Motivated by the work of T.E. Govindan in [5,8,9], this paper is concerned with a more general semilinear stochastic evolution equation. The difference between the equations considered in this paper and the previous one is that it makes some changes to the nonlinear function in random integral, which also depends on the probability distribution of stochastic process at that time. First, this paper considers the existence and uniqueness of mild solutions for such equations. Furthermore, Trotter-Kato approximation system is introduced for the mild solutions, and the weak convergence of induced probability measures and zeroth-order approximations are obtained. Then we consider the classical limit theorem about the parameter dependence of this kind of equations. Finally, an example of stochastic partial differential equation is given to illustrate our results.


title: Convective core sizes in rotating massive stars: I. Constraints from solar metallicity OB field stars

abstract: Spectroscopic studies of Galactic O and B stars show that many stars with masses above 8 M$_{\odot}$ are observed in the HR diagram just beyond the Main-Sequence (MS) band predicted by stellar models computed with a moderate overshooting. This may be an indication that the convective core sizes in stars in the upper part of the HR diagram are larger than predicted by these models. Combining stellar evolution models and spectroscopic parameters derived for a large sample of Galactic O and B stars, including brand new information about their projected rotational velocities, we reexamine the question of the convective core size in MS massive stars. We confirm that for stars more massive than about 8 M$_{\odot}$, the convective core size at the end of the MS phase increases more rapidly with the mass than in models computed with a constant step overshoot chosen to reproduce the main sequence width in the low mass range (around 2 M$_{\odot}$). This conclusion is valid for both the cases of non-rotating models and rotating models either with a moderate or a strong angular momentum transport. The increase of the convective core mass with the mass obtained from the TAMS position is, however, larger than the one deduced from the surface velocity drop for masses above about 15 M$_{\odot}$. Although observations available at the moment cannot decide what is the best choice between the core sizes given by the TAMS and the velocity drop, we discuss different methods to get out of this dilemma. At the moment, comparisons with eclipsing binaries seem to favor the solution given by the velocity drop. While we confirm the need for larger convective cores at higher masses, we find tensions in-between different methods for stars more massive than 15 M$_{\odot}$. The use of single-aged stellar populations (non-interacting binaries or stellar clusters) would be a great asset to resolve this tension.


title: An unfitted radial basis function generated finite difference method applied to thoracic diaphragm simulations

abstract: The thoracic diaphragm is the muscle that drives the respiratory cycle of a human being. Using a system of partial differential equations (PDEs) that models linear elasticity we compute displacements and stresses in a two-dimensional cross section of the diaphragm in its contracted state. The boundary data consists of a mix of displacement and traction conditions. If these are imposed as they are, and the conditions are not compatible, this leads to reduced smoothness of the solution. Therefore, the boundary data is first smoothed using the least-squares radial basis function generated finite difference (RBF-FD) framework. Then the boundary conditions are reformulated as a Robin boundary condition with smooth coefficients. The same framework is also used to approximate the boundary curve of the diaphragm cross section based on data obtained from a slice of a computed tomography (CT) scan. To solve the PDE we employ the unfitted least-squares RBF-FD method. This makes it easier to handle the geometry of the diaphragm, which is thin and non-convex. We show numerically that our solution converges with high-order towards a finite element solution evaluated on a fine grid. Through this simplified numerical model we also gain an insight into the challenges associated with the diaphragm geometry and the boundary conditions before approaching a more complex three-dimensional model.


title: High bandwidth temporal RF photonic signal processing with Kerr micro-combs: integration, fractional differentiation and Hilbert transforms

abstract: Integrated Kerr micro-combs, a powerful source of many wavelengths for photonic RF and microwave signal processing, are particularly useful for transversal filter systems. They have many advantages including a compact footprint, high versatility, large numbers of wavelengths, and wide bandwidths. We review recent progress on photonic RF and microwave high bandwidth temporal signal processing based on Kerr micro-combs with spacings from 49-200GHz. We cover integral and fractional Hilbert transforms, differentiators as well as integrators. The potential of optical micro-combs for RF photonic applications in functionality and ability to realize integrated solutions is also discussed.


title: Broadband super-resolution Terahertz Time domain spectroscopy applied to Gas analysis

abstract: Terahertz (THz) Time domain spectroscopy (THz-TDS) is a broadband spectroscopic technique spreading its uses in multiple fields: in science from material science to biology, in industry where it measures the thickness of a paint layer during the painting operation. Using such practical commercial apparatus with broad spectrum for gas spectroscopy could be a major asset for air quality monitoring and tracking of atmospheric composition. However, gas spectroscopy needs high resolution and the usual approach in THz-TDS, where the recorded time trace is Fourier transform, suffers from resolution limitation due to the size of the delay line in the system. In this letter, we introduce the concept of constraint reconstruction for super-resolution spectroscopy based on the modeling of the spectroscopic lines in a sparse spectrum. Light molecule gas typically shows sparse and narrow lines on a broad spectrum and we propose an algorithm reconstructing these lines with a resolution improvement of 10 the ultimate resolution reachable by the apparatus. We envision the proposed technique to lead to broadband, selective, rapid and cheap gas monitoring applications.


title: The Design of Dual Band Stacked Metasurfaces Using Integral Equations

abstract: An integral equation-based approach for the design of dual band stacked metasurfaces which are invariant in one-dimension is presented. The stacked metasurface will generate collimated beams at desired angles in each band upon reflection. The conductor-backed stacked metasurface consists of two metasurfaces (a patterned metallic cladding supported by a dielectric spacer) stacked one upon the other. The stacked metasurface is designed in three phases. First the patterned metallic cladding of each metasurface is homogenized and modeled as an inhomogeneous impedance sheet. An Electric Field Integral Equation (EFIE) is written to model the mutual coupling between the homogenized elements within each metasurface, and from metasurface to metasurface. The EFIE is transformed into matrix equations by the method of moments. The nonlinear matrix equations are solved at both bands iteratively resulting in dual band complex-valued impedance sheets. In the second phase, optimization is applied to transform these complex-valued impedance sheets into purely reactive sheets suitable for printed circuit board fabrication by introducing surface waves. In the third phase, the metallic claddings of each metasurface are patterned for full-wave simulation of the dual band stacked metasurface. Using this approach, two dual band stacked metasurfaces are designed

categories: physics.app-ph cs.SY eess.SY physics.optics


title: Control Barrier Functions in Sampled-Data Systems

abstract: This paper presents conditions for ensuring forward invariance of safe sets under sampled-data system dynamics with piecewise-constant controllers and fixed time-steps. First, we introduce two different metrics to compare the conservativeness of sufficient conditions on forward invariance under piecewise-constant controllers. Then, we propose three approaches for guaranteeing forward invariance, two motivated by continuous-time barrier functions, and one motivated by discrete-time barrier functions. All proposed conditions are control affine, and thus can be incorporated into quadratic programs for control synthesis. We show that the proposed conditions are less conservative than those in earlier studies, and show via simulation how this enables the use of barrier functions that are impossible to implement with the desired time-step using existing methods.


title: An Application-Driven Conceptualization of Corner Cases for Perception in Highly Automated Driving

abstract: Systems and functions that rely on machine learning (ML) are the basis of highly automated driving. An essential task of such ML models is to reliably detect and interpret unusual, new, and potentially dangerous situations. The detection of those situations, which we refer to as corner cases, is highly relevant for successfully developing, applying, and validating automotive perception functions in future vehicles where multiple sensor modalities will be used. A complication for the development of corner case detectors is the lack of consistent definitions, terms, and corner case descriptions, especially when taking into account various automotive sensors. In this work, we provide an application-driven view of corner cases in highly automated driving. To achieve this goal, we first consider existing definitions from the general outlier, novelty, anomaly, and out-of-distribution detection to show relations and differences to corner cases. Moreover, we extend an existing camera-focused systematization of corner cases by adding RADAR (radio detection and ranging) and LiDAR (light detection and ranging) sensors. For this, we describe an exemplary toolchain for data acquisition and processing, highlighting the interfaces of the corner case detection. We also define a novel level of corner cases, the method layer corner cases, which appear due to uncertainty inherent in the methodology or the data distribution.


title: Episodic memory governs choices: An RNN-based reinforcement learning model for decision-making task

abstract: Typical methods to study cognitive function are to record the electrical activities of animal neurons during the training of animals performing behavioral tasks. A key problem is that they fail to record all the relevant neurons in the animal brain. To alleviate this problem, we develop an RNN-based Actor-Critic framework, which is trained through reinforcement learning (RL) to solve two tasks analogous to the monkeys' decision-making tasks. The trained model is capable of reproducing some features of neural activities recorded from animal brain, or some behavior properties exhibited in animal experiments, suggesting that it can serve as a computational platform to explore other cognitive functions. Furthermore, we conduct behavioral experiments on our framework, trying to explore an open question in neuroscience: which episodic memory in the hippocampus should be selected to ultimately govern future decisions. We find that the retrieval of salient events sampled from episodic memories can effectively shorten deliberation time than common events in the decision-making process. The results indicate that salient events stored in the hippocampus could be prioritized to propagate reward information, and thus allow decision-makers to learn a strategy faster.


title: Uncertainty Quantification and Software Risk Analysis for Digital Twins in the Nearly Autonomous Management and Control Systems: A Review

abstract: A nearly autonomous management and control (NAMAC) system is designed to furnish recommendations to operators for achieving particular goals based on NAMAC's knowledge base. As a critical component in a NAMAC system, digital twins (DTs) are used to extract information from the knowledge base to support decision-making in reactor control and management during all modes of plant operations. With the advancement of artificial intelligence and data-driven methods, machine learning algorithms are used to build DTs of various functions in the NAMAC system. To evaluate the uncertainty of DTs and its impacts on the reactor digital instrumentation and control systems, uncertainty quantification (UQ) and software risk analysis is needed. As a comprehensive overview of prior research and a starting point for new investigations, this study selects and reviews relevant UQ techniques and software hazard and software risk analysis methods that may be suitable for DTs in the NAMAC system.


title: Conformal Transformation Electromagnetics Based on Schwarz-Christoffel Mapping for the Synthesis of Doubly Connected Metalenses

abstract: An innovative transformation electromagnetics (TE) paradigm, which leverages on the Schwarz-Christoffel (SC) theorem, is proposed to design effective and realistic field manipulation devices (FMDs). Thanks to the conformal property, such a TE design method allows one to considerably mitigate the anisotropy of the synthesized metalenses (i.e., devices with artificially engineered materials covering an antenna to modify its radiation features) with respect to those yielded by the competitive state-of-the-art TE techniques. Moreover, devices with doubly connected contours, thus including masts with arbitrary sections and lenses with holes/forbidden regions in which the material properties cannot be controlled, can be handled. A set of numerical experiments is presented to assess the features of the proposed method in terms of field-manipulation capabilities and complexity of the lens material in a comparative fashion.


title: Optical Frequency Comb Fourier Transform Spectroscopy of $^{14}$N$_2$$^{16}$O at 7.8 {\mu}m

abstract: We use a Fourier transform spectrometer based on a compact mid-infrared difference frequency generation comb source to perform broadband high-resolution measurements of nitrous oxide, $^{14}$N$_2$$^{16}$O, and retrieve line center frequencies of the $\nu$$_1$ fundamental band and the $\nu$$_1$ + $\nu$$_2$ - $\nu$$_2$ hot band. The spectrum spans 90 cm$^{-1}$ around 1285 cm$^{-1}$ with a sample point spacing of 3 ${\times}$ 10$^{-4}$ cm$^{-1}$ (9 MHz). We report line positions of 72 lines in the $\nu$$_1$ fundamental band between P(37) and R(38), and 112 lines in the $\nu$$_1$ + $\nu$$_2$ - $\nu$$_2$ hot band (split into two components with e/f rotationless parity) between P(34) and R(33), with uncertainties in the range of 90-600 kHz. We derive upper state constants of both bands from a fit of the effective ro-vibrational Hamiltonian to the line center positions. For the fundamental band, we observe excellent agreement in the retrieved line positions and upper state constants with those reported in a recent study by AlSaif et al. using a comb-referenced quantum cascade laser [J Quant Spectrosc Radiat Transf, 2018;211:172-178]. We determine the origin of the hot band with precision one order of magnitude better than previous work based on FTIR measurements by Toth [http://mark4sun.jpl.nasa.gov/n2o.html], which is the source of the HITRAN2016 data for these bands.


title: Electromagnetic Isolation Induced by Time-Varying Metasurfaces: Non-Reciprocal Bragg Grating

abstract: In this letter, we propose a magnet-less non-reciprocal isolating system based on time-varying metasurfaces. Two parallel time-varying metasurfaces, one for frequency up-conversion and one for down-conversion by the same amount, are used for realizing a region of space where incident waves from opposite directions experience an opposite Doppler frequency shift. As a result, any device within this region becomes sensitive to the illumination direction, exhibiting a different scattering response from opposite directions and thus breaking reciprocity. Very importantly, thanks to the opposite frequency shift of the metasurfaces, the frequency of the transmitted electromagnetic field is the same as for the incident one. Here, we demonstrate this general approach by using a Bragg grating as the device between the time-varying metasurfaces. The combined structure of the metasurfaces and the grating exhibits different transmission and reflection properties for opposite illumination direction, thereby realizing an isolator. More broadly, this letter presents a strategy for converting any conventional electromagnetic device to a non-reciprocal one by placing it between two time-varying metasurfaces. This approach opens the door to several new non-reciprocal components based on thin and lightweight metasurfaces, which are simpler to realize compared to their volumetric counterparts.


title: Low half-wave-voltage, ultra-high bandwidth thin-film LiNbO3 modulator based on hybrid waveguide and periodic capacitively loaded electrodes

abstract: A novel thin-film LiNbO3 (TFLN) electro-optic modulator is proposed and demonstrated. LiNbO3-silica hybrid waveguide is adopted to maintain low optical loss for an electrode spacing as narrow as 3 {\mu}m, resulting in a record low half-wave-voltage length product of only 1.7 V*cm. Capacitively loaded traveling-wave electrodes (CL-TWEs) are employed to reduce the microwave loss, while quartz substrate is used in place of silicon substrate to achieve velocity matching. The fabricated TFLN modulator with a 5-mm-long modulation region exhibits a half-wave-voltage of 3.4 V and merely 1.3 dB roll-off in electro-optic response up to 67 GHz, and a 3-dB modulation bandwidth over 110 GHz is predicted.


title: Boundary-Induced Embedded Eigenstate in a Single Resonator for Advanced Sensing

abstract: Electromagnetic embedded eigenstates, also known as bound states in the continuum (BICs), hold a great potential for applications in sensing, lasing, enhanced nonlinearities and energy harvesting. However, their demonstrations so far have been limited to large-area periodic arrays of suitably tailored elements, with fundamental restrictions on the overall footprint and performance in presence of inevitable disorder. In this work, we demonstrate a BIC localized in a single subwavelength resonator obtained by suitably tailoring the boundaries around it, enabling a new degree of control for on-demand symmetry breaking. We experimentally demonstrate how boundary-induced BICs open exciting opportunities for sensing by tracing the dissolution of NaCl in water, determining evaporation rates of distilled and saltwater with a resolution of less than 1e-6 L using a tabletop experimental setup.


title: Bilayer Photonic Graphene

abstract: Drawing inspiration from bilayer graphene, this paper introduces its photonic analog comprising two stacked graphene-like photonic crystals, that are coupled in the near-field through spoof surface plasmons. Beyond the twist degree of freedom that can radically alter the band structure of the bilayer photonic graphene, the photonic dispersion can be also tailored via the interlayer coupling which exhibits an exponential dependence on the distance between the two photonic crystals. We theoretically, numerically, and experimentally characterize the band structures of AA- and AB-stacked bilayer photonic graphene, as well as for twisted bilayer photonic graphene with even and odd sublattice exchange symmetries. Furthermore, we numerically predict the existence of magic angles in bilayer photonic graphene, which are associated with ultra-flat bands resulted from interlayer hybridization. Finally, we demonstrate that the bilayer photonic graphene at a particular twist angle satisfying even sublattice exchange symmetry is a high-order photonic topological insulator. The proposed bilayer photonic graphene could constitute a useful platform for identifying new quantum materials and inspiring next-generation photonic devices with new degrees of freedom and emerging functionality.


title: Functional surface formation by efficient laser ablation using single-pulse and burst-modes

abstract: Surfaces inspired by nature and their replication find great interest in science, technology, and medicine due to their unique functional properties. This research aimed to develop an efficient laser milling technology using single-pulse- and burst-modes of irradiation to replicate bio-inspired structures over large areas at high speed. The ability to form the trapezoidal-riblets inspired by shark skin at high production speeds while maintaining the lowest possible surface roughness was demonstrated.


title: On the Occasional Exactness of the Distributional Transform Approximation for Direct Gaussian Copula Models with Discrete Margins

abstract: The direct Gaussian copula model with discrete marginal distributions is an appealing data-analytic tool but poses difficult computational challenges due to its intractable likelihood. A number of approximations/surrogates for the likelihood have been proposed, including the continuous extension-based approximation (CE) and the distributional transform-based approximation (DT). The continuous extension approach is exact up to Monte Carlo error but does not scale well computationally. The distributional transform approach permits efficient computation but offers no theoretical guarantee that it is exact. In practice, though, the distributional transform-based approximate likelihood is so very nearly exact for some variants of the model as to permit genuine maximum likelihood or Bayesian inference. We demonstrate the exactness of the distributional transform-based objective function for two interesting variants of the model, and propose a quantity that can be used to assess exactness for experimentally observed datasets. Said diagnostic will permit practitioners to determine whether genuine Bayesian inference or ordinary maximum likelihood inference using the DT-based likelihood is possible for a given dataset.


title: Optimal Stationary State Estimation Over Multiple Markovian Packet Drop Channels

abstract: In this paper, we investigate the state estimation problem over multiple Markovian packet drop channels. In this problem setup, a remote estimator receives measurement data transmitted from multiple sensors over individual channels. By the method of Markovian jump linear systems, an optimal stationary estimator that minimizes the error variance in the steady state is obtained, based on the mean-square (MS) stabilizing solution to the coupled algebraic Riccati equations. An explicit necessary and sufficient condition is derived for the existence of the MS stabilizing solution, which coincides with that of the standard Kalman filter. More importantly, we provide a sufficient condition under which the MS detectability with multiple Markovian packet drop channels can be decoupled, and propose a locally optimal stationary estimator but computationally more tractable. Analytic sufficient and necessary MS detectability conditions are presented for the decoupled subsystems subsequently. Finally, numerical simulations are conducted to illustrate the results on the MS stabilizing solution, the MS detectability, and the performance of the optimal and locally optimal stationary estimators.


title: Addressing Research Software Sustainability via Institutes

abstract: Research software is essential to modern research, but it requires ongoing human effort to sustain: to continually adapt to changes in dependencies, to fix bugs, and to add new features. Software sustainability institutes, amongst others, develop, maintain, and disseminate best practices for research software sustainability, and build community around them. These practices can both reduce the amount of effort that is needed and create an environment where the effort is appreciated and rewarded. The UK SSI is such an institute, and the US URSSI and the Australian AuSSI are planning to become institutes, and this extended abstract discusses them and the strengths and weaknesses of this approach.


title: Experimental hierarchy and optimal robustness of quantum correlations of two-qubit states with controllable white noise

abstract: We demonstrate a hierarchy of various classes of quantum correlations on experimentally prepared two-qubit Werner-like states with controllable white noise. Werner states, which are white-noise-affected Bell states, are prototypal examples for studying such a hierarchy as a function of the amount of white noise. We experimentally generated Werner states and their generalizations, i.e., partially entangled pure states affected by white noise. These states enabled us to study the hierarchy of the following classes of correlations: separability, entanglement, steering in three- and two-measurement scenarios, and Bell nonlocality. We show that the generalized Werner states (GWSs) reveal fundamentally different aspects of the hierarchy compared to the Werner states. In particular, we find five different parameter regimes of the GWSs, including those steerable in a two-measurement scenario but not violating Bell inequalities. This regime cannot be observed for the usual Werner states. Furthermore, we find threshold curves separating different regimes of the quantum correlations and find the optimal states which allow for the largest amount of white noise, which does not destroy their specific quantum correlations (e.g., unsteerable entanglement). Thus, we could identify the optimal Bell-nondiagonal GWSs, which are, for this specific meaning, more robust against white noise compared to the Bell-diagonal GWSs (i.e., Werner states).


title: Signal-level Fusion for Indexing and Retrieval of Facial Biometric Data

abstract: The growing scope, scale, and number of biometric deployments around the world emphasise the need for research into technologies facilitating efficient and reliable biometric identification queries. This work presents a method of indexing biometric databases, which relies on signal-level fusion of facial images (morphing) to create a multi-stage data-structure and retrieval protocol. By successively pre-filtering the list of potential candidate identities, the proposed method makes it possible to reduce the necessary number of biometric template comparisons to complete a biometric identification transaction. The proposed method is extensively evaluated on publicly available databases using open-source and commercial off-the-shelf recognition systems. The results show that using the proposed method, the computational workload can be reduced down to around 30%, while the biometric performance of a baseline exhaustive search-based retrieval is fully maintained, both in closed-set and open-set identification scenarios.


title: Differential invariants of Kundt spacetimes

abstract: We find generators for the algebra of rational differential invariants for general and degenerate Kundt spacetimes and relate this to other approaches to the equivalence problem for Lorentzian metrics. Special attention is given to dimensions three and four.


title: Local heat emission due to unidirectional spin-wave heat conveyer effect observed by lock-in thermography

abstract: Lock-in thermography measurements were performed to reveal heat source distribution induced by the unidirectional spin-wave heat conveyer effect (USHCE) of magnetostatic surface spin waves. When the magnetostatic surface spin waves are excited in an yttrium iron garnet slab, the lock-in thermography images show spatially biased sharp and complicated heating patterns, indicating the importance of edge spin-wave dynamics for USHCE. The accessibility to the local heat emission properties allows us to clarify a capability of remote heating realized by USHCE; it can transfer energy for heating even through a macro-scale air gap between two magnetic materials owing to the long-range dipole-dipole coupling.


title: Multi-rate Control Design under Input Constraints via Fixed-Time Barrier Functions

abstract: In this paper, we introduce the notion of periodic safety, which requires that the system trajectories periodically visit a subset of a forward-invariant safe set, and utilize it in a multi-rate framework where a high-level planner generates a reference trajectory that is tracked by a low-level controller under input constraints. We introduce the notion of fixed-time barrier functions which is leveraged by the proposed low-level controller in a quadratic programming framework. Then, we design a model predictive control policy for high-level planning with a bound on the rate of change for the reference trajectory to guarantee that periodic safety is achieved. We demonstrate the effectiveness of the proposed strategy on a simulation example, where the proposed fixed-time stabilizing low-level controller shows successful satisfaction of control objectives, whereas an exponentially stabilizing low-level controller fails.


title: Measure equivalence rigidity of $\mathrm{Out}(F_N)$

abstract: We prove that for every $N\ge 3$, the group $\mathrm{Out}(F_N)$ of outer automorphisms of a free group of rank $N$ is superrigid from the point of view of measure equivalence: any countable group that is measure equivalent to $\mathrm{Out}(F_N)$, is in fact virtually isomorphic to $\mathrm{Out}(F_N)$. We introduce three new constructions of canonical splittings associated to a subgroup of $\mathrm{Out}(F_N)$ of independent interest. They encode respectively the collection of invariant free splittings, invariant cyclic splittings, and maximal invariant free factor systems. Our proof also relies on the following improvement of an amenability result by Bestvina and the authors: given a free factor system $\mathcal{F}$ of $F_N$, the action of $\mathrm{Out}(F_N,\mathcal{F})$ (the subgroup of $\mathrm{Out}(F_N)$ that preserves $\mathcal{F}$) on the space of relatively arational trees with amenable stabilizer is a Borel amenable action.


title: Bayesian Meta-Learning for Few-Shot Policy Adaptation Across Robotic Platforms

abstract: Reinforcement learning methods can achieve significant performance but require a large amount of training data collected on the same robotic platform. A policy trained with expensive data is rendered useless after making even a minor change to the robot hardware. In this paper, we address the challenging problem of adapting a policy, trained to perform a task, to a novel robotic hardware platform given only few demonstrations of robot motion trajectories on the target robot. We formulate it as a few-shot meta-learning problem where the goal is to find a meta-model that captures the common structure shared across different robotic platforms such that data-efficient adaptation can be performed. We achieve such adaptation by introducing a learning framework consisting of a probabilistic gradient-based meta-learning algorithm that models the uncertainty arising from the few-shot setting with a low-dimensional latent variable. We experimentally evaluate our framework on a simulated reaching and a real-robot picking task using 400 simulated robots generated by varying the physical parameters of an existing set of robotic platforms. Our results show that the proposed method can successfully adapt a trained policy to different robotic platforms with novel physical parameters and the superiority of our meta-learning algorithm compared to state-of-the-art methods for the introduced few-shot policy adaptation problem.


title: Quantum Zeno and Anti-Zeno probes of noise correlations in photon polarisation

abstract: We experimentally demonstrate, for the first time, noise diagnostics by repeated quantum measurements. Specifically, we establish the ability of a single photon, subjected to random polarisation noise, to diagnose non-Markovian temporal correlations of such a noise process. In the frequency domain, these noise correlations correspond to colored noise spectra, as opposed to the ones related to Markovian, white noise. Both the noise spectrum and its corresponding temporal correlations are diagnosed by probing the photon by means of frequent, (partially-)selective polarisation measurements. Our main result is the experimental demonstration that noise with positive temporal correlations corresponds to our single photon undergoing a dynamical regime enabled by the quantum Zeno effect (QZE), while noise characterized by negative (anti-) correlations corresponds to regimes associated with the anti-Zeno effect (AZE). This demonstration opens the way to a new kind of noise spectroscopy based on QZE and AZE in photon (or other single-particle) state probing.


title: Implementing Automated Market Makers with Constant Circle

abstract: This paper describe the implementation details of constant ellipse based automated market makers (CoinSwap). A CoinSwap prototype has been implemented at http://coinswapapp.io/ and the source codes are available at https://github.com/coinswapapp/


title: Analyzing the Influence of Dataset Composition for Emotion Recognition

abstract: Recognizing emotions from text in multimodal architectures has yielded promising results, surpassing video and audio modalities under certain circumstances. However, the method by which multimodal data is collected can be significant for recognizing emotional features in language. In this paper, we address the influence data collection methodology has on two multimodal emotion recognition datasets, the IEMOCAP dataset and the OMG-Emotion Behavior dataset, by analyzing textual dataset compositions and emotion recognition accuracy. Experiments with the full IEMOCAP dataset indicate that the composition negatively influences generalization performance when compared to the OMG-Emotion Behavior dataset. We conclude by discussing the impact this may have on HRI experiments.


title: Don't Forget to Sign the Gradients!

abstract: Engineering a top-notch deep learning model is an expensive procedure that involves collecting data, hiring human resources with expertise in machine learning, and providing high computational resources. For that reason, deep learning models are considered as valuable Intellectual Properties (IPs) of the model vendors. To ensure reliable commercialization of deep learning models, it is crucial to develop techniques to protect model vendors against IP infringements. One of such techniques that recently has shown great promise is digital watermarking. However, current watermarking approaches can embed very limited amount of information and are vulnerable against watermark removal attacks. In this paper, we present GradSigns, a novel watermarking framework for deep neural networks (DNNs). GradSigns embeds the owner's signature into the gradient of the cross-entropy cost function with respect to inputs to the model. Our approach has a negligible impact on the performance of the protected model and it allows model vendors to remotely verify the watermark through prediction APIs. We evaluate GradSigns on DNNs trained for different image classification tasks using CIFAR-10, SVHN, and YTF datasets. Experimental results show that GradSigns is robust against all known counter-watermark attacks and can embed a large amount of information into DNNs.


title: Some Properties and Applications of Burr III-Weibull Distribution

abstract: In this paper, we introduce a new distribution called Burr III-Weibull(BW) distribution using the concept of competing risk. We derive moments, conditional moments, mean deviation and quantiles of the proposed distribution. Also the Renyi's entropy and order statistics of the distribution are obtained. Estimation of parameters of the distribution is performed via maximum likelihood method. A simulation study is performed to validate the maximum likelihood estimator (MLE). A real practical data set is analyzed for illustration.


title: Semi-Supervised Federated Peer Learning for Skin Lesion Classification

abstract: Globally, Skin carcinoma is among the most lethal diseases. Millions of people are diagnosed with this cancer every year. Sill, early detection can decrease the medication cost and mortality rate substantially. The recent improvement in automated cancer classification using deep learning methods has reached a human-level performance requiring a large amount of annotated data assembled in one location, yet, finding such conditions usually is not feasible. Recently, federated learning (FL) has been proposed to train decentralized models in a privacy-preserved fashion depending on labeled data at the client-side, which is usually not available and costly. To address this, we propose \verb!FedPerl!, a semi-supervised federated learning method. Our method is inspired by peer learning from educational psychology and ensemble averaging from committee machines. FedPerl builds communities based on clients' similarities. Then it encourages communities members to learn from each other to generate more accurate pseudo labels for the unlabeled data. We also proposed the peer anonymization (PA) technique to anonymize clients. As a core component of our method, PA is orthogonal to other methods without additional complexity and reduces the communication cost while enhancing performance. Finally, we propose a dynamic peer-learning policy that controls the learning stream to avoid any degradation in the performance, especially for individual clients. Our experimental setup consists of 71,000 skin lesion images collected from 5 publicly available datasets. We test our method in four different scenarios in SSFL. With few annotated data, FedPerl is on par with a state-of-the-art method in skin lesion classification in the standard setup while outperforming SSFLs and the baselines by 1.8% and 15.8%, respectively. Also, it generalizes better to unseen clients while being less sensitive to noisy ones.


title: Abstraction and Symbolic Execution of Deep Neural Networks with Bayesian Approximation of Hidden Features

abstract: Intensive research has been conducted on the verification and validation of deep neural networks (DNNs), aiming to understand if, and how, DNNs can be applied to safety critical applications. However, existing verification and validation techniques are limited by their scalability, over both the size of the DNN and the size of the dataset. In this paper, we propose a novel abstraction method which abstracts a DNN and a dataset into a Bayesian network (BN). We make use of dimensionality reduction techniques to identify hidden features that have been learned by hidden layers of the DNN, and associate each hidden feature with a node of the BN. On this BN, we can conduct probabilistic inference to understand the behaviours of the DNN processing data. More importantly, we can derive a runtime monitoring approach to detect in operational time rare inputs and covariate shift of the input data. We can also adapt existing structural coverage-guided testing techniques (i.e., based on low-level elements of the DNN such as neurons), in order to generate test cases that better exercise hidden features. We implement and evaluate the BN abstraction technique using our DeepConcolic tool available at https://github.com/TrustAI/DeepConcolic.


title: FedDis: Disentangled Federated Learning for Unsupervised Brain Pathology Segmentation

abstract: In recent years, data-driven machine learning (ML) methods have revolutionized the computer vision community by providing novel efficient solutions to many unsolved (medical) image analysis problems. However, due to the increasing privacy concerns and data fragmentation on many different sites, existing medical data are not fully utilized, thus limiting the potential of ML. Federated learning (FL) enables multiple parties to collaboratively train a ML model without exchanging local data. However, data heterogeneity (non-IID) among the distributed clients is yet a challenge. To this end, we propose a novel federated method, denoted Federated Disentanglement (FedDis), to disentangle the parameter space into shape and appearance, and only share the shape parameter with the clients. FedDis is based on the assumption that the anatomical structure in brain MRI images is similar across multiple institutions, and sharing the shape knowledge would be beneficial in anomaly detection. In this paper, we leverage healthy brain scans of 623 subjects from multiple sites with real data (OASIS, ADNI) in a privacy-preserving fashion to learn a model of normal anatomy, that allows to segment abnormal structures. We demonstrate a superior performance of FedDis on real pathological databases containing 109 subjects; two publicly available MS Lesions (MSLUB, MSISBI), and an in-house database with MS and Glioblastoma (MSI and GBI). FedDis achieved an average dice performance of 0.38, outperforming the state-of-the-art (SOTA) auto-encoder by 42% and the SOTA federated method by 11%. Further, we illustrate that FedDis learns a shape embedding that is orthogonal to the appearance and consistent under different intensity augmentations.


title: DOPE: D-Optimal Pooling Experimental design with application for SARS-CoV-2 screening

abstract: Testing individuals for the presence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the pathogen causing the coronavirus disease 2019 (COVID-19), is crucial for curtailing transmission chains. Moreover, rapidly testing many potentially infected individuals is often a limiting factor in controlling COVID-19 outbreaks. Hence, pooling strategies, wherein individuals are grouped and tested simultaneously, are employed. We present a novel pooling strategy that implements D-Optimal Pooling Experimental design (DOPE). DOPE defines optimal pooled tests as those maximizing the mutual information between data and infection states. We estimate said mutual information via Monte-Carlo sampling and employ a discrete optimization heuristic for maximizing it. DOPE outperforms common pooling strategies both in terms of lower error rates and fewer tests utilized. DOPE holds several additional advantages: it provides posterior distributions of the probability of infection, rather than only binary classification outcomes; it naturally incorporates prior information of infection probabilities and test error rates; and finally, it can be easily extended to include other, newly discovered information regarding COVID-19. Hence, we believe that implementation of Bayesian D-optimal experimental design holds a great promise for the efforts of combating COVID-19 and other future pandemics.


title: SALT observations of the supernova remnant MCSNR J0127-7332 and its associated Be X-ray binary SXP 1062 in the SMC

abstract: We report the results of optical spectroscopy of the Small Magellanic Cloud supernova remnant (SNR) MCSNR J0127-7332 and the mass donor Be star, 2dFS 3831, in its associated high-mass X-ray binary SXP 1062 carried out with the Southern African Large Telescope (SALT). Using high-resolution long-slit spectra, we measured the expansion velocity of the SNR shell of \approx 140 km/s, indicating that MCSNR J0127-7332 is in the radiative phase. We found that the observed line ratios in the SNR spectrum can be understood if the local interstellar medium is ionized by 2dFS 3831 and/or OB stars around the SNR. We propose that MCSNR J0127-7332 is the result of supernova explosion within a bubble produced by the stellar wind of the supernova progenitor and that the bubble was surrounded by a massive shell at the moment of supernova explosion. We estimated the age of MCSNR J0127-7332 to be \la 10 000 yr. We found that the spectrum of 2dFS 3831 changes with orbital phase. Namely, the equivalent width of the Halpha emission line decreased by \approx 40 per cent in \approx 130 d after periastron passage of the neutron star and then almost returned to its original value in the next \approx 100 d. Also, the spectrum of 2dFS 3831 obtained closest to the periastron epoch (about three weeks after the periastron) shows a noticeable emission line of He II \lambda 4686, which disappeared in the next about two weeks. We interpret these changes as a result of the temporary perturbation and heating of the disk as the neutron star passes through it.


title: Energy-efficient Task Offloading for Relay Aided Mobile Edge Computing under Sequential Task Dependency

abstract: In this paper, we study a mobile edge computing (MEC) system in which the mobile device is assisted by a base station (BS) and a cooperative node. The mobile device has sequential tasks to complete, whereas the cooperative node assists the mobile device on both task offloading and task computation. In specific, two cases are investigated, which are 1) the cooperative node has no tasks to complete itself, and 2) the cooperative node has tasks to complete itself. Our target is to minimize the total energy consumption of the mobile device and the cooperative node through optimizing the transmit duration in task offloading, CPU frequency in task computing along with the task index to offload in the sequential tasks. In the first case, we decompose the mixed-integer non-convex problem into two levels. In the lower level problem, thanks to the convexity, Karush-Kuhn-Tucker (KKT) conditions are utilized to simplify the problem, which is then solved with bisection search. In the upper level problem, to find solution of the task index to offload, rather than utilizing traversal method, we develop a monotonic condition to simplify the searching process. In the second case, in order to guarantee the successful computation of the mobile device and cooperative node, the uploading transmission is classified into three schemes. Within each scheme, the non-convex problem is decomposed. In the lower level problem, semi-closed solution is found by Lagrangian dual method. In the upper level problem, traversal method is applied to find the optimal offloading index.


title: Social contagion on higher-order structures

abstract: In this Chapter, we discuss the effects of higher-order structures on SIS-like processes of social contagion. After a brief motivational introduction where we illustrate the standard SIS process on networks and the difference between simple and complex contagions, we introduce spreading processes on higher-order structures starting from the most general formulation on hypergraphs and then moving to several mean-field and heterogeneous mean-field approaches. The results highlight the rich phenomenology brought by taking into account higher-order contagion effects: both continuous and discontinuous transitions are observed, and critical mass effects emerge. We conclude with a short discussion on the theoretical results regarding the nature of the epidemic transition and the general need for data to validate these models.


title: Characterising different communities of Twitter users: Migrants and natives

abstract: Today, many users are actively using Twitter to express their opinions and to share information. Thanks to the availability of the data, researchers have studied behaviours and social networks of these users. International migration studies have also benefited from this social media platform to improve migration statistics. Although diverse types of social networks have been studied so far on Twitter, social networks of migrants and natives have not been studied before. This paper aims to fill this gap by studying characteristics and behaviours of migrants and natives on Twitter. To do so, we perform a general assessment of features including profiles and tweets, and an extensive network analysis on the network. We find that migrants have more followers than friends. They have also tweeted more despite that both of the groups have similar account ages. More interestingly, the assortativity scores showed that users tend to connect based on nationality more than country of residence, and this is more the case for migrants than natives. Furthermore, both natives and migrants tend to connect mostly with natives.


title: Destructive Controlled-Phase Gate Using Linear Optics

abstract: Knill, Laflamme, and Milburn [Nature 409, 46 (2001)] showed that linear optics techniques could be used to implement a nonlinear sign gate. They also showed that two of their nonlinear sign gates could be combined to implement a controlled-phase gate, which has a number of practical applications. Here we describe an alternative implementation of a controlled-phase gate that only requires the use of a single nonlinear sign gate. This gives a much higher average probability of success when the required ancilla photons are generated using heralding techniques. This implementation of a controlled-phase gate destroys the control qubit, which is acceptable in a number of applications where the control qubit would have been destroyed in any event, such as in a postselection process.


title: Reconfigurable and approximate computing for video coding

abstract: The Chapter begins with a discussion of the constraints and needs of video coding systems. The lack in flexibility of traditional monolithic codec specifications, not suitable to model commonalities among codecs and foster reusability among successive codec generations/updates, was the main trigger for the development of a new standard initiative within the ISO/IEC MPEG committee, called reconfigurable video coding (RVC). The MPEG-RVC framework exploits the dataflow nature behind video coding to foster flexible and reconfigurable codec design, as well as to support dynamic reconfiguration. The Chapter goes on to consider that the inherent resiliency of various functional blocks (like motion estimation in the high-efficiency video coding, HEVC) and the varying levels of user perception make video coding suitable to apply approximate computing techniques. Approximate computing, if properly supported at design time, allows achieving run-time trade-offs, representing a new direction in hardware-software codesign research. The main assumption behind approximate computing, exploited within video coding, is that the degree of accuracy (in this case during codec execution) is not required to be the same all the time. The final part of the Chapter attempts to put together the concepts addressed and remarks on which are, in the authors' opinion, some interesting research directions.


title: Ground-SLAM: Ground Constrained LiDAR SLAM for Structured Multi-Floor Environments

abstract: This paper proposes a 3D LiDAR SLAM algorithm named Ground-SLAM, which exploits grounds in structured multi-floor environments to compress the pose drift mainly caused by LiDAR measurement bias. Ground-SLAM is developed based on the well-known pose graph optimization framework. In the front-end, motion estimation is conducted using LiDAR Odometry (LO) with a novel sensor-centric sliding map introduced, which is maintained by filtering out expired features based on the model of error propagation. At each key-frame, the sliding map is recorded as a local map. The ground nearby is extracted and modelled as an infinite planar landmark in the form of Closest Point (CP) parameterization. Then, ground planes observed at different key-frames are associated, and the ground constraints are fused into the pose graph optimization framework to compress the pose drift of LO. Finally, loop-closure detection is carried out, and the residual error is jointly minimized, which could lead to a globally consistent map. Experimental results demonstrate superior performances in the accuracy of the proposed approach.


title: Quantum (dis)charge of black hole interiors

abstract: We analyze the "vacuum" polarization induced by a quantum charged scalar field near the inner horizon of a charged (Reissner-Nordstr\"om-de Sitter) black hole in quantum states that start out as regular states near an initial Cauchy surface. Contrary to the outer (i.e. event-) horizon, where polarization effects lead to a discharge, we find that near an inner horizon, the transversal component of the expected current density can have either sign depending on the black hole and field parameters. Thus, the inner horizon can be charged or discharged. But we find that it is always discharged close to extremality thus driving the black hole interior away from this critical point. Furthermore, we find that quantum effects dominate in that the strength of the blow up of the quantum current at the inner horizon is state-independent and stronger than that of the current of a classical solution.


title: Bruhat intervals, subword complexes and brick polyhedra for finite Coxeter groups

abstract: We study the interplay between the discrete geometry of Bruhat poset intervals and subword complexes of finite Coxeter systems. We establish connections between the cones generated by cover labels for Bruhat intervals and of root configurations for subword complexes, culminating in the notion of brick polyhedra for general subword complexes.


title: Golem: An algorithm for robust experiment and process optimization

abstract: Numerous challenges in science and engineering can be framed as optimization tasks, including the maximization of reaction yields, the optimization of molecular and materials properties, and the fine-tuning of automated hardware protocols. Design of experiment and optimization algorithms are often adopted to solve these tasks efficiently. Increasingly, these experiment planning strategies are coupled with automated hardware to enable autonomous experimental platforms. The vast majority of the strategies used, however, do not consider robustness against the variability of experiment and process conditions. In fact, it is generally assumed that these parameters are exact and reproducible. Yet some experiments may have considerable noise associated with some of their conditions, and process parameters optimized under precise control may be applied in the future under variable operating conditions. In either scenario, the optimal solutions found might not be robust against input variability, affecting the reproducibility of results and returning suboptimal performance in practice. Here, we introduce Golem, an algorithm that is agnostic to the choice of experiment planning strategy and that enables robust experiment and process optimization. Golem identifies optimal solutions that are robust to input uncertainty, thus ensuring the reproducible performance of optimized experimental protocols and processes. It can be used to analyze the robustness of past experiments, or to guide experiment planning algorithms toward robust solutions on the fly. We assess the performance and domain of applicability of Golem through extensive benchmark studies and demonstrate its practical relevance by optimizing an analytical chemistry protocol under the presence of significant noise in its experimental conditions.