title: Developing Biomarker Combinations in Multicenter Studies via Direct Maximization and Penalization

abstract: Motivated by a study of acute kidney injury, we consider the setting of biomarker studies involving patients at multiple centers where the goal is to develop a biomarker combination for diagnosis, prognosis, or screening. As biomarker studies become larger, this type of data structure will be encountered more frequently. In the presence of multiple centers, one way to assess the predictive capacity of a given combination is to consider the center-adjusted AUC (aAUC), a summary of the ability of the combination to discriminate between cases and controls in each center. Rather than using a general method, such as logistic regression, to construct the biomarker combination, we propose directly maximizing the aAUC. Furthermore, it may be desirable to have a biomarker combination with similar performance across centers. To that end, we allow for penalization of the variability in the center-specific AUCs. We demonstrate desirable asymptotic properties of the resulting combinations. Simulations provide small-sample evidence that maximizing the aAUC can lead to combinations with improved performance. We also use simulated data to illustrate the utility of constructing combinations by maximizing the aAUC while penalizing variability. Finally, we apply these methods to data from the study of acute kidney injury.


title: Direct Visual-Inertial Odometry with Semi-Dense Mapping

abstract: The paper presents a direct visual-inertial odometry system. In particular, a tightly coupled nonlinear optimization based method is proposed by integrating the recent advances in direct dense tracking and Inertial Measurement Unit (IMU) pre-integration, and a factor graph optimization is adapted to estimate the pose of the camera and rebuild a semi-dense map. Two sliding windows are maintained in the proposed approach. The first one, based on Direct Sparse Odometry (DSO), is to estimate the depths of candidate points for mapping and dense visual tracking. In the second one, measurements from the IMU pre-integration and dense visual tracking are fused probabilistically using a tightly-coupled, optimization-based sensor fusion framework. As a result, the IMU pre-integration provides additional constraints to suppress the scale drift induced by the visual odometry. Evaluations on real-world benchmark datasets show that the proposed method achieves competitive results in indoor scenes.


title: GENN: Predicting Correlated Drug-drug Interactions with Graph Energy Neural Networks

abstract: Gaining more comprehensive knowledge about drug-drug interactions (DDIs) is one of the most important tasks in drug development and medical practice. Recently graph neural networks have achieved great success in this task by modeling drugs as nodes and drug-drug interactions as links and casting DDI predictions as link prediction problems. However, correlations between link labels (e.g., DDI types) were rarely considered in existing works. We propose the graph energy neural network (GENN) to explicitly model link type correlations. We formulate the DDI prediction task as a structure prediction problem and introduce a new energy-based model where the energy function is defined by graph neural networks. Experiments on two real-world DDI datasets demonstrated that GENN is superior to many baselines without consideration of link type correlations and achieved $13.77\%$ and $5.01\%$ PR-AUC improvement on the two datasets, respectively. We also present a case study in which \mname can better capture meaningful DDI correlations compared with baseline models.


title: Propagation of Alfv\'en waves in the expanding solar wind with the fast-slow stream interaction

abstract: We carry out two-dimensional magnetohydrodynamic (MHD) simulations of an ensemble of Alfv\'enic fluctuations propagating in a structured, expanding solar wind including the presence of fast and slow solar wind streams. Using an appropriate expanding box model, the simulations incorporate the effects of fast-slow stream shear and compression and rarefaction self-consistently. We investigate the radial and longitudinal evolution of the cross-helicity, the total and residual energies and the power spectra of outward and inward Alfv\'enic fluctuations. The stream interaction is found to strongly affect the radial evolution of Alfv\'enic turbulence. The total energy in the Alfv\'en waves is depleted within the velocity shear regions, accompanied by the decrease of the normalized cross-helicity. The presence of stream-compression facilitates this process. Residual energy fluctuates around zero due to the correlation and de-correlation between the inward/outward waves but no net growth or decrease of the residual energy is observed. The radial power spectra of the inward/outward Alfv\'en waves show significant longitudinal variations. Kolmogorov-like spectra are developed only inside the fast and slow streams and when both the compression and shear are present. On the other hand, the spectra along the longitudinal direction show clear Kolmogorov-like inertial ranges in all cases.


title: Confederated Machine Learning on Horizontally and Vertically Separated Medical Data for Large-Scale Health System Intelligence

abstract: Health information is generally fragmented across silos. Though it is technically feasible to unite data for analysis in a manner that underpins a rapid learning healthcare system, privacy concerns and regulatory barriers limit data centralization. Machine learning can be conducted in a federated manner on patient datasets with the same set of variables, but separated across sites of care. But federated learning cannot handle the situation where different data types for a given patient are separated vertically across different organizations and when patient ID matching across different institutions is difficult. We call methods that enable machine learning model training on data separated by two or more degrees confederated machine learning. We proposed and evaluated a confederated learning to training machine learning model to stratify the risk of several diseases among when data are horizontally separated by individual, vertically separated by data type, and separated by identity without patient ID matching.


title: Entropy Stable h/p-Nonconforming Discretization with the Summation-by-Parts Property for the Compressible Euler and Navier-Stokes Equations

abstract: In this paper, the entropy conservative/stable algorithms presented by Del Rey Fernandez and coauthors [18,16,17] for the compressible Euler and Navier-Stokes equations on nonconforming p-refined/coarsened curvilinear grids is extended to h/p refinement/coarsening. The main difficulty in developing nonconforming algorithms is the construction of appropriate coupling procedures across nonconforming interfaces. Here, a computationally simple and efficient approach based upon using decoupled interpolation operators is utilized. The resulting scheme is entropy conservative/stable and element-wise conservative. Numerical simulations of the isentropic vortex and viscous shock propagation confirm the entropy conservation/stability and accuracy properties of the method (achieving ~ p + 1 convergence) which are comparable to those of the original conforming scheme [4,35]. Simulations of the Taylor-Green vortex at Re = 1,600 and turbulent flow past a sphere at Re = 2,000 show the robustness and stability properties of the overall spatial discretization for unstructured grids. Finally, to demonstrate the entropy conservation property of a fully-discrete explicit entropy stable algorithm with h/p refinement/coarsening, we present the time evolution of the entropy function obtained by simulating the propagation of the isentropic vortex using a relaxation Runge-Kutta scheme.


title: Tame fundamental groups of pure pairs and Abhyankar's lemma

abstract: Let $(R,\mathfrak{m}, k)$ be a strictly local normal $k$-domain of positive characteristic and $P$ be a prime divisor on $X=\text{Spec } R$. We study the Galois category of finite covers over $X$ that are at worst tamely ramified over $P$ in the sense of Grothendieck--Murre. Assuming that $(X,P)$ is a purely $F$-regular pair, our main result is that every Galois cover $f \: Y \to X$ in that Galois category satisfies that $\bigl(f^{-1}(P)\bigr)_{\text{red}}$ is a prime divisor. We shall explain why this should be thought as a (partial) generalization of a classical theorem due to S.S.~Abhyankar regarding the \'etale-local structure of tamely ramified covers between normal schemes with respect to a divisor with normal crossings. Additionally, we investigate the formal consequences this result has on the structure of the fundamental group representing the Galois category. We also obtain a characteristic zero analog by reduction to positive characteristics following Bhatt--Gabber--Olsson's methods.


title: A Convolution Bound Implies Tolerance to Time-variations and Unmodelled Dynamics

abstract: Recently it has been shown, in several settings, how to carry out adaptive control for an LTI plant so that a convolution bound holds on the closed-loop behavior; this, in turn, has been leveraged to prove robustness of the closed-loop system to time-varying parameters and unmodelled dynamics. The goal of this paper is to show that the same is true for a large class of finite-dimensional, nonlinear plant and controller combinations.


title: On the relation between Transversal and Longitudinal Scaling in Cities

abstract: Given that a group of cities follows a scaling law connecting urban population with socio-economic or infrastructural metrics (transversal scaling), should we expect that each city would follow the same behavior over time (longitudinal scaling)? This assumption has important policy implications, although rigorous empirical tests have been so far hindered by the lack of suitable data. Here, we advance the debate by looking into the temporal evolution of the scaling laws for 5507 municipalities in Brazil. We focus on the relationship between population size and two urban variables, GDP and water network length, analyzing the time evolution of the system of cities as well as their individual trajectory. We find that longitudinal (individual) scaling exponents are city-specific, but they are distributed around an average value that approaches to the transversal scaling exponent when the data are decomposed to eliminate external factors, and when we only consider cities with a sufficiently large growth rate. Such results give support to the idea that the longitudinal dynamics is a micro-scaling version of the transversal dynamics of the entire urban system. Finally, we propose a mathematical framework that connects the microscopic level to global behavior, and, in all analyzed cases, we find good agreement between theoretical prediction and empirical evidence.


title: A Comparison Study on Nonlinear Dimension Reduction Methods with Kernel Variations: Visualization, Optimization and Classification

abstract: Because of high dimensionality, correlation among covariates, and noise contained in data, dimension reduction (DR) techniques are often employed to the application of machine learning algorithms. Principal Component Analysis (PCA), Linear Discriminant Analysis (LDA), and their kernel variants (KPCA, KLDA) are among the most popular DR methods. Recently, Supervised Kernel Principal Component Analysis (SKPCA) has been shown as another successful alternative. In this paper, brief reviews of these popular techniques are presented first. We then conduct a comparative performance study based on three simulated datasets, after which the performance of the techniques are evaluated through application to a pattern recognition problem in face image analysis. The gender classification problem is considered on MORPH-II and FG-NET, two popular longitudinal face aging databases. Several feature extraction methods are used, including biologically-inspired features (BIF), local binary patterns (LBP), histogram of oriented gradients (HOG), and the Active Appearance Model (AAM). After applications of DR methods, a linear support vector machine (SVM) is deployed with gender classification accuracy rates exceeding 95% on MORPH-II, competitive with benchmark results. A parallel computational approach is also proposed, attaining faster processing speeds and similar recognition rates on MORPH-II. Our computational approach can be applied to practical gender classification systems and generalized to other face analysis tasks, such as race classification and age prediction.


