title: Time-dependent stochastic basis adaptation for uncertainty quantification

abstract: We extend stochastic basis adaptation and spatial domain decomposition methods to solve time varying stochastic partial differential equations (SPDEs) with a large number of input random parameters. Stochastic basis adaptation allows the determination of a low dimensional stochastic basis representation of a quantity of interest (QoI). Extending basis adaptation to time-dependent problems is challenging because small errors introduced in the previous time steps of the low dimensional approximate solution accumulate over time and cause divergence from the true solution. To address this issue we have introduced an approach where the basis adaptation varies at every time step so that the low dimensional basis is adapted to the QoI at that time step. We have coupled the time-dependent basis adaptation with domain decomposition to further increase the accuracy in the representation of the QoI. To illustrate the construction, we present numerical results for one-dimensional time varying linear and nonlinear diffusion equations with random space-dependent diffusion coefficients. Stochastic dimension reduction techniques proposed in the literature have mainly focused on quantifying the uncertainty in time independent and scalar QoI. To the best of our knowledge, this is the first attempt to extend dimensional reduction techniques to time varying and spatially dependent quantities such as the solution of SPDEs.


title: Technical Leverage in a Software Ecosystem: Development Opportunities and Security Risks

abstract: In finance, leverage is the ratio between assets borrowed from others and one's own assets. A matching situation is present in software: by using free open-source software (FOSS) libraries a developer leverages on other people's code to multiply the offered functionalities with a much smaller own codebase. In finance as in software, leverage magnifies profits when returns from borrowing exceed costs of integration, but it may also magnify losses, in particular in the presence of security vulnerabilities. We aim to understand the level of technical leverage in the FOSS ecosystem and whether it can be a potential source of security vulnerabilities. Also, we introduce two metrics change distance and change direction to capture the amount and the evolution of the dependency on third-party libraries. The application of the proposed metrics on 8494 distinct library versions from the FOSS Maven-based Java libraries shows that small and medium libraries (less than 100KLoC) have disproportionately more leverage on FOSS dependencies in comparison to large libraries. We show that leverage pays off as leveraged libraries only add a 4% delay in the time interval between library releases while providing four times more code than their own. However, libraries with such leverage (i.e., 75% of libraries in our sample) also have 1.6 higher odds of being vulnerable in comparison to the libraries with lower leverage. We provide an online demo for computing the proposed metrics for real-world software libraries available under the following URL: https://techleverage.eu/.


title: Non-minimizing connecting orbits for multi-well systems

abstract: Given a nonnegative, smooth potential $V: \mathbb{R}^k \to \mathbb{R}$ ($k \geq 2$) with multiple zeros, we say that a curve $\mathfrak{q}: \mathbb{R} \to \mathbb{R}^k$ is a connecting orbit if it solves the autonomous system of ordinary differential equations \begin{equation} \mathfrak{q}''= \nabla_{\mathbf{u}} V(\mathfrak{q}) , \hspace{2mm} \mbox{ in } \mathbb{R} \end{equation} and tends to a zero of $V$ at $\pm \infty$. Broadly, our goal is to study the existence of connecting orbits for the problem above using variational methods. Despite the rich previous literature concerning the existence of connecting orbits for other types of second order systems, to our knowledge only connecting orbits which minimize the associated energy functional in a suitable function space were proven to exist for autonomous multi-well potentials. The contribution of this paper is to provide, for a class of such potentials, some existence results regarding non-minimizing connecting orbits. Our results are closely related to the ones in the same spirit obtained by J. Bisgard in his PhD thesis (University of Wisconsin-Madison, 2005), where non-autonomous periodic multi-well potentials (ultimately excluding autonomous potentials) are considered. Our approach is based on several refined versions of the classical Mountain Pass Lemma.


title: Self-supervised 3D Representation Learning of Dressed Humans from Social Media Videos

abstract: A key challenge of learning a visual representation for the 3D high fidelity geometry of dressed humans lies in the limited availability of the ground truth data (e.g., 3D scanned models), which results in the performance degradation of 3D human reconstruction when applying to real-world imagery. We address this challenge by leveraging a new data resource: a number of social media dance videos that span diverse appearance, clothing styles, performances, and identities. Each video depicts dynamic movements of the body and clothes of a single person while lacking the 3D ground truth geometry. To learn a visual representation from these videos, we present a new self-supervised learning method to use the local transformation that warps the predicted local geometry of the person from an image to that of another image at a different time instant. This allows self-supervision by enforcing a temporal coherence over the predictions. In addition, we jointly learn the depths along with the surface normals that are highly responsive to local texture, wrinkle, and shade by maximizing their geometric consistency. Our method is end-to-end trainable, resulting in high fidelity depth estimation that predicts fine geometry faithful to the input real image. We further provide a theoretical bound of self-supervised learning via an uncertainty analysis that characterizes the performance of the self-supervised learning without training. We demonstrate that our method outperforms the state-of-the-art human depth estimation and human shape recovery approaches on both real and rendered images.


title: Heat flux in general quasifree fermionic right mover/left mover systems

abstract: With the help of time-dependent scattering theory on the observable algebra of infinitely extended quasifree fermionic chains, we introduce a general class of so-called right mover/left mover states which are inspired by the nonequilibrium steady states for the prototypical nonequilibrium configuration of a finite sample coupled to two thermal reservoirs at different temperatures. Under the assumption of spatial translation invariance, we relate the 2-point operator of such a right mover/left mover state to the asymptotic velocity of the system and prove that the system is thermodynamically nontrivial in the sense that its entropy production rate is strictly positive. Our study of these not necessarily gauge-invariant systems covers and substantially generalizes well-known quasifree fermionic chains and opens the way for a more systematic analysis of the heat flux in such systems.


title: On MCMC for variationally sparse Gaussian processes: A pseudo-marginal approach

abstract: Gaussian processes (GPs) are frequently used in machine learning and statistics to construct powerful models. However, when employing GPs in practice, important considerations must be made, regarding the high computational burden, approximation of the posterior, choice of the covariance function and inference of its hyperparmeters. To address these issues, Hensman et al. (2015) combine variationally sparse GPs with Markov chain Monte Carlo (MCMC) to derive a scalable, flexible and general framework for GP models. Nevertheless, the resulting approach requires intractable likelihood evaluations for many observation models. To bypass this problem, we propose a pseudo-marginal (PM) scheme that offers asymptotically exact inference as well as computational gains through doubly stochastic estimators for the intractable likelihood and large datasets. In complex models, the advantages of the PM scheme are particularly evident, and we demonstrate this on a two-level GP regression model with a nonparametric covariance function to capture non-stationarity.


title: Calculation of the nonlinear response functions of intra-exciton transitions in two-dimensional transition metal dichalcogenides

abstract: In this paper, we study the third-order nonlinear optical response due to transitions between excitonic levels in two-dimensional transition metal dichalcogeniedes. To accomplish this, we use methods not applied to the description of excitons in two-dimensional materials so far and combined with a variational approach to describe the $1s$ excitonic state. The aforementioned transitions allow to probe dark states which are not revealed in absorption experiments. We present general formulas capable of describing any third-order process. The specific case of two-photon absorption in WSe2 is studied. The case of the circular well is also studied as a benchmark of the theory.


title: Distribution-free uncertainty quantification for classification under label shift

abstract: Trustworthy deployment of ML models requires a proper measure of uncertainty, especially in safety-critical applications. We focus on uncertainty quantification (UQ) for classification problems via two avenues -- prediction sets using conformal prediction and calibration of probabilistic predictors by post-hoc binning -- since these possess distribution-free guarantees for i.i.d. data. Two common ways of generalizing beyond the i.i.d. setting include handling covariate and label shift. Within the context of distribution-free UQ, the former has already received attention, but not the latter. It is known that label shift hurts prediction, and we first argue that it also hurts UQ, by showing degradation in coverage and calibration. Piggybacking on recent progress in addressing label shift (for better prediction), we examine the right way to achieve UQ by reweighting the aforementioned conformal and calibration procedures whenever some unlabeled data from the target distribution is available. We examine these techniques theoretically in a distribution-free framework and demonstrate their excellent practical performance.


title: Fitting an Analytic Magnetic Field to a Prestellar Core

abstract: We deploy and demonstrate the capabilities of the magnetic field model developed by Ewertowski & Basu (2013) by fitting observed polarimetry data of the prestellar core FeSt 1-457. The analytic hourglass magnetic field function derived directly from Maxwell's equations yields a central-to-surface magnetic field strength ratio in the equatorial plane, as well as magnetic field directions with relative magnitudes throughout the core. This fit emerges from a comparison of a single plane of the model with the polarization map that results from the integrated properties of the magnetic field and dust throughout the core. Importantly, our fit is independent of any assumed density profile of the core. We check the robustness of the fit by using the POLARIS code to create synthetic polarization maps that result from the integrated scattering and emission properties of the dust grains and their radiative transfer, employing an observationally-motivated density profile. We find that the synthetic polarization maps obtained from the model also provides a good fit to the observed polarimetry. Our model fits the striking feature of significant curvature of magnetic field lines in the outer part of FeSt 1-457. Combined with independent column density estimates, we infer that the core of size $R_{\rm gas}$ has a mildly supercritical mass-to-flux ratio and may have formed through dynamical motions starting from a significantly larger radius $R$. A breakdown of flux-freezing through neutral-ion slip (ambipolar diffusion) could be responsible for effecting such a transition from a large-scale magnetic field structure to a more compact gas structure.


title: Hard-label Manifolds: Unexpected Advantages of Query Efficiency for Finding On-manifold Adversarial Examples

