title: Approximate decoherence free subspaces for distributed sensing

abstract: We consider the sensing of scalar valued fields with specific spatial dependence using a network of sensors, e.g. multiple atoms located at different positions within a trap. We show how to harness the spatial correlations to sense only a specific signal, and be insensitive to others at different positions or with unequal spatial dependence by constructing a decoherence-free subspace for noise sources at fixed, known positions. This can be extended to noise sources lying on certain surfaces, where we encounter a connection to mirror charges and equipotential surfaces in classical electrostatics. For general situations, we introduce the notion of an approximate decoherence-free subspace, where noise for all sources within some volume is significantly suppressed, at the cost of reducing the signal strength in a controlled way. We show that one can use this approach to maintain Heisenberg-scaling over long times and for a large number of sensors, despite the presence of multiple noise sources in large volumes. We introduce an efficient formalism to construct internal states and sensor configurations, and apply it to several examples to demonstrate the usefulness and wide applicability of our approach.


title: Resonance refraction and neutrino oscillations

abstract: The refraction index and matter potential depend on neutrino energy and this dependence has a resonance character associated to the production of the mediator in the $s-$channel. For light mediators and light particles of medium (background) the resonance can be realized at energies accessible to laboratory experiments. We study properties of the energy dependence of the potential for different C-asymmetries of background. Interplay of the background potential and the vacuum term leads to (i) bump in the oscillation probability in the resonance region, (ii) dip related to the MSW resonance in the background, (iii) substantial deviation of the effective $\Delta m^2$ above the resonance from the low energy value, etc. We considered generation of mixing in the background. Interactions with background shifts the energy of usual MSW resonance and produces new MSW resonances. Searches of the background effects allow us to put bounds on new interactions of neutrinos and properties of the background. We show that explanation of the MiniBooNE excess, as the bump due to resonance refraction, is excluded.


title: Bayesian Mechanics for Stationary Processes

abstract: This paper develops a Bayesian mechanics for adaptive systems. Firstly, we model the interface between a system and its environment with a Markov blanket. This affords conditions under which states internal to the blanket encode information about external states. Second, we introduce dynamics and represent adaptive systems as Markov blankets at steady-state. This allows us to identify a wide class of systems whose internal states appear to infer external states, consistent with variational inference in Bayesian statistics and theoretical neuroscience. Finally, we partition the blanket into sensory and active states. It follows that active states can be seen as performing active inference and well-known forms of stochastic control (such as PID control), which are prominent formulations of adaptive behaviour in theoretical biology and engineering.

categories: math-ph math.MP math.OC nlin.AO q-bio.NC


title: Electrohydrodynamic instabilities in freely suspended viscous films under normal electric fields

abstract: Electrohydrodynamic instabilities of fluid-fluid interfaces can be exploited in various microfluidic applications in order to enhance mixing, replicate well-controlled patterns or generate drops of a particular size. In this work, we study the stability and dynamics of a system of three superimposed layers of two immiscible fluids subject to a normal electric field. Following the Taylor-Melcher leaky dielectric model, the bulk remains electroneutral while a net charge accumulates on the interfaces. The interfacial charge dynamics is captured by a conservation equation accounting for Ohmic conduction, advection by the flow and finite charge relaxation. Using this model, we perform a linear stability analysis and identify different modes of instability, and we characterize the behavior of the system as a function of the relevant dimensionless groups in each mode. Further, we perform numerical simulations using the boundary element method in order to study the effect of nonlinearities on long-time interfacial dynamics. We demonstrate how the coupling of flow and surface charge transport in different modes of instability can give rise to nonlinear phenomena such as tip streaming or pinching of the film into droplets.


title: Three-dimensional de Sitter horizon thermodynamics

abstract: We explore thermodynamic contributions to the three-dimensional de Sitter horizon originating from metric and Chern-Simons gauge field fluctuations. In Euclidean signature these are computed by the partition function of gravity coupled to matter semi-classically expanded about the round three-sphere saddle. We investigate a corresponding Lorentzian picture - drawing inspiration from the topological entanglement entropy literature - in the form of an edge-mode theory residing at the de Sitter horizon. We extend the discussion to three-dimensional gravity with positive cosmological constant, viewed (semi-classically) as a complexified Chern-Simons theory. The putative gravitational edge-mode theory is a complexified version of the chiral Wess-Zumino-Witten model associated to the edge-modes of ordinary Chern-Simons theory. We introduce and solve a family of complexified Abelian Chern-Simons theories as a way to elucidate some of the more salient features of the gravitational edge-mode theories. We comment on the relation to the AdS$_4$/CFT$_3$ correspondence.


title: Persian Rhetorical Structure Theory

abstract: Over the past years, interest in discourse analysis and discourse parsing has steadily grown, and many discourse-annotated corpora and, as a result, discourse parsers have been built. In this paper, we present a discourse-annotated corpus for the Persian language built in the framework of Rhetorical Structure Theory as well as a discourse parser built upon the DPLP parser, an open-source discourse parser. Our corpus consists of 150 journalistic texts, each text having an average of around 400 words. Corpus texts were annotated using 18 discourse relations and based on the annotation guideline of the English RST Discourse Treebank corpus. Our text-level discourse parser is trained using gold segmentation and is built upon the DPLP discourse parser, which uses a large-margin transition-based approach to solve the problem of discourse parsing. The performance of our discourse parser in span (S), nuclearity (N) and relation (R) detection is around 78%, 64%, 44% respectively, in terms of F1 measure.


title: Ladder Polynomial Neural Networks

abstract: Polynomial functions have plenty of useful analytical properties, but they are rarely used as learning models because their function class is considered to be restricted. This work shows that when trained properly polynomial functions can be strong learning models. Particularly this work constructs polynomial feedforward neural networks using the product activation, a new activation function constructed from multiplications. The new neural network is a polynomial function and provides accurate control of its polynomial order. It can be trained by standard training techniques such as batch normalization and dropout. This new feedforward network covers several previous polynomial models as special cases. Compared with common feedforward neural networks, the polynomial feedforward network has closed-form calculations of a few interesting quantities, which are very useful in Bayesian learning. In a series of regression and classification tasks in the empirical study, the proposed model outperforms previous polynomial models.


title: Experimental Quantum Embedding for Machine Learning

abstract: The classification of big data usually requires a mapping onto new data clusters which can then be processed by machine learning algorithms by means of more efficient and feasible linear separators. Recently, Lloyd et al. have advanced the proposal to embed classical data into quantum ones: these live in the more complex Hilbert space where they can get split into linearly separable clusters. Here, we implement these ideas by engineering two different experimental platforms, based on quantum optics and ultra-cold atoms respectively, where we adapt and numerically optimize the quantum embedding protocol by deep learning methods, and test it for some trial classical data. We perform also a similar analysis on the Rigetti superconducting quantum computer. Therefore, we find that the quantum embedding approach successfully works also at the experimental level and, in particular, we show how different platforms could work in a complementary fashion to achieve this task. These studies might pave the way for future investigations on quantum machine learning techniques especially based on hybrid quantum technologies.

categories: quant-ph cond-mat.dis-nn cond-mat.quant-gas physics.optics


title: Spatio-Temporal Economic Properties of Multi-Product Supply Chains

abstract: In this work, we analyze the spatio-temporal economic properties of multi-product supply chains. Specifically, we interpret the supply chain as a coordinated market in which stakeholders (suppliers, consumers, and providers of transportation, storage, and transformation services) bid into a market that is cleared by an independent entity to obtain allocations and prices. The proposed model provides a general graph representation of spatio-temporal product transport that helps capture geographical transport, time delays, and storage (temporal transport) in a unified and compact manner. This representation allows us to establish fundamental economic properties for the supply chain (revenue adequacy, cost recovery, and competitiveness) and to establish bounds for space-time prices. To illustrate the concepts, we consider a case study in which organic waste is used for producing biogas and electricity. Our market model shows that incentives for waste storage emerge from electricity demand dynamics and illustrates how space-time price dynamics for waste and derived products emerge from geographical transport and storage.


title: Massive Scalar Perturbation of Extremal Rotating Braneworld Black hole: Superradiant Stability Analysis

abstract: We analyse the superradiant stability of braneworld extremal Kerr and Kerr-Newman black holes under massive scalar perturbation. These black hole solutions differ from their four dimensional counterpart by the presence of a tidal charge, which unlike electric charge, can take both positive and negative values and carries the signature of extra dimensions. We consider the perturbation of the brane geometry by a massive scalar field. From the radial equation of motion of the scalar field, we have found the effective potential felt by the field. For superradiant stability, there should not be any trapping well of the effective potential outside of the horizon. Using this condition we have specified the parameter space of superradiant stability which depends on parameters of both black hole and perturbation. In the case of Kerr-like extremal braneworld black holes, we have obtained a bound on tidal charge for superradiant stability under a massive scalar perturbation. Similarly, for slowly rotating Kerr-Newman like extremal braneworld black holes, we have obtained a lower bound of the product of mass to charge ratio of the black hole and scalar field, for which the system is superradiantly stable.


title: Room-temperature continuous-wave Dirac-vortex topological lasers on silicon

abstract: Robust laser sources are a fundamental building block for contemporary information technologies. Originating from condensed-matter physics, the concept of topology has recently entered the realm of optics, offering fundamentally new design principles for lasers with enhanced robustness. In analogy to the well-known Majorana fermions in topological superconductors, Dirac-vortex states have recently been investigated in passive photonic systems and are now considered as a promising candidate for single-mode large-area lasers. Here, we experimentally realize the first Dirac-vortex topological lasers in InAs/InGaAs quantum-dot materials monolithically grown on a silicon substrate. We observe room-temperature continuous-wave single-mode linearly polarized vertical laser emission at a telecom wavelength. Most importantly, we confirm that the wavelength of the Dirac-vortex laser is topologically robust against variations in the cavity size, and its free spectral range defies the universal inverse scaling law with the cavity size. These lasers will play an important role in CMOS-compatible photonic and optoelectronic systems on a chip.