title: Privacy Preserving Stochastic Channel-Based Federated Learning with Neural Network Pruning

abstract: Artificial neural network has achieved unprecedented success in a wide variety of domains such as classifying, predicting and recognizing objects. This success depends on the availability of big data since the training process requires massive and representative data sets. However, data collection is often prevented by privacy concerns and people want to take control over their sensitive information during both training and using processes. To address this problem, we propose a privacy-preserving method for the distributed system, Stochastic Channel-Based Federated Learning (SCBF), which enables the participants to train a high-performance model cooperatively without sharing their inputs. We design, implement and evaluate a channel-based update algorithm for the central server in a distributed system, which selects the channels with regard to the most active features in a training loop and uploads them as learned information from local datasets. A pruning process is applied to the algorithm based on the validation set, which serves as a model accelerator. In the experiment, our model presents equal performances and higher saturating speed than the Federated Averaging method which reveals all the parameters of local models to the server when updating. We also demonstrate that the converging rates could be increased by introducing a pruning process.


title: The Bayesian Inversion Problem for Thermal Average Sampling of Quantum Systems

abstract: In this article, we propose a novel method for sampling potential functions based on noisy observation data of a finite number of observables in quantum canonical ensembles, which leads to the accurate sampling of a wide class of test observables. The method is based on the Bayesian inversion framework, which provides a platform for analyzing the posterior distribution and naturally leads to an efficient numerical sampling algorithm. We highlight that, the stability estimate is obtained by treating the potential functions as intermediate variables in the following way: the discrepancy between two sets of observation data of training observables can bound the distance between corresponding posterior distributions of potential functions, while the latter naturally leads to a bound of the discrepancies between corresponding thermal averages of test observables. Besides, the training observables can be more flexible than finite samples of the local density function, which are mostly used in previous researches. The method also applies to the multi-level quantum systems in the non-adiabatic regime. In addition, we provide extensive numerical tests to verify the accuracy and efficiency of the proposed algorithm.


title: Commensurability of Baumslag-Solitar groups

abstract: In this paper we classify Baumslag-Solitar groups up to commensurability. In order to prove our main result we give a solution to the isomorphism problem for a subclass of Generalised Baumslag-Solitar groups.


title: Stack Pressure Considerations for Room Temperature All-Solid-State Lithium Metal Batteries

abstract: All-solid-state batteries are expected to enable batteries with high energy density with the use of lithium metal anodes. Although solid electrolytes are believed to be mechanically strong enough to prevent lithium dendrites from propagating, various reports today still show cell failure due to lithium dendritic growth at room temperature. While cell parameters such as current density, electrolyte porosity and interfacial properties have been investigated, mechanical properties of lithium metal and the role of applied stack pressure on the shorting behavior is still poorly understood. Here, we investigated failure mechanisms of lithium metal in all-solid-state batteries as a function of stack pressure, and conducted in situ characterization of the interfacial and morphological properties of the buried lithium in solid electrolytes. We found that a low stack pressure of 5 MPa allows reliable plating and stripping in a lithium symmetric cell for more than 1000 hours, and a Li | Li6PS5Cl | LiNi0.80Co0.15Al0.05O2 full cell, plating more than 4 um of lithium per charge, is able to cycle over 200 cycles at room temperature. These results suggest the possibility of enabling the lithium metal anode in all-solid-state batteries at reasonable stack pressures.


title: AKM$^2$D : An Adaptive Framework for Online Sensing and Anomaly Quantification

abstract: In point-based sensing systems such as coordinate measuring machines (CMM) and laser ultrasonics where complete sensing is impractical due to the high sensing time and cost, adaptive sensing through a systematic exploration is vital for online inspection and anomaly quantification. Most of the existing sequential sampling methodologies focus on reducing the overall fitting error for the entire sampling space. However, in many anomaly quantification applications, the main goal is to estimate sparse anomalous regions in the pixel-level accurately. In this paper, we develop a novel framework named Adaptive Kernelized Maximum-Minimum Distance AKM$^2$D to speed up the inspection and anomaly detection process through an intelligent sequential sampling scheme integrated with fast estimation and detection. The proposed method balances the sampling efforts between the space-filling sampling (exploration) and focused sampling near the anomalous region (exploitation). The proposed methodology is validated by conducting simulations and a case study of anomaly detection in composite sheets using a guided wave test.


title: Distributed Learning of Deep Neural Networks using Independent Subnet Training

abstract: Distributed machine learning (ML) can bring more computational resources to bear than single-machine learning, thus enabling reductions in training time. Distributed learning partitions models and data over many machines, allowing model and dataset sizes beyond the available compute power and memory of a single machine. In practice though, distributed ML is challenging when distribution is mandatory, rather than chosen by the practitioner. In such scenarios, data could unavoidably be separated among workers due to limited memory capacity per worker or even because of data privacy issues. There, existing distributed methods will utterly fail due to dominant transfer costs across workers, or do not even apply. We propose a new approach to distributed fully connected neural network learning, called independent subnet training (IST), to handle these cases. In IST, the original network is decomposed into a set of narrow subnetworks with the same depth. These subnetworks are then trained locally before parameters are exchanged to produce new subnets and the training cycle repeats. Such a naturally "model parallel" approach limits memory usage by storing only a portion of network parameters on each device. Additionally, no requirements exist for sharing data between workers (i.e., subnet training is local and independent) and communication volume and frequency are reduced by decomposing the original network into independent subnets. These properties of IST can cope with issues due to distributed data, slow interconnects, or limited device memory, making IST a suitable approach for cases of mandatory distribution. We show experimentally that IST results in training times that are much lower than common distributed learning approaches.


title: Particle-like ultracompact objects in Einstein-scalar-Gauss-Bonnet theories

abstract: We present a new type of ultracompact objects, featuring lightrings and echoes in the gravitational-wave spectrum. These particle-like solutions arise in Einstein-scalar-Gauss-Bonnet theories in four spacetime dimensions, representing globally regular spacetime manifolds. The scalar field diverges at the center, but the effective stress-energy tensor is free from pathologies. We determine their domain of existence and compare with wormhole solutions, black holes and the Fisher solution.


title: Upper bounds on polynomials with small Galois group

abstract: When monic integral polynomials of degree $n \geq 2$ are ordered by the maximum of the absolute value of their coefficients, the Hilbert irreducibility theorem implies that asymptotically 100% are irreducible and have Galois group isomorphic to $S_n$. In particular, amongst such polynomials whose coefficients are bounded by $B$ in absolute value, asymptotically $(1+o(1))(2B+1)^n$ are irreducible and have Galois group $S_n$. When $G$ is a proper transitive subgroup of $S_n$, however, the asymptotic count of polynomials with Galois group $G$ has been determined only in very few cases. Here, we show that if there are strong upper bounds on the number of degree $n$ fields with Galois group $G$, then there are also strong bounds on the number of polynomials with Galois group $G$. For example, for any prime $p$, we show that there are at most $O(B^{3 - \frac{2}{p}} (\log B)^{p - 1})$ polynomials with Galois group $C_p$ and coefficients bounded by $B$.


title: Maximum Matchings in Geometric Intersection Graphs

abstract: Let $G$ be an intersection graph of $n$ geometric objects in the plane. We show that a maximum matching in $G$ can be found in $O(\rho^{3\omega/2}n^{\omega/2})$ time with high probability, where $\rho$ is the density of the geometric objects and $\omega>2$ is a constant such that $n \times n$ matrices can be multiplied in $O(n^\omega)$ time. The same result holds for any subgraph of $G$, as long as a geometric representation is at hand. For this, we combine algebraic methods, namely computing the rank of a matrix via Gaussian elimination, with the fact that geometric intersection graphs have small separators. We also show that in many interesting cases, the maximum matching problem in a general geometric intersection graph can be reduced to the case of bounded density. In particular, a maximum matching in the intersection graph of any family of translates of a convex object in the plane can be found in $O(n^{\omega/2})$ time with high probability, and a maximum matching in the intersection graph of a family of planar disks with radii in $[1, \Psi]$ can be found in $O(\Psi^6\log^{11} n + \Psi^{12 \omega} n^{\omega/2})$ time with high probability.


title: LRP2020: Probing Diverse Phenomena through Data-Intensive Astronomy

abstract: The era of data-intensive astronomy is being ushered in with the increasing size and complexity of observational data across wavelength and time domains, the development of algorithms to extract information from this complexity, and the computational power to apply these algorithms to the growing repositories of data. Data-intensive approaches are pushing the boundaries of nearly all fields of astronomy, from exoplanet science to cosmology, and they are becoming a critical modality for how we understand the universe. The success of these approaches range from the discovery of rare or unexpected phenomena, to characterizing processes that are now accessible with precision astrophysics and a deep statistical understanding of the datasets, to developing algorithms that maximize the science that can be extracted from any set of observations. In this white paper, we propose a number of initiatives to maximize Canada's ability to compete in this data-intensive era. We propose joining international collaborations and leveraging Canadian facilities for legacy data potential. We propose continuing to build a more agile computing infrastructure that's responsive to the needs of tackling larger and more complex data, as well as enabling quick prototyping and scaling of algorithms. We recognize that developing the fundamental skills of the field will be critical for Canadian astronomers, and discuss avenues through with the appropriate computational and statistical training could occur. Finally, we note that the transition to data-intensive techniques is not limited to astronomy, and we should coordinate with other disciplines to develop and make use of best practises in methods, infrastructure, and education.


title: Requirements for Developing Robust Neural Networks