abstract: Designing deep networks robust to adversarial examples remains an open problem. Likewise, recent zeroth order hard-label attacks on image classification models have shown comparable performance to their first-order, gradient-level alternatives. It was recently shown in the gradient-level setting that regular adversarial examples leave the data manifold, while their on-manifold counterparts are in fact generalization errors. In this paper, we argue that query efficiency in the zeroth-order setting is connected to an adversary's traversal through the data manifold. To explain this behavior, we propose an information-theoretic argument based on a noisy manifold distance oracle, which leaks manifold information through the adversary's gradient estimate. Through numerical experiments of manifold-gradient mutual information, we show this behavior acts as a function of the effective problem dimensionality and number of training points. On real-world datasets and multiple zeroth-order attacks using dimension-reduction, we observe the same universal behavior to produce samples closer to the data manifold. This results in up to two-fold decrease in the manifold distance measure, regardless of the model robustness. Our results suggest that taking the manifold-gradient mutual information into account can thus inform better robust model design in the future, and avoid leakage of the sensitive data manifold.


title: A case study of non-Fourier heat conduction using Internal Variables and GENERIC

abstract: Applying simultaneously the methodology of Non-Equilibrium Thermodynamics with Internal Variables (NET-IV) and the framework of General Equation for the Non-Equilibrium Reversible-Irreversible Coupling (GENERIC), we demonstrate that, in heat conduction theories, entropy current multipliers can be interpreted as relaxed state variables. Fourier's law and its various extensions -- the Maxwell-Cattaneo-Vernotte, Guyer-Krumhansl, Jeffreys type, Ginzburg-Landau (Allen-Cahn) type and ballistic-diffusive -- heat conduction equations are derived in both formulations. Along these lines, a comparison of NET-IV and GENERIC is also performed. Our results may pave the way for microscopic/multiscale understanding of beyond-Fourier heat conduction, and open new ways for numerical simulations of heat-conduction problems.


title: Null expectations and null hypothesis testing for the species abundance distribution

abstract: The number of elements (N) and types (S) sampled from an ecological system are among the most powerful constraints on observations of abundance, distribution, and diversity. Together, N and S determine sets of possible forms (i.e., feasible sets) for the species abundance distribution (SAD). There are three approaches to the description of the null SAD (= the average feasible SAD). The first approach is based on the random uniform sampling of surjections. I calculate the probability of a given SAD, given N and S, under this approach (Eq. 4). The second approach is based on the random sampling of compositions. I calculate the probability of a given SAD, given N and S, under this approach (Eq. 8). The third approach is based on the random uniform sampling of partitions. I review the approach, which was developed by Locey & White (2013), and provide some asymptotic results useful for ecologists. The center of a feasible set is a null expectation, which should deviate enough from the alternative models before invoking for biological or ecological mechanisms underlying the SAD. Here, I integrate the feasible set approach with the typical framework of inference in ecology (goodness-of-fit, null hypothesis testing, model comparison, null modelling). I describe how to perform numerical simulations to describe expectations under different approaches to the feasible set. I develop objective or fitness functions to allow the estimation of the most likely SAD using numerical optimization. I provide tools to compare null expectations based on the feasible set approach with the observations, in the context of model comparison and null hypothesis testing.


title: Evaluation of Complexity Measures for Deep Learning Generalization in Medical Image Analysis

abstract: The generalization performance of deep learning models for medical image analysis often decreases on images collected with different devices for data acquisition, device settings, or patient population. A better understanding of the generalization capacity on new images is crucial for clinicians' trustworthiness in deep learning. Although significant research efforts have been recently directed toward establishing generalization bounds and complexity measures, still, there is often a significant discrepancy between the predicted and actual generalization performance. As well, related large empirical studies have been primarily based on validation with general-purpose image datasets. This paper presents an empirical study that investigates the correlation between 25 complexity measures and the generalization abilities of supervised deep learning classifiers for breast ultrasound images. The results indicate that PAC-Bayes flatness-based and path norm-based measures produce the most consistent explanation for the combination of models and data. We also investigate the use of multi-task classification and segmentation approach for breast images, and report that such learning approach acts as an implicit regularizer and is conducive toward improved generalization.


title: Smooth matching of $\hat{q}$ from hadronic to quark and gluon degrees of freedom

abstract: One of the key signatures of the Quark Gluon Plasma (QGP) is the energy loss of high momentum particles as they traverse the strongly interacting medium. The energy loss of these particles is governed by the jet transport coefficient $\hat{q}/T^3$, wherein it has been thought that there is a large jump as the system transitions between the hadron gas and Quark Gluon Plasma phases. Here we calculate $\hat{q}/T^3$ within the Hadron Resonance Gas (HRG) model with the particle list PDG16+ and find that, if one incorporates the experimental error in the hadronic calculation of $\hat{q}/T^3$ and assumes a higher pseudo-critical temperature, then a smooth transition from the hadron gas phase into the Quark Gluon Plasma phase is possible. We also find a significant enhancement in $\hat{q}/T^3$ at finite baryon chemical potential and find issues with the relationship between the shear viscosity and the jet transport coefficient within a hadron gas phase.


title: Large Graph Convolutional Network Training with GPU-Oriented Data Communication Architecture

abstract: Graph Convolutional Networks (GCNs) are increasingly adopted in large-scale graph-based recommender systems. Training GCN requires the minibatch generator traversing graphs and sampling the sparsely located neighboring nodes to obtain their features. Since real-world graphs often exceed the capacity of GPU memory, current GCN training systems keep the feature table in host memory and rely on the CPU to collect sparse features before sending them to the GPUs. This approach, however, puts tremendous pressure on host memory bandwidth and the CPU. This is because the CPU needs to (1) read sparse features from memory, (2) write features into memory as a dense format, and (3) transfer the features from memory to the GPUs. In this work, we propose a novel GPU-oriented data communication approach for GCN training, where GPU threads directly access sparse features in host memory through zero-copy accesses without much CPU help. By removing the CPU gathering stage, our method significantly reduces the consumption of the host resources and data access latency. We further present two important techniques to achieve high host memory access efficiency by the GPU: (1) automatic data access address alignment to maximize PCIe packet efficiency, and (2) asynchronous zero-copy access and kernel execution to fully overlap data transfer with training. We incorporate our method into PyTorch and evaluate its effectiveness using several graphs with sizes up to 111 million nodes and 1.6 billion edges. In a multi-GPU training setup, our method is 65-92% faster than the conventional data transfer method, and can even match the performance of all-in-GPU-memory training for some graphs that fit in GPU memory.


title: Secure Software Development in the Era of Fluid Multi-party Open Software and Services

abstract: Pushed by market forces, software development has become fast-paced. As a consequence, modern development projects are assembled from 3rd-party components. Security & privacy assurance techniques once designed for large, controlled updates over months or years, must now cope with small, continuous changes taking place within a week, and happening in sub-components that are controlled by third-party developers one might not even know they existed. In this paper, we aim to provide an overview of the current software security approaches and evaluate their appropriateness in the face of the changed nature in software development. Software security assurance could benefit by switching from a process-based to an artefact-based approach. Further, security evaluation might need to be more incremental, automated and decentralized. We believe this can be achieved by supporting mechanisms for lightweight and scalable screenings that are applicable to the entire population of software components albeit there might be a price to pay.


title: An Agent-based Model to Evaluate Interventions on Online Dating Platforms to Decrease Racial Homogamy

abstract: Perhaps the most controversial questions in the study of online platforms today surround the extent to which platforms can intervene to reduce the societal ills perpetrated on them. Up for debate is whether there exist any effective and lasting interventions a platform can adopt to address, e.g., online bullying, or if other, more far-reaching change is necessary to address such problems. Empirical work is critical to addressing such questions. But it is also challenging, because it is time-consuming, expensive, and sometimes limited to the questions companies are willing to ask. To help focus and inform this empirical work, we here propose an agent-based modeling (ABM) approach. As an application, we analyze the impact of a set of interventions on a simulated online dating platform on the lack of long-term interracial relationships in an artificial society. In the real world, a lack of interracial relationships are a critical vehicle through which inequality is maintained. Our work shows that many previously hypothesized interventions online dating platforms could take to increase the number of interracial relationships from their website have limited effects, and that the effectiveness of any intervention is subject to assumptions about sociocultural structure. Further, interventions that are effective in increasing diversity in long-term relationships are at odds with platforms' profit-oriented goals. At a general level, the present work shows the value of using an ABM approach to help understand the potential effects and side effects of different interventions that a platform could take.


title: Measurement of the generalized spin polarizabilities of the neutron in the low $Q^2$ region

abstract: Understanding the nucleon spin structure in the regime where the strong interaction becomes truly strong poses a challenge to both experiment and theory. At energy scales below the nucleon mass of about 1 GeV, the intense interaction among the quarks and gluons inside the nucleon makes them highly correlated. Their coherent behaviour causes the emergence of effective degrees of freedom, requiring the application of non-perturbative techniques, such as chiral effective field theory. Here, we present measurements of the neutron's generalized spin-polarizabilities that quantify the neutron's spin precession under electromagnetic fields at very low energy-momentum transfer squared down to 0.035 GeV$^2$. In this regime, chiral effective field theory calculations are expected to be applicable. Our data, however, show a strong discrepancy with these predictions, presenting a challenge to the current description of the neutron's spin properties.


title: On uncertainties in the reconstruction of nanostructures in EUV scatterometry and grazing incidence small-angle X-ray scattering

