title: Shared Control Between Pilots and Autopilots: Illustration of a Cyber-Physical Human System

abstract: Although increased automation has made it easier to control aircraft, ensuring a safe interaction between the pilots and the autopilots is still a challenging problem, especially in the presence of severe anomalies. Current approach consists of autopilot solutions that disengage themselves when they become ineffective. This may cause reengagement of the pilot at the worst possible time, which can result in undesired consequences. In this paper, a series of research studies that propose pilot-autopilot interaction schemes based on the Capacity for Maneuver (CfM) concept, are covered. CfM refers to the remaining capacity of the actuators that can be used for bringing the aircraft to safety. It is demonstrated that CfM-based pilot-autopilot interaction schemes, or Shared Control Architectures (SCA), can be promising alternatives to the existing schemes. Two SCAs are tested in the experiments. In SCA1, the pilot takes over the control from the autopilot using the monitored CfM information. In SCA2, the pilot takes on the role of a supervisor helping the autopilot whenever a need arises, based on the CfM information. Whenever the aircraft experiences a severe anomaly, the pilot assesses the situation based on his/her CfM monitoring and intervenes by providing two control system parameter estimates to the autopilot. This increases the effectiveness of the autopilot. Using human-in-the-loop simulations, it is shown that CfM based interactions provides smaller tracking errors and larger overall CfM. The subjects including a commercial airline pilot and several university students are trained using a systematic procedure. Creation of new cyber physical & human systems is inevitable along with deeper engagement with the social science community so as to get better insight into human decision making. The results reported here should be viewed as a first of several steps in this research direction.


title: A proof of Kolyvagin's Conjecture via the BDP main conjecture

abstract: We adapt Wei Zhang's proof of Kolyvagin's conjecture for modular abelian varieties over $\mathbb{Q}$ to rely on the BDP main conjecture instead of on the cyclotomic main conjecture. The main ingredient is a reduction to a case that is tractable by the BDP main conjecture, in a similar spirit to Zhang's reduction to the rank one case. By using the BDP main conjecture instead of the cyclotomic main conjecture, our approach is more suitable than Zhang's to extend to modular abelian varieties over totally real fields.


title: Two-Sample Test Based on Classification Probability

abstract: Robust classification algorithms have been developed in recent years with great success. We take advantage of this development and recast the classical two-sample test problem in the framework of classification. Based on the estimates of classification probabilities from a classifier trained from the samples, a test statistic is proposed. We explain why such a test can be a powerful test and compare its performance in terms of the power and efficiency with those of some other recently proposed tests with simulation and real-life data. The test proposed is nonparametric and can be applied to complex and high dimensional data wherever there is a classifier that provides consistent estimate of the classification probability for such data.


title: The value of knowing the market price of risk

abstract: This paper presents an optimal allocation problem in a financial market with one risk-free and one risky asset, when the market is driven by a stochastic market price of risk. We solve the problem in continuous time, for an investor with a Constant Relative Risk Aversion (CRRA) utility, under two scenarios: when the market price of risk is observable (the {\em full information case}), and when it is not (the {\em partial information case}). The corresponding market models are complete in the partial information case and incomplete in the other case, hence the two scenarios exhibit rather different features. We study how the access to more accurate information on the market price of risk affects the optimal strategies and we determine the maximal price that the investor would be willing to pay to get such information. In particular, we examine two cases of additional information, when an exact observation of the market price of risk is available either at time $0$ only (the {\em initial information case}), or during the whole investment period (the {\em dynamic information case}).


title: Spatial Scales and Locality of Magnetic Helicity: Part 1

abstract: Magnetic helicity is approximately conserved in resistive MHD models. It quantifies the entanglement of the magnetic field within the plasma. The transport and removal of helicity is crucial in both the dynamo in the solar interior and active region evolution in the solar corona. This transport typically leads to highly inhomogeneous distributions of entanglement.There exists no consistent systematic means of decomposing helicity over varying spatial scales and in localised regions. We apply a multiresolution wavelet decomposition to the magnetic field and demonstrate how it can be applied to various quantities associated with magnetic helicity, including the field line helicity. We use a geometrical definition of helicity which allows these quantities to be calculated for fields with arbitrary boundary conditions. It is shown that the multiresolution decomposition of helicity has the crucial property of local additivity and demonstrate a general linear energy-topology conservation law which is a significant generalisation of the two point correlation decomposition used in the analysis of homogeneous turbulence and periodic fields. The localisation property of the wavelet representation is shown to characterise inhomogeneous distributions which a Fourier representation cannot. Using an analytic representation of a resistive braided field relaxation we demonstrate a clear correlation between the variations in energy at various length scales and the variations in helicity at the same spatial scales. Its application to helicity flows in a surface flux transport model show how various contributions to the global helicity input from active region field evolution and polar field development are naturally separated by this representation.


title: Entirety of Quantum Uncertainty and Its Experimental Verification

abstract: As a foundation of modern physics, uncertainty relations describe an ultimate limit for the measurement uncertainty of incompatible observables. Traditionally, uncertain relations are formulated by mathematical bounds for a specific state. Here we present a method for geometrically characterizing uncertainty relations as an entire area of variances of the observables, ranging over all possible input states. We find that for the pair of position $x$ and momentum $p$ operators, Heisenberg's uncertainty principle points exactly to the area of the variances of $x$ and $p$. Moreover, for finite-dimensional systems, we prove that the corresponding area is necessarily semialgebraic; in other words, this set can be represented via finite polynomial equations and inequalities, or any finite union of such sets. In particular, we give the analytical characterization of the areas of variances of (a) a pair of one-qubit observables, (b) a pair of projective observables for arbitrary dimension, and give the first experimental observation of such areas in a photonic system.


title: Realization of Universal Quantum Gates with Spin-Qudits in Colloidal Quantum Dots

abstract: We exploit hyperfine interactions in a single Mn-ion confined in a quantum dot (QD) to create a qudit, i.e. a multi-level quantum-bit system, with well defined, addressable and robust set of spin states for the realization of universal quantum gates. We generate and probe an arbitrary superposition of states between selected hyperfine energy level pairs by using electron double resonance detected nuclear magnetic resonance (EDNMR). This enables the observation of Rabi oscillations and the experimental realization of NOT and SWAP universal quantum gates that are robust against decoherence. Our protocol for cyclical preparation, manipulation and read-out of logic gates offers opportunities for integration of qudits in scalable quantum circuit architectures beyond solid state electron spin qubits.


title: A descriptive Main Gap Theorem

abstract: Answering one of the main questions of [FHK14, Chapter 7], we show that there is a tight connection between the depth of a classifiable shallow theory $T$ and the Borel rank of the isomorphism relation $\cong^\kappa_T$ on its models of size $\kappa$, for $\kappa$ any cardinal satisfying $\kappa^{< \kappa} = \kappa > 2^{\aleph_0}$. This is achieved by establishing a link between said rank and the $\mathcal{L}_{\infty \kappa}$-Scott height of the $\kappa$-sized models of $T$, and yields to the following descriptive set-theoretical analogue of Shelah's Main Gap Theorem: Given a countable complete first-order theory $T$, either $\cong^\kappa_T$ is Borel with a countable Borel rank (i.e. very simple, given that the length of the relevant Borel hierarchy is $\kappa^+ > \aleph_1$), or it is not Borel at all. The dividing line between the two situations is the same as in Shelah's theorem, namely that of classifiable shallow theories. We also provide a Borel reducibility version of the above theorem, discuss some limitations to the possible (Borel) complexities of $\cong^\kappa_T$, and provide a characterization of categoricity of $T$ in terms of the descriptive set-theoretical complexity of $\cong^\kappa_T$.


title: Spatial and spectral dynamics in STEM hyperspectral imaging using random scan patterns

abstract: The evolution of the scanning modules for scanning transmission electron microscopes (STEM) has realized the possibility to generate arbitrary scan pathways, an approach currently explored to improve acquisition speed and to reduce electron dose effects. In this work, we present the implementation of a random scan operating mode in STEM achieved at the hardware level via a custom scan control module. A pre-defined pattern with fully shuffled raster order is used to sample the entire region of interest. Subsampled random sparse images can then be extracted at successive time frames, to which suitable image reconstruction techniques can be applied. With respect to the conventional raster scan mode, this method permits to limit dose accumulation effects, but also to decouple the spatial and temporal information in hyperspectral images. We provide some proofs of concept of the flexibility of the random scan operating mode, presenting examples of its applications in different spectro-microscopy contexts: atomically-resolved elemental maps with electron energy loss spectroscopy and nanoscale-cathodoluminescence spectrum images. By employing adapted post-processing tools, it is demonstrated that the method allows to precisely track and correct for sample instabilities and to follow spectral diffusion with a high spatial localization.


title: Active Learning for Risk-Sensitive Inverse Reinforcement Learning

abstract: One typical assumption in inverse reinforcement learning (IRL) is that human experts act to optimize the expected utility of a stochastic cost with a fixed distribution. This assumption deviates from actual human behaviors under ambiguity. Risk-sensitive inverse reinforcement learning (RS-IRL) bridges such gap by assuming that humans act according to a random cost with respect to a set of subjectively distorted distributions instead of a fixed one. Such assumption provides the additional flexibility to model human's risk preferences, represented by a risk envelope, in safe-critical tasks. However, like other learning from demonstration techniques, RS-IRL could also suffer inefficient learning due to redundant demonstrations. Inspired by the concept of active learning, this research derives a probabilistic disturbance sampling scheme to enable an RS-IRL agent to query expert support that is likely to expose unrevealed boundaries of the expert's risk envelope. Experimental results confirm that our approach accelerates the convergence of RS-IRL algorithms with lower variance while still guaranteeing unbiased convergence.