abstract: Validation accuracy is a necessary, but not sufficient, measure of a neural network classifier's quality. High validation accuracy during development does not guarantee that a model is free of serious flaws, such as vulnerability to adversarial attacks or a tendency to misclassify (with high confidence) data it was not trained on. The model may also be incomprehensible to a human or base its decisions on unreasonable criteria. These problems, which are not unique to classifiers, have been the focus of a substantial amount of recent research. However, they are not prioritized during model development, which almost always optimizes on validation accuracy to the exclusion of everything else. The product of this approach is likely to fail in unexpected ways outside of the training environment. We believe that, in addition to validation accuracy, the model development process must give added weight to other performance metrics such as explainability, resistance to adversarial attacks, and overconfidence on out-of-distribution data.


title: Quantum Physical Unclonable Functions: Possibilities and Impossibilities

abstract: A Physical Unclonable Function (PUF) is a device with unique behaviour that is hard to clone hence providing a secure fingerprint. A variety of PUF structures and PUF-based applications have been explored theoretically as well as being implemented in practical settings. Recently, the inherent unclonability of quantum states has been exploited to derive the quantum analogue of PUF as well as new proposals for the implementation of PUF. We present the first comprehensive study of quantum Physical Unclonable Functions (qPUFs) with quantum cryptographic tools. We formally define qPUFs, encapsulating all requirements of classical PUFs as well as introducing a new testability feature inherent to the quantum setting only. We use a quantum game-based framework to define different levels of security for qPUFs: quantum exponential unforgeability, quantum existential unforgeability and quantum selective unforgeability. We introduce a new quantum attack technique based on the universal quantum emulator algorithm of Marvin and Lloyd to prove no qPUF can provide quantum existential unforgeability. On the other hand, we prove that a large family of qPUFs (called unitary PUFs) can provide quantum selective unforgeability which is the desired level of security for most PUF-based applications.


title: Modeling the Comb Filter Effect and Interaural Coherence for Binaural Source Separation

abstract: Typical methods for binaural source separation consider only the direct sound as the target signal in a mixture. However, in most scenarios, this assumption limits the source separation performance. It is well known that the early reflections interact with the direct sound, producing acoustic effects at the listening position, e.g. the so-called comb filter effect. In this article, we propose a novel source separation model, that utilizes both the direct sound and the first early reflection information to model the comb filter effect. This is done by observing the interaural phase difference obtained from the time-frequency representation of binaural mixtures. Furthermore, a method is proposed to model the interaural coherence of the signals. Including information related to the sound multipath propagation, the performance of the proposed separation method is improved with respect to the baselines that did not use such information, as illustrated by using binaural recordings made in four rooms, having different sizes and reverberation times.


title: Optical turbulence at Ali, China -- Results from the first year of lunar scintillometer observations

abstract: The location of an astronomical observatory is a key factor that affects its scientific productivity. The best astronomical sites are generally those found at high altitudes. Several such sites in western China and the Tibetan plateau are presently under development for astronomy. One of these is Ali, which at over 5000 m is one of the highest astronomical sites in the world. In order to further investigate the astronomical potential of Ali, we have installed a lunar scintillometer, for the primary purpose of profiling atmospheric turbulence. This paper describes the instrument and technique, and reports results from the first year of observations. We find that ground layer (GL) turbulence at Ali is remarkably weak and relatively thin. The median seeing, from turbulence in the range 11- 500 m above ground is 0.34 arcsec, with seeing better than 0.26 arcsec occurring 25 per cent of the time. Under median conditions, half of the GL turbulence lies below a height of 62 m. These initial results, and the high altitude and relatively low temperatures, suggest that Ali could prove to be an outstanding site for ground-based astronomy.


title: Do non-dipolar magnetic fields contribute to spin-down torques?

abstract: Main sequence low-mass stars are known to spin-down as a consequence of their magnetised stellar winds. However, estimating the precise rate of this spin-down is an open problem. The mass-loss rate, angular momentum-loss rate and the magnetic field properties of low-mass stars are fundamentally linked making this a challenging task. Of particular interest is the stellar magnetic field geometry. In this work, we consider whether non-dipolar field modes contribute significantly to the spin-down of low-mass stars. We do this using a sample of stars that have all been previously mapped with Zeeman-Doppler imaging. For a given star, as long as its mass-loss rate is below some critical mass-loss rate, only the dipolar fields contribute to its spin-down torque. However, if it has a larger mass-loss rate, higher order modes need to be considered. For each star, we calculate this critical mass-loss rate, which is a simple function of the field geometry. Additionally, we use two methods of estimating mass-loss rates for our sample of stars. In the majority of cases, we find that the estimated mass-loss rates do not exceed the critical mass-loss rate and hence, the dipolar magnetic field alone is sufficient to determine the spin-down torque. However, we find some evidence that, at large Rossby numbers, non-dipolar modes may start to contribute.


title: Online Active Perception for Partially Observable Markov Decision Processes with Limited Budget

abstract: Active perception strategies enable an agent to selectively gather information in a way to improve its performance. In applications in which the agent does not have prior knowledge about the available information sources, it is crucial to synthesize active perception strategies at runtime. We consider a setting in which at runtime an agent is capable of gathering information under a limited budget. We pose the problem in the context of partially observable Markov decision processes. We propose a generalized greedy strategy that selects a subset of information sources with near-optimality guarantees on uncertainty reduction. Our theoretical analysis establishes that the proposed active perception strategy achieves near-optimal performance in terms of expected cumulative reward. We demonstrate the resulting strategies in simulations on a robotic navigation problem.


title: A Silicon MEMS EM vibration energy harvester

abstract: This paper presents an optimized silicon-MEMS electromagnetic vibration energy harvester suitable for applications such as machine health monitoring. The harvester comprises a DRIE-etched silicon suspension, and pick-and-place N42 NdBFe magnets and copper coils, housed in a 3D-printed package. The harvester is designed to operate near 50 Hz with 0.5-1 g vibrations using a long-stroke suspension. Multi-domain harvester optimization results in an open-circuit voltage of 1.7 V, a matched-load power output of 2.2 mW, and a matched-load power-output density of 1.23 mW/cm3 at 1.1 g with a resonance frequency of 76.3 Hz.


title: How to measure CMB spectral distortions with an imaging telescope

abstract: We propose the implementation of an imaging telescope in combination with an inter-frequency calibrator to measure the spectral shape of the microwave sky by exploiting the differences in the sky intensity between multiple pairs of frequency channels. By jointly sampling the cosmological and foreground parameters in a Bayesian framework for $600$ instrument configurations, we determine the minimum calibration accuracy required in order to obtain measurements of spectral distortions and simultaneously measure spectral and spatial fluctuations of the CMB. We demonstrate the feasibility of this technique for a CMB mission like PICO (Probe of Inflation and Cosmic Origins), and show that a $10$-$\sigma$ measurement of the $y$-distortion along with a two orders of magnitude improvement on the FIRAS (Far-Infrared Absolute Spectrophotometer) $\mu$-distortion limit is feasible from this technique. We argue that longer term applications may be envisaged at even higher sensitivity, capable of attaining the $\mu$-distortion that provide a robust prediction of the $\Lambda$CDM model and even primordial recombination lines of hydrogen and helium, in the context of the ESA (European Space Agency) Voyage 2035-2050 program.


title: A Conditional Generative Model for Predicting Material Microstructures from Processing Methods

abstract: Microstructures of a material form the bridge linking processing conditions - which can be controlled, to the material property - which is the primary interest in engineering applications. Thus a critical task in material design is establishing the processing-structure relationship, which requires domain expertise and techniques that can model the high-dimensional material microstructure. This work proposes a deep learning based approach that models the processing-structure relationship as a conditional image synthesis problem. In particular, we develop an auxiliary classifier Wasserstein GAN with gradient penalty (ACWGAN-GP) to synthesize microstructures under a given processing condition. This approach is free of feature engineering, requires modest domain knowledge and is applicable to a wide range of material systems. We demonstrate this approach using the ultra high carbon steel (UHCS) database, where each microstructure is annotated with a label describing the cooling method it was subjected to. Our results show that ACWGAN-GP can synthesize high-quality multiphase microstructures for a given cooling method.


title: Eraser morphisms and membership problem in groups and monoids

abstract: We develop the theory of fragile words by introducing the concept of eraser morphism and extending the concept to more general contexts such as (free) inverse monoids. We characterize the image of the eraser morphism in the free group case, and show that it has decidable membership problem. We establish several algorithmic properties of the class of finite-${\cal{J}}$-above (inverse) monoids. We prove that the image of the eraser morphism in the free inverse monoid case (and more generally, in the finite-${\cal{J}}$-above case) has decidable membership problem, and relate its kernel to the free group fragile words.


title: Energy localization in interacting atomic chains with topological solitons

abstract: Topological defects in low-dimensional non-linear systems feature a sliding-to-pinning transition of relevance for a variety of research fields, ranging from biophysics to nano- and solid-state physics. We find that the dynamics after a local excitation results in a highly-non-trivial energy transport in the presence of a topological soliton, characterized by a strongly enhanced energy localization in the pinning regime. Moreover, we show that the energy flux in ion crystals with a topological defect can be sensitively regulated by experimentally accessible environmental parameters. Whereas, third-order non-linear resonances can cause an enhanced long-time energy delocalization, robust energy localization persists for distinct parameter ranges even for long evolution times and large local excitations.


title: Risks of Using Non-verified Open Data: A case study on using Machine Learning techniques for predicting Pregnancy Outcomes in India

abstract: Artificial intelligence (AI) has evolved considerably in the last few years. While applications of AI is now becoming more common in fields like retail and marketing, application of AI in solving problems related to developing countries is still an emerging topic. Specially, AI applications in resource-poor settings remains relatively nascent. There is a huge scope of AI being used in such settings. For example, researchers have started exploring AI applications to reduce poverty and deliver a broad range of critical public services. However, despite many promising use cases, there are many dataset related challenges that one has to overcome in such projects. These challenges often take the form of missing data, incorrectly collected data and improperly labeled variables, among other factors. As a result, we can often end up using data that is not representative of the problem we are trying to solve. In this case study, we explore the challenges of using such an open dataset from India, to predict an important health outcome. We highlight how the use of AI without proper understanding of reporting metrics can lead to erroneous conclusions.


title: Microfoundations of Discounting