abstract: Increasing miniaturization and complexity of nanostructures require innovative metrology solutions with high throughput that can assess complex 3D structures in a non-destructive manner. EUV scatterometry is investigated for the characterization of nanostructured surfaces. The reconstruction is based on a rigorous simulation using a Maxwell solver based on finite-elements and is statistically validated with a Markov-Chain Monte Carlo sampling method. Here it is shown that this method is suitable for the dimensional characterization of the nanostructures and the investigation of oxide or contamination layers. In comparison to grazing-incidence small-angle X-rayscattering (GISAXS) EUV allows to probe smaller areas. The influence of the divergence on the diffracted intensities in EUV is much lower than in GISAXS, which also reduces the computational effort of the reconstruction.


title: A Systematic Evaluation of Transfer Learning and Pseudo-labeling with BERT-based Ranking Models

abstract: Due to high annotation costs making the best use of existing human-created training data is an important research direction. We, therefore, carry out a systematic evaluation of transferability of BERT-based neural ranking models across five English datasets. Previous studies focused primarily on zero-shot and few-shot transfer from a large dataset to a dataset with a small number of queries. In contrast, each of our collections has a substantial number of queries, which enables a full-shot evaluation mode and improves reliability of our results. Furthermore, since source datasets licences often prohibit commercial use, we compare transfer learning to training on pseudo-labels generated by a BM25 scorer. We find that training on pseudo-labels -- possibly with subsequent fine-tuning using a modest number of annotated queries -- can produce a competitive or better model compared to transfer learning. Yet, it is necessary to improve the stability and/or effectiveness of the few-shot training, which, sometimes, can degrade performance of a pretrained model.


title: Triangles and triple products of Laplace eigenfunctions

abstract: Consider an $L^2$-normalized Laplace-Beltrami eigenfunction $e_\lambda$ on a compact, boundary-less Riemannian manifold with $\Delta e_\lambda = -\lambda^2 e_\lambda$. We study eigenfunction triple products \[ \langle e_\lambda e_\mu, e_\nu \rangle = \int e_\lambda e_\mu \overline{e_\nu} \, dV. \] We show the overall $\ell^2$-concentration of these triple products is determined by the measure of some set of configurations of triangles with side lengths equal to the frequencies $\lambda,\mu,$ and $\nu$. A rapidly vanishing proportion of this mass lies in the `classically forbidden' regime where $\lambda, \mu,$ and $\nu$ fail to satisfy the triangle inequality. As a consequence, we improve a result by Lu, Sogge, and Steinerberger.


title: Revisiting Priority $k$-Center: Fairness and Outliers

abstract: In the Priority $k$-Center problem, the input consists of a metric space $(X,d)$, an integer $k$, and for each point $v \in X$ a priority radius $r(v)$. The goal is to choose $k$-centers $S \subseteq X$ to minimize $\max_{v \in X} \frac{1}{r(v)} d(v,S)$. If all $r(v)$'s are uniform, one obtains the $k$-Center problem. Plesn\'ik [Plesn\'ik, Disc. Appl. Math. 1987] introduced the Priority $k$-Center problem and gave a $2$-approximation algorithm matching the best possible algorithm for $k$-Center. We show how the problem is related to two different notions of fair clustering [Harris et al., NeurIPS 2018; Jung et al., FORC 2020]. Motivated by these developments we revisit the problem and, in our main technical contribution, develop a framework that yields constant factor approximation algorithms for Priority $k$-Center with outliers. Our framework extends to generalizations of Priority $k$-Center to matroid and knapsack constraints, and as a corollary, also yields algorithms with fairness guarantees in the lottery model of Harris et al [Harris et al, JMLR 2019].


title: Stabilization of periodic sweeping processes and asymptotic average velocity for soft locomotors with dry friction

abstract: We study the asymptotic stability of periodic solutions for sweeping processes defined by a polyhedron with translationally moving faces. Previous results are improved by obtaining a stronger $W^{1,2}$ convergence. Then we present an application to a model of crawling locomotion. Our stronger convergence allows us to prove the stabilization of the system to a running-periodic (or derivo-periodic, or relative-periodic) solution and the well-posedness of an average asymptotic velocity depending only on the gait adopted by the crawler. Finally, we discuss some examples of finite-time versus asymptotic-only convergence.


title: Optimal Control of Differentially Flat Systems is Surprisingly Easy

abstract: As we move to increasingly complex cyber-physical systems (CPS), new approaches are needed to plan efficient state trajectories in real-time. In this paper, we propose an approach to significantly reduce the complexity of solving optimal control problems for a class of CPS with nonlinear dynamics. We exploit the property of differential flatness to simplify the Euler-Lagrange equations that arise during optimization, and this simplification eliminates the numerical instabilities that plague optimal control in general. We also present an explicit differential equation that describes the evolution of the optimal state trajectory, and we extend our results to consider both the unconstrained and constrained cases. Furthermore, we demonstrate the performance of our approach by generating the optimal trajectory for a planar manipulator with two revolute joints. We show in simulation that our approach is able to generate the constrained optimal trajectory in $4.5$ ms while respecting workspace constraints and switching between a `left' and `right' bend in the elbow joint.


title: Legendrian Satellites and Decomposable Cobordisms

abstract: We investigate the interactions between the Legendrian satellite construction and the existence of exact, orientable Lagrangian cobordisms between Legendrian knots. Given Lagrangian cobordisms between two Legendrian knots and between two Legendrian tangles, we construct a Lagrangian cobordism between Legendrian satellites of the knots by the closures of the tangles, with extra twists on both the top and the bottom satellite to compensate for the genus of the cobordism. If the original cobordisms were decomposable, then a decomposable cobordism between satellites exists as well, again with extra twists.


title: PSFs for mapping artificial night sky luminance over large territories

abstract: Knowledge of the night sky radiance over a large territory may be valuable informationto identify sites appropriate to astronomical observations or for assessing the impacts ofartificial light at night on ecosystems. Measuring the sky radiance can be a complex endeavourdepending on the desired temporal and spatial resolution. Similarly, modelling of artificialnight sky radiance for multiple points of a territory can represent a significant amount ofcomputing time depending on the complexity of the model used. The use of the convolutionof a point spread function with the light sources geographical distribution has been suggestedin order to model the sky radiance over large territories of hundreds of kilometres in size.We determine how the point spread function is sensitive to the main driving parameters ofthe artificial night sky radiance such as the wavelength, the ground reflectance, the obstaclesproperties, the Upward Light Output Ratio and the Aerosol Optical Depth using the Illuminav2 model. The obtained functions were used to model the artificial night sky brightness ofthe Mont-M\'egantic International Dark Sky Reserve for winter and summer conditions. Theresults were compared to the New world atlas of artificial night sky brightness, the Illuminav2 model and in situ Sky Quality Camera measurements. We found that the New world atlasoverestimates the artificial sky brightness by 55% whereas the Illumina model underestimatesit by 48%. This may be due to varying atmospherical conditions and the fact that the modelonly accounts for public light sources.


title: Performance Analysis and Optimization of Uplink Cellular Networks with Flexible Frame Structure

abstract: Future wireless cellular networks must support both enhanced mobile broadband (eMBB) and ultra reliable low latency communication (URLLC) to manage heterogeneous data traffic for emerging wireless services. To achieve this goal, a promising technique is to enable flexible frame structure by dynamically changing the data frame's numerology according to the channel information as well as traffic quality of service requirements. However, due to nonorthogonal subcarriers, this technique can result in an interference, known as inter numerology interference (INI), thus, degrading the network performance. In this work, a novel framework is proposed to analyze the INI in the uplink cellular communications. In particular, a closed form expression is derived for the INI power in the uplink with a flexible frame structure, and a new resource allocation problem is formulated to maximize the network spectral efficiency (SE) by jointly optimizing the power allocation and numerology selection in a multi user uplink scenario. The simulation results validate the derived theoretical INI analyses and provide guidelines for power allocation and numerology selection.


title: A converse of the Banach contraction principle for partial metric spaces and the continuum hypothesis

abstract: A version of the Bessaga inverse of the Banach contraction principle for partial metric spaces is presented. Equivalence of that version and the continuum hypothesis is shown as well.


title: WaveGuard: Understanding and Mitigating Audio Adversarial Examples

abstract: There has been a recent surge in adversarial attacks on deep learning based automatic speech recognition (ASR) systems. These attacks pose new challenges to deep learning security and have raised significant concerns in deploying ASR systems in safety-critical applications. In this work, we introduce WaveGuard: a framework for detecting adversarial inputs that are crafted to attack ASR systems. Our framework incorporates audio transformation functions and analyses the ASR transcriptions of the original and transformed audio to detect adversarial inputs. We demonstrate that our defense framework is able to reliably detect adversarial examples constructed by four recent audio adversarial attacks, with a variety of audio transformation functions. With careful regard for best practices in defense evaluations, we analyze our proposed defense and its strength to withstand adaptive and robust attacks in the audio domain. We empirically demonstrate that audio transformations that recover audio from perceptually informed representations can lead to a strong defense that is robust against an adaptive adversary even in a complete white-box setting. Furthermore, WaveGuard can be used out-of-the box and integrated directly with any ASR model to efficiently detect audio adversarial examples, without the need for model retraining.


title: $^{16}\mathrm{O}^{16}\mathrm{O}$ at RHIC and the LHC comparing $\alpha$ clustering vs substructure

abstract: Collisions of light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed $^{16}\mathrm{O}^{16}\mathrm{O}$ run at the LHC and RHIC we study the potential for finding $\alpha$ clustering in $^{16}$O. Here we use the state-of-the-art iEBE-VISHNU package with $^{16}$O nucleonic configurations from {\rm ab initio} nuclear lattice simulations. This setup was tuned using a Bayesian analysis on pPb and PbPb systems. We find that the $^{16}\mathrm{O}^{16}\mathrm{O}$ system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between $\alpha$-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.


title: Maximising the total weight of on-time jobs on parallel machines subject to a conflict graph