title: Universal quantum circuits for quantum chemistry

abstract: Universal gate sets for quantum computing have been known for decades, yet no universal gate set has been proposed for particle-conserving unitaries, which are the operations of interest in quantum chemistry. In this work, we show that controlled single-excitation gates in the form of Givens rotations are universal for particle-conserving unitaries. Single-excitation gates describe an arbitrary $U(2)$ rotation on the two-qubit subspace spanned by the states $|01\rangle, |10\rangle$, while leaving other states unchanged -- a transformation that is analogous to a single-qubit rotation on a dual-rail qubit. The proof is constructive, so our result also provides an explicit method for compiling arbitrary particle-conserving unitaries. Additionally, we describe a method for using controlled single-excitation gates to prepare an arbitrary state of a fixed number of particles. We derive analytical gradient formulas for Givens rotations as well as decompositions into single-qubit and CNOT gates. Our results offer a unifying framework for quantum computational chemistry where every algorithm is a unique recipe built from the same universal ingredients: Givens rotations.


title: Variational quantum algorithm for molecular geometry optimization

abstract: Classical algorithms for predicting the equilibrium geometry of strongly correlated molecules require expensive wave function methods that become impractical already for few-atom systems. In this work, we introduce a variational quantum algorithm for finding the most stable structure of a molecule by explicitly considering the parametric dependence of the electronic Hamiltonian on the nuclear coordinates. The equilibrium geometry of the molecule is obtained by minimizing a more general cost function that depends on both the quantum circuit and the Hamiltonian parameters, which are simultaneously optimized at each step. The algorithm is applied to find the equilibrium geometries of the $\mathrm{H}_2$, $\mathrm{H}_3^+$, $\mathrm{BeH}_2$ and $\mathrm{H}_2\mathrm{O}$ molecules. The quantum circuits used to prepare the electronic ground state for each molecule were designed using an adaptive algorithm where excitation gates in the form of Givens rotations are selected according to the norm of their gradient. All quantum simulations are performed using the PennyLane library for quantum differentiable programming. The optimized geometrical parameters for the simulated molecules show an excellent agreement with their counterparts computed using classical quantum chemistry methods.


title: Many-body localization in the interpolating Aubry-Andr\'e-Fibonacci model

abstract: We investigate the localization properties of a spin chain with an antiferromagnetic nearest-neighbour coupling, subject to an external quasiperiodic on-site magnetic field. The quasiperiodic modulation interpolates between two paradigmatic models, namely the Aubry-Andr\'e and the Fibonacci models. We find that stronger many-body interactions extend the ergodic phase in the former, whereas they shrink it in the latter. Furthermore, the many-body localization transition points at the two limits of the interpolation appear to be continuously connected along the deformation. As a result, the position of the many-body localization transition depends on the interaction strength for an intermediate degree of deformation of the quasiperiodic modulation. Moreover, in the region of parameter space where the single-particle spectrum contains both localized and extended states, many-body interactions induce an anomalous effect: weak interactions localize the system, whereas stronger interactions enhance ergodicity. We map the model's localization phase diagram using the decay of the quenched spin imbalance in relatively long chains. This is accomplished employing a time-dependent variational approach applied to a matrix product state decomposition of the many-body state. Our model serves as a rich playground for testing many-body localization under tunable potentials.


title: EARLIN: Early Out-of-Distribution Detection for Resource-efficient Collaborative Inference

abstract: Collaborative inference enables resource-constrained edge devices to make inferences by uploading inputs (e.g., images) to a server (i.e., cloud) where the heavy deep learning models run. While this setup works cost-effectively for successful inferences, it severely underperforms when the model faces input samples on which the model was not trained (known as Out-of-Distribution (OOD) samples). If the edge devices could, at least, detect that an input sample is an OOD, that could potentially save communication and computation resources by not uploading those inputs to the server for inference workload. In this paper, we propose a novel lightweight OOD detection approach that mines important features from the shallow layers of a pretrained CNN model and detects an input sample as ID (In-Distribution) or OOD based on a distance function defined on the reduced feature space. Our technique (a) works on pretrained models without any retraining of those models, and (b) does not expose itself to any OOD dataset (all detection parameters are obtained from the ID training dataset). To this end, we develop EARLIN (EARLy OOD detection for Collaborative INference) that takes a pretrained model and partitions the model at the OOD detection layer and deploys the considerably small OOD part on an edge device and the rest on the cloud. By experimenting using real datasets and a prototype implementation, we show that our technique achieves better results than other approaches in terms of overall accuracy and cost when tested against popular OOD datasets on top of popular deep learning models pretrained on benchmark datasets.


title: A logical framework with a graph meta-language

abstract: We conjecture that the relative unpopularity of logical frameworks among practitioners is partly due to their complex meta-languages, which often demand both programming skills and theoretical knowledge of the meta-language in question for them to be fruitfully used. We present ongoing work on a logical framework with a meta-language based on graphs. A simpler meta-language leads to a shallower embedding of the object language, but hopefully leads to easier implementation and usage. A graph-based logical framework also opens up interesting possibilities in time and space performance by using heavily optimized graph databases as backends and by proof compression algorithms geared towards graphs. Deductive systems can be specified using simple domain-specific languages built on top of the graph database's query language. There is support for interactive (through a web-based interface) and semiautomatic (following user-defined tactics specified by a domain-specific language) proof modes. We have so far implemented nine systems for propositional logic, including Fitch-style and backward-directed Natural Deduction systems for intuitionistic and classical logic (with the classical systems reusing the rules of the intuitionistic ones), and a Hilbert-style system for the K modal logic.


title: Extensions to Multifidelity Monte Carlo Methods for Simulations of Chaotic Systems

abstract: Multifidelity Monte Carlo methods often rely on a preprocessing phase consisting of standard Monte Carlo sampling to estimate correlation coefficients between models of different fidelity to determine the weights and number of samples for each level. For computationally intensive models, as are often encountered in simulations of chaotic systems, this up-front cost can be prohibitive. In this work, a correlation estimation procedure is developed for the case in which the highest and next highest fidelity models are generated via discretizing the same mathematical model using different resolution. The procedure uses discretization error estimates to estimate the required correlation coefficient without the need to sample the highest fidelity model, which can dramatically decrease the cost of the preprocessing phase. The method is extended to chaotic problems by using discretization error estimates that account for the statistical nature of common quantities of interest and the accompanying finite sampling errors that pollute estimates of such quantities of interest. The methodology is then demonstrated on a model problem based on the Kuramoto-Sivashinsky equation.


title: Focusing Atom Laser Beams

abstract: We theoretically study the focusing of a quasi-continuous atom laser beam of rubidium-85 ($^{85}$Rb). A two-sate model analysis based on the Gross-Pitaevskii equation is used which comprises the effects of two-body atom-atom interactions and three-body recombination losses. Utilizing optical focusing potentials such as harmonic potentials, the essential factors such as the width, peak density and atom loss rate of the focused atom laser beam profile are investigated. Our analysis predicts that using an atom laser offers a dramatic improvement in resolution of up to $8$ nm.


title: On the Nature of AGN and Star Formation Enhancement in the $z = 3.1$ SSA22 Protocluster: The HST WFC3 IR View