abstract: An important question in economics is how people choose between different payments in the future. The classical normative model predicts that a decision maker discounts a later payment relative to an earlier one by an exponential function of the time between them. Descriptive models use non-exponential functions to fit observed behavioral phenomena, such as preference reversal. Here we propose a model of discounting, consistent with standard axioms of choice, in which decision makers maximize the growth rate of their wealth. Four specifications of the model produce four forms of discounting -- no discounting, exponential, hyperbolic, and a hybrid of exponential and hyperbolic -- two of which predict preference reversal. Our model requires no assumption of behavioral bias or payment risk.


title: A Method of EV Detour-to-Recharge Behavior Modeling and Charging Station Deployment

abstract: Electric vehicles (EVs) are increasingly used in transportation. Worldwide use of EVs, for their limited battery capacity, calls for effective planning of EVs charging stations to enhance the efficiency of using EVs. This paper provides a methodology of describing EV detouring behavior for recharging, and based on this, we adopt the extra driving length caused by detouring and the length of uncompleted route as the indicators of evaluating an EV charging station deployment plan. In this way, we can simulate EV behavior based on travel data (demand). Then, a genetic algorithm (GA) based EV charging station sitting optimization method is developed to obtain an effective plan. A detailed case study based on a 100-node 203-branch transportation network within a 30 km * 30 km region is included to test the effectiveness of our method. Insights from our method may be applicable for charging station planning in various transportation networks.


title: HDMI-Walk: Attacking HDMI Distribution Networks via Consumer Electronic Control Protocol

abstract: The High Definition Multimedia Interface (HDMI) is the de-facto standard for Audio/Video interfacing between video-enabled devices. Today, almost tens of billions of HDMI devices exist worldwide and are widely used to distribute A/V signals in smart homes, offices, concert halls, and sporting events making HDMI one of the most highly deployed systems in the world. An important component in HDMI is the Consumer Electronics Control (CEC) protocol, which allows for the interaction between devices within an HDMI distribution network. Nonetheless, existing network security mechanisms only protect traditional networking components, leaving CEC outside of their scope. In this work, we identify and tap into CEC protocol vulnerabilities, using them to implement realistic proof-of-work attacks on HDMI distribution networks. We study, how current insecure CEC protocol practices and HDMI distributions may grant an adversary a novel attack surface for HDMI devices otherwise thought to be unreachable. To introduce this novel attack surface, we present HDMI-Walk, which opens a realm of remote and local CEC attacks to HDMI devices. Specifically, with HDMI-Walk, an attacker can perform malicious analysis of devices, eavesdropping, Denial of Service attacks, targeted device attacks, and even facilitate well-known existing attacks through HDMI. With HDMI-Walk, we prove it is feasible for an attacker to gain arbitrary control of HDMI devices. We demonstrate the implementations of both local and remote attacks with commodity HDMI devices. Finally, we discuss security mechanisms to provide impactful and comprehensive security evaluation to these real-world systems while guaranteeing deployability and providing minimal overhead considering the current limitations of the CEC protocol. To the best of our knowledge, this is the first work solely investigating the security of HDMI device distribution networks.


title: Discounted Reinforcement Learning Is Not an Optimization Problem

abstract: Discounted reinforcement learning is fundamentally incompatible with function approximation for control in continuing tasks. It is not an optimization problem in its usual formulation, so when using function approximation there is no optimal policy. We substantiate these claims, then go on to address some misconceptions about discounting and its connection to the average reward formulation. We encourage researchers to adopt rigorous optimization approaches, such as maximizing average reward, for reinforcement learning in continuing tasks.


title: An Efficient and Accurate Method for Modeling Nonlinear Fractional Viscoelastic Biomaterials

abstract: Computational biomechanics plays an important role in biomedical engineering: using modeling to understand pathophysiology, treatment and device design. While experimental evidence indicates that the mechanical response of most tissues is viscoelasticity, current biomechanical models in the computation community often assume only hyperelasticity. Fractional viscoelastic constitutive models have been successfully used in literature to capture the material response. However, the translation of these models into computational platforms remains limited. Many experimentally derived viscoelastic constitutive models are not suitable for three-dimensional simulations. Furthermore, the use of fractional derivatives can be computationally prohibitive, with a number of current numerical approximations having a computational cost that is $ \mathcal{O} ( N_T^2) $ and a storage cost that is $ \mathcal{O}(N_T) $ ($ N_T $ denotes the number of time steps). In this paper, we present a novel numerical approximation to the Caputo derivative which exploits a recurrence relation similar to those used to discretize classic temporal derivatives, giving a computational cost that is $ \mathcal{O} (N) $ and a storage cost that is fixed over time. The approximation is optimized for numerical applications, and the error estimate is presented to demonstrate efficacy of the method. The method is shown to be unconditionally stable in the linear viscoelastic case. It was then integrated into a computational biomechanical framework, with several numerical examples verifying accuracy and computational efficiency of the method, including in an analytic test, in an analytic fractional differential equation, as well as in a computational biomechanical model problem.


title: Lipschitz Learning for Signal Recovery

abstract: We consider the recovery of signals from their observations, which are samples of a transform of the signals rather than the signals themselves, by using machine learning (ML). We will develop a theoretical framework to characterize the signals that can be robustly recovered from their observations by an ML algorithm, and establish a Lipschitz condition on signals and observations that is both necessary and sufficient for the existence of a robust recovery. We will compare the Lipschitz condition with the well-known restricted isometry property of the sparse recovery of compressive sensing, and show the former is more general and less restrictive. For linear observations, our work also suggests an ML method in which the output space is reduced to the lowest possible dimension.


title: A DECam View of the Diffuse Dwarf Galaxy Crater II: The Colour-Magnitude Diagram

abstract: We present a deep Blanco/DECam colour-magnitude diagram (CMD) for the large but very diffuse Milky Way satellite dwarf galaxy Crater II. The CMD shows only old stars with a clearly bifurcated subgiant branch (SGB) that feeds a narrow red giant branch. The horizontal branch (HB) shows many RR Lyrae and red HB stars. Comparing the CMD with [Fe/H] = -2.0 and [$\alpha$/Fe] = +0.3 alpha-enhanced BaSTI isochrones indicates a mean age of 12.5 Gyr for the main event and a mean age of 10.5 Gyr for the brighter SGB. With such multiple star formation events Crater II shows similarity to more massive dwarfs that have intermediate age populations, however for Crater II there was early quenching of the star formation and no intermediate age or younger stars are present. The spatial distribution of Crater II stars overall is elliptical in the plane of the sky, the detailed distribution shows a lack of strong central concentration, and some inhomogeneities. The 10.5 Gyr subgiant and upper main sequence stars show a slightly higher central concentration when compared to the 12.5 Gyr population. Matching to Gaia DR2 we find the proper motion of Crater II: $\mu_{\alpha}\cos \delta$=-0.14 $\pm$ 0.07 , $\mu_{\delta}$=-0.10 $\pm$ 0.04 mas yr$^{-1}$, approximately perpendicular to the semi-major axis of Crater II. Our results provide constraints on the star formation and chemical enrichment history of Crater II, but cannot definitively determine whether or not substantial mass has been lost over its lifetime.


title: Concepts, Components and Collections of Trading Strategies and Market Color

abstract: This paper acts as a collection of various trading strategies and useful pieces of market information that might help to implement such strategies. This list is meant to be comprehensive (though by no means exhaustive) and hence we only provide pointers and give further sources to explore each strategy further. To set the stage for this exploration, we consider the factors that determine good and bad trades, the notions of market efficiency, the real prospect amidst the seemingly high expectations of homogeneous expectations from human beings and the catch-22 connotations that arise while comprehending the true meaning of rational investing. We can broadly classify trading ideas and client market color material into Delta-One and Derivative strategies since this acts as a natural categorization that depends on the expertise of the various trading desks that will implement these strategies. For each strategy, we will have a core idea and we will present different flavors of this central theme to demonstrate that we can easily cater to the varying risk appetites, regional preferences, asset management styles, investment philosophies, liability constraints, investment horizons, notional trading size, trading frequency and other preferences of different market participants.


title: Types for Parallel Complexity in the Pi-calculus

abstract: Type systems as a way to control or analyze programs have been largely studied in the context of functional programming languages. Some of those work allow to extract from a typing derivation for a program a complexity bound on this program. We present how to adapt this result for parallel complexity in the pi-calculus, as a model of concurrency and parallel communication. We study two notions of time complexity: the total computation time without parallelism (the work) and the computation time under maximal parallelism (the span). We define reduction relations in the pi-calculus to capture those two notions, and we present two type systems from which one can extract a complexity bound on a process. The type systems are inspired by input/output types and size types, with temporal information about communications.


title: RecordFlux: Formal Message Specification and Generation of Verifiable Binary Parsers

abstract: Various vulnerabilities have been found in message parsers of protocol implementations in the past. Even highly sensitive software components like TLS libraries are affected regularly. Resulting issues range from denial-of-service attacks to the extraction of sensitive information. The complexity of protocols and imprecise specifications in natural language are the core reasons for subtle bugs in implementations, which are hard to find. The lack of precise specifications impedes formal verification. In this paper, we propose a model and a corresponding domain-specific language to formally specify message formats of existing real-world binary protocols. A unique feature of the model is the capability to define invariants, which specify relations and dependencies between message fields. Furthermore, the model allows defining the relation of messages between different protocol layers and thus ensures correct interpretation of payload data. We present a technique to derive verifiable parsers based on the model, generate efficient code for their implementation, and automatically prove the absence of runtime errors. Examples of parser specifications for Ethernet and TLS demonstrate the applicability of our approach.


title: Distance, de Sitter and Trans-Planckian Censorship conjectures: the status quo of Warm Inflation