title: Proper Kasparov Cycles and the Baum-Connes Conjecture

abstract: We introduce the notion of proper Kasparov cycles for Kasparov's G-equivariant KK-theory for a general locally compact, second countable topological group G. We show that for any proper Kasparov cycle, its induced map on K-theory factors through the left-hand side of the Baum-Connes conjecture. This allows us to upgrade the direct splitting method, a recent new approach to the Baum-Connes conjecture which, in contrast to the standard gamma element method (the Dirac dual-Dirac method), avoids the need of constructing proper algebras and the Dirac and the dual-Dirac elements. We introduce the notion of Kasparov cycles with Property (gamma) removing the G-compact assumption on the universal space EG in the previous paper "Direct Splitting Method for the Baum-Connes Conjecture". We show that the existence of a cycle with Property (gamma) implies the split-injectivity of the Baum-Connes assembly map for all coefficients. We also obtain results concerning the surjectivity of the assembly map.


title: First Integrals vs Limit Cycles

abstract: This paper applies a recent result determining periodic orbits on the basis of first integrals, for Li\'enard systems. By solving a first order ODE with singularities, a crucial result is proved to locate intervals of single and isolated maximum amplitudes periodic orbits (limit cycles). With this result an upper bound for the number of limit cycles is provided. Some examples are presented along with conclusions and future work


title: Multimodal Multitask Representation Learning for Pathology Biobank Metadata Prediction

abstract: Metadata are general characteristics of the data in a well-curated and condensed format, and have been proven to be useful for decision making, knowledge discovery, and also heterogeneous data organization of biobank. Among all data types in the biobank, pathology is the key component of the biobank and also serves as the gold standard of diagnosis. To maximize the utility of biobank and allow the rapid progress of biomedical science, it is essential to organize the data with well-populated pathology metadata. However, manual annotation of such information is tedious and time-consuming. In the study, we develop a multimodal multitask learning framework to predict four major slide-level metadata of pathology images. The framework learns generalizable representations across tissue slides, pathology reports, and case-level structured data. We demonstrate improved performance across all four tasks with the proposed method compared to a single modal single task baseline on two test sets, one external test set from a distinct data source (TCGA) and one internal held-out test set (TTH). In the test sets, the performance improvements on the averaged area under receiver operating characteristic curve across the four tasks are 16.48% and 9.05% on TCGA and TTH, respectively. Such pathology metadata prediction system may be adopted to mitigate the effort of expert annotation and ultimately accelerate the data-driven research by better utilization of the pathology biobank.


title: Impact of magnetic field on the gas mass fraction of galaxy clusters

abstract: Magnetic fields have been observed in galaxy clusters with strengths of the order of $\sim \mu$G. The non-thermal pressure exerted by magnetic fields also contributes to the total pressure in galaxy clusters and can in turn affect the estimates of the gas mass fraction, $f_{gas}$. In this paper, we have considered a central magnetic field strength of $5\mu$G, motivated by observations and simulations of galaxy clusters. The profile of the magnetic field has also been taken from the results obtained from simulations and observations. The role of magnetic field has been taken into account in inferring the gas density distribution through the hydrostatic equilibrium condition (HSE) by including the magnetic pressure. We have found that the resultant gas mass fraction is smaller with magnetic field as compared to that without magnetic field. However, this decrease is dependent on the strength and the profile of the magnetic field. We have also determined the total mass using the NFW profile to check for the dependency of $f_{gas}$ estimates on total mass estimators. From our analysis, we conclude that for the magnetic field strength that galaxy clusters seem to possess, the non-thermal pressure from magnetic fields has an impact of $\approx 1~\%$ on the gas mass fraction of galaxy clusters. However, with upcoming facilities like Square Kilometre Array (SKA), it can be further expected to improve with more precise observations of the magnetic field strength and profile in galaxy clusters, particularly in the interior region.


title: Ultra-subwavelength focusing and giant magnetic-field enhancement in a low-loss one-way waveguide based on remanence

abstract: The subwavelength focusing based on surface plasmon polaritons (SPP) has been widely explored in tapered metallic structures. However, the efficiency of energy localization is relatively weak, largely due to high propagation loss and strong back reflection. Here, we propose a straight-tapered 3-dimensional (3D) one-way surface magnetoplasmon (SMP) waveguide with the ending surface of $\sim 10^{-4}\lambda_0 \times 10^{-4}\lambda_0$ to achieve energy focusing in the microwave regime. Due to low propagation loss of SMP, we achieve huge magnetic field enhancement in such an ultra-subwavelength area, by five orders of magnitude. Instead of using an external static magnetic field, our proposed SMP waveguide relies on remanence, which is very convenient for operating practical 3D applications. These results show promising applications in magnetic-field enhancing or quenching fluorescence, luminescence or nonlinearity of 2D materials, novel scanning near-field microwave microscopy and energy storage.


title: Convergence of the random Abelian sandpile

abstract: We prove that Abelian sandpiles with random initial states converge almost surely to unique scaling limits. The proof follows the Armstrong-Smart program for stochastic homogenization of uniformly elliptic equations. Using simple random walk estimates, we prove an analogous result for the divisible sandpile and identify its scaling limit as exactly that of the averaged divisible sandpile. As a corollary, this gives a new quantitative proof of known results on the stabilizability of Abelian sandpiles.


title: Gallium phosphide as a piezoelectric platform for quantum optomechanics

abstract: Recent years have seen extraordinary progress in creating quantum states of mechanical oscillators, leading to great interest in potential applications for such systems in both fundamental as well as applied quantum science. One example is the use of these devices as transducers between otherwise disparate quantum systems. In this regard, a promising approach is to build integrated piezoelectric optomechanical devices, that are then coupled to microwave circuits. Optical absorption, low quality factors and other challenges have up to now prevented operation in the quantum regime, however. Here, we design and characterize such a piezoelectric optomechanical device fabricated from gallium phosphide in which a 2.9~GHz mechanical mode is coupled to a high quality factor optical resonator in the telecom band. The large electronic bandgap and the resulting low optical absorption of this new material, on par with devices fabricated from silicon, allows us to demonstrate quantum behavior of the structure. This not only opens the way for realizing noise-free quantum transduction between microwaves and optics, but in principle also from various color centers with optical transitions in the near visible to the telecom band.


title: Adaptive Leader-Following Consensus for Multiple Euler-Lagrange Systems with an Uncertain Leader System

abstract: In this paper, we study the leader-following consensus problem of multiple Euler-Lagrange systems subject to an uncertain leader system. We first establish an adaptive distributed observer for a neutrally stable linear leader system whose system matrix is not known exactly. Under standard assumptions, this adaptive distributed observer can estimate and pass the leader's state to each follower through the communication network of the system without knowing the leader's system matrix exactly. Under the additional assumption that the leader's state is persistently exciting, this adaptive distributed observer can also asymptotically learn the parameters of the leader's system matrix. On the basis of this adaptive distributed observer, we further synthesize an adaptive distributed control law to solve our problem via the certainty equivalence principle. Our result allows the leader-following consensus problem of multiple Euler-Lagrange systems to be solved even if none of the followers knows the system matrix of the leader system exactly.


title: Dynamical Transition from a Naked Singularity to a Black Hole

abstract: We show that a Reissner-Nordstr\"{o}m (RN) black hole can be formed by dropping a charged thin dust shell onto a RN naked singularity. This is in contrast to the fact that a RN naked singularity is prohibited from forming by dropping a charged thin dust shell onto a RN black hole. This implies the strong tendency of the RN singularity to be covered by a horizon in favour of cosmic censorship. We show that an extreme RN black hole can also be formed from a RN naked singularity by the same process in a finite advanced time. We also discuss the evolution of the charged thin dust shells and the causal structure of the resultant spacetimes.


title: Dynamic description of ternary and quaternary splits of heavy nuclear systems in the deep-inelastic regime

abstract: Colliding heavy nuclear systems in the deep-inelastic regime may undergo partitioning into multiple fragments when fusion can not be achieved. While multiple breakups are common at Fermi energy, they are rather exotic in the deep-inelastic regime, where density, excitation and, in general, transport conditions, are expected to be different. Abundant ternary and quaternary splits have been observed in recent experiments, for instance in symmetric semi-central and semi-peripheral collisions with heavy systems, like $^{197}\mathrm{Au}$ + $^{197}\mathrm{Au}$ at 15 MeV per nucleon. In these conditions, we undertook a microscopic description of the reaction dynamics. Relying on the full solution of the Boltzmann-Langevin equation implemented in the BLOB approach, we could follow in time the development of instabilities along deformation.


title: Trade-offs in dynamic coloring for bipartite and general graphs

abstract: We present trade-offs in the incremental and fully dynamic settings to maintian a proper coloring. For any fully dynamic $2$-coloring algorithm, the maximum of the update time, number of recolorings, and query time is $\Omega(\log n)$. We present a deterministic fully dynamic $2$-coloring algorithm with $O(\log^2 n)$ amortized update time, $O(\log n)$ amortized query time, and one recoloring in the worst case. For any incremental $2$-coloring algorithm which explicitly maintains the color of every vertex after each update, the amortized update time and the amortized number of recolorings is $\Omega(\log n)$. For such an algorithm, in the worst case the update time and the number of recolorings is $\Omega(n)$. We then design a deterministic incremental $2$-coloring algorithm which explicitly maintains the color of every vertex after each update, with amortized $O(\log n)$ update time and amortized $O(\log n)$ many recolorings. Further, in the worst case the update time and the number of recolorings is $O(n)$. Further, we present a deterministic incremental $(1+2 \log n)$-coloring algorithm which explicitly maintains the color of every vertex after each update, with $O(\alpha(n))$ amortized update time, at most one recoloring and $O(1)$ query time. We then show a deterministic incremental $2$-coloring algorithm which does not maintain color of every vertex after each update, with amortized $O(\alpha(n))$ update time, amortized $O(\alpha(n))$ recolorings, and amortized $O(\alpha(n))$ query time. For general graphs and graphs of bounded arboricity $\gamma$ and maximum degree $\Delta$ we present a deterministic $(\Delta+1)$-coloring algorithm with $O(\sqrt{m})$ update time, $O(1)$ query time, and one recoloring. Finally, we show a deterministic $(\Delta+1)$-coloring algorithm with amortized $O(\gamma + \log{n})$ update time, $O(1)$ query time, and one recoloring.