abstract: The paper considers scheduling on parallel machines under the constraint that some pairs of jobs cannot be processed concurrently. Each job has an associated weight, and all jobs have the same deadline. The objective is to maximise the total weight of on-time jobs. The problem is known to be strongly NP-hard in general. A polynomial-time algorithm for scheduling unit execution time jobs on two machines is proposed. The performance of a broad family of approximation algorithms for scheduling unit execution time jobs on more than two machines is analysed. For the case of arbitrary job processing times, two integer linear programming formulations are proposed and compared with two formulations known from the earlier literature. An iterated variable neighborhood search algorithm is also proposed and evaluated by means of computational experiments.


title: Broadband Imaging to Study the Spectral Distribution of Meteor Radio Afterglows

abstract: We present observations of 86 meteor radio afterglows (MRAs) using the new broadband imager at the Long Wavelength Array Sevilleta (LWA-SV) station. The MRAs were detected using the all-sky images with a bandwidth up to 20 MHz. We fit the spectra with both a power law and a log-normal function. When fit with a power law, the spectra varied from flat to steep and the derived spectral index distribution from the fit peaked at -1.73. When fit with a log-normal function, the spectra exhibits turnovers at frequencies between 30-40 MHz, and appear to be a better functional fit to the spectra. We compared the spectral parameters from the two fitting methods with the physical properties of MRAs. We observe a weak correlation between the log-normal turnover frequency and the altitude of MRAs. The spectral indices from the power law fit do not show any strong correlations with the physical properties of MRAs. However, the full width half maximum (FWHM) duration of MRAs is correlated with the local time, incidence angle, luminosity and optically derived kinetic energy of parent meteoroid. Also, the average luminosity of MRAs seems to be correlated with the kinetic energy of parent meteoroid and the altitude at which they occur.

categories: astro-ph.EP astro-ph.IM physics.plasm-ph physics.space-ph


title: From viral evolution to spatial contagion: a biologically modulated Hawkes model

abstract: Mutations sometimes increase contagiousness for evolving pathogens. During an epidemic, scientists use viral genome data to infer a shared evolutionary history and connect this history to geographic spread. We propose a model that directly relates a pathogen's evolution to its spatial contagion dynamics -- effectively combining the two epidemiological paradigms of phylogenetic inference and self-exciting process modeling -- and apply this \emph{phylogenetic Hawkes process} to a Bayesian analysis of 23,422 viral cases from the 2014-2016 Ebola outbreak in West Africa. The proposed model is able to detect individual viruses with significantly elevated rates of spatiotemporal propagation for a subset of 1,610 samples that provide genome data. Finally, to facilitate model application in big data settings, we develop massively parallel implementations for the gradient and Hessian of the log-likelihood and apply our high performance computing framework within an adaptively preconditioned Hamiltonian Monte Carlo routine.


title: Bound states of a quartic and sextic inverse-powerlaw potential for all angular momenta

abstract: We use the tridiagonal representation approach to solve the radial Schr\"odinger equation for an inverse power-law potential of a combined quartic and sextic degrees and for all angular momenta. The amplitude of the quartic singularity is larger than that of the sextic but the signs are negative and positive, respectively. It turns out that the system has a finite number of bound states, which is determined by the larger ratio of the two singularity amplitudes. The solution is written as a finite series of square integrable functions written in terms of the Bessel polynomial.


title: Avoiding Monochromatic Solutions to 3-term Equations

abstract: Given an equation, the integers $[n] = \{1, 2, \dots, n\}$ as inputs, and the colors red and blue, how can we color $[n]$ in order to minimize the number of monochromatic solutions to the equation, and what is the minimum? The answer is only known for a handful of equations, but much progress has been made on improving upper and lower bounds on minima for various equations. A well-studied characteristic an equation, which has its roots in graph Ramsey theory, is to determine if the minimum number of monochromatic solutions can be achieved (asymptotically) by uniformly random colorings. Such equations are called common. We prove that no 3-term equations are common and provide a lower bound for a specific class of 3-term equations.


title: So you want to be a Super Researcher?

abstract: Publishing original scientific research is inherent to the work of a researcher. However, the pressure to maintain productivity and scientific impact can lead to research group publishing excessively, negatively affecting the mental health of a researcher. Ph.D. students and early career researchers are particularly susceptible to this pressure due to the inherent vulnerability of their positions. At present, there are no resources that concisely summarise the publication culture of a research group to help the researcher make an informed decision before joining. In this article, we present the 'Super Researcher' app, an R Shiny application(app) with a user-friendly interface. Using text-mining methodology to extract publicly available author data from Scopus, this pilot app has four fundamental functions to provide snapshot information that will help researchers grasp the publication culture of a research group within minutes. The 'Super Researcher' app provides information on: 1) institution data, 2) author's publication, 3) co-author network plots and 4) publication journals. The 'Super Researcher' app is built on R shiny which provides an interactive interface to users. This app utilizes the Big Data framework Apache Spark to mine relevant information from a huge author information database. The author's information is stored and manipulated using both SQL(SQLite) and NoSQL(HBase) databases. Hbase is used for local data storage and manipulation while SQLite feeds data to the R Shiny interface. In this paper, we introduce these functionalities and illustrate how this information can help guide a researcher to select a new Principle Investigator (PI) with better compatibility in terms of publication attitude using a case study. Available: https://researchmind.co.uk/super-researcher/


title: Manipulation of Orbital Angular Momentum Spectrum Using Shape-Tailored Metasurfaces

abstract: Vortex beams carrying orbital angular momentum (OAM) have been widely applied in various electromagnetic, optical, and quantum systems. A tailored OAM spectrum composed of several specific modes as expected holds a promise for expanding the degrees of freedom of the systems. However, such a broadband high-purity tailored spectrum is difficult to be achieved by the present devices, where the broadband amplitude manipulation has not been explored yet. In this work, inspired by the envelope-modulation theory, an elegant and universal way to manipulate the OAM spectrum in wide bandwidth is proposed by using a shape-tailored metasurface. Firstly, the rotating meta-atoms on a triangular lattice are proved to have smaller coupling distortion than that on a square lattice, and this behavior is critical for high-purity vortex spectrum generation by the Pancharatnam-Berry-based metasurfaces. Secondly, a universal modulation relation is established between the spatial arrangement of metasurfaces and the generated vortex beams. Finally, the broadband modulated OAM spectra and the comb-like OAM spectra are theoretically and experimentally demonstrated by the shape-tailored metasurfaces. The proposed amplitude-modulation scheme offers a novel concept and engineering route to manipulate the OAM spectrum in wide bandwidth, which could promote the development of OAM-based applications.


title: Broadband photonic RF channelizer with a Kerr soliton crystal micro-comb

abstract: We report a 92 channel RF channelizer based on a 48.9 GHz integrated micro-comb that operates via soliton crystals, together with a passive high-Q ring resonator that acts as a periodic filter with an optical 3dB bandwidth of 121.4 MHz. We obtain an instant RF bandwidth of 8.08 GHz and 17.55 GHz achieved through temperature tuning. These results represent a major advance to achieving fully integrated photonic RF spectrum channelizers with reduced low complexity, size, and high performance for digital-compatible signal detection and broadband analog signal processing.


title: Optical data transmission field trial @ 44Tb/s with a 49GHz Kerr soliton crystal microcomb

abstract: We report world record high data transmission over standard optical fiber from a single optical source. We achieve a line rate of 44.2 Terabits per second (Tb/s) employing only the C-band at 1550nm, resulting in a spectral efficiency of 10.4 bits/s/Hz. We use a new and powerful class of micro-comb called soliton crystals that exhibit robust operation and stable generation as well as a high intrinsic efficiency that, together with an extremely low spacing of 48.9 GHz enables a very high coherent data modulation format of 64 QAM. We achieve error free transmission across 75 km of standard optical fiber in the lab and over a field trial with a metropolitan optical fiber network. This work demonstrates the ability of optical micro-combs to exceed other approaches in performance for the most demanding practical optical communications applications.


title: An integral's journey over the real line

abstract: In 1826 Cauchy presented an Integral over the real line. Al and I thought a derivation would be mighty fine. So we packed our contour integral bags that day, and we now present an analytic continuation this time.


title: Conservation laws for a $q$-deformed nonrelativistic particle

abstract: We derive $q$-versions of Green's theorem from the Leibniz rules of partial derivatives for the $q$-deformed Euclidean space. Using these results and the Schr\"{o}dinger equations for a $q$-deformed nonrelativistic particle, we derive continuity equations for the probability density, the energy density, and the momentum density of a $q$-deformed nonrelativistic particle.


title: Note on the Application of Divergent Series for Finding a Particular Solution to a Nonhomogeneous Linear Ordinary Differential Equation with Constant Coefficients

abstract: There are many methods for finding a particular solution to a nonhomogeneous linear ordinary differential equation (ODE) with constant coefficients. The method of undetermined coefficients, Laplace transform method and differential operator method are generally known. The latter mentioned method sometimes uses the Maclaurin expansion of an inverse differential operator but only in the case when the obtained series is convergent. The present work deals also with how to find a particular solution if the corresponding infinite series is divergent, using only the terms of that series and the method of summation of divergent series.


title: Non-invasive Current Density Imaging of Lithium-Ion Batteries

abstract: The rapid pace of replacing fossil fuel propelled transport by electric vehicles is critically dependent on high-performing, high energy density and efficient batteries. Optimal and safe use of existing battery cells and development of much-needed novel battery chemistries and geometries require a large range of diagnostic and monitoring tools. While structural and chemical information is readily extracted through a host of imaging techniques, non-invasive functional detection of interior battery processes remains limited. Here we introduce sensitive magnetometry performed outside the battery that reveals a battery cell's internal current distribution. As a key application, we use an array of sensors to image the magnetic field present under cycling of a pouch cell between charge states. We find good agreement between measured and modelled fields with sufficient resolution to detect percent-level deviations around areas of high current density. This opens the path towards rapid and reliable assessment throughout the battery life cycle, from battery development and manufacturing quality assurance to operational safety and optimised use.