abstract: In this short note, we pointed out that Warm inflationary scenario remains to be favoured over its cold counterpart by the recently proposed conjectures, which aim to overrule de Sitter like constructions in String Landscapes. On the other hand, the canonical cold (single-field slow-roll) inflationary models, which were in tension with the previously proposed de Sitter conjecture, have now become even more unlikely to realise in String Landscapes as these scenarios fail to cope up with both the conjectures, de Sitter and trans-Planckian censorship, at one go.


title: Involutive latin solutions of the Yang-Baxter equation

abstract: Wolfgang Rump showed that there is a one-to-one correspondence between nondegenerate involutive set-theoretic solutions of the Yang-Baxter equation and binary algebras in which all left translations $L_x$ are bijections, the squaring map is a bijection, and the identity $(xy)(xz) = (yx)(yz)$ holds. We call these algebras \emph{rumples} in analogy with quandles, another class of binary algebras giving solutions of the Yang-Baxter equation. We focus on latin rumples, that is, on rumples in which all right translations are bijections as well. We prove that an affine latin rumple of order $n$ exists if and only if $n=p_1^{p_1 k_1}\cdots p_m^{p_m k_m}$ for some distinct primes $p_i$ and positive integers $k_i$. A large class of affine solutions is obtained from nonsingular near-circulant matrices $A$, $B$ satisfying $[A,B]=A^2$. We characterize affine latin rumples as those latin rumples for which the displacement group generated by $L_x L_y\inv$ is abelian and normal in the group generated by all translations. We develop the extension theory of rumples sufficiently to obtain examples of latin rumples that are not affine, not even isotopic to a group. Finally, we investigate latin rumples in which the dual identity $(zx)(yx) = (zy)(xy)$ holds as well, and we show, among other results, that the generators $L_x L_y\inv$ of their displacement group have order dividing four.


title: Effect of nucleation sites on the growth and quality of single-crystal boron arsenide

abstract: Boron arsenide (BAs) has been the least investigated cubic III-V compound, but it has recently attracted significant attention since the confirmation of its unusually high thermal conductivity above 1000 W/m-K. However, determining how to achieve growth of a BAs single crystal on the centimeter scale remains unsolved, which strongly limits further research into, and potential applications of, this interesting material. Here we report our technique to grow a 7-mm-long BAs single crystal via the chemical vapor transport method by applying an additional nucleation site. The different thermal conductivity values obtained from BAs single crystals grown on nucleation sites of different compositions show the importance of choosing the proper nucleation-site material. We believe these findings will inspire further research into the growth of this unique semiconductor.


title: Tensor-based algorithms for image classification

abstract: The interest in machine learning with tensor networks has been growing rapidly in recent years. We show that tensor-based methods developed for learning the governing equations of dynamical systems from data can, in the same way, be used for supervised learning problems and propose two novel approaches for image classification. One is a kernel-based reformulation of the previously introduced MANDy (multidimensional approximation of nonlinear dynamics), the other an alternating ridge regression in the tensor-train format. We apply both methods to the MNIST and fashion MNIST data set and show that the approaches are competitive with state-of-the-art neural network-based classifiers.


title: Towards a Definitive Compressibility Measure for Repetitive Sequences

abstract: Unlike in statistical compression, where Shannon's entropy is a definitive lower bound, no such clear measure exists for the compressibility of repetitive sequences. Since statistical entropy does not capture repetitiveness, ad-hoc measures like the size $z$ of the Lempel--Ziv parse are frequently used to estimate it. The size $b \le z$ of the smallest bidirectional macro scheme captures better what can be achieved via copy-paste processes, though it is NP-complete to compute and it is not monotonic upon symbol appends. Recently, a more principled measure, the size $\gamma$ of the smallest string \emph{attractor}, was introduced. The measure $\gamma \le b$ lower bounds all the previous relevant ones, yet length-$n$ strings can be represented and efficiently indexed within space $O(\gamma\log\frac{n}{\gamma})$, which also upper bounds most measures. While $\gamma$ is certainly a better measure of repetitiveness than $b$, it is also NP-complete to compute and not monotonic, and it is unknown if one can always represent a string in $o(\gamma\log n)$ space. In this paper, we study an even smaller measure, $\delta \le \gamma$, which can be computed in linear time, is monotonic, and allows encoding every string in $O(\delta\log\frac{n}{\delta})$ space because $z = O(\delta\log\frac{n}{\delta})$. We show that $\delta$ better captures the compressibility of repetitive strings. Concretely, we show that (1) $\delta$ can be strictly smaller than $\gamma$, by up to a logarithmic factor; (2) there are string families needing $\Omega(\delta\log\frac{n}{\delta})$ space to be encoded, so this space is optimal for every $n$ and $\delta$; (3) one can build run-length context-free grammars of size $O(\delta\log\frac{n}{\delta})$, whereas the smallest (non-run-length) grammar can be up to $\Theta(\log n/\log\log n)$ times larger; and (4) within $O(\delta\log\frac{n}{\delta})$ space we can not only...


title: Boron isotope effect on the thermal conductivity of boron arsenide single crystals

abstract: Boron arsenide (BAs) with a zinc blende structure has recently been discovered to exhibit unusual and ultrahigh thermal conductivity (k), providing a new outlook for research on BAs and other high thermal conductivity materials. Technology for BAs crystal growth has been continuously improving, however, the influence of boron isotopes, pure or mixed, on the thermal conductivity in BAs is still not completely clear. Here we report detailed studies on the growth of single crystals of BAs with different isotopic ratios and demonstrate that the k of isotopically pure BAs is at least 10% higher than that of BAs grown from natural B. Raman spectroscopy characterization shows differences in scattering among various BAs samples. The presented results will be helpful in guiding further studies on the influence of isotopes on optimizing k in BAs.


title: Solid-state laser refrigeration of a semiconductor optomechanical resonator

abstract: Photothermal heating represents a major constraint that limits the performance of many nanoscale optoelectronic and optomechanical devices including nanolasers, quantum optomechanical resonators, and integrated photonic circuits. Although radiation-pressure damping has been reported to cool an individual vibrational mode of an optomechanical resonator to its quantum ground state, to date the internal material temperature within an optomechanical resonator has not been reported to cool via laser excitation. Here we demonstrate the direct laser refrigeration of a semiconductor optomechanical resonator >20K below room temperature based on the emission of upconverted, anti-Stokes photoluminescence of trivalent ytterbium ions doped within a yttrium-lithium-fluoride (YLF) host crystal. Optically-refrigerating the lattice of a dielectric resonator has the potential to impact several fields including scanning probe microscopy, the sensing of weak forces, the measurement of atomic masses, and the development of radiation-balanced solid-state lasers. In addition, optically refrigerated resonators may be used in the future as a promising starting point to perform motional cooling for exploration of quantum effects at mesoscopic length scales,temperature control within integrated photonic devices, and solid-state laser refrigeration of quantum materials


title: Local rigidity of manifolds with hyperbolic cusps II. Nonlinear theory

abstract: This article is the second in a series of two whose aim is to extend a recent result of Guillarmou-Lefeuvre [arXiv:1806.04218] on the local rigidity of the marked length spectrum from the case of compact negatively-curved Riemannian manifolds to the case of manifolds with hyperbolic cusps. We deal with the nonlinear version of the problem and prove that such manifolds are locally rigid for nonlinear perturbations of the metric that slightly decrease at infinity. Our proof relies on the linear theory addressed in [arXiv:1907.01809] and on a careful analytic study of the generalized X-ray transform operator $\Pi_2$. In particular, we prove that the latter fits in the microlocal theory for cusp manifolds developed in [arXiv:1411.5083, arXiv:1712.07832, arXiv:1907.01809].


title: The Sparse Reverse of Principal Component Analysis for Fast Low-Rank Matrix Completion

abstract: Matrix completion constantly receives tremendous attention from many research fields. It is commonly applied for recommender systems such as movie ratings, computer vision such as image reconstruction or completion, multi-task learning such as collaboratively modeling time-series trends of multiple sensors, and many other applications. Matrix completion techniques are usually computationally exhaustive and/or fail to capture the heterogeneity in the data. For example, images usually contain a heterogeneous set of objects, and thus it is a challenging task to reconstruct images with high levels of missing data. In this paper, we propose the sparse reverse of principal component analysis for matrix completion. The proposed approach maintains smoothness across the matrix, produces accurate estimates of the missing data, converges iteratively, and it is computationally tractable with a controllable upper bound on the number of iterations until convergence. The accuracy of the proposed technique is validated on natural images, movie ratings, and multisensor data. It is also compared with common benchmark methods used for matrix completion.


title: Continuous Comagnetometry using Transversely Polarized Xe Isotopes

abstract: We demonstrate a transversely polarized spin-exchange pumped noble gas comagnetometer which suppresses systematic errors from longitudinal polarization. Rb atoms as well as $^{131}$Xe and $^{129}$Xe nuclei are simultaneously polarized perpendicular to a pulsed bias field. Both Xe isotopes' nuclear magnetic resonance conditions are simultaneously satisfied by frequency modulation of the pulse repetition rate. The Rb atoms detect the Xe precession. We highlight the importance of magnetometer phase shifts when performing comagnetometry. For detection of non-magnetic spin-dependent interactions the sensing bandwidth is 1 Hz, the white-noise level is 7 $\mu$Hz /$\sqrt{\text{Hz}}$, and the bias stability is $\approx1$ $\mu$Hz.


title: Energy Resource Control via Privacy Preserving Data

abstract: Although the frequent monitoring of smart meters enables granular control over energy resources, it also increases the risk of leakage of private information such as income, home occupancy, and power consumption behavior that can be inferred from the data by an adversary. We propose a method of releasing modified smart meter data so specific private attributes are obscured while the utility of the data for use in an energy resource controller is preserved. The method achieves privatization by injecting noise conditional on the private attribute through a linear filter learned via a minimax optimization. The optimization contains the loss function of a classifier for the private attribute, which we maximize, and the energy resource controller's objective formulated as a canonical form optimization, which we minimize. We perform our experiment on a dataset of household consumption with solar generation and another from the Commission for Energy Regulation that contains household smart meter data with sensitive attributes such as income and home occupancy. We demonstrate that our method is able to significantly reduce the ability of an adversary to classify the private attribute while maintaining a similar objective value for an energy storage controller.