title: Stability and Metastability of Skyrmions in Thin Lamellae of Cu$_2$OSeO$_3$

abstract: We report small angle X-ray scattering (SAXS) measurements of the skyrmion lattice in two 200~nm thick Cu$_2$OSeO$_3$ lamellae aligned with the applied magnetic field parallel to the out of plane [110] or [100] crystallographic directions. Our measurements show that the equilibrium skyrmion phase in both samples is expanded significantly compared to bulk crystals, existing between approximately 30 and 50~K over a wide region of magnetic field. This skyrmion state is elliptically distorted at low fields for the [110] sample, and symmetric for the [100] sample, possibly due to crystalline anisotropy becoming more important at this sample thickness than it is in bulk samples. Furthermore, we find that a metastable skyrmion state can be observed at low temperature by field cooling through the equilibrium skyrmion pocket in both samples. In contrast to the behavior in bulk samples, the volume fraction of metastable skyrmions does not significantly depend on cooling rate. We show that a possible explanation for this is the change in the lowest temperature of the skyrmion state in this lamellae compared to bulk, without requiring different energetics of the skyrmion state.


title: An existence theory for nonlinear equations on metric graphs via energy methods

abstract: The purpose of this paper is to develop a general existence theory for constrained minimization problems for functionals defined on function spaces on metric measure spaces $(\mathcal M, d, \mu)$. We apply this theory to functionals defined on metric graphs $\mathcal G$, in particular $L^2$-constrained minimization problems for functionals of the form $$E(u) = \frac{1}{2} a(u,u) - \frac{1}{q}\int_{\mathcal K} |u|^q \, \mathrm dx,$$ where $q>2$, $a(\cdot, \cdot)$ is a suitable symmetric sesquilinear form on some function space on $\mathcal G$ and $\mathcal K \subseteq \mathcal G$ is given. We show how the existence of solutions can be obtained via decomposition methods using spectral properties of the operator $A$ associated with the form $a(\cdot, \cdot)$ and discuss the spectral quantities involved. An example that we consider is the higher-order variant of the stationary NLS (nonlinear Schr\"odinger) energy functional with potential $V\in L^2+ L^\infty(\mathcal G)$ $$ E^{(k)}(u)= \frac{1}{2} \int_{\mathcal G} |u^{(k)}|^2+ V(x) |u|^2 \, \mathrm dx - \frac{1}{p} \int_{\mathcal K} |u|^q \, \mathrm dx$$ defined on a class of higher-order Sobolev spaces $H^k(\mathcal G)$ that we introduce. When $\mathcal K$ is a bounded subgraph, one has localized nonlinearities, which we treat as a special case. When $k=1$ we also consider metric graphs with infinite edge set as well as magnetic potentials. Then the operator $A$ associated to the linear form is a Schr\"odinger operator, and in the $L^2$-subcritical case $2<q<6$, we obtain generalizations of existence results for the NLS functional as for instance obtained by Adami, Serra and Tilli [JFA 271 (2016), 201-223], and Cacciapuoti, Finco and Noja [Nonlinearity 30 (2017), 3271-3303], among others.


title: A Control Theorem for Primitive ideals in Iwasawa algebras

abstract: Let p be a prime, K a p-adic field, G a nilpotent, uniform pro-p group. We prove that all faithful, primitive ideals in the Iwasawa algebra KG are controlled by the centraliser of the second term in the upper central series for G.


title: Learning to Search for MIMO Detection

abstract: This paper proposes a novel learning to learn method, called learning to learn iterative search algorithm (LISA), for signal detection in a multi-input multi-output (MIMO) system. The idea is to regard the signal detection problem as a decision making problem over tree. The goal is to learn the optimal decision policy. In LISA, deep neural networks are used as parameterized policy function. Through training, optimal parameters of the neural networks are learned and thus optimal policy can be approximated. Different neural network based architectures are used for fixed and varying channel models, respectively. LISA provides soft decisions and does not require any information about the additive white Gaussian noise. Simulation results show that LISA 1) obtains near maximum likelihood detection performance in both fixed and varying channel models under QPSK modulation; 2) achieves significantly better bit error rate (BER) performance than classical detectors and recently proposed deep/machine learning based detectors at various modulations and signal to noise (SNR) ratios both under i.i.d and correlated Rayleigh fading channels in the simulation experiments; 3) is robust to MIMO detection problems with imperfect channel state information; and 4) generalizes very well against channel correlation and SNRs.


title: Steady-State Optimal Frequency Control for Lossy Power Grids with Distributed Communication

abstract: We present a distributed and price-based control approach for frequency regulation in power grids with nonzero line conductances. Both grid and controller are modeled as a port-Hamiltonian system, where the grid model consists of differential as well as algebraic equations. Simulations show that the resulting controller asymptotically stabilizes the frequency while maintaining minimum overall generation costs in steady state and being robust in terms of clock drifts and uncontrollable loads. Moreover, it is shown that active power sharing can be achieved by an appropriate choice of the cost function.


title: Theoretical analysis of quantum random walks with stress-engineered optics

abstract: Quantum random walks (QRWs) are random processes in which the resulting probability density of the "walker" state, whose movement is governed by a "coin" state, is described in a non-classical manner. Previously, Q-plates have been used to demonstrate QRWs with polarization and orbital angular momentum playing the role of coin and walker states, respectively. In this theoretical analysis, we show how stress-engineered optics can be used to develop new platforms for complex QRWs through relative simple optical elements. Our work opens up new paths to speed up classical-to-quantum transitions in robust photonic networks.


title: Markov Processes with Jumps on Manifolds and Lie Groups