abstract: We examine possible environmental sources of the enhanced star formation and active galactic nucleus (AGN) activity in the $z = 3.09$ SSA22 protocluster using Hubble WFC3 F160W ($\sim1.6\ \rm \mu m$) observations of the SSA22 field, including new observations centered on eight X-ray selected protocluster AGN. To investigate the role of mergers in the observed AGN and star formation enhancement, we apply both quantitative (S\'ersic-fit and Gini-$M_{20}$) and visual morphological classifications to F160W images of protocluster Lyman break galaxies (LBGs) in the fields of the X-ray AGN and $z \sim 3$ field LBGs in SSA22 and GOODS-N. We find no statistically significant differences between the morphologies and merger fractions of protocluster and field LBGs, though we are limited by small number statistics in the protocluster. We also fit the UV-to-near-IR spectral energy distributions (SED) of F160W-detected protocluster and field LBGs to characterize their stellar masses and star formation histories (SFH). We find that the mean protocluster LBG is by a factor of $\sim2$ times more massive and more attenuated than the mean $z \sim 3$ field LBG. We take our results to suggest that ongoing mergers are not more common among protocluster LBGs than field LBGs, though protocluster LBGs appear to be more massive. We speculate that the larger mass of the protocluster LBGs contributes to the enhancement of SMBH mass and accretion rate in the protocluster, which in turn drives the observed protocluster AGN enhancement.


title: Probing protoplanetary disk evolution in the Chamaeleon II region

abstract: Context. Characterizing the evolution of protoplanetary disks is necessary to improve our understanding of planet formation. Constraints on both dust and gas are needed to determine the dominant disk dissipation mechanisms. Aims. We aim to compare the disk dust masses in the Chamaeleon II (Cha II) star-forming region with other regions with ages between 1 and 10Myr. Methods. We use ALMA band 6 observations (1.3 mm) to survey 29 protoplanetary disks in Cha II. Dust mass estimates are derived from the continuum data. Results. Out of our initial sample of 29 disks, we detect 22 sources in the continuum, 10 in 12CO, 3 in 13CO, and none in C18O (J=2-1). Additionally, we detect two companion candidates in the continuum and 12CO emission. Most disk dust masses are lower than 10Mearth, assuming thermal emission from optically thin dust. We compare consistent estimations of the distributions of the disk dust mass and the disk-to-stellar mass ratios in Cha II with six other low mass and isolated star-forming regions in the age range of 1-10Myr: Upper Sco, CrA, IC 348, Cha I, Lupus, and Taurus. When comparing the dust-to-stellar mass ratio, we find that the masses of disks in Cha II are statistically different from those in Upper Sco and Taurus, and we confirm that disks in Upper Sco, the oldest region of the sample, are statistically less massive than in all other regions. Performing a second statistical test of the dust mass distributions from similar mass bins, we find no statistical differences between these regions and Cha II. Conclusions. We interpret these trends, most simply, as a sign of decline in the disk dust masses with time or dust evolution. Different global initial conditions in star-forming regions may also play a role, but their impact on the properties of a disk population is difficult to isolate in star-forming regions lacking nearby massive stars.


title: A Generalization of Thorpe's Inequality

abstract: We present a generalization of the topological inequality of Thorpe between the Euler characteristic and $k^{th}$-Pontryagin number of a $4k$-manifold. We also correct and complete some of the arguments from the work of Thorpe in which this inequality originally appeared.


title: A CNN Segmentation-Based Approach to Object Detection and Tracking in Ultrasound Scans with Application to the Vagus Nerve Detection

abstract: Ultrasound scanning is essential in several medical diagnostic and therapeutic applications. It is used to visualize and analyze anatomical features and structures that influence treatment plans. However, it is both labor intensive, and its effectiveness is operator dependent. Real-time accurate and robust automatic detection and tracking of anatomical structures while scanning would significantly impact diagnostic and therapeutic procedures to be consistent and efficient. In this paper, we propose a deep learning framework to automatically detect and track a specific anatomical target structure in ultrasound scans. Our framework is designed to be accurate and robust across subjects and imaging devices, to operate in real-time, and to not require a large training set. It maintains a localization precision and recall higher than 90% when trained on training sets that are as small as 20% in size of the original training set. The framework backbone is a weakly trained segmentation neural network based on U-Net. We tested the framework on two different ultrasound datasets with the aim to detect and track the Vagus nerve, where it outperformed current state-of-the-art real-time object detection networks.


title: Long-Term Environmental Stability of Nitrogen-Healed Black Phosphorus

abstract: The unique optoelectronic properties of black phosphorus (BP) have triggered great interest in its applications in areas not fulfilled by other layered materials (LMs). However, its poor stability (fast degradation, i.e. <<1 h for monolayers) under ambient conditions restricts its practical application. We demonstrate here, by an experimental-theoretical approach, that the incorporation of nitrogen molecules (N2) into the BP structure results in a relevant improvement of its stability in air, up to 8 days without optical degradation signs. Our strategy involves the generation of defects (phosphorus vacancies) by electron-beam irradiation, followed by their healing with N2 molecules. As an additional route, N2 plasma treatment is presented as an alternative for large area application. Our first principles calculations elucidate the mechanisms involved in the nitrogen incorporation as well as on the stabilization of the modified BP, which corroborates with our experimental observations. This stabilization approach can be applied in the processing of BP, allowing for its use in environmentally stable van der Waals heterostructures with other LMs as well as in optoelectronic and wearable devices.


title: Approximate Maximum Halfspace Discrepancy

abstract: Consider the geometric range space $(X, \mathcal{H}_d)$ where $X \subset \mathbb{R}^d$ and $\mathcal{H}_d$ is the set of ranges defined by $d$-dimensional halfspaces. In this setting we consider that $X$ is the disjoint union of a red and blue set. For each halfspace $h \in \mathcal{H}_d$ define a function $\Phi(h)$ that measures the "difference" between the fraction of red and fraction of blue points which fall in the range $h$. In this context the maximum discrepancy problem is to find the $h^* = \arg \max_{h \in (X, \mathcal{H}_d)} \Phi(h)$. We aim to instead find an $\hat{h}$ such that $\Phi(h^*) - \Phi(\hat{h}) \le \varepsilon$. This is the central problem in linear classification for machine learning, in spatial scan statistics for spatial anomaly detection, and shows up in many other areas. We provide a solution for this problem in $O(|X| + (1/\varepsilon^d) \log^4 (1/\varepsilon))$ time, which improves polynomially over the previous best solutions. For $d=2$ we show that this is nearly tight through conditional lower bounds. For different classes of $\Phi$ we can either provide a $\Omega(|X|^{3/2 - o(1)})$ time lower bound for the exact solution with a reduction to APSP, or an $\Omega(|X| + 1/\varepsilon^{2-o(1)})$ lower bound for the approximate solution with a reduction to 3SUM. A key technical result is a $\varepsilon$-approximate halfspace range counting data structure of size $O(1/\varepsilon^d)$ with $O(\log (1/\varepsilon))$ query time, which we can build in $O(|X| + (1/\varepsilon^d) \log^4 (1/\varepsilon))$ time.


title: Decomposition of transition systems into sets of synchronizing state machines

abstract: Transition systems (TS) and Petri nets (PN) are important models of computation ubiquitous in formal methods for modeling systems. An important problem is how to extract from a given TS a PN whose reachability graph is equivalent (with a suitable notion of equivalence) to the original TS. This paper addresses the decomposition of transition systems into synchronizing state machines (SMs), which are a class of Petri nets where each transition has one incoming and one outgoing arc and all markings have exactly one token. This is an important case of the general problem of extracting a PN from a TS. The decomposition is based on the theory of regions, and it is shown that a property of regions called excitation-closure is a sufficient condition to guarantee the equivalence between the original TS and a decomposition into SMs. An efficient algorithm is provided which solves the problem by reducing its critical steps to the maximal independent set problem (to compute a minimal set of irredundant SMs) or to satisfiability (to merge the SMs). We report experimental results that show a good trade-off between quality of results vs. computation time.


title: Hierarchical Online Convex Optimization

abstract: We consider online convex optimization (OCO) over a heterogeneous network with communication delay, where multiple workers together with a master execute a sequence of decisions to minimize the accumulation of time-varying global costs. The local data may not be independent or identically distributed, and the global cost functions may not be locally separable. Due to communication delay, neither the master nor the workers have in-time information about the current global cost function. We propose a new algorithm, termed Hierarchical OCO (HiOCO), which takes full advantage of the network heterogeneity in information timeliness and computation capacity to enable multi-step gradient descent at both the workers and the master. We analyze the impacts of the unique hierarchical architecture, multi-slot delay, and gradient estimation error to derive upper bounds on the dynamic regret of HiOCO, which measures the gap of costs between HiOCO and an offline globally optimal performance benchmark.


title: Interplay between polarization, strain and defect-pairs in Fe-doped SrMnO$_{3-\delta}$

abstract: Defect chemistry, strain, and structural, magnetic and electronic degrees of freedom constitute a rich space for the design of functional properties in transition metal oxides. Here, we show that it is possible to engineer polarity and ferroelectricity in non-polar perovskite oxides via polar defect pairs formed by anion vacancies coupled to substitutional cations. We use a self-consistent site-dependent DFT+$U$ approach that accounts for local structural and chemical changes upon defect creation and which is crucial to reconcile predictions with the available experimental data. Our results for Fe-doped oxygen-deficient SrMnO$_3$ show that substitutional Fe and oxygen vacancies can promote polarity due to an off-center displacement of the defect charge resulting in a net electric dipole moment, which polarizes the lattice in the defect neighborhood. The formation of these defects and the resulting polarization can be tuned by epitaxial strain, resulting in enhanced polarization also for strain values lower than the ones necessary to induce a polar phase transition in undoped SrMnO$_3$. For high enough defect concentrations, these defect dipoles couple in a parallel fashion, thus enabling defect- and strain-based engineering of ferroelectricity in SrMnO$_3$.


title: Unbinned Angular Analysis of $B\to D^*\ell \nu_\ell$ and the Right-handed Current

abstract: In this article, we perform a sensitivity study of an unbinned angular analysis of the $B\to D^*\ell \nu_\ell$ decay, including the contributions from the right-handed current. We show that the angular observable can constrain very strongly the right-handed current without the intervention of the yet unsolved $V_{cb}$ puzzle.


title: Nonparametric inference on counterfactuals in first-price auctions

abstract: In a classical model of the first-price sealed-bid auction with independent private values, we develop nonparametric estimation and inference procedures for a class of policy-relevant metrics, such as total expected surplus and expected revenue under counterfactual reserve prices. Motivated by the linearity of these metrics in the quantile function of bidders' values, we propose a bid spacings-based estimator of the latter and derive its Bahadur-Kiefer expansion. This makes it possible to construct exact uniform confidence bands and assess the optimality of a given auction rule. Using the data on U.S. Forest Service timber auctions, we test whether setting zero reserve prices in these auctions was revenue maximizing.


title: Majorana Dark Matter and Neutrino mass in a singlet-doublet extension of the Standard Model

abstract: A minimal extension of the Standard Model (SM) by a vector-like fermion doublet and three right handed (RH) singlet neutrinos is proposed in order to explain dark matter and tiny neutrino mass simultaneously. The DM arises as a mixture of the neutral component of the fermion doublet and one of the RH neutrinos, both assumed to be odd under an imposed $\mathcal{Z}_2$ symmetry. Being Majorana in nature, the DM escapes from $Z$-mediated direct search constraints to mark a significant difference from singlet-doublet Dirac DM. The other two $\mathcal{Z}_2$ even heavy RH neutrinos give rise masses and mixing of light neutrinos via Type-I Seesaw mechanism. Relic density and direct search allowed parameter space for the model is investigated through detailed numerical scan.


title: Semantic Parsing Natural Language into Relational Algebra

abstract: Natural interface to database (NLIDB) has been researched a lot during the past decades. In the core of NLIDB, is a semantic parser used to convert natural language into SQL. Solutions from traditional NLP methodology focuses on grammar rule pattern learning and pairing via intermediate logic forms. Although those methods give an acceptable performance on certain specific database and parsing tasks, they are hard to generalize and scale. On the other hand, recent progress in neural deep learning seems to provide a promising direction towards building a general NLIDB system. Unlike the traditional approach, those neural methodologies treat the parsing problem as a sequence-to-sequence learning problem. In this paper, we experimented on several sequence-to-sequence learning models and evaluate their performance on general database parsing task.


title: rFaaS: Enabling High Performance Serverless with RDMA and Leases

abstract: High performance is needed in many computing systems, from batch-managed supercomputers to general-purpose cloud platforms. However, scientific clusters lack elastic parallelism, while clouds cannot offer competitive costs for high-performance applications. In this work, we investigate how modern cloud programming paradigms can bring the elasticity needed to allocate idle resources, decreasing computation costs and improving overall data center efficiency. Function-as-a-Service (FaaS) brings the pay-as-you-go execution of stateless functions, but its performance characteristics cannot match coarse-grained cloud and cluster allocations. To make serverless computing viable for high-performance and latency-sensitive applications, we present rFaaS, an RDMA-accelerated FaaS platform. We identify critical limitations of serverless - centralized scheduling and inefficient network transport - and improve the FaaS architecture with allocation leases and microsecond invocations. We show that our remote functions add only negligible overhead on top of the fastest available networks, and we decrease the execution latency by orders of magnitude compared to contemporary FaaS systems. Furthermore, we demonstrate the performance of rFaaS by evaluating real-world FaaS benchmarks and parallel applications. Overall, our results show that new allocation policies and remote memory access help FaaS applications achieve high performance and bring serverless computing to HPC.


title: Threshold-Based Quantum Optimization

abstract: We propose and study Th-QAOA (pronounced Threshold QAOA), a variation of the Quantum Alternating Operator Ansatz (QAOA) that replaces the standard phase separator operator, which encodes the objective function, with a threshold function that returns a value $1$ for solutions with an objective value above the threshold and a $0$ otherwise. We vary the threshold value to arrive at a quantum optimization algorithm. We focus on a combination with the Grover Mixer operator; the resulting GM-Th-QAOA can be viewed as a generalization of Grover's quantum search algorithm and its minimum/maximum finding cousin to approximate optimization. Our main findings include: (i) we provide intuitive arguments and show empirically that the optimum parameter values of GM-Th-QAOA (angles and threshold value) can be found with $O(\log(p) \times \log M)$ iterations of the classical outer loop, where $p$ is the number of QAOA rounds and $M$ is an upper bound on the solution value (often the number of vertices or edges in an input graph), thus eliminating the notorious outer-loop parameter finding issue of other QAOA algorithms; (ii) GM-Th-QAOA can be simulated classically with little effort up to 100 qubits through a set of tricks that cut down memory requirements; (iii) somewhat surprisingly, GM-Th-QAOA outperforms non-thresholded GM-QAOA in terms of approximation ratios achieved. This third result holds across a range of optimization problems (MaxCut, Max k-VertexCover, Max k-DensestSubgraph, MaxBisection) and various experimental design parameters, such as different input edge densities and constraint sizes.


title: Auto-Pipeline: Synthesizing Complex Data Pipelines By-Target Using Reinforcement Learning and Search

abstract: Recent work has made significant progress in helping users to automate single data preparation steps, such as string-transformations and table-manipulation operators (e.g., Join, GroupBy, Pivot, etc.). We in this work propose to automate multiple such steps end-to-end, by synthesizing complex data pipelines with both string transformations and table-manipulation operators. We propose a novel "by-target" paradigm that allows users to easily specify the desired pipeline, which is a significant departure from the traditional by-example paradigm. Using by-target, users would provide input tables (e.g., csv or json files), and point us to a "target table" (e.g., an existing database table or BI dashboard) to demonstrate how the output from the desired pipeline would schematically "look like". While the problem is seemingly underspecified, our unique insight is that implicit table constraints such as FDs and keys can be exploited to significantly constrain the space to make the problem tractable. We develop an Auto-Pipeline system that learns to synthesize pipelines using reinforcement learning and search. Experiments on large numbers of real pipelines crawled from GitHub suggest that Auto-Pipeline can successfully synthesize 60-70% of these complex pipelines with up to 10 steps.


title: Accessing Intermediate-Mass Black Holes in 728 Globular Star Clusters in NGC\,4472

abstract: Intermediate-mass black holes (IMBHs) by definition have masses of $M_{\rm IMBH} \sim 10^{2-5}~M_\odot$, a range with few observational constraints. Finding IMBHs in globular star clusters (GCs) would validate a formation channel for massive black-hole seeds in the early universe. Here, we simulate a 60-hour observation with the next-generation Very Large Array (ngVLA) of 728 GC candidates in the Virgo Cluster galaxy NGC\,4472. Interpreting the radio detection thresholds as signatures of accretion onto IMBHs, we benchmark IMBH mass thresholds in three scenarios and find the following: (1) Radio analogs of ESO\,243-49 HLX-1, a strong IMBH candidate with $M_{\rm IMBH}^{\rm HLX} \sim 10^{4-5}~M_\odot$ in a star cluster, are easy to access in all 728 GC candidates. (2) For the 30 GC candidates with extant X-ray detections, the empirical fundamental-plane relation involving black hole mass plus X-ray and radio luminosities suggests access to $M_{\rm IMBH}^{\rm FP} \sim 10^{1.7-3.6}~M_\odot$, with an uncertainty of 0.44 dex. (3) A fiducial Bondi accretion model was applied to all 728 GC candidates and to radio stacks of GC candidates. This model suggests access to IMBH masses, with a statistical uncertainty of 0.39 dex, of $M_{\rm IMBH}^{\rm B} \sim 10^{4.9-5.1}~M_\odot$ for individual GC candidates and $M_{\rm IMBH}^{\rm B,stack} \sim 10^{4.5}~M_\odot$ for radio stacks of about 100-200 GC candidates. The fiducial Bondi model offers initial guidance, but is subject to additional systematic uncertainties and should be superseded by hydrodynamical simulations of gas flows in GCs.


title: Steerable 3D Spherical Neurons

abstract: Emerging from low-level vision theory, steerable filters found their counterpart in prior work on steerable convolutional neural networks equivariant to rigid transformations. In our work, we propose a steerable feed-forward learning-based approach that consists of neurons with spherical decision surfaces and operates on point clouds. Such spherical neurons are obtained by conformal embedding of Euclidean space and have recently been revisited in the context of learning representations of point sets. Focusing on 3D geometry, we exploit the isometry property of spherical neurons and derive a 3D steerability constraint. After training spherical neurons to classify point clouds in a canonical orientation, we use a tetrahedron basis to quadruplicate the neurons and construct rotation-equivariant spherical filter banks. We then apply the derived constraint to interpolate the filter bank outputs and, thus, obtain a rotation-invariant network. Finally, we use a synthetic point set and real-world 3D skeleton data to verify our theoretical findings. The code is available at https://github.com/pavlo-melnyk/steerable-3d-neurons.


title: Nonuniform Defocus Removal for Image Classification

abstract: We propose and study the single-frame anisoplanatic deconvolution problem associated with image classification using machine learning algorithms, named the nonuniform defocus removal (NDR) problem. Mathematical analysis of the NDR problem is done and the so-called defocus removal (DR) algorithm for solving it is proposed. Global convergence of the DR algorithm is established without imposing any unverifiable assumption. Numerical results on simulation data show significant features of DR including solvability, noise robustness, convergence, model insensitivity and computational efficiency. Physical relevance of the NDR problem and practicability of the DR algorithm are tested on experimental data. Back to the application that originally motivated the investigation of the NDR problem, we show that the DR algorithm can improve the accuracy of classifying distorted images using convolutional neural networks. The key difference of this paper compared to most existing works on single-frame anisoplanatic deconvolution is that the new method does not require the data image to be decomposable into isoplanatic subregions. Therefore, solution approaches partitioning the image into isoplanatic zones are not applicable to the NDR problem and those handling the entire image such as the DR algorithm need to be developed and analyzed.


title: A Photonic-Circuits-Inspired Compact Network: Toward Real-Time Wireless Signal Classification at the Edge

abstract: Machine learning (ML) methods are ubiquitous in wireless communication systems and have proven powerful for applications including radio-frequency (RF) fingerprinting, automatic modulation classification, and cognitive radio. However, the large size of ML models can make them difficult to implement on edge devices for latency-sensitive downstream tasks. In wireless communication systems, ML data processing at a sub-millisecond scale will enable real-time network monitoring to improve security and prevent infiltration. In addition, compact and integratable hardware platforms which can implement ML models at the chip scale will find much broader application to wireless communication networks. Toward real-time wireless signal classification at the edge, we propose a novel compact deep network that consists of a photonic-hardware-inspired recurrent neural network model in combination with a simplified convolutional classifier, and we demonstrate its application to the identification of RF emitters by their random transmissions. With the proposed model, we achieve 96.32% classification accuracy over a set of 30 identical ZigBee devices when using 50 times fewer training parameters than an existing state-of-the-art CNN classifier. Thanks to the large reduction in network size, we demonstrate real-time RF fingerprinting with 0.219 ms latency using a small-scale FPGA board, the PYNQ-Z1.


title: Certain residual properties of HNN-extensions with central associated subgroups

abstract: Suppose that $G$ is a group, $H$ and $K$ are proper isomorphic central subgroups of $G$, and $\mathfrak{G}$ is an HNN-extension of $G$ with the associated subgroups $H$ and $K$. We prove necessary and sufficient conditions for $\mathfrak{G}$ to be residually a $\mathcal{C}$-group, where $\mathcal{C}$ is a class of groups closed under taking subgroups, extensions, homomorphic images, and Cartesian products of the form $\prod_{y \in Y}X_{y}$, where $X, Y \in \mathcal{C}$ and $X_{y}$ is an isomorphic copy of $X$ for each $y \in Y$.


title: POLAR: A Polynomial Arithmetic Framework for Verifying Neural-Network Controlled Systems

abstract: We present POLAR, a polynomial arithmetic-based framework for efficient bounded-time reachability analysis of neural-network controlled systems (NNCSs). Existing approaches that leverage the standard Taylor Model (TM) arithmetic for approximating the neural-network controller cannot deal with non-differentiable activation functions and suffer from rapid explosion of the remainder when propagating the TMs. POLAR overcomes these shortcomings by integrating TM arithmetic with \textbf{Bernstein B{\'e}zier Form} and \textbf{symbolic remainder}. The former enables TM propagation across non-differentiable activation functions and local refinement of TMs, and the latter reduces error accumulation in the TM remainder for linear mappings in the network. Experimental results show that POLAR significantly outperforms the current state-of-the-art tools in terms of both efficiency and tightness of the reachable set overapproximation. The source code can be found in https://github.com/ChaoHuang2018/POLAR_Tool


title: Raman Linewidth Contributions from Four-Phonon and Electron-Phonon Interactions in Graphene

abstract: The Raman peak position and linewidth provide insight into phonon anharmonicity and electron-phonon interactions (EPI) in materials. For monolayer graphene, prior first-principles calculations have yielded decreasing linewidth with increasing temperature, which is opposite to measurement results. Here, we explicitly consider four-phonon anharmonicity, phonon renormalization, and electron-phonon coupling, and find all to be important to successfully explain both the $G$ peak frequency shift and linewidths in our suspended graphene sample at a wide temperature range. Four-phonon scattering contributes a prominent linewidth that increases with temperature, while temperature dependence from EPI is found to be reversed above a doping threshold ($\hbar\omega_G/2$, with $\omega_G$ being the frequency of the $G$ phonon).


title: A multi-stage machine learning model on diagnosis of esophageal manometry

abstract: High-resolution manometry (HRM) is the primary procedure used to diagnose esophageal motility disorders. Its interpretation and classification includes an initial evaluation of swallow-level outcomes and then derivation of a study-level diagnosis based on Chicago Classification (CC), using a tree-like algorithm. This diagnostic approach on motility disordered using HRM was mirrored using a multi-stage modeling framework developed using a combination of various machine learning approaches. Specifically, the framework includes deep-learning models at the swallow-level stage and feature-based machine learning models at the study-level stage. In the swallow-level stage, three models based on convolutional neural networks (CNNs) were developed to predict swallow type, swallow pressurization, and integrated relaxation pressure (IRP). At the study-level stage, model selection from families of the expert-knowledge-based rule models, xgboost models and artificial neural network(ANN) models were conducted, with the latter two model designed and augmented with motivation from the export knowledge. A simple model-agnostic strategy of model balancing motivated by Bayesian principles was utilized, which gave rise to model averaging weighted by precision scores. The averaged (blended) models and individual models were compared and evaluated, of which the best performance on test dataset is 0.81 in top-1 prediction, 0.92 in top-2 predictions. This is the first artificial-intelligence-style model to automatically predict CC diagnosis of HRM study from raw multi-swallow data. Moreover, the proposed modeling framework could be easily extended to multi-modal tasks, such as diagnosis of esophageal patients based on clinical data from both HRM and functional luminal imaging probe panometry (FLIP).


title: Scene Uncertainty and the Wellington Posterior of Deterministic Image Classifiers

abstract: We propose a method to estimate the uncertainty of the outcome of an image classifier on a given input datum. Deep neural networks commonly used for image classification are deterministic maps from an input image to an output class. As such, their outcome on a given datum involves no uncertainty, so we must specify what variability we are referring to when defining, measuring and interpreting uncertainty, and attributing "confidence" to the outcome. To this end, we introduce the Wellington Posterior, which is the distribution of outcomes that would have been obtained in response to data that could have been generated by the same scene that produced the given image. Since there are infinitely many scenes that could have generated any given image, the Wellington Posterior involves inductive transfer from scenes other than the one portrayed. We explore the use of data augmentation, dropout, ensembling, single-view reconstruction, and model linearization to compute a Wellington Posterior. Additional methods include the use of conditional generative models such as generative adversarial networks, neural radiance fields, and conditional prior networks. We test these methods against the empirical posterior obtained by performing inference on multiple images of the same underlying scene. These developments are only a small step towards assessing the reliability of deep network classifiers in a manner that is compatible with safety-critical applications and human interpretation.


title: Transflower: probabilistic autoregressive dance generation with multimodal attention

abstract: Dance requires skillful composition of complex movements that follow rhythmic, tonal and timbral features of music. Formally, generating dance conditioned on a piece of music can be expressed as a problem of modelling a high-dimensional continuous motion signal, conditioned on an audio signal. In this work we make two contributions to tackle this problem. First, we present a novel probabilistic autoregressive architecture that models the distribution over future poses with a normalizing flow conditioned on previous poses as well as music context, using a multimodal transformer encoder. Second, we introduce the currently largest 3D dance-motion dataset, obtained with a variety of motion-capture technologies, and including both professional and casual dancers. Using this dataset, we compare our new model against two baselines, via objective metrics and a user study, and show that both the ability to model a probability distribution, as well as being able to attend over a large motion and music context are necessary to produce interesting, diverse, and realistic dance that matches the music.


title: Faraday Waves in strongly interacting superfluids

abstract: We report on the observation of Faraday waves in a cigar-shaped Fermi superfluid of $^6$Li parametrically excited by modulating the radial trap frequency. We characterize the phenomenon as a function of the interaction parameter by means of a Feshbach resonance. Starting from the BEC side of the resonance we observe a drop on the visibility of the Faraday pattern as we approach to unitarity, possibly due to the increased incompressibility of the system. We probe the superfluid excitation spectrum by extracting an effective 1D speed of sound for different values of the interaction parameter, in good agreement with numerical simulations. Finally, we perform a stability analysis in the parameter space using a simplified model and we show the emergence of the Faraday waves as unstable solutions to a Mathieu-like equation.


title: On classical inequalities for autocorrelations and autoconvolutions

abstract: In this paper we study an autocorrelation inequality proposed by Barnard and Steinerberger. The study of these problems is motivated by a classical problem in additive combinatorics. We establish the existence of extremizers to this inequality, for a general class of weights, including Gaussian functions (as studied by the second author and Ramos) and characteristic function (as originally studied by Barnard and Steinerberger). Moreover, via a discretization argument and numerical analysis, we find some almost optimal approximation for the best constant allowed in this inequality. We also discuss some other related problem about autoconvolutions.


title: Statistical Methods for the meta-analysis paper by Itzhaky et al

abstract: This document describes the statistical methods used in Itzhaky et al ("Systematic Review and Meta-analysis: Twenty-six Years of Randomized Clinical Trials of Psychosocial Interventions to Reduce Suicide Risk in Adolescents"). That paper is a meta-analysis of randomized controlled clinical trials testing methods for preventing suicidal behavior and/or ideation in youth. Particularly on the behavior side the meta-data are challenging to analyze. This paper has two parts. The first is an informal discussion of the statistical methods used. The second gives detailed mathematical derivations of some formulas and methods.


title: Geometry, Number Theory and the Butterfly Spectrum of Two-Dimensional Bloch Electrons

abstract: We take a deeper dive into the geometry and the number theory that underlay the butterfly graphs of the Harper and the generalized Harper models of Bloch electrons in a magnetic field. Root of the number theoretical characteristics of the fractal spectrum is traced to a close relationship between the Farey tree -- the hierarchical tree that generates all rationals and the Wannier diagram -- a graph that labels all the gaps of the butterfly graph. The resulting Farey-Wannier hierarchical lattice of trapezoids provides geometrical representation of the nested pattern of butterflies in the butterfly graph. Some features of the energy spectrum such as absence of some of the Wannier trajectories in the butterfly graph fall outside the number theoretical framework, can be stated as a simple rule of "minimal violation of mirror symmetry". In a generalized Harper model, Farey-Wannier representation prevails as the lattice regroups to form some hexagonal unit cells creating new {\it species} of butterflies


title: Knowledge-Grounded Self-Rationalization via Extractive and Natural Language Explanations

abstract: Models that generate extractive rationales (i.e., subsets of features) or natural language explanations (NLEs) for their predictions are important for explainable AI. While an extractive rationale provides a quick view of the features most responsible for a prediction, an NLE allows for a comprehensive description of the decision-making process behind a prediction. However, current models that generate the best extractive rationales or NLEs often fall behind the state-of-the-art (SOTA) in terms of task performance. In this work, we bridge this gap by introducing RExC, a self-rationalizing framework that grounds its predictions and two complementary types of explanations (NLEs and extractive rationales) in background knowledge. Our framework improves over previous methods by: (i) reaching SOTA task performance while also providing explanations, (ii) providing two types of explanations, while existing models usually provide only one type, and (iii) beating by a large margin the previous SOTA in terms of quality of both types of explanations. Furthermore, a perturbation analysis in RExC shows a high degree of association between explanations and predictions, a necessary property of faithful explanations.


title: Numerical analysis of the LDG method for large deformations of prestrained plates

abstract: A local discontinuous Galerkin (LDG) method for approximating large deformations of prestrained plates is introduced and tested on several insightful numerical examples in our previous computational work. This paper presents a numerical analysis of this LDG method, focusing on the free boundary case. The problem consists of minimizing a fourth order bending energy subject to a nonlinear and nonconvex metric constraint. The energy is discretized using LDG and a discrete gradient flow is used for computing discrete minimizers. We first show $\Gamma$-convergence of the discrete energy to the continuous one. Then we prove that the discrete gradient flow decreases the energy at each step and computes discrete minimizers with control of the metric constraint defect. We also present a numerical scheme for initialization of the gradient flow, and discuss the conditional stability of it.


title: Convergence Analysis and Numerical Studies for Linearly Elastic Peridynamics with Dirichlet-Type Boundary Conditions

abstract: The nonlocal models of peridynamics have successfully predicted fractures and deformations for a variety of materials. In contrast to local mechanics, peridynamic boundary conditions must be defined on a finite volume region outside the body. Therefore, theoretical and numerical challenges arise in order to properly formulate Dirichlet-type nonlocal boundary conditions, while connecting them to the local counterparts. While a careless imposition of local boundary conditions leads to a smaller effective material stiffness close to the boundary and an artificial softening of the material, several strategies were proposed to avoid this unphysical surface effect. In this work, we study convergence of solutions to nonlocal state-based linear elastic model to their local counterparts as the interaction horizon vanishes, under different formulations and smoothness assumptions for nonlocal Dirichlet-type boundary conditions. Our results provide explicit rates of convergence that are sensitive to the compatibility of the nonlocal boundary data and the extension of the solution for the local model. In particular, under appropriate assumptions, constant extensions yield $\frac{1}{2}$ order convergence rates and linear extensions yield $\frac{3}{2}$ order convergence rates. With smooth extensions, these rates are improved to quadratic convergence. We illustrate the theory for any dimension $d\geq 2$ and numerically verify the convergence rates with a number of two dimensional benchmarks, including linear patch tests, manufactured solutions, and domains with curvilinear surfaces. Numerical results show a first order convergence for constant extensions and second order convergence for linear extensions, which suggests a possible room of improvement in the future convergence analysis.


title: Optimal Checkpointing for Adjoint Multistage Time-Stepping Schemes

abstract: We consider checkpointing strategies that minimize the number of recomputations needed when performing discrete adjoint computations using multistage time-stepping schemes, which requires computing several substeps within one complete time step. In this case we propose two algorithms that can generate optimal checkpointing schedules under weak assumptions. The first is an extension of the seminal Revolve algorithm adapted to multistage schemes. The second algorithm, named CAMS, is developed based on dynamic programming, and it requires the least number of recomputations when compared with other algorithms. The CAMS algorithm is made publicly available in a library with bindings to C and Python. Numerical results illustrate that the proposed algorithms can deliver up to two times the speedup compared with that of classical Revolve. Moreover, we discuss a tailored implementation of an adjoint computation that is arguably better suited for mature scientific computing libraries by avoiding the central control assumed by the original checkpointing strategy. The proposed algorithms have been adopted by the PETSc TSAdjoint library. Their performance has been demonstrated with a large-scale PDE-constrained optimization problem on a leadership-class supercomputer.


title: Self-paced Principal Component Analysis

abstract: Principal Component Analysis (PCA) has been widely used for dimensionality reduction and feature extraction. Robust PCA (RPCA), under different robust distance metrics, such as l1-norm and l2, p-norm, can deal with noise or outliers to some extent. However, real-world data may display structures that can not be fully captured by these simple functions. In addition, existing methods treat complex and simple samples equally. By contrast, a learning pattern typically adopted by human beings is to learn from simple to complex and less to more. Based on this principle, we propose a novel method called Self-paced PCA (SPCA) to further reduce the effect of noise and outliers. Notably, the complexity of each sample is calculated at the beginning of each iteration in order to integrate samples from simple to more complex into training. Based on an alternating optimization, SPCA finds an optimal projection matrix and filters out outliers iteratively. Theoretical analysis is presented to show the rationality of SPCA. Extensive experiments on popular data sets demonstrate that the proposed method can improve the state of-the-art results considerably.


title: Pastprop-RNN: improved predictions of the future by correcting the past

abstract: Forecasting accuracy is reliant on the quality of available past data. Data disruptions can adversely affect the quality of the generated model (e.g. unexpected events such as out-of-stock products when forecasting demand). We address this problem by pastcasting: predicting how data should have been in the past to explain the future better. We propose Pastprop-LSTM, a data-centric backpropagation algorithm that assigns part of the responsibility for errors to the training data and changes it accordingly. We test three variants of Pastprop-LSTM on forecasting competition datasets, M4 and M5, plus the Numenta Anomaly Benchmark. Empirical evaluation indicates that the proposed method can improve forecasting accuracy, especially when the prediction errors of standard LSTM are high. It also demonstrates the potential of the algorithm on datasets containing anomalies.


title: Can Buyers Reveal for a Better Deal?

abstract: We study market interactions in which buyers are allowed to credibly reveal partial information about their types to the seller. Previous recent work has studied the special case of one buyer and one good, showing that such communication can simultaneously improve social welfare and ex ante buyer utility. However, with multiple buyers, we find that the buyer-optimal signalling schemes from the one-buyer case are actually harmful to buyer welfare. Moreover, we prove several impossibility results showing that, with either multiple i.i.d. buyers or multiple i.i.d. goods, maximizing buyer utility can be at odds with social efficiency, which is surprising in contrast with the one-buyer, one-good case. Finally, we investigate the computational tractability of implementing desirable equilibrium outcomes. We find that, even with one buyer and one good, optimizing buyer utility is generally NP-hard but tractable in a practical restricted setting.


title: Semi-Supervised Raw-to-Raw Mapping

abstract: The raw-RGB colors of a camera sensor vary due to the spectral sensitivity differences across different sensor makes and models. This paper focuses on the task of mapping between different sensor raw-RGB color spaces. Prior work addressed this problem using a pairwise calibration to achieve accurate color mapping. Although being accurate, this approach is less practical as it requires: (1) capturing pair of images by both camera devices with a color calibration object placed in each new scene; (2) accurate image alignment or manual annotation of the color calibration object. This paper aims to tackle color mapping in the raw space through a more practical setup. Specifically, we present a semi-supervised raw-to-raw mapping method trained on a small set of paired images alongside an unpaired set of images captured by each camera device. Through extensive experiments, we show that our method achieves better results compared to other domain adaptation alternatives in addition to the single-calibration solution. We have generated a new dataset of raw images from two different smartphone cameras as part of this effort. Our dataset includes unpaired and paired sets for our semi-supervised training and evaluation.


title: Multimodal Few-Shot Learning with Frozen Language Models

abstract: When trained at sufficient scale, auto-regressive language models exhibit the notable ability to learn a new language task after being prompted with just a few examples. Here, we present a simple, yet effective, approach for transferring this few-shot learning ability to a multimodal setting (vision and language). Using aligned image and caption data, we train a vision encoder to represent each image as a sequence of continuous embeddings, such that a pre-trained, frozen language model prompted with this prefix generates the appropriate caption. The resulting system is a multimodal few-shot learner, with the surprising ability to learn a variety of new tasks when conditioned on examples, represented as a sequence of multiple interleaved image and text embeddings. We demonstrate that it can rapidly learn words for new objects and novel visual categories, do visual question-answering with only a handful of examples, and make use of outside knowledge, by measuring a single model on a variety of established and new benchmarks.


title: FCC-ee overview: new opportunities create new challenges

abstract: With its high luminosity, its clean experimental conditions, and a range of energies that cover the four heaviest particles known today, FCC-ee offers a wealth of physics possibilities, with high potential for discoveries. The FCC-ee is an essential and complementary step towards a 100 TeV hadron collider, and as such offers a uniquely powerful combined physics program. This vision is the backbone of the 2020 European Strategy for Particle Physics. One of the main challenges is now to design experimental systems that can, demonstrably, fully exploit these extraordinary opportunities.


title: The Bulk Flow Velocity and Acceleration of the Inner Jet in M\,87

abstract: A high sensitivity, 7mm Very Long Baseline Array image of M\,87 is analyzed in order to estimate the jet velocity within 0.65 mas of the point of origin. The image captured a high signal to noise, double-ridged, counter-jet extending $\sim 1$ mas from the nucleus. After defining conditions and requirements that justify approximate time averaged bilateral symmetry, a continuous set of Lorentz transformations are found that map the double-ridged counter-jet intensity profile into the double-ridged jet intensity profile. The mapping is realized by a uniformly accelerating flow with intrinsic velocity of $\sim 0.27$c at 0.4 mas (a de-projected distance of 0.38 lt-yrs) to $0.38$c at 0.65 mas (a de-projected distance of 0.61 lt-yrs) from the nucleus. Since the velocity field is derived from the global surface brightness profile and does not depend on the motion of enhanced features, it is most likely a bulk flow velocity as opposed to a pattern velocity. This interpretation is corroborated by the fact that the distribution of the apparent velocities of previously identified individual features (from the literature) within 0.65 mas of the nucleus are consistent with local hydrodynamic shocks being advected with the local bulk flow velocity. The bulk flow velocity of the visible inner jet is a constraint that can potentially break degeneracies between numerical simulations that are designed to replicate both the annulus that was imaged by the Event Horizon Telescope as well as the base of the inner jet.


title: On a Novel Effective Equation of the Reduced Hartree-Fock Theory

abstract: We show that there is an one-to-one correspondence between solutions to the Poisson-Landscape equations and the reduced Hartree-Fock equations in the semi-classical limit at low temperature. Moreover, we prove that the difference between the two corresponding solutions is small by providing explicit estimates.


title: Optimal investment and proportional reinsurance in a regime-switching market model under forward preferences

abstract: In this paper we study the optimal investment and reinsurance problem of an insurance company whose investment preferences are described via a forward dynamic exponential utility in a regime-switching market model. Financial and actuarial frameworks are dependent since stock prices and insurance claims vary according to a common factor given by a continuous time finite state Markov chain. We construct the value function and we prove that it is a forward dynamic utility. Then, we characterize the investment strategy and the optimal proportional level of reinsurance. We also perform numerical experiments and provide sensitivity analyses with respect to some model parameters.


title: Charged Defects and Phonon Hall Effects in Ionic Crystals

abstract: It has been known for decades that a magnetic field can deflect phonons as they flow in response to a thermal gradient, producing a thermal Hall effect. Several recent experiments have revealed ratios of the phonon Hall conductivity to the phonon longitudinal conductivity in oxide dielectrics that are larger than $10^{-3}$ when phonon mean-free-paths exceed phonon wavelengths. At the same time $\kappa_{H}/\kappa_{L}$ is not strongly temperature dependent. We argue that these two properties together imply a mechanism related to phonon scattering from defects that break time-reversal symmetry, and we show that Lorentz forces acting on charged defects produce substantial skew-scattering amplitudes, and related thermal Hall effects that are consistent with recent observations.


title: Deconstruction and reconstruction of image-degrading effects in the human abdomen using Fullwave: phase aberration, multiple reverberation, and trailing reverberation

abstract: Ultrasound image degradation in the human body is complex and occurs due to the distortion of the wave as it propagates to and from the target. Here, we establish a simulation based framework that deconstructs the sources of image degradation into a separable parameter space that includes phase aberration from speed variation, multiple reverberations, and trailing reverberation. These separable parameters are then used to reconstruct images with known and independently modulable amounts of degradation using methods that depend on the additive or multiplicative nature of the degradation. Experimental measurements and Fullwave simulations in the human abdomen demonstrate this calibrated process in abdominal imaging by matching relevant imaging metrics such as phase aberration, reverberation strength, speckle brightness and coherence length. Applications of the reconstruction technique are illustrated for beamforming strategies (phase aberration correction, spatial coherence imaging), in a standard abdominal environment, as well as in impedance ranges much higher than those naturally occurring in the body.


title: Test Case Selection and Prioritization Using Machine Learning: A Systematic Literature Review

abstract: Regression testing is an essential activity to assure that software code changes do not adversely affect existing functionalities. With the wide adoption of Continuous Integration (CI) in software projects, which increases the frequency of running software builds, running all tests can be time-consuming and resource-intensive. To alleviate that problem, Test case Selection and Prioritization (TSP) techniques have been proposed to improve regression testing by selecting and prioritizing test cases in order to provide early feedback to developers. In recent years, researchers have relied on Machine Learning (ML) techniques to achieve effective TSP (ML-based TSP). Such techniques help combine information about test cases, from partial and imperfect sources, into accurate prediction models. This work conducts a systematic literature review focused on ML-based TSP techniques, aiming to perform an in-depth analysis of the state of the art, thus gaining insights regarding future avenues of research. To that end, we analyze 29 primary studies published from 2006 to 2020, which have been identified through a systematic and documented process. This paper addresses five research questions addressing variations in ML-based TSP techniques and feature sets for training and testing ML models, alternative metrics used for evaluating the techniques, the performance of techniques, and the reproducibility of the published studies. We summarize the results related to our research questions in a high-level summary that can be used as a taxonomy for classifying future TSP studies.


title: Light in correlated disordered media

abstract: The optics of correlated disordered media is a fascinating research topic emerging at the interface between the physics of waves in complex media and nanophotonics. Inspired by photonic structures in nature and enabled by advances in nanofabrication processes, recent investigations have unveiled how the design of structural correlations down to the subwavelength scale could be exploited to control the scattering, transport and localization of light in matter. From optical transparency to superdiffusive light transport to photonic gaps, the optics of correlated disordered media challenges our physical intuition and offers new perspectives for applications. This article reviews the theoretical foundations, state-of-the-art experimental techniques and major achievements in the study of light interaction with correlated disorder, covering a wide range of systems -- from short-range correlated photonic liquids, to L\'evy glasses containing fractal heterogeneities, to hyperuniform disordered photonic materials. The mechanisms underlying light scattering and transport phenomena are elucidated on the basis of rigorous theoretical arguments. We overview the exciting ongoing research on mesoscopic phenomena, such as transport phase transitions and speckle statistics, and the current development of disorder engineering for applications such as light-energy management and visual appearance design. Special efforts are finally made to identify the main theoretical and experimental challenges to address in the near future.

categories: physics.optics cond-mat.dis-nn cond-mat.mes-hall cond-mat.soft


title: Observational optical constraints of regular black holes

abstract: In this work, we consider two recently introduced novel regular black hole solutions and investigate the circular null geodesics to find the connection between the photon sphere, the horizon and the black hole shadow radii. We also study the energy emission rate for these solutions and discuss how the parameters of models affect the emission of particles around the black holes. Furthermore, we compare the resulting shadow of these regular black holes with observational data of the Event Horizon Telescope and find the allowed regions of the model parameters for which the obtained shadow is consistent with the data. Finally, we employ the correspondence between the quasinormal modes in the eikonal limit and shadow radius to study the scalar field perturbations in these backgrounds.


title: Extreme Nonnegative Quadratics over Stanley Reisner Varieties

abstract: We consider the convex geometry of the cone of nonnegative quadratics over Stanley-Reisner varieties. Stanley-Reisner varieties (which are unions of coordinate planes) are amongst the simplest real projective varieties, so this is potentially a starting point that can generalize to more complicated real projective varieties. This subject has some suprising connections to algebraic topology and category theory, which we exploit heavily in our work. These questions are also valuable in applied math, because they directly translate to questions about positive semidefinite (PSD) matrices. In particular, this relates to a long line of work concerning the extent to which it is possible to approximately check that a matrix is PSD by checking that some principle submatrices are PSD, or to check if a partial matrix can be approximately completed to full PSD matrix. We systematize both these practical and theoretical questions using a framework based on algebraic topology, category theory, and convex geometry. As applications of this framework we are able to classify the extreme nonnegative quadratics over many Stanley-Reisner varieties. We plan to follow these structural results with a paper that is more focused on quantitative questions about PSD matrix completion, which have applications in sparse semidefinite programming.


title: Knowledge Infused Policy Gradients with Upper Confidence Bound for Relational Bandits

abstract: Contextual Bandits find important use cases in various real-life scenarios such as online advertising, recommendation systems, healthcare, etc. However, most of the algorithms use flat feature vectors to represent context whereas, in the real world, there is a varying number of objects and relations among them to model in the context. For example, in a music recommendation system, the user context contains what music they listen to, which artists create this music, the artist albums, etc. Adding richer relational context representations also introduces a much larger context space making exploration-exploitation harder. To improve the efficiency of exploration-exploitation knowledge about the context can be infused to guide the exploration-exploitation strategy. Relational context representations allow a natural way for humans to specify knowledge owing to their descriptive nature. We propose an adaptation of Knowledge Infused Policy Gradients to the Contextual Bandit setting and a novel Knowledge Infused Policy Gradients Upper Confidence Bound algorithm and perform an experimental analysis of a simulated music recommendation dataset and various real-life datasets where expert knowledge can drastically reduce the total regret and where it cannot.


title: Optical MIMO Communication Using Holographic Spectral Multiplexing of Pulsed Ultrashort Laser

abstract: In this paper, we introduce Holographic Spectral Multiplexing (HSM) as a novel technique to enable multiple-input multiple-output (MIMO) communication in optical networks. HSM uses the spectral space of ultrashort laser pulses to create line codes in the form of 2D holograms. The pulse processing is performed in the temporal Fourier domain by spatially dispersing the pulse frequency components in a spectral processing device (SPD). The 2D holograms are composed of the patterns of intensity disparities that an SLM inscribes on the spectrally-decomposed Fourier plane of the pulse. The holographic line codes defined in this way transform the ultrashort laser pulses into high-entropy data symbols, hence, enhance the communication's spectral efficiency. Unlike conventional optical multiplexing schemes (e.g., TDM, WDM, or SDM), HSM does not physically or abstractly separate the communication propagation space into subchannels. Rather, HSM realizes a MIMO communication paradigm by allowing the photonic waves under the pulse envelope to propagate in the same space so they scatter and interfere by chromatic dispersion. This allows HSM to form beams between the pixels of SLM at the sender and receiver sides and optimize the beam to adapt to channel scattering situations. In this way, HSM delivers a rate gain that in the best case exponentially increases the information rate of communication.


title: Implicit Gradient Alignment in Distributed and Federated Learning

abstract: A major obstacle to achieving global convergence in distributed and federated learning is the misalignment of gradients across clients, or mini-batches due to heterogeneity and stochasticity of the distributed data. In this work, we show that data heterogeneity can in fact be exploited to improve generalization performance through implicit regularization. One way to alleviate the effects of heterogeneity is to encourage the alignment of gradients across different clients throughout training. Our analysis reveals that this goal can be accomplished by utilizing the right optimization method that replicates the implicit regularization effect of SGD, leading to gradient alignment as well as improvements in test accuracies. Since the existence of this regularization in SGD completely relies on the sequential use of different mini-batches during training, it is inherently absent when training with large mini-batches. To obtain the generalization benefits of this regularization while increasing parallelism, we propose a novel GradAlign algorithm that induces the same implicit regularization while allowing the use of arbitrarily large batches in each update. We experimentally validate the benefits of our algorithm in different distributed and federated learning settings.


title: Closed-form Continuous-time Neural Models

abstract: Continuous-time neural processes are performant sequential decision-makers that are built by differential equations (DE). However, their expressive power when they are deployed on computers is bottlenecked by numerical DE solvers. This limitation has significantly slowed down the scaling and understanding of numerous natural physical phenomena such as the dynamics of nervous systems. Ideally, we would circumvent this bottleneck by solving the given dynamical system in closed form. This is known to be intractable in general. Here, we show it is possible to closely approximate the interaction between neurons and synapses -- the building blocks of natural and artificial neural networks -- constructed by liquid time-constant networks (LTCs) efficiently in closed-form. To this end, we compute a tightly-bounded approximation of the solution of an integral appearing in LTCs' dynamics, that has had no known closed-form solution so far. This closed-form solution substantially impacts the design of continuous-time and continuous-depth neural models; for instance, since time appears explicitly in closed-form, the formulation relaxes the need for complex numerical solvers. Consequently, we obtain models that are between one and five orders of magnitude faster in training and inference compared to differential equation-based counterparts. More importantly, in contrast to ODE-based continuous networks, closed-form networks can scale remarkably well compared to other deep learning instances. Lastly, as these models are derived from liquid networks, they show remarkable performance in time series modeling, compared to advanced recurrent models.


title: Domain Conditional Predictors for Domain Adaptation

abstract: Learning guarantees often rely on assumptions of i.i.d. data, which will likely be violated in practice once predictors are deployed to perform real-world tasks. Domain adaptation approaches thus appeared as a useful framework yielding extra flexibility in that distinct train and test data distributions are supported, provided that other assumptions are satisfied such as covariate shift, which expects the conditional distributions over labels to be independent of the underlying data distribution. Several approaches were introduced in order to induce generalization across varying train and test data sources, and those often rely on the general idea of domain-invariance, in such a way that the data-generating distributions are to be disregarded by the prediction model. In this contribution, we tackle the problem of generalizing across data sources by approaching it from the opposite direction: we consider a conditional modeling approach in which predictions, in addition to being dependent on the input data, use information relative to the underlying data-generating distribution. For instance, the model has an explicit mechanism to adapt to changing environments and/or new data sources. We argue that such an approach is more generally applicable than current domain adaptation methods since it does not require extra assumptions such as covariate shift and further yields simpler training algorithms that avoid a common source of training instabilities caused by minimax formulations, often employed in domain-invariant methods.


title: On Preserving the Behavior in Software Refactoring: A Systematic Mapping Study

abstract: Context: Refactoring is the art of modifying the design of a system without altering its behavior. The idea is to reorganize variables, classes and methods to facilitate their future adaptations and comprehension. As the concept of behavior preservation is fundamental for refactoring, several studies, using formal verification, language transformation and dynamic analysis, have been proposed to monitor the execution of refactoring operations and their impact on the program semantics. However, there is no existing study that examines the available behavior preservation strategies for each refactoring operation. Objective: This paper identifies behavior preservation approaches in the research literature. Method: We conduct, in this paper, a systematic mapping study, to capture all existing behavior preservation approaches that we classify based on several criteria including their methodology, applicability, and their degree of automation. Results: The results indicate that several behavior preservation approaches have been proposed in the literature. The approaches vary between using formalisms and techniques, developing automatic refactoring safety tools, and performing a manual analysis of the source code. Conclusion: Our taxonomy reveals that there exist some types of refactoring operations whose behavior preservation is under-researched. Our classification also indicates that several possible strategies can be combined to better detect any violation of the program semantics.


title: Building Bridges: Generative Artworks to Explore AI Ethics

abstract: In recent years, there has been an increased emphasis on understanding and mitigating adverse impacts of artificial intelligence (AI) technologies on society. Across academia, industry, and government bodies, a variety of endeavours are being pursued towards enhancing AI ethics. A significant challenge in the design of ethical AI systems is that there are multiple stakeholders in the AI pipeline, each with their own set of constraints and interests. These different perspectives are often not understood, due in part to communication gaps.For example, AI researchers who design and develop AI models are not necessarily aware of the instability induced in consumers' lives by the compounded effects of AI decisions. Educating different stakeholders about their roles and responsibilities in the broader context becomes necessary. In this position paper, we outline some potential ways in which generative artworks can play this role by serving as accessible and powerful educational tools for surfacing different perspectives. We hope to spark interdisciplinary discussions about computational creativity broadly as a tool for enhancing AI ethics.


title: The conjugate gradient algorithm on a general class of spiked covariance matrices

abstract: We consider the conjugate gradient algorithm applied to a general class of spiked sample covariance matrices. The main result of the paper is that the norms of the error and residual vectors at any finite step concentrate on deterministic values determined by orthogonal polynomials with respect to a deformed Marchenko--Pastur law. The first-order limits and fluctuations are shown to be universal. Additionally, for the case where the bulk eigenvalues lie in a single interval we show a stronger universality result in that the asymptotic rate of convergence of the conjugate gradient algorithm only depends on the support of the bulk, provided the spikes are well-separated from the bulk. In particular, this shows that the classical condition number bound for the conjugate gradient algorithm is pessimistic for spiked matrices.


title: Eigenvalue Estimates for $p$-Laplace Problems on Domains Expressed in Fermi Coordinates

abstract: We prove explicit and sharp eigenvalue estimates for Neumann $p$-Laplace eigenvalues in domains that admit a representation in Fermi coordinates. More precisely, if $\gamma$ denotes a non-closed curve in $\mathbb{R}^2$ symmetric with respect to the $y$-axis, let $D\subset \mathbb{R}^2$ denote the domain of points that lie on one side of $\gamma$ and within a prescribed distance $\delta(s)$ from $\gamma(s)$ (here $s$ denotes the arc length parameter for $\gamma$). Write $\mu_1^{odd}(D)$ for the lowest nonzero eigenvalue of the Neumann $p$-Laplacian with an eigenfunction that is odd with respect to the $y$-axis. For all $p>1$, we provide a lower bound on $\mu_1^{odd}(D)$ when the distance function $\delta$ and the signed curvature $k$ of $\gamma$ satisfy certain geometric constraints. In the linear case ($p=2$), we establish sufficient conditions to guarantee $\mu_1^{odd}(D)=\mu_1(D)$. We finally study the asymptotics of $\mu_1(D)$ as the distance function tends to zero. We show that in the limit, the eigenvalues converge to the lowest nonzero eigenvalue of a weighted one-dimensional Neumann $p$-Laplace problem.


title: An Intriguing Globular Cluster in the Galactic Bulge from the VVV Survey

abstract: Recent near-IR Surveys have discovered a number of new bulge globular cluster (GC) candidates that need to be further investigated. Our main objective is to use public data from the Gaia Mission, VVV, 2MASS and WISE in order to measure the physical parameters of Minni48, a new candidate GC located in the inner bulge of the Galaxy at l=359.35 deg, b=2.79 deg. Even though there is a bright foreground star contaminating the field, the cluster appears quite bright in near- and mid-IR images. We obtain deep decontaminated optical and near-IR colour-magnitude diagrams (CMDs) for this cluster. The heliocentric cluster distance is determined from the red clump (RC) and the red giant branch (RGB) tip magnitudes in the near-IR CMD, while the cluster metallicity is estimated from the RGB slope and the fit to theoretical isochrones. The GC size is found to be r = 6' +/- 1', while reddening and extinction values are E(J-Ks)=0.60 +/- 0.05 mag, A_G=3.23 +/- 0.10 mag, A_Ks=0.45 +/- 0.05 mag. The resulting mean Gaia proper motions are PMRA=-3.5 +/- 0.5 mas/yr, PMDEC=-6.0 +/- 0.5 mas/yr. The IR magnitude of the RC yields an accurate distance modulus estimate of (m-M)_0=14.61 mag, equivalent to a distance D=8.4 +/- 1.0 kpc. This is consistent with the optical distance estimate: (m-M)_0=14.67 mag, D=8.6 +/- 1.0 kpc, and with the RGB tip distance: (m-M)_0=14.45 mag, D=7.8 +/- 1.0 kpc. The derived metallicity is [Fe/H]=-0.20 +/- 0.30 dex. A good fit to the PARSEC stellar isochrones is obtained in all CMDs using Age = 10 +/- 2 Gyr. The total absolute magnitude of this GC is estimated to be M_Ks= -9.04 +/- 0.66 mag. Based on its position, kinematics, metallicity and age, we conclude that Minni48 is a genuine GC, similar to other well known metal-rich bulge GCs. It is located at a projected Galactocentric angular distance of 2.9 deg, equivalent to 0.4 kpc, being one of the closest GCs to the Galactic centre.


title: Approximation Schemes for Path Integration on Riemannian Manifolds

abstract: In this paper, we prove a finite dimensional approximation scheme for the Wiener measure on closed Riemannian manifolds, establishing a generalization for $L^{1}$-functionals, of the approach followed by Andersson and Driver on [1]. We follow a new approach motived by categorical concept of colimit.


title: Compositional Reinforcement Learning from Logical Specifications

abstract: We study the problem of learning control policies for complex tasks given by logical specifications. Recent approaches automatically generate a reward function from a given specification and use a suitable reinforcement learning algorithm to learn a policy that maximizes the expected reward. These approaches, however, scale poorly to complex tasks that require high-level planning. In this work, we develop a compositional learning approach, called DiRL, that interleaves high-level planning and reinforcement learning. First, DiRL encodes the specification as an abstract graph; intuitively, vertices and edges of the graph correspond to regions of the state space and simpler sub-tasks, respectively. Our approach then incorporates reinforcement learning to learn neural network policies for each edge (sub-task) within a Dijkstra-style planning algorithm to compute a high-level plan in the graph. An evaluation of the proposed approach on a set of challenging control benchmarks with continuous state and action spaces demonstrates that it outperforms state-of-the-art baselines.


title: Discovery of 13 New Orbital Periods for Classical Novae

abstract: I report new orbital periods (P) for 13 classical novae, based on light curves from TESS, AAVSO, and other public archives. These new nova periods now constitute nearly one-seventh of all known nova periods. Five of my systems have P>1 day, which doubles the number of such systems that must have evolved companion stars. (This is simply because ground-based time series have neither the coverage nor the stability required to discover these small-amplitude long periods.) V1016 Sgr has a rare P below the period gap, and suddenly becomes useful for current debates on evolution of novae. Five of the novae (FM Cir, V399 Del, V407 Lup, YZ Ret, and V549 Vel) have the orbital modulations in the tail of the eruption after the transition phase. Soon after the transition, YZ Ret shows a unique set of aperiodic diminishing oscillations, plus YZ Ret shows two highly-significant transient periods, 1.1% and 4.5% longer than the orbital period, much like for the superhump phenomenon. I also report an optical 591.27465 second periodicity for V407 Lup, which is coherent and must be tied to the white dwarf spin period. The new orbital periods in days are 0.1883907 +- 0.0000048 for V1405 Cas, 3.4898 +- 0.0072 for FM Cir, 0.162941 +- 0.000060 for V339 Del, 3.513 +- 0.020 for V407 Lup, 1.32379 +- 0.00048 for V2109 Oph, 3.21997 +- 0.00039 for V392 Per, 0.1628714 +- 0.0000110 for V598 Pup, 0.1324539 +- 0.0000098 for YZ Ret, 0.07579635 +- 0.00000017 for V1016 Sgr, 7.101 +- 0.016 for V5583 Sgr, 0.61075 +- 0.00071 for V1534 Sco, 0.40319 +- 0.00005 for V549 Vel, and 0.146501 +- 0.000058 for NQ Vul.


title: Towards auto-completion on software requirements statements

abstract: As software systems become more complex, modern software development requires more attention to human perspectives, and active participation of development teams in requirements elicitation tasks. In this context, incomplete or ambiguous requirements descriptions do not guide the development of good software products. We hypothesize that the text auto-completion feature improves the quality of the software requirements artifacts. We present the motivation for this study, related works, our approach and future research efforts.