title: RDMA vs. RPC for Implementing Distributed Data Structures

abstract: Distributed data structures are key to implementing scalable applications for scientific simulations and data analysis. In this paper we look at two implementation styles for distributed data structures: remote direct memory access (RDMA) and remote procedure call (RPC). We focus on operations that require individual accesses to remote portions of a distributed data structure, e.g., accessing a hash table bucket or distributed queue, rather than global operations in which all processors collectively exchange information. We look at the trade-offs between the two styles through microbenchmarks and a performance model that approximates the cost of each. The RDMA operations have direct hardware support in the network and therefore lower latency and overhead, while the RPC operations are more expressive but higher cost and can suffer from lack of attentiveness from the remote side. We also run experiments to compare the real-world performance of RDMA- and RPC-based data structure operations with the predicted performance to evaluate the accuracy of our model, and show that while the model does not always precisely predict running time, it allows us to choose the best implementation in the examples shown. We believe this analysis will assist developers in designing data structures that will perform well on current network architectures, as well as network architects in providing better support for this class of distributed data structures.


title: On distinct consecutive differences

abstract: We show that if $A=\{a_1 < a_2 < \ldots < a_k\}$ is a set of real numbers such that the differences of the consecutive elements are distinct, then for and finite $B \subset \mathbb{R}$, $$|A+B|\gg |A|^{1/2}|B|.$$ The bound is tight up to the constant.


title: Variable Selection with Random Survival Forest and Bayesian Additive Regression Tree for Survival Data

abstract: In this paper we utilize a survival analysis methodology incorporating Bayesian additive regression trees to account for nonlinear and additive covariate effects. We compare the performance of Bayesian additive regression trees, Cox proportional hazards and random survival forests models for censored survival data, using simulation studies and survival analysis for breast cancer with U.S. SEER database for the year 2005. In simulation studies, we compare the three models across varying sample sizes and censoring rates on the basis of bias and prediction accuracy. In survival analysis for breast cancer, we retrospectively analyze a subset of 1500 patients having invasive ductal carcinoma that is a common form of breast cancer mostly affecting older woman. Predictive potential of the three models are then compared using some widely used performance assessment measures in survival literature.


title: Existence of traveling wave solutions of a deterministic vector-host epidemic model with direct transmission

abstract: We consider an epidemic model with direct transmission given by a system of nonlinear partial differential equations and study the existence of traveling wave solutions. When the basic reproductive number of the considered model is less than one, we show that there is no nontrivial traveling wave solution. On the other hand, when the basic reproductive number is greater than one, we prove that there is a minimum wave speed $c^*$ such that the system has a traveling wave solution with speed $c$ connecting both equilibrium points for any $c\ge c^*$. Moreover, under suitable assumption on the diffusion rates, we show that there is no traveling wave solution with speed less than $c^*$. We conclude with numerical simulations to illustrate our findings. The numerical experiments supports the validity of our theoretical results.


title: Measuring the Impact Ionization and Charge Trapping Probabilities in SuperCDMS HVeV Phonon Sensing Detectors

abstract: A 0.93 gram $1{\times}1{\times}0.4$ cm$^3$ SuperCDMS silicon HVeV detector operated at 30 mK was illuminated by 1.91 eV photons using a room temperature pulsed laser coupled to the cryostat via fiber optic. The detector's response under a variety of specific operating conditions was used to study the detector leakage current, charge trapping and impact ionization in the high-purity Si substrate. The measured probabilities for a charge carrier in the detector to undergo charge trapping (0.713 $\pm$ 0.093%) or cause impact ionization (1.576 $\pm$ 0.110%) were found to be nearly independent of bias polarity and charge-carrier type (electron or hole) for substrate biases of $\pm$ 140 V.


title: 4f Conduction in the Magnetic Semiconductor NdN

abstract: We report the growth of films of the intrinsic ferromagnetic semiconductor NdN, and an investigation of their optical and transport properties. There is clear evidence of a strong anomalous Hall effect as expected from a 4f conduction channel, supported by an optical absorption into a 4f or 4f/5d hybridized tail at the base of the conduction band. The results reveal a heavy-fermion 4f/5d band lying where it can be occupied at controllable levels with nitrogen-vacancy donors.


title: Unambiguous separators for tropical tree automata

abstract: In this paper we show that given a max-plus automaton (over trees, and with real weights) computing a function $f$ and a min-plus automaton (similar) computing a function $g$ such that $f\leqslant g$, there exists effectively an unambiguous tropical automaton computing $h$ such that $f\leqslant h\leqslant g$. This generalizes a result of Lombardy and Mairesse of 2006 stating that series which are both max-plus and min-plus rational are unambiguous. This generalization goes in two directions: trees are considered instead of words, and separation is established instead of characterization (separation implies characterization). The techniques in the two proofs are very different.


title: SLAM-based Integrity Monitoring Using GPS and Fish-eye Camera

abstract: Urban navigation using GPS and fish-eye camera suffers from multipath effects in GPS measurements and data association errors in pixel intensities across image frames. We propose a Simultaneous Localization and Mapping (SLAM)-based Integrity Monitoring (IM) algorithm to compute the position protection levels while accounting for multiple faults in both GPS and vision. We perform graph optimization using the sequential data of GPS pseudoranges, pixel intensities, vehicle dynamics, and satellite ephemeris to simultaneously localize the vehicle as well as the landmarks, namely GPS satellites and key image pixels in the world frame. We estimate the fault mode vector by analyzing the temporal correlation across the GPS measurement residuals and spatial correlation across the vision intensity residuals. In particular, to detect and isolate the vision faults, we developed a superpixel-based piecewise Random Sample Consensus (RANSAC) technique to perform spatial voting across image pixels. For an estimated fault mode, we compute the protection levels by applying worst-case failure slope analysis to the linearized Graph-SLAM framework. We perform ground vehicle experiments in the semi-urban area of Champaign, IL and have demonstrated the successful detection and isolation of multiple faults. We also validate tighter protection levels and lower localization errors achieved via the proposed algorithm as compared to SLAM-based IM that utilizes only GPS measurements.


title: Reconstruction and Measurement of $\mathcal{O}$(100) MeV Energy Electromagnetic Activity from $\pi^0 \rightarrow \gamma\gamma$ Decays in the MicroBooNE LArTPC

abstract: We present results on the reconstruction of electromagnetic (EM) activity from photons produced in charged current $\nu_{\mu}$ interactions with final state $\pi^0$s. We employ a fully-automated reconstruction chain capable of identifying EM showers of $\mathcal{O}$(100) MeV energy, relying on a combination of traditional reconstruction techniques together with novel machine-learning approaches. These studies demonstrate good energy resolution, and good agreement between data and simulation, relying on the reconstructed invariant $\pi^0$ mass and other photon distributions for validation. The reconstruction techniques developed are applied to a selection of $\nu_{\mu} + {\rm Ar} \rightarrow \mu + \pi^0 + X$ candidate events to demonstrate the potential for calorimetric separation of photons from electrons and reconstruction of $\pi^0$ kinematics.


title: A characteristic map for the holonomy groupoid of a foliation

abstract: We prove a generalisation of Bott's vanishing theorem for the full transverse frame holonomy groupoid of any transversely orientable foliated manifold. As a consequence we obtain a characteristic map encoding both primary and secondary characteristic classes. Previous descriptions of this characteristic map are formulated for the Morita equivalent \'{e}tale groupoid obtained via a choice of complete transversal. By working with the full holonomy groupoid we obtain novel geometric representatives of characteristic classes. In particular we give a geometric, non-\'{e}tale analogue of the codimension 1 Godbillon-Vey cyclic cocycle of Connes and Moscovici in terms of path integrals of the curvature form of a Bott connection.


title: Cross lingual transfer learning for zero-resource domain adaptation

abstract: We propose a method for zero-resource domain adaptation of DNN acoustic models, for use in low-resource situations where the only in-language training data available may be poorly matched to the intended target domain. Our method uses a multi-lingual model in which several DNN layers are shared between languages. This architecture enables domain adaptation transforms learned for one well-resourced language to be applied to an entirely different low-resource language. First, to develop the technique we use English as a well-resourced language and take Spanish to mimic a low-resource language. Experiments in domain adaptation between the conversational telephone speech (CTS) domain and broadcast news (BN) domain demonstrate a 29% relative WER improvement on Spanish BN test data by using only English adaptation data. Second, we demonstrate the effectiveness of the method for low-resource languages with a poor match to the well-resourced language. Even in this scenario, the proposed method achieves relative WER improvements of 18-27% by using solely English data for domain adaptation. Compared to other related approaches based on multi-task and multi-condition training, the proposed method is able to better exploit well-resource language data for improved acoustic modelling of the low-resource target domain.


title: High-field Spatial Imaging of Charge Transport in Silicon at Low Temperature

abstract: We present direct imaging measurements of charge transport across a 1 cm x 1 cm x 4 mm-thick crystal of high purity silicon ($\sim$15 k$\Omega$-cm) at temperatures of 5 K and 500 mK. We use these data to measure the lateral diffusion of electrons and holes as a function of the electric field applied along the [111] crystal axis, and to verify our low-temperature Monte Carlo software. The range of field strengths in this paper exceed those used in the previous study (DOI: 10.1063/1.5049691) by a factor of 10, and now encompasses the region in which some recent silicon dark matter detectors operate (DOI: 10.1103/PhysRevLett.121.051301). We also report on a phenomenon of surface charge trapping which can reduce expected charge collection.


title: Donor's Deferral and Return Behavior: Partial Identification from a Regression Discontinuity Design with Manipulation