title: Cognitive Homeostatic Agents

abstract: Human brain has been used as an inspiration for building autonomous agents, but it is not obvious what level of computational description of the brain one should use. This has led to overly opinionated symbolic approaches and overly unstructured connectionist approaches. We propose that using homeostasis as the computational description provides a good compromise. Similar to how physiological homeostasis is the regulation of certain homeostatic variables, cognition can be interpreted as the regulation of certain 'cognitive homeostatic variables'. We present an outline of a Cognitive Homeostatic Agent, built as a hierarchy of physiological and cognitive homeostatic subsystems and describe structures and processes to guide future exploration. We expect this to be a fruitful line of investigation towards building sophisticated artificial agents that can act flexibly in complex environments, and produce behaviors indicating planning, thinking and feelings.


title: Pure and Linear Frequency Converter Temporal Metasurface

abstract: Metasurfaces are ultrathin structures which are constituted by an array of subwavelength scatterers with designable scattering responses. They have opened up unprecedented exciting opportunities for extraordinary wave engineering processes. On the other hand, frequency converters have drawn wide attention due to their vital applications in telecommunication systems, health care devices, radio astronomy, military radars and biological sensing systems. Here, we show that a spurious-free and linear frequency converter metasurface can be realized by leveraging unique properties of engineered transmissive temporal supercells. Such a metasurface is formed by time-modulated supercells; themselves are composed of temporal and static patch resonators and phase shifters. This represents the first frequency converter metasurface possessing large frequency conversion ratio with controllable frequency bands and transmission magnitude. In contrast to conventional nonlinear mixers, the proposed temporal frequency converter offers a linear response. In addition, by taking advantage of the proposed surface-interconnector-phaser-surface (SIPS) architecture, a spurious-free and linear frequency conversion is achievable, where all undesired mixing products are strongly suppressed. The proposed metasurface may be digitally controlled and programmed through a field programmable gate array. This makes the spurious-free and linear frequency converter metasurface a prominent solution for wireless and satellite telecommunication systems, as well as invisibility cloaks and radars. This study opens a way to realize more complicated and enhanced-efficiency spectrum-changing metasurface.


title: From Quantifying Vagueness To Pan-niftyism