abstract: We review some developments concerning Markov and Feller processes with jumps in geometric settings. These include stochastic differential equations in Markus canonical form, the Courr\`{e}ge theorem on Lie groups, and invariant Markov processes on manifolds under both transitive and more general Lie group actions.


title: Minimax Confidence Intervals for the Sliced Wasserstein Distance

abstract: Motivated by the growing popularity of variants of the Wasserstein distance in statistics and machine learning, we study statistical inference for the Sliced Wasserstein distance--an easily computable variant of the Wasserstein distance. Specifically, we construct confidence intervals for the Sliced Wasserstein distance which have finite-sample validity under no assumptions or under mild moment assumptions. These intervals are adaptive in length to the regularity of the underlying distributions. We also bound the minimax risk of estimating the Sliced Wasserstein distance, and as a consequence establish that the lengths of our proposed confidence intervals are minimax optimal over appropriate distribution classes. To motivate the choice of these classes, we also study minimax rates of estimating a distribution under the Sliced Wasserstein distance. These theoretical findings are complemented with a simulation study demonstrating the deficiencies of the classical bootstrap, and the advantages of our proposed methods. We also show strong correspondences between our theoretical predictions and the adaptivity of our confidence interval lengths in simulations. We conclude by demonstrating the use of our confidence intervals in the setting of simulator-based likelihood-free inference. In this setting, contrasting popular approximate Bayesian computation methods, we develop uncertainty quantification methods with rigorous frequentist coverage guarantees.


title: Character-Centric Storytelling

abstract: Sequential vision-to-language or visual storytelling has recently been one of the areas of focus in computer vision and language modeling domains. Though existing models generate narratives that read subjectively well, there could be cases when these models miss out on generating stories that account and address all prospective human and animal characters in the image sequences. Considering this scenario, we propose a model that implicitly learns relationships between provided characters and thereby generates stories with respective characters in scope. We use the VIST dataset for this purpose and report numerous statistics on the dataset. Eventually, we describe the model, explain the experiment and discuss our current status and future work.


title: Time scale design for network resilience

abstract: In this paper we consider the $\mathcal{H}_2$-norm of networked systems with multi-time scale consensus dynamics. We develop a general framework for such systems that allows for edge weighting, independent agent-based time scales, as well as measurement and process noise. From this general system description, we highlight an interesting case where the influences of the weighting and scaling can be separated in the design problem. We then consider the design of the time scale parameters for minimizing the $\mathcal{H}_2$-norm for the purpose of network resilience.


title: Mitigating Network Noise on Dragonfly Networks through Application-Aware Routing

abstract: System noise can negatively impact the performance of HPC systems, and the interconnection network is one of the main factors contributing to this problem. To mitigate this effect, adaptive routing sends packets on non-minimal paths if they are less congested. However, while this may mitigate interference caused by congestion, it also generates more traffic since packets traverse additional hops, causing in turn congestion on other applications and on the application itself. In this paper, we first describe how to estimate network noise. By following these guidelines, we show how noise can be reduced by using routing algorithms which select minimal paths with a higher probability. We exploit this knowledge to design an algorithm which changes the probability of selecting minimal paths according to the application characteristics. We validate our solution on microbenchmarks and real-world applications on two systems relying on a Dragonfly interconnection network, showing noise reduction and performance improvement.


title: Walling up Backdoors in Intrusion Detection Systems

abstract: Interest in poisoning attacks and backdoors recently resurfaced for Deep Learning (DL) applications. Several successful defense mechanisms have been recently proposed for Convolutional Neural Networks (CNNs), for example in the context of autonomous driving. We show that visualization approaches can aid in identifying a backdoor independent of the used classifier. Surprisingly, we find that common defense mechanisms fail utterly to remove backdoors in DL for Intrusion Detection Systems (IDSs). Finally, we devise pruning-based approaches to remove backdoors for Decision Trees (DTs) and Random Forests (RFs) and demonstrate their effectiveness for two different network security datasets.


title: TruPercept: Trust Modelling for Autonomous Vehicle Cooperative Perception from Synthetic Data

abstract: Inter-vehicle communication for autonomous vehicles (AVs) stands to provide significant benefits in terms of perception robustness. We propose a novel approach for AVs to communicate perceptual observations, tempered by trust modelling of peers providing reports. Based on the accuracy of reported object detections as verified locally, communicated messages can be fused to augment perception performance beyond line of sight and at great distance from the ego vehicle. Also presented is a new synthetic dataset which can be used to test cooperative perception. The TruPercept dataset includes unreliable and malicious behaviour scenarios to experiment with some challenges cooperative perception introduces. The TruPercept runtime and evaluation framework allows modular component replacement to facilitate ablation studies as well as the creation of new trust scenarios we are able to show.


title: Simple Model for Estimation of the Influence of Velocity on Advancing Dynamic Contact Angles

abstract: Given the importance of dynamic contact angle, its numerous applications and the complexity and difficulty of use of available approaches, here we present a new simple semi-empirical models for estimation of dynamic contact angle's dependence on the wetting line velocity. These models should be applicable to any geometry, a very large range of capillary numbers and static contact angles and all solid-liquid-gas systems without requiring further experiments. Two simple equations are intuitively derived from the most promising theoretical dynamic contact angle models, the hydrodynamic and the molecular-kinetic models. Then the models, along with the basic form of the Hoffman model, are fitted to a large pool of data. The data are extracted from numerous studies and cover over 5 decades of capillary number, include static contact angles up to the superhydrophobic region and comprise of various geometries. The resulting models are compared to each other where the hydrodynamic model's predictions are found to be superior to the other two models by all statistical measures. Then, noting that the molecular forces become dominant at lower capillary numbers we separate the data into low and high capillary regions and repeat the process. There was only a minuscule difference between the results obtained for general models and high capillary models ($\textrm{Ca}>10^{-4}$) but the empirical approach resulted in the most accurate model at low capillary numbers ($\textrm{Ca}<10^{-4}$).


title: Visualizing Movement Control Optimization Landscapes

abstract: A large body of animation research focuses on optimization of movement control, either as action sequences or policy parameters. However, as closed-form expressions of the objective functions are often not available, our understanding of the optimization problems is limited. Building on recent work on analyzing neural network training, we contribute novel visualizations of high-dimensional control optimization landscapes; this yields insights into why control optimization is hard and why common practices like early termination and spline-based action parameterizations make optimization easier. For example, our experiments show how trajectory optimization can become increasingly ill-conditioned with longer trajectories, but parameterizing control as partial target states---e.g., target angles converted to torques using a PD-controller---can act as an efficient preconditioner. Both our visualizations and quantitative empirical data also indicate that neural network policy optimization scales better than trajectory optimization for long planning horizons. Our work advances the understanding of movement optimization and our visualizations should also provide value in educational use.


title: The Feigin-Odesskii Wheel Conditions and Sheaves on Surfaces

abstract: In this paper, we start the study of Feigin-Odesskii wheel conditions from a geometric viewpoint, and generalize it to the K-theory Hall algebra of any surface.


title: Field line winding of braided vector fields in tubular subdomains

abstract: Braided vector fields on spatial subdomains homeomorphic to the cylinder play a crucial role in applications such as solar and plasma physics, relativistic astrophysics, fluid and vortex dynamics, elasticity, and bio-elasticity. Often the vector field's topology -- the entanglement of its field lines -- is non-trivial, and can play a significant role in the vector field's evolution. We present a complete topological characterisation of such vector fields (up to isotopy) using a quantity called field line winding. This measures the entanglement of each field line with all other field lines of the vector field, and may be defined for an arbitrary tubular subdomain by prescribing a minimally distorted coordinate system. We propose how to define such coordinates, and prove that the resulting field line winding distribution uniquely classifies the topology of a braided vector field. The field line winding is similar to the field line helicity considered previously for magnetic (solenoidal) fields, but is a more fundamental measure of the field line topology because it does not conflate linking information with field strength.


title: sktime: A Unified Interface for Machine Learning with Time Series

abstract: We present sktime -- a new scikit-learn compatible Python library with a unified interface for machine learning with time series. Time series data gives rise to various distinct but closely related learning tasks, such as forecasting and time series classification, many of which can be solved by reducing them to related simpler tasks. We discuss the main rationale for creating a unified interface, including reduction, as well as the design of sktime's core API, supported by a clear overview of common time series tasks and reduction approaches.


title: Generating Black-Box Adversarial Examples for Text Classifiers Using a Deep Reinforced Model

abstract: Recently, generating adversarial examples has become an important means of measuring robustness of a deep learning model. Adversarial examples help us identify the susceptibilities of the model and further counter those vulnerabilities by applying adversarial training techniques. In natural language domain, small perturbations in the form of misspellings or paraphrases can drastically change the semantics of the text. We propose a reinforcement learning based approach towards generating adversarial examples in black-box settings. We demonstrate that our method is able to fool well-trained models for (a) IMDB sentiment classification task and (b) AG's news corpus news categorization task with significantly high success rates. We find that the adversarial examples generated are semantics-preserving perturbations to the original text.


title: Efficiency limits of electronically-coupled upconverter and quantum ratchet solar cells using detailed balance

abstract: The intermediate band solar cell (IBSC) and quantum ratchet solar cell (QRSC) have the potential to surpass the efficiency of standard single-junction solar cells by allowing sub-gap photon absorption through states deep inside the band gap. High efficiency IBSC and QRSC devices have not yet been achieved, however, since introducing mid-gap states also increases recombination, which can harm the device. We consider the electronically coupled upconverter (ECUC) solar cell and show that it can achieve the same efficiencies as the QRSC. Although they are equivalent in the detailed balance limit, the ECUC is less sensitive to nonradiative processes, which makes it a more practical implementation for IB devices. We perform a case study of crystalline-silicon based ECUC cells, focusing on hydrogenated amorphous silicon as the upconverter material and highlighting potential dopants for the ECUC. These results illustrate a new path for the development of IB-based devices.


title: Simulation of nonlinear Compton scattering from bound electrons

abstract: Recent investigations by Fuchs et al. [Nat. Phys. 11, 964 (2015)] revealed an anomalous frequency shift in non-linear Compton scattering of high-intensity X-rays by electrons in solid beryllium. This frequency shift was at least 800 eV to the red of the values predicted by analytical free-electron models for the same process. In this paper, we describe a method for simulating non-linear Compton scattering. The method is applied to the case of bound electrons in a local, spherical potential to explore the role of binding energy in the frequency shift of scattered X-rays for different scattered angles. The results of the calculation do not exhibit an additional redshift for the scattered X-rays beyond the non-linear Compton shift predicted by the free-electron model. However, it does reveal a small blue-shift relative to the free electron prediction for non-linear Compton scattering. The effect of electron-electron correlation effects is calculated and determined to be unlikely to be the source of the redshift. The case of linear Compton scattering from a photoionized electron followed by electron recapture is examined as a possible source of the redshift and ruled out.


title: Learning to Manipulate Object Collections Using Grounded State Representations

abstract: We propose a method for sim-to-real robot learning which exploits simulator state information in a way that scales to many objects. We first train a pair of encoder networks to capture multi-object state information in a latent space. One of these encoders is a CNN, which enables our system to operate on RGB images in the real world; the other is a graph neural network (GNN) state encoder, which directly consumes a set of raw object poses and enables more accurate reward calculation and value estimation. Once trained, we use these encoders in a reinforcement learning algorithm to train image-based policies that can manipulate many objects. We evaluate our method on the task of pushing a collection of objects to desired tabletop regions. Compared to methods which rely only on images or use fixed-length state encodings, our method achieves higher success rates, performs well in the real world without fine tuning, and generalizes to different numbers and types of objects not seen during training.


title: Multi-mapping Image-to-Image Translation via Learning Disentanglement

abstract: Recent advances of image-to-image translation focus on learning the one-to-many mapping from two aspects: multi-modal translation and multi-domain translation. However, the existing methods only consider one of the two perspectives, which makes them unable to solve each other's problem. To address this issue, we propose a novel unified model, which bridges these two objectives. First, we disentangle the input images into the latent representations by an encoder-decoder architecture with a conditional adversarial training in the feature space. Then, we encourage the generator to learn multi-mappings by a random cross-domain translation. As a result, we can manipulate different parts of the latent representations to perform multi-modal and multi-domain translations simultaneously. Experiments demonstrate that our method outperforms state-of-the-art methods.


title: Fractional matching preclusion number of graphs

abstract: The \emph{fractional matching preclusion number} of a graph $G$, denoted by $fmp(G)$, is the minimum number of edges whose deletion results in a graph that has no fractional perfect matchings. In this paper, we first give some sharp upper and lower bounds of fractional matching preclusion number. Next, graphs with large and small fractional matching preclusion number are characterized, respectively. In the end, we investigate some extremal problems on fractional matching preclusion number.


title: Another approach to track reconstruction: cluster analysis

abstract: A novel combination of data analysis techniques is proposed for the reconstruction of all tracks of primary charged particles, as well as of daughters of displaced vertices (decays, photon conversions, nuclear interactions), created in high energy collisions. Instead of performing a classical trajectory building or an image transformation, an efficient use of both local and global information is undertaken while keeping competing choices open. The measured hits of adjacent tracking layers are clustered first with the help of a mutual nearest neighbor search in the angular distance. The resulted chains of connected hits are used as initial clusters and as input for cluster analysis algorithms, such as the robust $k$-medians clustering. This latter proceeds by alternating between the hit-to-track assignment and the track-fit update steps, until convergence. The calculation of the hit-to-track distance and that of the track-fit $\chi^2$ is performed through the global covariance of the measured hits. The clustering is complemented with elements from a more sophisticated Metropolis-Hastings MCMC algorithm, with the possibility of adding new track hypotheses or removing unnecessary ones. Simplified but realistic models of today's silicon trackers, including the relevant physics processes, are employed to test and study the performance (efficiency, purity) of the proposed method as a function of the particle multiplicity in the collision event.


title: Katugampola Fractional Calculus With Generalized $k-$Wright Function

abstract: In this article, we presented some properties of the Katugampola fractional integrals and derivatives. Also we studied the fractional calculus properties involving Katugampola Fractional integrals and derivatives of generalized $k-$Wright function $_{n}\Phi_{m}^{k}(z).$\\[2mm]


title: Estimating Glycemic Impact of Cooking Recipes via Online Crowdsourcing and Machine Learning

abstract: Consumption of diets with low glycemic impact is highly recommended for diabetics and pre-diabetics as it helps maintain their blood glucose levels. However, laboratory analysis of dietary glycemic potency is time-consuming and expensive. In this paper, we explore a data-driven approach utilizing online crowdsourcing and machine learning to estimate the glycemic impact of cooking recipes. We show that a commonly used healthiness metric may not always be effective in determining recipes suitable for diabetics, thus emphasizing the importance of the glycemic-impact estimation task. Our best classification model, trained on nutritional and crowdsourced data obtained from Amazon Mechanical Turk (AMT), can accurately identify recipes which are unhealthful for diabetics.


title: A Note on Explicit Milstein-Type Scheme for Stochastic Differential Equation with Markovian Switching

abstract: An explicit Milstein-type scheme for stochastic differential equation with Markovian switching is derived and its strong convergence in $\mathcal{L}^2$-sense is established without using It\^o-Taylor expansion formula. Rate of strong convergence is shown to be equal to $1.0$ under the assumptions that coefficients satisfy mild regularity conditions. More precisely, coefficients are assumed to be only once differentiable which are more relaxed conditions than those made in existing literature.


title: Uncertainty Quantification of Tunable Elastic Metamaterials using Polynomial Chaos

abstract: Owing to their periodic and intricate configurations, metamaterials engineered for acoustic and elastic wave control inevitably suffer from manufacturing anomalies and deviate from theoretical dispersion predictions. This work exploits the Polynomial Chaos theory to quantify the magnitude and extent of these deviations and assess their impact on the desired behavior. It is shown that uncertainties stemming from surface roughness, tolerances, and other inconsistencies in a metamaterial's unit cell parameters alter the targetted band gap width, location and the confidence level with which it is guaranteed. The effect of uncertainties are projected from a Bloch-wave dispersion analysis of three distinct phononic and resonant cellular configurations and are further confirmed in the frequency response the finite structures. The analysis concludes with a unique algorithm intended to guide the design of metamaterials in the presence of system uncertainties.


title: Enhanced hydrodynamic transport in near magic angle twisted bilayer graphene

abstract: Using the semiclassical quantum Boltzmann theory and employing the Dirac model with twist angle-dependent Fermi velocity we obtain results for the electrical resistivity, the electronic thermal resistivity, the Seebeck coefficient, and the Wiedemann-Franz ratio in near magic angle twisted bilayer graphene, as functions of doping density (around the charge-neutrality-point) and modified Fermi velocity $\tilde v$. The $\tilde v$-dependence of the relevant scattering mechanisms, i.e. electron-hole Coulomb, long-ranged impurities, and acoustic gauge phonons, is considered in detail. We find a range of twist angles and temperatures, where the combined effect of momentum-non-conserving collisions (long-ranged impurities and phonons) is minimal, opening a window for the observation of strong hydrodynamic transport. Several experimental signatures are identified, such as a sharp dependence of the electric resistivity on doping density and a large enhancement of the Wiedemann-Franz ratio and the Seebeck coefficient.


title: Microscopic study of strange hadron productions at the LHC with$\sqrt{s_{NN}}$=2.76 TeV

abstract: Ratio of the yield of strange hadrons to pions is considered as an important observable in studying the properties of the system produced in relativistic heavy ion collisions. Production of strange hadrons $K, \bar{K}, \Lambda, \Sigma, \Xi$ and $\Omega$ have been evaluated microscopically using rate equations by considering their hadronic interaction cross sections in an expanding medium. The yields obtained from rate equations are normalized with thermal pions and compared with the measurements from Pb-Pb collisions at various multiplicities at LHC energy. The calculation has been done for various initial and freeze out conditions. At 2760 GeV, LHC energy, $\Lambda, \Xi, \Omega$ freeze out close to $T_C$ and $K^0_s$ freezes out little later. But there is a subtle difference in freeze out temperatures of different species which may be distinguishable at lower colliding energies.


title: On Explicit Tamed Milstein-type scheme for Stochastic Differential Equation with Markovian Switching

abstract: We propose a new tamed Milstein-type scheme for stochastic differential equation with Markovian switching when drift coefficient is assumed to grow super-linearly. The strong rate of convergence is shown to be equal to $1.0$ under mild regularity (e.g. once differentiability) requirements on drift and diffusion coefficients. Novel techniques are developed to tackle two-fold difficulties arising due to jumps of the Markov chain and the reduction of regularity requirements on the coefficients.


title: Inferring and Learning Multi-Robot Policies by Observing an Expert

abstract: We present a technique for learning how to solve a multi-robot mission that requires interaction with an external environment by observing an expert system executing the same mission. We define the expert system as a team of robots equipped with a library of controllers, each designed to solve a specific task, supervised by an expert policy that appropriately selects controllers based on the states of robots and environment. The objective is for an un-trained team of robots (i.e., imitator system) equipped with the same library of controllers, but agnostic to the expert policy, to execute the mission, with performances comparable to those of the expert system. From un-annotated observations of the expert system, a multi-hypothesis filtering technique is used to estimate individual controllers executed by the expert policy. Then, the history of estimated controllers and environmental states is used to train a neural network policy for the imitator system. Considering a perimeter protection scenario on a team of differential-drive robots, we show that the learned policy endows the imitator system with performances comparable to those of the expert system.


title: The most massive, passive, and oldest galaxies at 0.5 < z < 2.1: Downsizing signature from galaxies selected from MgUV index

abstract: Aims. We seek is to identify old and massive galaxies at 0.5<z<2.1 on the basis of the magnesium index MgUV and then study their physical properties. We computed the MgUV index based on the best spectral fitting template of $\sim$3700 galaxies using data from the VLT VIMOS Deep Survey (VVDS) and VIMOS Ultra Deep Survey (VUDS) galaxy redshift surveys. Based on galaxies with the largest signal to noise and the best fit spectra we selected 103 objects with the highest spectral MgUV signature. We performed an independent fit of the photometric data of these galaxies and computed their stellar masses, star formation rates, extinction by dust and age, and we related these quantities to the MgUV index. We find that the MgUV index is a suitable tracer of early-type galaxies at an advanced stage of evolution. Selecting galaxies with the highest MgUV index allows us to choose the most massive, passive, and oldest galaxies at any epoch. The formation epoch t_f computed from the fitted age as a function of the total mass in stars supports the downsizing formation paradigm in which galaxies with the highest mass formed most of their stars at an earlier epoch.


title: Distributional conformal prediction

abstract: We propose a robust method for constructing conditionally valid prediction intervals based on models for conditional distributions such as quantile and distribution regression. Our approach can be applied to important prediction problems including cross-sectional prediction, k-step-ahead forecasts, synthetic controls and counterfactual prediction, and individual treatment effects prediction. Our method exploits the probability integral transform and relies on permuting estimated ranks. Unlike regression residuals, ranks are independent of the predictors, allowing us to construct conditionally valid prediction intervals under heteroskedasticity. We establish approximate conditional validity under consistent estimation and provide approximate unconditional validity under model misspecification, overfitting, and with time series data. We also propose a simple "shape" adjustment of our baseline method that yields optimal prediction intervals.


title: More trouble for unitary quantum theory

abstract: Single-world unitary quantum theories imply that some measurements have results whose probabilities can not be calculated by the Born rule.


title: Machine Learning based IoT Edge Node Security Attack and Countermeasures

abstract: Advances in technology have enabled tremendous progress in the development of a highly connected ecosystem of ubiquitous computing devices collectively called the Internet of Things (IoT). Ensuring the security of IoT devices is a high priority due to the sensitive nature of the collected data. Physically Unclonable Functions (PUFs) have emerged as critical hardware primitive for ensuring the security of IoT nodes. Malicious modeling of PUF architectures has proven to be difficult due to the inherently stochastic nature of PUF architectures. Extant approaches to malicious PUF modeling assume that a priori knowledge and physical access to the PUF architecture is available for malicious attack on the IoT node. However, many IoT networks make the underlying assumption that the PUF architecture is sufficiently tamper-proof, both physically and mathematically. In this work, we show that knowledge of the underlying PUF structure is not necessary to clone a PUF. We present a novel non-invasive, architecture independent, machine learning attack for strong PUF designs with a cloning accuracy of 93.5% and improvements of up to 48.31% over an alternative, two-stage brute force attack model. We also propose a machine-learning based countermeasure, discriminator, which can distinguish cloned PUF devices and authentic PUFs with an average accuracy of 96.01%. The proposed discriminator can be used for rapidly authenticating millions of IoT nodes remotely from the cloud server.


title: Infinitesimal symmetries in Contact Hamiltonian systems

abstract: In this paper, we extend the well-known Noether theorem for Lagrangian systems to contact Lagrangian systems. We introduce a classification of infinitesimal symmetries and obtain the corresponding dissipated quantities. We notice that in contact dynamics, the existence of infinitesimal symmetries does not produce conserved quantities, but functions that dissipate at the same rate than the energy; so, the corresponding quotients are true conserved quantities.


title: Vehicle routing by learning from historical solutions

abstract: The goal of this paper is to investigate a decision support system for vehicle routing, where the routing engine learns from the subjective decisions that human planners have made in the past, rather than optimizing a distance-based objective criterion. This is an alternative to the practice of formulating a custom VRP for every company with its own routing requirements. Instead, we assume the presence of past vehicle routing solutions over similar sets of customers, and learn to make similar choices. The approach is based on the concept of learning a first-order Markov model, which corresponds to a probabilistic transition matrix, rather than a deterministic distance matrix. This nevertheless allows us to use existing arc routing VRP software in creating the actual route plans. For the learning, we explore different schemes to construct the probabilistic transition matrix. Our results on a use-case with a small transportation company show that our method is able to generate results that are close to the manually created solutions, without needing to characterize all constraints and sub-objectives explicitly. Even in the case of changes in the client sets, our method is able to find solutions that are closer to the actual route plans than when using distances, and hence, solutions that would require fewer manual changes to transform into the actual route plan.


title: Composite Dynamics in the Early Universe

abstract: We study the occurrence of a strong first-order electroweak phase transition in composite Higgs models. Minimal constructions realising this scenario are based on the coset SO(6)/SO(5) which delivers an extended Higgs sector with an additional scalar. In such models, a two-step phase transition can be obtained with the scalar singlet acquiring a vacuum expectation value at intermediate temperatures. A bonus of the Nambu-Goldstone boson nature of the scalar-sector dynamics is the presence of non-renormalisable Higgs interactions that can trigger additional sources of CP violation needed to realise baryogenesis at the electroweak scale. Another interesting aspect of this scenario is the generation of gravitational wave signatures that can be observed at future space-based interferometers.


title: On the Optimality of the Greedy Policy for Battery Limited Energy Harvesting Communications

abstract: Consider a battery limited energy harvesting communication system with online power control. Assuming independent and identically distributed (i.i.d.) energy arrivals and the harvest-store-use architecture, it is shown that the greedy policy achieves the maximum throughput if and only if the battery capacity is below a certain positive threshold that admits a precise characterization. Simple lower and upper bounds on this threshold are established. The asymptotic relationship between the threshold and the mean of the energy arrival process is analyzed for several examples.


title: No-Arbitrage Commodity Option Pricing with Market Manipulation

abstract: We design three continuous--time models in finite horizon of a commodity price, whose dynamics can be affected by the actions of a representative risk--neutral producer and a representative risk--neutral trader. Depending on the model, the producer can control the drift and/or the volatility of the price whereas the trader can at most affect the volatility. The producer can affect the volatility in two ways: either by randomizing her production rate or, as the trader, using other means such as spreading false information. Moreover, the producer contracts at time zero a fixed position in a European convex derivative with the trader. The trader can be price-taker, as in the first two models, or she can also affect the volatility of the commodity price, as in the third model. We solve all three models semi--explicitly and give closed--form expressions of the derivative price over a small time horizon, preventing arbitrage opportunities to arise. We find that when the trader is price-taker, the producer can always compensate the loss in expected production profit generated by an increase of volatility by a gain in the derivative position by driving the price at maturity to a suitable level. Finally, in case the trader is active, the model takes the form of a nonzero-sum linear-quadratic stochastic differential game and we find that when the production rate is already at its optimal stationary level, there is an amount of derivative position that makes both players better off when entering the game.


title: Diffusion-controlled formation and collapse of a d-dimensional A-particle island in the B-particle sea

abstract: We consider diffusion-controlled evolution of a $d$-dimensional $A$-particle island in the $B$-particle sea at propagation of the sharp reaction front $A+B\to 0$ at equal species diffusivities. The $A$-particle island is formed by a localized (point)$A$-source with a strength $\lambda$ that acts for a finite time $T$. We reveal the conditions under which the island collapse time $t_{c}$ becomes much longer than the injection period $T$ (long-living island) and demonstrate that regardless of $d$ the evolution of the long-living island radius $r_{f}(t)$ is described by the universal law $\zeta_{f}=r_{f}/r_{f}^{M}=\sqrt{e\tau|\ln\tau|}$ where $\tau=t/t_{c}$ and $r_{f}^{M}$ is the maximal island expansion radius at the front turning point $t_{M}=t_{c}/e$. We find that in the long-living island regime the ratio $t_{c}/T$ changes with the increase of the injection period $T$ by the law $\propto (\lambda^{2}T^{2-d})^{1/d}$ i.e. increases with the increase of $T$ in the one-dimensional (1D) case, does not change with the increase of $T$ in the 2D case and decreases with the increase of $T$ in the 3D case. We derive the scaling laws for particles death in the long-living island and determine the limits of their applicability. We demonstrate also that these laws describe asymptotically the evolution of the $d$-dimensional spherical island with a uniform initial particle distribution generalizing the results obtained earlier for the quasi-one-dimensional geometry. As striking results we present a systematic analysis of the front relative width evolution for fluctuation, logarithmically modified and mean-field regimes and demonstrate that in a wide range of parameters the front remains sharp up to a narrow vicinity of the collapse point.


title: An approach for security evaluation and certification of a complete quantum communication system

abstract: Although quantum communication systems are being deployed on a global scale, their realistic security certification is not yet available. Here we present a security evaluation and improvement protocol for complete quantum communication systems. The protocol subdivides a system by defining seven system implementation sub-layers based on a hierarchical order of information flow; then it categorises the known system implementation imperfections by hardness of protection and practical risk. Next, an initial analysis report lists all potential loopholes in its quantum-optical part. It is followed by interactions with the system manufacturer, testing and patching most loopholes, and re-assessing their status. Our protocol has been applied on multiple commercial quantum key distribution systems to improve their security. A detailed description of our methodology is presented with the example of a subcarrier-wave system. Our protocol is a step towards future security evaluation and security certification standards.


title: Fast Search with Poor OCR

abstract: The indexing and searching of historical documents have garnered attention in recent years due to massive digitization efforts of important collections worldwide. Pure textual search in these corpora is a problem since optical character recognition (OCR) is infamous for performing poorly on such historical material, which often suffer from poor preservation. We propose a novel text-based method for searching through noisy text. Our system represents words as vectors, projects queries and candidates obtained from the OCR into a common space, and ranks the candidates using a metric suited to nearest-neighbor search. We demonstrate the practicality of our method on typewritten German documents from the WWII era.


title: A Near-Ideal Molecule-Based Haldane Spin-Chain

abstract: The molecular coordination complex NiI$_2$(3,5-lut)$_4$ [where (3,5-lut) $=$ (3,5-lutidine) $=$ (C$_7$H$_9$N)] has been synthesized and characterized by several techniques including synchrotron X-ray diffraction, ESR, SQUID magnetometry, pulsed-field magnetization, inelastic neutron scattering and muon spin relaxation. Templated by the configuration of 3,5-lut ligands the molecules pack in-registry with the Ni--I$\cdots$I--Ni chains aligned along the $c$--axis. This arrangement leads to through-space I$\cdots$I magnetic coupling which is directly measured for the first time in this work. The net result is a near-ideal realization of the $S = 1$ Haldane chain with $J = 17.5~\rm{K}$ and energy gaps of $\Delta^{\parallel} = 5.3~{\rm K}$ $\Delta^{\perp} =7.7~{\rm K}$, split by the easy-axis single-ion anisotropy $D=-1.2~{\rm K}$. The ratio $D/J = -0.07$ affords one of the most isotropic Haldane systems yet discovered, while the ratio $\Delta_0/J = 0.40(1)$ (where $\Delta_0$ is the average gap size) is close to its ideal theoretical value, suggesting a very high degree of magnetic isolation of the spin chains in this material. The Haldane gap is closed by orientation-dependent critical fields $\mu_0H_{\rm c}^{\parallel} = 5.3~\rm{T}$ and $\mu_0H_{\rm c}^{\perp} = 4.3~\rm{T}$, which are readily accessible experimentally and permit investigations across the entirety of the Haldane phase, with the fully polarized state occurring at $\mu_0 H_{\rm s}^{\parallel}=46.0~\rm{T}$ and $\mu_0 H_{\rm s}^{\perp}=50.7~\rm{T}$. The results are explicable within the so-called fermion model, in contrast to other reported easy-axis Haldane systems. Zero-field magnetic order is absent down to $20~{\rm mK}$ and emergent end-chain effects are observed in the gapped state, as evidenced by detailed low-temperature measurements.


title: Asymptotic variational analysis of incompressible elastic strings

abstract: Starting from three-dimensional nonlinear elasticity under the restriction of incompressibility, we derive reduced models to capture the behavior of strings in response to external forces. Our $\Gamma$-convergence analysis of the constrained energy functionals in the limit of shrinking cross sections gives rise to explicit one-dimensional limit energies. The latter depend on the scaling of the applied forces. The effect of local volume preservation is reflected either in their energy densities through a constrained minimization over the cross-section variables or in the class of admissible deformations. Interestingly, all scaling regimes allow for compression and/or stretching of the string. The main difficulty in the proof of the $\Gamma$-limit is to establish recovery sequences that accommodate the nonlinear differential constraint imposed by the incompressibility. To this end, we modify classical constructions in the unconstrained case with the help of an inner perturbation argument tailored for $3$d-$1$d dimension reduction problems.


title: Top quark measurements by the D0 collaboration

abstract: The 2019 European Physical Society Prize for High Energy and Particle Physics was awarded to the CDF and D0 collaborations for "the discovery of the top quark and the detailed measurement of its properties". This paper is based upon the talk accepting the prize on behalf of the D0 Collaboration.


title: Graph Nets for Partial Charge Prediction

abstract: Atomic partial charges are crucial parameters for Molecular Dynamics (MD) simulations, molecular mechanics calculations, and virtual screening, as they determine the electrostatic contributions to interaction energies. Current methods for calculating partial charges, however, are either slow and scale poorly with molecular size (quantum chemical methods) or unreliable (empirical methods). Here, we present a new charge derivation method based on Graph Nets---a set of update and aggregate functions that operate on molecular topologies and propagate information thereon---that could approximate charges derived from Density Functional Theory (DFT) calculations with high accuracy and an over 500-fold speed up.


title: Strong cosmic censorship for the Dirac field in the higher dimensional Reissner-Nordstrom--de Sitter black hole

abstract: We investigate the strong cosmic censorship for the Dirac field in the higher dimensional Reissner-Norstrom--de Sitter black hole. To achieve this, we first use the conformal transformation trick to massage the Dirac equation to a pair of coupled equations in a meticulously chosen orthonormal basis and derive the criterion on the quasinormal modes for the violation of the strong cosmic censorship, which turns out to be independent of the spacetime dimension. Then we apply the Crank-Nicolson method to evolve our Dirac equation in the double null coordinates and extract the low-lying quasinormal modes from the evolution data by the Prony method. It is shown for the spacetime dimension $D=4, 5, 6$ under consideration that although the strong cosmic censorship is violated by the perturbation from the neutral Dirac field in the near-extremal black hole, the strong cosmic censorship can be restored when the charge of the Dirac field is increased beyond a critical value. The closer to the extremal limit the black hole is, the larger the critical charge of the Dirac field is.


title: Angular measures and Birkhoff orthogonality in Minkowski planes

abstract: Let $x$ and $y$ be two unit vectors in a normed plane $\mathbb{R}^2$. We say that $x$ is Birkhoff orthogonal to $y$ if the line through $x$ in the direction $y$ supports the unit disc. A B-measure (Fankh\"anel 2011) is an angular measure $\mu$ on the unit circle for which $\mu(C)=\pi/2$ whenever $C$ is a shorter arc of the unit circle connecting two Birkhoff orthogonal points. We present a characterization of the normed planes that admit a B-measure.


title: Strain-driven InAs island growth on top of GaAs(111) nanopillars

abstract: We analyze the shape and position of heteroepitaxial InAs islands on the top face of cylindrical GaAs(111)A nanopillars experimentally and theoretically. Catalyst-free molecular beam epitaxial growth of InAs at low temperatures on GaAs nanopillars results in InAs islands with diameters < 30 nm exhibiting predominantly rounded triangular in-plane shapes. The islands show a tendency to grow at positions displaced from the center towards the pillar edge. Atomistic molecular statics simulations evidence that triangular-prismatic islands centered to the pillar axis with diameters smaller than that of the nanopillars are energetically favored. Moreover, we reveal the existence of minimum-energy states for off-axis island positions, in agreement with the experiment. These findings are interpreted by evaluating the spatial strain distributions and the number of broken bonds of surface atoms as a measure for the surface energy. The preferred off-axis island positions can be understood in terms of an increased compliancy of the GaAs nanopillar beneath the island because of the vicinity of free surfaces, leading to a reduction of strain energy. The influence of surface steps on the energy of the system is addressed as well.


title: Pointer-based Fusion of Bilingual Lexicons into Neural Machine Translation

abstract: Neural machine translation (NMT) systems require large amounts of high quality in-domain parallel corpora for training. State-of-the-art NMT systems still face challenges related to out-of-vocabulary words and dealing with low-resource language pairs. In this paper, we propose and compare several models for fusion of bilingual lexicons with an end-to-end trained sequence-to-sequence model for machine translation. The result is a fusion model with two information sources for the decoder: a neural conditional language model and a bilingual lexicon. This fusion model learns how to combine both sources of information in order to produce higher quality translation output. Our experiments show that our proposed models work well in relatively low-resource scenarios, and also effectively reduce the parameter size and training cost for NMT without sacrificing performance.


title: Algorithm for Training Neural Networks on Resistive Device Arrays

abstract: Hardware architectures composed of resistive cross-point device arrays can provide significant power and speed benefits for deep neural network training workloads using stochastic gradient descent (SGD) and backpropagation (BP) algorithm. The training accuracy on this imminent analog hardware however strongly depends on the switching characteristics of the cross-point elements. One of the key requirements is that these resistive devices must change conductance in a symmetrical fashion when subjected to positive or negative pulse stimuli. Here, we present a new training algorithm, so-called the "Tiki-Taka" algorithm, that eliminates this stringent symmetry requirement. We show that device asymmetry introduces an unintentional implicit cost term into the SGD algorithm, whereas in the "Tiki-Taka" algorithm a coupled dynamical system simultaneously minimizes the original objective function of the neural network and the unintentional cost term due to device asymmetry in a self-consistent fashion. We tested the validity of this new algorithm on a range of network architectures such as fully connected, convolutional and LSTM networks. Simulation results on these various networks show that whatever accuracy is achieved using the conventional SGD algorithm with symmetric (ideal) device switching characteristics the same accuracy is also achieved using the "Tiki-Taka" algorithm with non-symmetric (non-ideal) device switching characteristics. Moreover, all the operations performed on the arrays are still parallel and therefore the implementation cost of this new algorithm on array architectures is minimal; and it maintains the aforementioned power and speed benefits. These algorithmic improvements are crucial to relax the material specification and to realize technologically viable resistive crossbar arrays that outperform digital accelerators for similar training tasks.


title: hm-toolbox: Matlab software for HODLR and HSS matrices

abstract: Matrices with hierarchical low-rank structure, including HODLR and HSS matrices, constitute a versatile tool to develop fast algorithms for addressing large-scale problems. While existing software packages for such matrices often focus on linear systems, their scope of applications is in fact much wider and includes, for example, matrix functions and eigenvalue problems. In this work, we present a new Matlab toolbox called hm-toolbox, which encompasses this versatility with a broad set of tools for HODLR and HSS matrices, unmatched by existing software. While mostly based on algorithms that can be found in the literature, our toolbox also contains a few new algorithms as well as novel auxiliary functions. Being entirely based on Matlab, our implementation does not strive for optimal performance. Nevertheless, it maintains the favorable complexity of hierarchical low-rank matrices and offers, at the same time, a convenient way of prototyping and experimenting with algorithms. A number of applications illustrate the use of the hm-toolbox.


title: Nontrivial topology in the continuous spectrum of a magnetized plasma

abstract: Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can coexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the Alfv\'{e}n continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear Alfv\'{e}n eigenmode (RSAE) is well known within magnetic fusion as it has been observed in several tokamaks. In examining how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology. More broadly, this work raises the possibility that fusion devices could be carefully tailored to produce other topological states with potentially useful behavior.

categories: physics.plasm-ph astro-ph.SR cond-mat.mes-hall physics.flu-dyn


title: Design of High-Performance Photon Number Resolving Photodetectors Based on Coherently Interacting Nanoscale Elements

abstract: A number of applications in basic science and technology would benefit from high fidelity photon number resolving photodetectors. While some recent experimental progress has been made in this direction, the requirements for true photon number resolution are stringent, and no design currently exists that achieves this goal. Here we employ techniques from fundamental quantum optics to demonstrate that detectors composed of subwavelength elements interacting collectively with the photon field can achieve high-performance photon number resolution. We propose a new design that simultaneously achieves photon number resolution, high efficiency, low jitter, low dark counts, and high count rate. We discuss specific systems that satisfy the design requirements, pointing to the important role of nanoscale device elements.


title: Representation stability for pure braid group Milnor fibers

abstract: We prove a representation stability result for the Milnor fiber associated to the pure braid group. Our result connects previous work of Simona Settepenella to representation stability in the sense of Church--Ellenberg--Farb, answering a question of Graham Denham. We also use our result to compute the stable integral homology of the Milnor fiber.


title: Learning to Deceive with Attention-Based Explanations

abstract: Attention mechanisms are ubiquitous components in neural architectures applied to natural language processing. In addition to yielding gains in predictive accuracy, attention weights are often claimed to confer interpretability, purportedly useful both for providing insights to practitioners and for explaining why a model makes its decisions to stakeholders. We call the latter use of attention mechanisms into question by demonstrating a simple method for training models to produce deceptive attention masks. Our method diminishes the total weight assigned to designated impermissible tokens, even when the models can be shown to nevertheless rely on these features to drive predictions. Across multiple models and tasks, our approach manipulates attention weights while paying surprisingly little cost in accuracy. Through a human study, we show that our manipulated attention-based explanations deceive people into thinking that predictions from a model biased against gender minorities do not rely on the gender. Consequently, our results cast doubt on attention's reliability as a tool for auditing algorithms in the context of fairness and accountability.


title: Blockchain of Signature Material Combining Cryptographic Hash Function and DNA Steganography

abstract: An ideal signature material and method, which can be used to prove the authenticity of a physical item and against forgery, should be immune to the fast developments in digital and engineering technologies. Herein, the design of signature material combining cryptographic hash function and DNA steganography is proposed. The encrypting materials are used to construct a series of time-stamped records (blockchain) associated with published hash values, while each DNA-encrypted block is associated with a set of DNA keys. The decrypted DNA information, as digital keys, can be validated through a hash function to compare with the published hash values. The blocks can also be cross-referenced among different related signatures. While both digital cryptography and DNA steganography can have large key size, automated brutal force search is far more labor intensive and expensive for DNA steganography with wet lab experiments, as compared to its digital counterpart. Moreover, the time-stamped blockchain structure permits the incorporation of new cryptographic functions and DNA steganographies over time, thus can evolve over time without losing the continuous history line.


title: The Maximum Binary Tree Problem

abstract: We introduce and investigate the approximability of the maximum binary tree problem (MBT) in directed and undirected graphs. The goal in MBT is to find a maximum-sized binary tree in a given graph. MBT is a natural variant of the well-studied longest path problem, since both can be viewed as finding a maximum-sized tree of bounded degree in a given graph. The connection to longest path motivates the study of MBT in directed acyclic graphs (DAGs), since the longest path problem is solvable efficiently in DAGs. In contrast, we show that MBT in DAGs is in fact hard: it has no efficient $\exp(-O(\log n/ \log \log n))$-approximation algorithm under the exponential time hypothesis, where $n$ is the number of vertices in the input graph. In undirected graphs, we show that MBT has no efficient $\exp(-O(\log^{0.63}{n}))$-approximation under the exponential time hypothesis. Our inapproximability results rely on self-improving reductions and structural properties of binary trees. We also show constant-factor inapproximability assuming $\text{P}\neq \text{NP}$. In addition to inapproximability results, we present algorithmic results along two different flavors: (1) We design a randomized algorithm to verify if a given directed graph on $n$ vertices contains a binary tree of size $k$ in $2^k \text{poly}(n)$ time. (2) Motivated by the longest heapable subsequence problem, introduced by Byers, Heeringa, Mitzenmacher, and Zervas (ANALCO 2011), which is equivalent to MBT in permutation DAGs, we design efficient algorithms for MBT in bipartite permutation graphs.


title: Safety-Critical Adaptive Control with Nonlinear Reference Model Systems

abstract: In this paper, a model reference adaptive control architecture is proposed for uncertain nonlinear systems to achieve prescribed performance guarantees. Specifically, a general nonlinear reference model system is considered that captures an ideal and safe system behavior. An adaptive control architecture is then proposed to suppress the effects of system uncertainties without any prior knowledge of their magnitude and rate upper bounds. More importantly, the proposed control architecture enforces the system state trajectories to evolve within a user-specified prescribed distance from the reference system trajectories, satisfying the safety constraints. This eliminates the ad-hoc tuning process for the adaptation rate that is conventionally required in model reference adaptive control to ensure safety. The efficacy of the proposed control architecture is also demonstrated through an illustrative numerical example.


title: Toward Efficient Evaluation of Logic Encryption Schemes: Models and Metrics

abstract: Research in logic encryption over the last decade has resulted in various techniques to prevent different security threats such as Trojan insertion, intellectual property leakage, and reverse engineering. However, there is little agreement on a uniform set of metrics and models to efficiently assess the achieved security level and the trade-offs between security and overhead. This paper addresses the above challenges by relying on a general logic encryption model that can encompass all the existing techniques, and a uniform set of metrics that can capture multiple, possibly conflicting, security concerns. We apply our modeling approach to four state-of-the-art encryption techniques, showing that it enables fast and accurate evaluation of design trade-offs, average prediction errors that are at least 2X smaller than previous approaches, and the evaluation of compound encryption methods.


title: A Programming Framework for Differential Privacy with Accuracy Concentration Bounds

abstract: Differential privacy offers a formal framework for reasoning about privacy and accuracy of computations on private data. It also offers a rich set of building blocks for constructing data analyses. When carefully calibrated, these analyses simultaneously guarantee privacy of the individuals contributing their data, and accuracy of their results for inferring useful properties about the population. The compositional nature of differential privacy has motivated the design and implementation of several programming languages aimed at helping a data analyst in programming differentially private analyses. However, most of the programming languages for differential privacy proposed so far provide support for reasoning about privacy but not for reasoning about the accuracy of data analyses. To overcome this limitation, in this work we present DPella, a programming framework providing data analysts with support for reasoning about privacy, accuracy and their trade-offs. The distinguishing feature of DPella is a novel component which statically tracks the accuracy of different data analyses. In order to make tighter accuracy estimations, this component leverages taint analysis for automatically inferring statistical independence of the different noise quantities added for guaranteeing privacy. We show the flexibility of our approach by not only implementing classical counting queries (e.g., CDFs) but also by analyzing hierarchical counting queries (like those done by Census Bureaus), where accuracy have different constraints per level and data analysts should figure out the best manner to calibrate privacy to meet the accuracy requirements.


title: Finding Maximum Edge-Disjoint Paths Between Multiple Terminals

abstract: Let $G=(V,E)$ be a multigraph with a set $T\subseteq V$ of terminals. A path in $G$ is called a $T$-path if its ends are distinct vertices in $T$ and no internal vertices belong to $T$. In 1978, Mader showed a characterization of the maximum number of edge-disjoint $T$-paths. In this paper, we provide a combinatorial, deterministic algorithm for finding the maximum number of edge-disjoint $T$-paths. The algorithm adopts an augmenting path approach. More specifically, we utilize a new concept of short augmenting walks in auxiliary labeled graphs to capture a possible augmentation of the number of edge-disjoint $T$-paths. To design a search procedure for a short augmenting walk, we introduce blossoms analogously to the matching algorithm of Edmonds (1965). When the search procedure terminates without finding a short augmenting walk, the algorithm provides a certificate for the optimality of the current edge-disjoint $T$-paths. From this certificate, one can obtain the Edmonds--Gallai type decomposition introduced by Seb\H{o} and Szeg\H{o} (2004). The algorithm runs in $O(|E|^2)$ time, which is much faster than the best known deterministic algorithm based on a reduction to linear matroid parity. We also present a strongly polynomial algorithm for the maximum integer free multiflow problem, which asks for a nonnegative integer combination of $T$-paths maximizing the sum of the coefficients subject to capacity constraints on the edges.


title: Calorimetry for low-energy electrons using charge and light in liquid argon

abstract: Precise calorimetric reconstruction of 5-50 MeV electrons in liquid argon time projection chambers (LArTPCs) will enable the study of astrophysical neutrinos in DUNE and could enhance the physics reach of oscillation analyses. Liquid argon scintillation light has the potential to improve energy reconstruction for low-energy electrons over charge-based measurements alone. Here we demonstrate light-augmented calorimetry for low-energy electrons in a single-phase LArTPC using a sample of Michel electrons from decays of stopping cosmic muons in the LArIAT experiment at Fermilab. Michel electron energy spectra are reconstructed using both a traditional charge-based approach as well as a more holistic approach that incorporates both charge and light. A maximum-likelihood fitter, using LArIAT's well-tuned simulation, is developed for combining these quantities to achieve optimal energy resolution. A sample of isolated electrons is simulated to better determine the energy resolution expected for astrophysical electron-neutrino charged-current interaction final states. In LArIAT, which has very low wire noise and an average light yield of 18 pe/MeV, an energy resolution of $\sigma/E \simeq 9.3\%/\sqrt{E} \oplus 1.3\%$ is achieved. Samples are then generated with varying wire noise levels and light yields to gauge the impact of light-augmented calorimetry in larger LArTPCs. At a charge-readout signal-to-noise of S/N $\simeq$ 30, for example, the energy resolution for electrons below 40 MeV is improved by $\approx$ 10%, $\approx$ 20%, and $\approx$ 40% over charge-only calorimetry for average light yields of 10 pe/MeV, 20 pe/MeV, and 100 pe/MeV, respectively.


title: Adaptive and Dynamically Constrained Process Noise Estimation for Orbit Determination

abstract: This paper introduces two new algorithms to accurately estimate the process noise covariance of a discrete-time Kalman filter online for robust orbit determination in the presence of dynamics model uncertainties. Common orbit determination process noise techniques, such as state noise compensation and dynamic model compensation, require offline tuning and a priori knowledge of the dynamical environment. Alternatively, the process noise covariance can be estimated through adaptive filtering. However, many adaptive filtering techniques are not applicable to onboard orbit determination due to computational cost or the assumption of a linear time-invariant system. Furthermore, existing adaptive filtering techniques do not constrain the process noise covariance according to the underlying continuous-time dynamical model, and there has been limited work on adaptive filtering with colored process noise. To overcome these limitations, a novel approach is developed which optimally fuses state noise compensation and dynamic model compensation with covariance matching adaptive filtering. This yields two adaptive and dynamically constrained process noise covariance estimation techniques. Unlike many adaptive filtering approaches, the new techniques accurately extrapolate over measurement outages and do not rely on ad hoc methods to ensure the process noise covariance is positive semi-definite. The benefits of the proposed algorithms are demonstrated through two case studies: an illustrative linear system and the autonomous navigation of two spacecraft orbiting an asteroid.


title: Distributed Function Minimization in Apache Spark

abstract: We report on an open-source implementation for distributed function minimization on top of Apache Spark by using gradient and quasi-Newton methods. We show-case it with an application to Optimal Transport and some scalability tests on classification and regression problems.


title: Lightfield Coordinates Adapted to Asgeirsson's Theorem

abstract: John's differential equation and its canonical form, the ultrahyperbolic equation, plays important role in lightfield imaging. The equation describes a local constraint on the lightfield, that was first observed as a "dimensionality gap" in the frequency representation. Related to the ultrahyperbolic equation, Asgeirsson's theorems describe global properties. These indicate new, global, constraints on the lightfield. In order to help validate those theorems on real captured images, we introduce a coordinate system for the lightfield, which suits better the Asgeirsson theorems, and analyze behaviour in terms of the new coordinates.