abstract: Volunteer labor can temporarily yield lower benefits to charities than its costs. In such instances, organizations may wish to defer volunteer donations to a later date. Exploiting a discontinuity in blood donations' eligibility criteria, we show that deferring donors reduces their future volunteerism. In our setting, medical staff manipulates donors' reported hemoglobin levels over a threshold to facilitate donation. Such manipulation invalidates standard regression discontinuity design. To circumvent this issue, we propose a procedure for obtaining partial identification bounds where manipulation is present. Our procedure is applicable in various regression discontinuity settings where the running variable is manipulated and discrete.


title: Light- and strange-quark mass dependence of the $\rho(770)$ meson properties

abstract: From an analysis of recent ($I=J=1$)-$\pi\pi$-phase-shift and pseudoscalar-meson decay-constant lattice data on two distinct chiral trajectories, where either the sum of the up, down and strange quark masses, or the mass of the strange quark is kept fixed, we extract the light and strange quark mass dependence of the rho meson parameters, and make predictions of those on chiral trajectories which involve lighter masses than the physical strange quark mass. We find that the mass of the rho meson can get as light as $700$ MeV for strange quark mass zero at physical pion masses. While the ratio of the couplings to the $\pi\pi$ and $K\bar{K}$ channels is equal to $\sqrt{2}$ at the SU(3) symmetric chiral trajectory.


title: Development of a multiwire proportional chamber with good tolerance to burst hits

abstract: The DeeMe experiment to search for muon-to-electron conversions with a sensitivity 10--100 times better than those achieved by previous experiments is in preparation at the Japan Proton Accelerator Research Complex. The magnetic spectrometer used by the DeeMe experiment consists of an electromagnet and four multiwire proportional chambers (MWPCs). The newly developed MWPCs are operated with a high voltage (HV) switching technique and have good burst-hit tolerance. In this article, the final designs of the MWPCs, amplifiers for readout, and HV switching modules are described. Additionally, some results of MWPC performance evaluation are presented.


title: HSC-XD 52: An X-ray detected AGN in a low-mass galaxy at $z\sim0.56$

abstract: The properties of low-mass galaxies hosting central black holes provide clues about the formation and evolution of the progenitors of supermassive black holes. In this letter, we present HSC-XD 52, a spectroscopically confirmed low-mass active galactic nucleus (AGN) at an intermediate redshift of $z\sim0.56$. We detect this object as a very luminous X-ray source coincident with a galaxy observed by the Hyper Suprime-Cam (HSC) as part of a broader search for low-mass AGN. We constrain its stellar mass through spectral energy distribution modeling to be LMC-like at $M_\star \approx 3 \times 10^9 M_\odot$, placing it in the dwarf regime. We estimate a central black hole mass of $M_\mathrm{BH} \sim 10^{6} M_\odot$. With an average X-ray luminosity of $L_X \approx 3.5 \times 10^{43}~\mathrm{erg}~\mathrm{s}^{-1}$, HSC-XD 52 is among the most luminous X-ray selected AGN in dwarf galaxies. The spectroscopic and photometric properties of HSC-XD 52 indicate that it is an intermediate redshift counterpart to local low-mass AGN.


title: Quantifying conjugacy separability in wreath products of groups

abstract: We study generalisations of conjugacy separability in restricted wreath products of groups. We provide an effective upper bound for $\mathcal{C}$-conjugacy separability of a wreath product $A \wr B$ in terms of the $\mathcal{C}$-conjugacy separability of $A$ and $B$, the growth of $\mathcal{C}$-cyclic subgroup separability of $B$, and the $\mathcal{C}$-residual girth of $B.$ As an application, we provide a characterisation of when $A \wr B$ is $p$-conjugacy separable. We use this characterisation to the provide for each prime $p$ an example of wreath products with infinite base group that are $p$-conjugacy separable. We also provide asymptotic upper bounds for conjugacy separability for wreath products of nilpotent groups which include the lamplighter groups and provide asymptotic upper bounds for conjugacy separability of the free metabelian groups. Along the way, we provide a polynomial upper bound for the shortest conjugator between two elements of length at most $n$ in a finitely generated nilpotent group.


title: Pi-PE: A Pipeline for Pulmonary Embolism Detection using Sparsely Annotated 3D CT Images

abstract: Pulmonary embolisms (PE) are known to be one of the leading causes for cardiac-related mortality. Due to inherent variabilities in how PE manifests and the cumbersome nature of manual diagnosis, there is growing interest in leveraging AI tools for detecting PE. In this paper, we build a two-stage detection pipeline that is accurate, computationally efficient, robust to variations in PE types and kernels used for CT reconstruction, and most importantly, does not require dense annotations. Given the challenges in acquiring expert annotations in large-scale datasets, our approach produces state-of-the-art results with very sparse emboli contours (at 10mm slice spacing), while using models with significantly lower number of parameters. We achieve AUC scores of 0.94 on the validation set and 0.85 on the test set of highly severe PEs. Using a large, real-world dataset characterized by complex PE types and patients from multiple hospitals, we present an elaborate empirical study and provide guidelines for designing highly generalizable pipelines.


title: Straight-Through Estimator as Projected Wasserstein Gradient Flow

abstract: The Straight-Through (ST) estimator is a widely used technique for back-propagating gradients through discrete random variables. However, this effective method lacks theoretical justification. In this paper, we show that ST can be interpreted as the simulation of the projected Wasserstein gradient flow (pWGF). Based on this understanding, a theoretical foundation is established to justify the convergence properties of ST. Further, another pWGF estimator variant is proposed, which exhibits superior performance on distributions with infinite support,e.g., Poisson distributions. Empirically, we show that ST and our proposed estimator, while applied to different types of discrete structures (including both Bernoulli and Poisson latent variables), exhibit comparable or even better performances relative to other state-of-the-art methods. Our results uncover the origin of the widespread adoption of the ST estimator and represent a helpful step towards exploring alternative gradient estimators for discrete variables.


title: Modelling Quantum Devices and the Reconstruction of Physics in Practical Systems

abstract: Modelling quantum devices is to find a model according to quantum theory that can explain the result of experiments in a quantum device. We find that usually we cannot correctly identify the model describing the actual physics of the device regardless of the experimental effort given a limited set of operations. According to sufficient conditions that we find, correctly reconstructing the model requires either a particular set of pure states and projective measurements or a set of evolution operators that can generate all unitary operators.


title: The general dual-polar Orlicz-Minkowski problem

abstract: This paper gives a systematic study to the general dual-polar Orlicz-Minkowski problem (e.g., Problem \ref{general-dual-polar}). This problem involves the general dual volume $\widetilde{V}_G(\cdot)$ recently proposed in \cite{GHWXY, GHXY} in order to study the general dual Orlicz-Minkowski problem. As $\widetilde{V}_G(\cdot)$ extends the volume and the $q$th dual volume, the general dual-polar Orlicz-Minkowski problem is "polar" to the recently initiated general dual Orlicz-Minkowski problem in \cite{GHWXY, GHXY} and "dual" to the newly proposed polar Orlicz-Minkowski problem in \cite{LuoYeZhu}. The existence, continuity and uniqueness, if applicable, for the solutions to the general dual-polar Orlicz-Minkowski problem are established. Polytopal solutions and/or counterexamples to the general dual-polar Orlicz-Minkowski problem for discrete measures are also provided. Several variations of the general dual-polar Orlicz-Minkowski problem are discussed as well, in particular the one leading to the general Orlicz-Petty bodies.


title: Template-based Minor Embedding for Adiabatic Quantum Optimization

abstract: Quantum Annealing (QA) can be used to quickly obtain near-optimal solutions for Quadratic Unconstrained Binary Optimization (QUBO) problems. In QA hardware, each decision variable of a QUBO should be mapped to one or more adjacent qubits in such a way that pairs of variables defining a quadratic term in the objective function are mapped to some pair of adjacent qubits. However, qubits have limited connectivity in existing QA hardware. This has spurred work on preprocessing algorithms for embedding the graph representing problem variables with quadratic terms into the hardware graph representing qubits adjacencies, such as the Chimera graph in hardware produced by D-Wave Systems. In this paper, we use integer linear programming to search for an embedding of the problem graph into certain classes of minors of the Chimera graph, which we call template embeddings. One of these classes corresponds to complete bipartite graphs, for which we show the limitation of the existing approach based on minimum Odd Cycle Transversals (OCTs). One of the formulations presented is exact, and thus can be used to certify the absence of a minor embedding using that template. On an extensive test set consisting of random graphs from five different classes of varying size and sparsity, we can embed more graphs than a state-of-the-art OCT-based approach, our approach scales better with the hardware size, and the runtime is generally orders of magnitude smaller.


title: Scaling Law for Cracking in Shrinkable Granular Packings

abstract: Hydrated granular packings often crack into discrete clusters of grains when dried. Despite its ubiquity, accurate prediction of cracking remains elusive. Here, we elucidate the previously overlooked role of individual grain shrinkage---a feature common to many materials---in determining crack patterning using both experiments and simulations. By extending the classical Griffith crack theory, we obtain a scaling law that quantifies how cluster size depends on the interplay between grain shrinkage, stiffness, and size---applicable to a diverse array of shrinkable, granular packings.

categories: cond-mat.soft cond-mat.dis-nn cond-mat.mes-hall cond-mat.mtrl-sci nlin.PS


title: To React or not to React: End-to-End Visual Pose Forecasting for Personalized Avatar during Dyadic Conversations

abstract: Non verbal behaviours such as gestures, facial expressions, body posture, and para-linguistic cues have been shown to complement or clarify verbal messages. Hence to improve telepresence, in form of an avatar, it is important to model these behaviours, especially in dyadic interactions. Creating such personalized avatars not only requires to model intrapersonal dynamics between a avatar's speech and their body pose, but it also needs to model interpersonal dynamics with the interlocutor present in the conversation. In this paper, we introduce a neural architecture named Dyadic Residual-Attention Model (DRAM), which integrates intrapersonal (monadic) and interpersonal (dyadic) dynamics using selective attention to generate sequences of body pose conditioned on audio and body pose of the interlocutor and audio of the human operating the avatar. We evaluate our proposed model on dyadic conversational data consisting of pose and audio of both participants, confirming the importance of adaptive attention between monadic and dyadic dynamics when predicting avatar pose. We also conduct a user study to analyze judgments of human observers. Our results confirm that the generated body pose is more natural, models intrapersonal dynamics and interpersonal dynamics better than non-adaptive monadic/dyadic models.