abstract: In this short paper, we will introduce a simple model for quantifying philosophical vagueness. There is growing interest in this endeavor to quantify vague concepts of consciousness, agency, etc. We will then discuss some of the implications of this model including the conditions under which the quantification of `nifty' leads to pan-nifty-ism. Understanding this leads to an interesting insight - the reason a framework to quantify consciousness like Integrated Information Theory implies (forms of) panpsychism is because there is favorable structure already implicitly encoded in the construction of the quantification metric.


title: Information Encoding/Decoding using the Memory Effect in Fractional-order Capacitive Devices

abstract: In this study, we show that the discharge voltage pattern of a fractional-order supercapacitor from the same initial steady-state voltage into a constant resistor is dependent on the past charging voltage profile. The charging voltage was designed to follow a power-law function, i.e. $v_c(t)=V_{cc} \left( {t}/{t_{ss}}\right)^p \;(0<t \leqslant t_{ss})$, in which $t_{ss}$ (charging time duration between zero voltage to the terminal voltage $V_{cc}$) and $p$ ($0<p<1$) act as two variable parameters. We used this history-dependence of the dynamic behavior of the device to uniquely retrieve information pre-coded in the charging waveform pattern. Furthermore, we provide an analytical model based on fractional calculus that explains phenomenologically the information storage mechanism. The use of this intrinsic material memory effect may lead to new types of methods for information storage and retrieval.


title: Koopman Operator Based Modeling for Quadrotor Control on $SE(3)$

abstract: In this paper, we propose a Koopman operator based approach to describe the nonlinear dynamics of a quadrotor on SE(3) in terms of an infinite-dimensional linear system which evolves in the space of observable functions (lifted space) and which is more appropriate for control design purposes. The major challenge when using the Koopman operator is the characterization of a set of observable functions that can span the lifted space. Recent methods either use tools from machine learning to learn the observable functions or guess a suitable set of observables that best describes the nonlinear dynamics. Instead of guessing or learning the observables, in this work we derive them in a systematic way for the quadrotor dynamics on SE(3). In addition, we prove that the proposed sequence of observable functions converges pointwise to the zero function, which allows us to select only a finite set of observable functions to form (an approximation of) the lifted space. Our theoretical analysis is also confirmed by numerical simulations which demonstrate that by increasing the dimension of the lifted space, the derived linear state space model can approximate the nonlinear quadrotor dynamics more accurately.


title: Spacetime Paths as a Whole

abstract: The mathematical similarities between non-relativistic wavefunction propagation in quantum mechanics and image propagation in scalar diffraction theory are used to develop a novel understanding of time and paths through spacetime as a whole. It is well known that Feynman's original derivation of the path integral formulation of non-relativistic quantum mechanics uses time-slicing to calculate amplitudes as sums over all possible paths through space, but along a definite curve through time. Here, a 3+1D spacetime wave distribution and its 4-momentum dual are formally developed which have no external time parameter and therefore cannot change or evolve in the usual sense. Time is thus seen "from the outside". A given 3+1D momentum representation of a system encodes complete dynamical information, describing the system's spacetime behavior as a whole. A comparison is made to the mathematics of holograms, and properties of motion for simple systems are derived.


title: Electro-optically modulated polarization mode conversion in lithium niobate ridge waveguides

abstract: Lithium niobate on insulator (LNOI) waveguides, as an emerging technology, have proven to offer a promising platform for integrated optics, due to their strong optical confinement comparable to silicon on insulator (SOI) waveguides, while possessing the versatile properties of lithium niobate, such as high electro-optic coefficients. In this paper, we show that mode hybridization, a phenomenon widely found in vertically asymmetric waveguides, can be efficiently modulated in an LNOI ridge waveguide by electro-optic effect, leading to a polarization mode converter with 97% efficiency. Moreover, the proposed device does not require tapering or periodic poling, thereby greatly simplifying the fabrication process. It can also be actively switched by external fields. Such a platform facilitates technological progress of photonic circuits and sensors.


title: Fukaya categories of higher-genus surfaces and pants decompositions

abstract: In this paper we prove a local-to-global principle for the Fukaya category of a closed Riemann surface $\Sigma$ of genus $g \geq 2$. We show that $\mathrm{Fuk}(\Sigma)$ can be glued from the Fukaya categories of the pairs-of-pants making up a pants decomposition of $\Sigma$. This extends our earlier results for the case of punctured Riemann surfaces. Our result has several interesting consequences: we obtain simple proofs of old and new HMS statements for Riemann surfaces, and establish a geometrization theorem for the objects of $\mathrm{Fuk}(\Sigma)$.


title: Training artificial neural networks for precision orientation and strain mapping using 4D electron diffraction datasets

abstract: Techniques for training artificial neural networks (ANNs) and convolutional neural networks (CNNs) using simulated dynamical electron diffraction patterns are described. The premise is based on the following facts. First, given a suitable crystal structure model and scattering potential, electron diffraction patterns can be simulated accurately using dynamical diffraction theory. Secondly, using simulated diffraction patterns as input, ANNs can be trained for the determination of crystal structural properties, such as crystal orientation and local strain. Further, by applying the trained ANNs to four-dimensional diffraction datasets (4D-DD) collected using the scanning electron nanodiffraction (SEND) or 4D scanning transmission electron microscopy (4D-STEM) techniques, the crystal structural properties can be mapped at high spatial resolution. Here, we demonstrate the ANN-enabled possibilities for the analysis of crystal orientation and strain at high precision and benchmark the performance of ANNs and CNNs by comparing with previous methods. A factor of thirty improvement in angular resolution at 0.009 degrees (0.16 mrad) for orientation mapping, sensitivity at 0.04% or less for strain mapping, and improvements in computational performance are demonstrated.


title: Casimir Effect Nearby and Through a Cosmological Wormhole

abstract: In this letter, we investigate the changes in the quantum vacuum energy density of a massless scalar field inside a Casimir cavity that orbits a wormhole, by considering the cosmological model with an isotropic form of the Morris-Thorne wormhole, embedded in the FLRW universe. In this sense, we examine the effects of its global curvature and scale factor in an instant of the cosmic history, besides the influences of the local geometry as well as of inertial forces, on the Casimirenergy density. We also study the behavior of this quantity when each plate is fixed without rotation at the opposite sides of the wormhole throat, at zero and finite temperatures, taking into account the effective distance between the plates through the wormhole throat.


title: X-Ray Scattering from Light-Driven Spin Fluctuations in a Doped Mott Insulator

abstract: Manipulating spin fluctuations with ultrafast laser pulses is a promising route to dynamically control collective phenomena in strongly correlated materials. However, understanding how photoexcited spin degrees of freedom evolve at a microscopic level requires a momentum- and energy-resolved characterization of their nonequilibrium dynamics. Here, we study the photoinduced dynamics of finite-momentum spin excitations in two-dimensional Mott insulators on a square lattice. By calculating the time-resolved resonant inelastic x-ray scattering cross-section, we show that an ultrafast pump above the Mott gap induces a prompt softening of the spin excitation energy, compatible with a transient renormalization of the exchange interaction. While spin fluctuations in a hole-doped system (paramagnons) are well described by Floquet theory, magnons at half filling are found to deviate from this picture. Furthermore, we show that the paramagnon softening is accompanied by an ultrafast suppression of $d$-wave pairing correlations, indicating a link between the transient spin excitation dynamics and superconducting pairing far from equilibrium.


title: Multi-task Learning with High-Dimensional Noisy Images

abstract: Recent medical imaging studies have given rise to distinct but inter-related datasets corresponding to multiple experimental tasks or longitudinal visits. Standard scalar-on-image regression models that fit each dataset separately are not equipped to leverage information across inter-related images, and existing multi-task learning approaches are compromised by the inability to account for the noise that is often observed in images. We propose a novel joint scalar-on-image regression framework involving wavelet-based image representations with grouped penalties that are designed to pool information across inter-related images for joint learning, and which explicitly accounts for noise in high-dimensional images via a projection-based approach. In the presence of non-convexity arising due to noisy images, we derive non-asymptotic error bounds under non-convex as well as convex grouped penalties, even when the number of voxels increases exponentially with sample size. A projected gradient descent algorithm is used for computation, which is shown to approximate the optimal solution via well-defined non-asymptotic optimization error bounds under noisy images. Extensive simulations and application to a motivating longitudinal Alzheimer's disease study illustrate significantly improved predictive ability and greater power to detect true signals, that are simply missed by existing methods without noise correction due to the attenuation to null phenomenon.


title: Atom-in-jellium equations of state and melt curves in the white dwarf regime

abstract: Atom-in-jellium calculations of the electron states, and perturbative calculations of the Einstein frequency, were used to construct equations of state (EOS) from around $10^{-5}$ to $10^7$g/cm$^3$ and $10^{-4}$ to $10^{6}$eV for elements relevant to white dwarf (WD) stars. This is the widest range reported for self-consistent electronic shell structure calculations. Elements of the same ratio of atomic weight to atomic number were predicted to asymptote to the same $T=0$ isotherm, suggesting that, contrary to recent studies of the crystallization of WDs, the amount of gravitational energy that could be released by separation of oxygen and carbon is small. A generalized Lindemann criterion based on the amplitude of the ion-thermal oscillations calculated using atom-in-jellium theory, previously used to extrapolate melt curves for metals, was found to reproduce previous thermodynamic studies of the melt curve of the one component plasma with a choice of vibration amplitude consistent with low pressure results. For elements for which low pressure melting satisfies the same amplitude criterion, such as Al, this melt model thus gives a likely estimate of the melt curve over the full range of normal electronic matter; for the other elements, it provides a useful constraint on the melt locus.

categories: astro-ph.SR cond-mat.mtrl-sci physics.comp-ph physics.plasm-ph


title: Optimization of graded filleted lattice structures subject to yield and buckling constraints

abstract: To reduce the stress concentration and ensure the structural safety for lattice structure designs, in this paper, a new optimization framework is developed for the optimal design of graded lattice structures, innovatively integrating fillet designs as well as yield and elastic buckling constraints. Both strut and fillet radii are defined as design variables. Homogenization method is employed to characterize the effective elastic constants and yield stresses of the lattice metamaterials. Metamaterial models are developed to represent the relationships between the metamaterial effective properties and lattice geometric variables. A yield constraint, based on the modified Hills yield criterion, is developed as a function of relative strut radii and fillet parameters. An elastic buckling constraint, based on the Euler buckling formula and the Johnson formula, is developed as a function of relative strut radii. Both yield and buckling constraints are integrated into an optimization problem formulation; a new optimization framework is proposed and a case study of minimizing the compliance of a Messerschmitt-Bolkow-Blohm beam is conducted. The yield and buckling constraints guarantee the safety of the optimized beams composed of BCC and PC lattices. Reductions in compliance and stress concentration are achieved by the optimized MBB beams.


title: Neural model robustness for skill routing in large-scale conversational AI systems: A design choice exploration

abstract: Current state-of-the-art large-scale conversational AI or intelligent digital assistant systems in industry comprises a set of components such as Automatic Speech Recognition (ASR) and Natural Language Understanding (NLU). For some of these systems that leverage a shared NLU ontology (e.g., a centralized intent/slot schema), there exists a separate skill routing component to correctly route a request to an appropriate skill, which is either a first-party or third-party application that actually executes on a user request. The skill routing component is needed as there are thousands of skills that can either subscribe to the same intent and/or subscribe to an intent under specific contextual conditions (e.g., device has a screen). Ensuring model robustness or resilience in the skill routing component is an important problem since skills may dynamically change their subscription in the ontology after the skill routing model has been deployed to production. We show how different modeling design choices impact the model robustness in the context of skill routing on a state-of-the-art commercial conversational AI system, specifically on the choices around data augmentation, model architecture, and optimization method. We show that applying data augmentation can be a very effective and practical way to drastically improve model robustness.


title: Motif-driven Dense Subgraph Discovery in Directed and Labeled Networks

abstract: Dense regions in networks are an indicator of interesting and unusual information. However, most existing methods only consider simple, undirected, unweighted networks. Complex networks in the real-world often have rich information though: edges are asymmetrical and nodes/edges have categorical and numerical attributes. Finding dense subgraphs in such networks in accordance with this rich information is an important problem with many applications. Furthermore, most existing algorithms ignore the higher-order relationships (i.e., motifs) among the nodes. Motifs are shown to be helpful for dense subgraph discovery but their wide spectrum in heterogeneous networks makes it challenging to utilize them effectively. In this work, we propose quark decomposition framework to locate dense subgraphs that are rich with a given motif. We focus on networks with directed edges and categorical attributes on nodes/edges. For a given motif, our framework builds subgraphs, called quarks, in varying quality and with hierarchical relations. Our framework is versatile, efficient, and extendible. We discuss the limitations and practical instantiations of our framework as well as the role confusion problem that needs to be considered in directed networks. We give an extensive evaluation of our framework in directed, signed-directed, and node-labeled networks. We consider various motifs and evaluate the quark decomposition using several real-world networks. Results show that quark decomposition performs better than the state-of-the-art techniques. Our framework is also practical and scalable to networks with up to 101M edges.


title: Nutrition5k: Towards Automatic Nutritional Understanding of Generic Food

abstract: Understanding the nutritional content of food from visual data is a challenging computer vision problem, with the potential to have a positive and widespread impact on public health. Studies in this area are limited to existing datasets in the field that lack sufficient diversity or labels required for training models with nutritional understanding capability. We introduce Nutrition5k, a novel dataset of 5k diverse, real world food dishes with corresponding video streams, depth images, component weights, and high accuracy nutritional content annotation. We demonstrate the potential of this dataset by training a computer vision algorithm capable of predicting the caloric and macronutrient values of a complex, real world dish at an accuracy that outperforms professional nutritionists. Further we present a baseline for incorporating depth sensor data to improve nutrition predictions. We will publicly release Nutrition5k in the hope that it will accelerate innovation in the space of nutritional understanding.


title: DeepLocalize: Fault Localization for Deep Neural Networks

abstract: Deep neural networks (DNNs) are becoming an integral part of most software systems. Previous work has shown that DNNs have bugs. Unfortunately, existing debugging techniques do not support localizing DNN bugs because of the lack of understanding of model behaviors. The entire DNN model appears as a black box. To address these problems, we propose an approach that automatically determines whether the model is buggy or not, and identifies the root causes. Our key insight is that historic trends in values propagated between layers can be analyzed to identify faults, and localize faults. To that end, we first enable dynamic analysis of deep learning applications: by converting it into an imperative representation and alternatively using a callback mechanism. Both mechanisms allows us to insert probes that enable dynamic analysis over the traces produced by the DNN while it is being trained on the training data. We then conduct dynamic analysis over the traces to identify the faulty layer that causes the error. We propose an algorithm for identifying root causes by capturing any numerical error and monitoring the model during training and finding the relevance of every layer on the DNN outcome. We have collected a benchmark containing 40 buggy models and patches that contain real errors in deep learning applications from Stack Overflow and GitHub. Our benchmark can be used to evaluate automated debugging tools and repair techniques. We have evaluated our approach using this DNN bug-and-patch benchmark, and the results showed that our approach is much more effective than the existing debugging approach used in the state of the practice Keras library. For 34 out of 40 cases, our approach was able to detect faults whereas the best debugging approach provided by Keras detected 32 out of 40 faults. Our approach was able to localize 21 out of 40 bugs whereas Keras did not localize any faults.


title: Linear Depth Deduction with Subformula Property for Intuitionistic Epistemic Logic

abstract: In their seminal paper Artemov and Protopopescu provide Hilbert formal systems, Brower-Heyting-Kolmogorov and Kripke semantics for the logics of intuitionistic belief and knowledge. Subsequently Krupski has proved that the logic of intuitionistic knowledge is PSPACE-complete and Su and Sano have provided calculi enjoying the subformula property. This paper continues the investigations around to sequent calculi for Intuitionistic Epistemic Logics by providing sequent calculi that have the subformula property and that are terminating in linear depth. Our calculi allow us to design a procedure that for invalid formulas returns a Kripke model of minimal depth. Finally we also discuss refutational sequent calculi, that is sequent calculi to prove the invalidity.


title: Electromagnetic Signatures from the Tidal Tail of a Black Hole - Neutron Star Merger

abstract: Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave (GW) observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation, and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle due to two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of $\sim 3$ with orientation while the peak in the optical bands varies by $\sim 3$ magnitudes. The spectrum is blue-shifted at viewing angles along the bulk motion, which increases the $V$-band peak magnitude to $\sim -14$ despite the lanthanide-rich composition.


title: Shapes and recession cones in mixed-integer convex representability

abstract: Mixed-integer convex representable (MICP-R) sets are those sets that can be represented exactly through a mixed-integer convex programming formulation. Following up on recent work by Lubin et al. (2017, 2020) we investigate structural geometric properties of MICP-R sets, which strongly differentiate them from the class of mixed-integer linear representable sets (MILP-R). First, we provide an example of an MICP-R set which is the countably infinite union of convex sets with countably infinitely many different recession cones. This is in sharp contrast with MILP-R sets which are at most infinite unions of polyhedra that share the same recession cone. Second, we provide an example of an MICP-R set which is the countably infinite union of polytopes all of which have different shapes (no pair is combinatorially equivalent, which implies they are not affine transformations of each other). Again, this is in sharp contrast with MILP-R sets which are at most infinite unions of polyhedra that are all translations of a finite subset of themselves. Interestingly, we show that a countably infinite union of convex sets sharing the same volume can be MICP-R only if the sets are all translations of a finite subset of themselves (i.e. the natural conceptual analogue to the MILP-R case).


title: A squeezed quantum microcomb on a chip

abstract: The optical microresonator-based frequency comb (microcomb) provides a versatile platform for nonlinear physics studies and has wide applications ranging from metrology to spectroscopy. Deterministic quantum regime is an unexplored aspect of microcombs, in which unconditional entanglements among hundreds of equidistant frequency modes can serve as critical ingredients to scalable universal quantum computing and quantum networking. Here, we demonstrate a deterministic quantum microcomb in a silica microresonator on a silicon chip. 40 continuous-variable quantum modes, in the form of 20 simultaneously two-mode squeezed comb pairs, are observed within 1 THz optical span at telecommunication wavelengths. A maximum raw squeezing of 1.6 dB is attained. A high-resolution spectroscopy measurement is developed to characterize the frequency equidistance of quantum microcombs. Our demonstration offers the possibility to leverage deterministically generated, frequency multiplexed quantum states and integrated photonics to open up new avenues in fields of spectroscopy, quantum metrology, and scalable quantum information processing.


title: Heat kernels for reflected diffusions with jumps on inner uniform domains

abstract: In this paper, we study sharp two-sided heat kernel estimates for a large class of symmetric reflected diffusions with jumps on the closure of an inner uniform domain $D$ in a length metric space. The length metric is the intrinsic metric of a strongly local Dirichlet form. When $D$ is an inner uniform domain in the Euclidean space, a prototype for a special case of the processes under consideration are symmetric reflected diffusions with jumps on $D$, whose infinitesimal generators are non-local (pseudo-differential) operators $L$ on $D$ of the form $$ L u(x) =\frac12 \sum_{i, j=1}^d \frac{\partial}{\partial x_i} \left(a_{ij}(x) \frac{\partial u(x)}{\partial x_j}\right) + \lim_{\eps \downarrow 0} \int_{\{y\in D: \, \rho_D(y, x)>\eps\}} (u(y)-u(x)) J(x, y)\, dy $$ satisfying "Neumann boundary condition". Here, $\rho_D(x,y)$ is the length metric on $D$, $A(x)=(a_{ij}(x))_{1\leq i,j\leq d}$ is a measurable $d\times d$ matrix-valued function on $D$ that is uniformly elliptic and bounded, and $$ J(x,y):= \frac{1}{\Phi(\rho_D(x,y))} \int_{[\alpha_1, \alpha_2]} \frac{c(\alpha, x,y)} {\rho_D(x,y)^{d+\alpha}} \,\nu(d\alpha) , $$ where $\nu$ is a finite measure on $[\alpha_1, \alpha_2] \subset (0, 2)$, $\Phi$ is an increasing function on $[ 0, \infty )$ with $c_1e^{c_2r^{\beta}} \le \Phi(r) \le c_3 e^{c_4r^{\beta}}$ for some $\beta \in [0,\infty]$, and $c(\alpha , x, y)$ is a jointly measurable function that is bounded between two positive constants and is symmetric in $(x, y)$.


title: The Merger Dynamics of the Galaxy Cluster Abell 1775: New Insights from Chandra and XMM-Newton for a Cluster Simultaneously Hosting a WAT and a NAT Radio Sources

abstract: We present a new study of the merger dynamics of Abell~1775 by analyzing the high-quality Chandra and XMM-Newton archival data. We confirm/identify an arc-shaped edge (i.e., the head) at $\sim48$~kpc west of the X-ray peak, a split cold gas tail that extends eastward to $\sim163$~kpc, and a plume of spiral-like X-ray excess (within about $81-324$~kpc northeast of the cluster core) that connects to the end of the tail. The head, across which the projected gas temperature rises outward from $3.39_{-0.18}^{+0.28}$~keV to $5.30_{-0.43}^{+0.54}$~keV, is found to be a cold front with a Mach number of $\mathcal{M}\sim0.79$. Along the surfaces of the cold front and tail, typical KHI features (noses and wings, etc.) are found and are used to constrain the upper limit of the magnetic field ($\sim11.2~\mu$G) and the viscosity suppression factor ($\sim0.01$). Combining optical and radio evidence we propose a two-body merger (instead of systematic motion in a large-scale gas environment) scenario and have carried out idealized hydrodynamic simulations to verify it. We find that the observed X-ray emission and temperature distributions can be best reproduced with a merger mass ratio of 5 after the first pericentric passage. The NAT radio galaxy is thus more likely to be a single galaxy falling into the cluster center at a relative velocity of 2800~$\rm km~s^{-1}$, a speed constrained by its radio morphology. The infalling subcluster is expected to have a relatively low gas content, because only a gas-poor subcluster can cause central-only disturbances as observed in such an off-axis merger.


title: Onshore to Near-Shore Outsourcing Transitions: Unpacking Tensions

abstract: This study is directed towards highlighting tensions of incoming and outgoing vendors during outsourcing in a near-shore context. Incoming-and-outgoing of vendors generate a complex form of relationship in which the participating organizations cooperate and compete simultaneously. It is of great importance to develop knowledge about this kind of relationship typically in the current GSE-related multi-sourcing environment. We carried out a longitudinal case study and utilized data from the 'Novopay' project, which is available in the public domain. This project involved an outgoing New Zealand based vendor and incoming Australian based vendor. The results show that the demand for the same human resources, dependency upon cooperation and collaboration between vendors, reliance on each other system's configurations and utilizing similar strategies by the client, which worked for the previous vendor, generated a set of tensions which needed to be continuously managed throughout the project.


title: Re-polarisation of macrophages within a multi-scale moving boundary tumour invasion model

abstract: Cancer invasion of the surrounding tissue is a multiscale process that involves not only tumour cells but also other immune cells in the environment, such as the tumour-associated macrophages (TAMs). The heterogeneity of these immune cells, with the two extremes being the pro-inflammatory and anti-tumour M1 cells, and the anti-inflammatory and pro-tumour M2 cells, has a significant impact on cancer invasion as these cell interact in different ways with the tumour cells and with the ExtraCellular Matrix (ECM). Experimental studies have shown that cancer cells co-migrate with TAMs, but the impact of these different TAM sub-populations (which can change their phenotype and re-polarise depending on the microenvironment) on this co-migration is not fully understood. In this study, we extend a previous multi-scale moving boundary mathematical model, by introducing the M1-like macrophages alongside with their exerted multi-scale effects on the tumour invasion process. With the help of this model we investigate numerically the impact of re-polarising the M2 TAMs into the anti-tumoral M1 phenotype and how such a strategy affects the overall tumour progression. In particular, we investigate numerically whether the M2-->M1 re-polarisation could depend on time and/or space, and what would be the macroscopic effects of this spatial- and temporal-dependent re-polarisation on tumour invasion.


title: Gaussian processes meet NeuralODEs: A Bayesian framework for learning the dynamics of partially observed systems from scarce and noisy data

abstract: This paper presents a machine learning framework (GP-NODE) for Bayesian systems identification from partial, noisy and irregular observations of nonlinear dynamical systems. The proposed method takes advantage of recent developments in differentiable programming to propagate gradient information through ordinary differential equation solvers and perform Bayesian inference with respect to unknown model parameters using Hamiltonian Monte Carlo sampling and Gaussian Process priors over the observed system states. This allows us to exploit temporal correlations in the observed data, and efficiently infer posterior distributions over plausible models with quantified uncertainty. Moreover, the use of sparsity-promoting priors such as the Finnish Horseshoe for free model parameters enables the discovery of interpretable and parsimonious representations for the underlying latent dynamics. A series of numerical studies is presented to demonstrate the effectiveness of the proposed GP-NODE method including predator-prey systems, systems biology, and a 50-dimensional human motion dynamical system. Taken together, our findings put forth a novel, flexible and robust workflow for data-driven model discovery under uncertainty. All code and data accompanying this manuscript are available online at \url{https://github.com/PredictiveIntelligenceLab/GP-NODEs}.


title: Clusterability in Neural Networks

abstract: The learned weights of a neural network have often been considered devoid of scrutable internal structure. In this paper, however, we look for structure in the form of clusterability: how well a network can be divided into groups of neurons with strong internal connectivity but weak external connectivity. We find that a trained neural network is typically more clusterable than randomly initialized networks, and often clusterable relative to random networks with the same distribution of weights. We also exhibit novel methods to promote clusterability in neural network training, and find that in multi-layer perceptrons they lead to more clusterable networks with little reduction in accuracy. Understanding and controlling the clusterability of neural networks will hopefully render their inner workings more interpretable to engineers by facilitating partitioning into meaningful clusters.


title: PolarNet: Accelerated Deep Open Space Segmentation Using Automotive Radar in Polar Domain

abstract: Camera and Lidar processing have been revolutionized with the rapid development of deep learning model architectures. Automotive radar is one of the crucial elements of automated driver assistance and autonomous driving systems. Radar still relies on traditional signal processing techniques, unlike camera and Lidar based methods. We believe this is the missing link to achieve the most robust perception system. Identifying drivable space and occupied space is the first step in any autonomous decision making task. Occupancy grid map representation of the environment is often used for this purpose. In this paper, we propose PolarNet, a deep neural model to process radar information in polar domain for open space segmentation. We explore various input-output representations. Our experiments show that PolarNet is a effective way to process radar data that achieves state-of-the-art performance and processing speeds while maintaining a compact size.


title: Limits of Probabilistic Safety Guarantees when Considering Human Uncertainty

abstract: When autonomous robots interact with humans, such as during autonomous driving, explicit safety guarantees are crucial in order to avoid potentially life-threatening accidents. Many data-driven methods have explored learning probabilistic bounds over human agents' trajectories (i.e. confidence tubes that contain trajectories with probability $\delta$), which can then be used to guarantee safety with probability $1-\delta$. However, almost all existing works consider $\delta \geq 0.001$. The purpose of this paper is to argue that (1) in safety-critical applications, it is necessary to provide safety guarantees with $\delta < 10^{-8}$, and (2) current learning-based methods are ill-equipped to compute accurate confidence bounds at such low $\delta$. Using human driving data (from the highD dataset), as well as synthetically generated data, we show that current uncertainty models use inaccurate distributional assumptions to describe human behavior and/or require infeasible amounts of data to accurately learn confidence bounds for $\delta \leq 10^{-8}$. These two issues result in unreliable confidence bounds, which can have dangerous implications if deployed on safety-critical systems.


title: An Analytical Solution to the IMU Initialization Problem for Visual-Inertial Systems

abstract: The fusion of visual and inertial measurements is becoming more and more popular in the robotics community since both sources of information complement well each other. However, in order to perform this fusion, the biases of the Inertial Measurement Unit (IMU) as well as the direction of gravity must be initialized first. Additionally, in case of a monocular camera, the metric scale is also needed. The most popular visual-inertial initialization approaches rely on accurate vision-only motion estimates to build a non-linear optimization problem that solves for these parameters in an iterative way. In this paper, we rely on the previous work in [1] and propose an analytical solution to estimate the accelerometer bias, the direction of gravity and the scale factor in a maximum-likelihood framework. This formulation results in a very efficient estimation approach and, due to the non-iterative nature of the solution, avoids the intrinsic issues of previous iterative solutions. We present an extensive validation of the proposed IMU initialization approach and a performance comparison against the state-of-the-art approach described in [2] with real data from the publicly available EuRoC dataset, achieving comparable accuracy at a fraction of its computational cost and without requiring an initial guess for the scale factor. We also provide a C++ open source reference implementation.


title: An Effective Loss Function for Generating 3D Models from Single 2D Image without Rendering

abstract: Differentiable rendering is a very successful technique that applies to a Single-View 3D Reconstruction. Current renderers use losses based on pixels between a rendered image of some 3D reconstructed object and ground-truth images from given matched viewpoints to optimise parameters of the 3D shape. These models require a rendering step, along with visibility handling and evaluation of the shading model. The main goal of this paper is to demonstrate that we can avoid these steps and still get reconstruction results as other state-of-the-art models that are equal or even better than existing category-specific reconstruction methods. First, we use the same CNN architecture for the prediction of a point cloud shape and pose prediction like the one used by Insafutdinov & Dosovitskiy. Secondly, we propose the novel effective loss function that evaluates how well the projections of reconstructed 3D point clouds cover the ground truth object's silhouette. Then we use Poisson Surface Reconstruction to transform the reconstructed point cloud into a 3D mesh. Finally, we perform a GAN-based texture mapping on a particular 3D mesh and produce a textured 3D mesh from a single 2D image. We evaluate our method on different datasets (including ShapeNet, CUB-200-2011, and Pascal3D+) and achieve state-of-the-art results, outperforming all the other supervised and unsupervised methods and 3D representations, all in terms of performance, accuracy, and training time.


title: Gemini: Dynamic Bias Correction for Autonomous Experimentation and Molecular Simulation

abstract: Bayesian optimization has emerged as a powerful strategy to accelerate scientific discovery by means of autonomous experimentation. However, expensive measurements are required to accurately estimate materials properties, and can quickly become a hindrance to exhaustive materials discovery campaigns. Here, we introduce Gemini: a data-driven model capable of using inexpensive measurements as proxies for expensive measurements by correcting systematic biases between property evaluation methods. We recommend using Gemini for regression tasks with sparse data and in an autonomous workflow setting where its predictions of expensive to evaluate objectives can be used to construct a more informative acquisition function, thus reducing the number of expensive evaluations an optimizer needs to achieve desired target values. In a regression setting, we showcase the ability of our method to make accurate predictions of DFT calculated bandgaps of hybrid organic-inorganic perovskite materials. We further demonstrate the benefits that Gemini provides to autonomous workflows by augmenting the Bayesian optimizer Phoenics to yeild a scalable optimization framework leveraging multiple sources of measurement. Finally, we simulate an autonomous materials discovery platform for optimizing the activity of electrocatalysts for the oxygen evolution reaction. Realizing autonomous workflows with Gemini, we show that the number of measurements of a composition space comprising expensive and rare metals needed to achieve a target overpotential is significantly reduced when measurements from a proxy composition system with less expensive metals are available.


title: The Effect of Behavioral Probability Weighting in a Simultaneous Multi-Target Attacker-Defender Game

abstract: We consider a security game in a setting consisting of two players (an attacker and a defender), each with a given budget to allocate towards attack and defense, respectively, of a set of nodes. Each node has a certain value to the attacker and the defender, along with a probability of being successfully compromised, which is a function of the investments in that node by both players. For such games, we characterize the optimal investment strategies by the players at the (unique) Nash Equilibrium. We then investigate the impacts of behavioral probability weighting on the investment strategies; such probability weighting, where humans overweight low probabilities and underweight high probabilities, has been identified by behavioral economists to be a common feature of human decision-making. We show via numerical experiments that behavioral decision-making by the defender causes the Nash Equilibrium investments in each node to change (where the defender overinvests in the high-value nodes and underinvests in the low-value nodes).


title: A Platoon Formation Framework in a Mixed Traffic Environment

abstract: Connected and automated vehicles (CAVs) provide the most intriguing opportunity to reduce pollution, energy consumption, and travel delays. In this paper, we address the problem of vehicle platoon formation in a traffic network with partial CAV penetration rates. We investigate the interaction between CAV and human-driven vehicle (HDV) dynamics, and provide a rigorous control framework that enables platoon formation with the HDVs by only controlling the CAVs within the network. We present a complete analytical solution of the CAV control policy and the conditions under which a platoon formation is feasible. We evaluate the solution and demonstrate the efficacy of the proposed framework using numerical simulation.


title: Point Cloud based Hierarchical Deep Odometry Estimation

abstract: Processing point clouds using deep neural networks is still a challenging task. Most existing models focus on object detection and registration with deep neural networks using point clouds. In this paper, we propose a deep model that learns to estimate odometry in driving scenarios using point cloud data. The proposed model consumes raw point clouds in order to extract frame-to-frame odometry estimation through a hierarchical model architecture. Also, a local bundle adjustment variation of this model using LSTM layers is implemented. These two approaches are comprehensively evaluated and are compared against the state-of-the-art.


title: Rover Relocalization for Mars Sample Return by Virtual Template Synthesis and Matching

abstract: We consider the problem of rover relocalization in the context of the notional Mars Sample Return campaign. In this campaign, a rover (R1) needs to be capable of autonomously navigating and localizing itself within an area of approximately 50 x 50 m using reference images collected years earlier by another rover (R0). We propose a visual localizer that exhibits robustness to the relatively barren terrain that we expect to find in relevant areas, and to large lighting and viewpoint differences between R0 and R1. The localizer synthesizes partial renderings of a mesh built from reference R0 images and matches those to R1 images. We evaluate our method on a dataset totaling 2160 images covering the range of expected environmental conditions (terrain, lighting, approach angle). Experimental results show the effectiveness of our approach. This work informs the Mars Sample Return campaign on the choice of a site where Perseverance (R0) will place a set of sample tubes for future retrieval by another rover (R1).


title: Robust Electric-field Input Circuits for Clocked Molecular Quantum-dot Cellular Automata

abstract: Quantum-dot cellular automata (QCA) is a paradigm for low-power, general-purpose, classical computing designed to overcome the challenges facing CMOS in the extreme limits of scaling. A molecular implementation of QCA offers nanometer-scale devices with device densities and operating speeds which may surpass CMOS device densities and speeds by several orders of magnitude, all at room temperature. Here, a proposal for electric field bit write-in to molecular QCA circuits is extended to synchronous QCA circuits clocked using an applied electric field, \(\vec{E}\). Input electrodes, which may be much larger than the cells themselves, immerse an input circuit in an input field \(E_y \hat{y}\), in addition to the applied clocking field \(E_z \hat{z}\). The input field selects the input bit on a field-sensitive portion of the circuit. Another portion of the circuit with reduced \(E_y\)-sensitivity functions as a shift register, transmitting the input bit to downstream QCA logic for processing. It is shown that a simple rotation of the molecules comprising the shift register makes them immune to unwanted effects from the input field or fringing fields in the direction of the input field. Furthermore, the circuits also tolerate a significant unwanted field component \(E_x \hat{x}\) in the third direction, which is neither the clocking nor input direction. The write-in of classical bits to molecular QCA circuits is a road-block that must be cleared in order to realize energy-efficient molecular computation using QCA. The results presented here show that interconnecting shift registers may be designed to function in the presence of significant unwanted fringing fields from large input electrodes. Furthermore, the techniques devloped here may also enable molecular QCA logic to tolerate these same unwanted fringing fields.