title: On Tractable Computation of Expected Predictions

abstract: Computing expected predictions of discriminative models is a fundamental task in machine learning that appears in many interesting applications such as fairness, handling missing values, and data analysis. Unfortunately, computing expectations of a discriminative model with respect to a probability distribution defined by an arbitrary generative model has been proven to be hard in general. In fact, the task is intractable even for simple models such as logistic regression and a naive Bayes distribution. In this paper, we identify a pair of generative and discriminative models that enables tractable computation of expectations, as well as moments of any order, of the latter with respect to the former in case of regression. Specifically, we consider expressive probabilistic circuits with certain structural constraints that support tractable probabilistic inference. Moreover, we exploit the tractable computation of high-order moments to derive an algorithm to approximate the expectations for classification scenarios in which exact computations are intractable. Our framework to compute expected predictions allows for handling of missing data during prediction time in a principled and accurate way and enables reasoning about the behavior of discriminative models. We empirically show our algorithm to consistently outperform standard imputation techniques on a variety of datasets. Finally, we illustrate how our framework can be used for exploratory data analysis.


title: The ${\cal N}=4$ Coset Model and the Higher Spin Algebra

abstract: By computing the operator product expansions between the first two ${\cal N}=4$ higher spin multiplets in the unitary coset model, the (anti)commutators of higher spin currents are obtained under the large $(N,k)$ 't Hooft-like limit. The free field realization with complex bosons and fermions is presented. The (anti)commutators for generic spins $s_1$ and $s_2$ with manifest $SO(4)$ symmetry at vanishing 't Hooft-like coupling constant are completely determined. The structure constants can be written in terms of the ones in the ${\cal N}=2$ ${\cal W}_{\infty}$ algebra found by Bergshoeff, Pope, Romans, Sezgin and Shen previously, in addition to the spin-dependent fractional coefficients and two $SO(4)$ invariant tensors. We also describe the ${\cal N}=4$ higher spin generators, by using the above coset construction results, for general super spin $s$ in terms of oscillators in the matrix generalization of $AdS_3$ Vasiliev higher spin theory at nonzero 't Hooft-like coupling constant. We obtain the ${\cal N}=4$ higher spin algebra for low spins and present how to determine the structure constants, which depend on the higher spin algebra parameter, in general, for fixed spins $s_1$ and $s_2$.


title: Effects of niche overlap on co-existence, fixation and invasion in a population of two interacting species

abstract: Synergistic and antagonistic interactions in multi-species populations - such as resource sharing and competition - result in remarkably diverse behaviors in populations of interacting cells, such as in soil or human microbiomes, or clonal competition in cancer. The degree of inter- and intra-specific interaction can often be quantified through the notion of an ecological "niche". Typically, weakly interacting species that occupy largely distinct niches result in stable mixed populations, while strong interactions and competition for the same niche results in rapid extinctions of some species and fixations of others. We investigate the transition of a deterministically stable mixed population to a stochasticity-induced fixation as a function of the niche overlap between the two species. We also investigate the effect of the niche overlap on the population stability with respect to external invasions. Our results have important implications for a number of experimental systems.


title: Prostate cancer inference via weakly-supervised learning using a large collection of negative MRI

abstract: Recent advances in medical imaging techniques have led to significant improvements in the management of prostate cancer (PCa). In particular, multi-parametric MRI (mp-MRI) continues to gain clinical acceptance as the preferred imaging technique for non-invasive detection and grading of PCa. However, the machine learning-based diagnosis systems for PCa are often constrained by the limited access to accurate lesion ground truth annotations for training. The performance of the machine learning system is highly dependable on both quality and quantity of lesion annotations associated with histopathologic findings, resulting in limited scalability and clinical validation. Here, we propose the baseline MRI model to alternatively learn the appearance of mp-MRI using radiology-confirmed negative MRI cases via weakly supervised learning. Since PCa lesions are case-specific and highly heterogeneous, it is assumed to be challenging to synthesize PCa lesions using the baseline MRI model, while it would be relatively easier to synthesize the normal appearance in mp-MRI. We then utilize the baseline MRI model to infer the pixel-wise suspiciousness of PCa by comparing the original and synthesized MRI with two distance functions. We trained and validated the baseline MRI model using 1,145 negative prostate mp-MRI scans. For evaluation, we used separated 232 mp-MRI scans, consisting of both positive and negative MRI cases. The 116 positive MRI scans were annotated by radiologists, confirmed with post-surgical whole-gland specimens. The suspiciousness map was evaluated by receiver operating characteristic (ROC) analysis for PCa lesions versus non-PCa regions classification and free-response receiver operating characteristic (FROC) analysis for PCa localization. Our proposed method achieved 0.84 area under the ROC curve and 77.0% sensitivity at one false positive per patient in FROC analysis.


title: Quantum Spacetime Instability and Breakdown of Semiclassical Gravity

abstract: The semiclassical gravity describes gravitational back-reactions of the classical spacetime interacting with quantum matter fields but the quantum effects on the background is formally defined as higher derivative curvatures. These induce catastrophic instabilities and classic solutions become unstable under small perturbations or their evolutions. In this paper we discuss validity of the semiclassical gravity from the perspective of the spacetime instabilities and consider cosmological dynamics of the Universe in this theory. We clearly show that the homogenous and isotropic flat Universe is unstable and the solutions either grow exponentially or oscillate even in Planckian time $t_{\rm I}=(\alpha_{ 1 }G_N)^{1/2}\approx \alpha_{ 1 }10^{-43}\ {\rm sec}$. The subsequent curvature evolution leads to Planck-scale spacetime curvature in a short time and causes a catastrophe of the Universe unless one takes extremely large values of the gravitational couplings. Furthermore, we confirm the above suggestion by comparing the semiclassical solutions and $\Lambda$CDM with the Planck data and it is found that the semiclassical solutions are not consistent with the cosmological observations. Thus, the standard semiclassical gravity using quantum energy momentum tensor $\left< { T }_{ \mu \nu } \right>$ is not appropriate to describe our Universe.


title: Dynamic Embedding on Textual Networks via a Gaussian Process

abstract: Textual network embedding aims to learn low-dimensional representations of text-annotated nodes in a graph. Prior work in this area has typically focused on fixed graph structures; however, real-world networks are often dynamic. We address this challenge with a novel end-to-end node-embedding model, called Dynamic Embedding for Textual Networks with a Gaussian Process (DetGP). After training, DetGP can be applied efficiently to dynamic graphs without re-training or backpropagation. The learned representation of each node is a combination of textual and structural embeddings. Because the structure is allowed to be dynamic, our method uses the Gaussian process to take advantage of its non-parametric properties. To use both local and global graph structures, diffusion is used to model multiple hops between neighbors. The relative importance of global versus local structure for the embeddings is learned automatically. With the non-parametric nature of the Gaussian process, updating the embeddings for a changed graph structure requires only a forward pass through the learned model. Considering link prediction and node classification, experiments demonstrate the empirical effectiveness of our method compared to baseline approaches. We further show that DetGP can be straightforwardly and efficiently applied to dynamic textual networks.


title: Observation of Electric-Dipole Transitions in the Laser-Cooling Candidate Th$^-$

abstract: Despite the fact that the laser cooling method is a well-established technique to obtain ultra-cold neutral atoms and atomic cations, it has so far never been applied to atomic anions due to the lack of suitable electric-dipole transitions. Efforts of more than a decade currently has La$^-$ as the only promising candidate for laser cooling. Our previous work [Tang et al., Phys. Rev. Lett. 123, 203002(2019)] showed that Th$^-$ is also a potential candidate. Here we report on a combination of experimental and theoretical studies to determine the relevant transition frequencies, transition rates, and branching ratios in Th$^-$. The resonant frequency of the laser cooling transition is determined to be $\nu/c$ = 4118.0 (10) cm$^{-1}$. The transition rate is calculated as A=1.17x10^4 s$^{-1}$. The branching fraction to dark states is very small, 1.47x10$^{-10}$, thus this represents an ideal closed cycle for laser cooling. Since Th has zero nuclear spin, it is an excellent candidate to be used to sympathetically cool antiprotons in a Penning trap.


title: Exponential Time Differencing for the Tracer Equations Appearing in Primitive Equation Ocean Models

abstract: The tracer equations are part of the primitive equations used in ocean modeling and describe the transport of tracers, such as temperature, salinity or chemicals, in the ocean. Depending on the number of tracers considered, several equations may be added to and coupled to the dynamics system. In many relevant situations, the time-step requirements of explicit methods imposed by the transport and mixing in the vertical direction are more restrictive than those for the horizontal, and this may cause the need to use very small time steps if a fully explicit method is employed. To overcome this issue, we propose an exponential time differencing (ETD) solver where the vertical terms (transport and diffusion) are treated with a matrix exponential, whereas the horizontal terms are dealt with in an explicit way. We investigate numerically the computational speed-ups that can be obtained over other semi-implicit methods, and we analyze the advantages of the method in the case of multiple tracers.


title: Kornia: an Open Source Differentiable Computer Vision Library for PyTorch

abstract: This work presents Kornia -- an open source computer vision library which consists of a set of differentiable routines and modules to solve generic computer vision problems. The package uses PyTorch as its main backend both for efficiency and to take advantage of the reverse-mode auto-differentiation to define and compute the gradient of complex functions. Inspired by OpenCV, Kornia is composed of a set of modules containing operators that can be inserted inside neural networks to train models to perform image transformations, camera calibration, epipolar geometry, and low level image processing techniques, such as filtering and edge detection that operate directly on high dimensional tensor representations. Examples of classical vision problems implemented using our framework are provided including a benchmark comparing to existing vision libraries.