title: Relating Structure and Power: Extended Version

abstract: Combinatorial games are widely used in finite model theory, constraint satisfaction, modal logic and concurrency theory to characterize logical equivalences between structures. In particular, Ehrenfeucht-Fraisse games, pebble games, and bisimulation games play a central role. We show how each of these types of games can be described in terms of an indexed family of comonads on the category of relational structures and homomorphisms. The index $k$ is a resource parameter which bounds the degree of access to the underlying structure. The coKleisli categories for these comonads can be used to give syntax-free characterizations of a wide range of important logical equivalences. Moreover, the coalgebras for these indexed comonads can be used to characterize key combinatorial parameters: tree-depth for the Ehrenfeucht-Fraisse comonad, tree-width for the pebbling comonad, and synchronization-tree depth for the modal unfolding comonad. These results pave the way for systematic connections between two major branches of the field of logic in computer science which hitherto have been almost disjoint: categorical semantics, and finite and algorithmic model theory.


title: Satellite Image Classification with Deep Learning

abstract: Satellite imagery is important for many applications including disaster response, law enforcement, and environmental monitoring. These applications require the manual identification of objects and facilities in the imagery. Because the geographic expanses to be covered are great and the analysts available to conduct the searches are few, automation is required. Yet traditional object detection and classification algorithms are too inaccurate and unreliable to solve the problem. Deep learning is a family of machine learning algorithms that have shown promise for the automation of such tasks. It has achieved success in image understanding by means of convolutional neural networks. In this paper we apply them to the problem of object and facility recognition in high-resolution, multi-spectral satellite imagery. We describe a deep learning system for classifying objects and facilities from the IARPA Functional Map of the World (fMoW) dataset into 63 different classes. The system consists of an ensemble of convolutional neural networks and additional neural networks that integrate satellite metadata with image features. It is implemented in Python using the Keras and TensorFlow deep learning libraries and runs on a Linux server with an NVIDIA Titan X graphics card. At the time of writing the system is in 2nd place in the fMoW TopCoder competition. Its total accuracy is 83%, the F1 score is 0.797, and it classifies 15 of the classes with accuracies of 95% or better.


title: Cross-Domain Few-Shot Learning by Representation Fusion

abstract: In order to quickly adapt to new data, few-shot learning aims at learning from few examples, often by using already acquired knowledge. The new data often differs from the previously seen data due to a domain shift, that is, a change of the input-target distribution. While several methods perform well on small domain shifts like new target classes with similar inputs, larger domain shifts are still challenging. Large domain shifts may result in high-level concepts that are not shared between the original and the new domain, whereas low-level concepts like edges in images might still be shared and useful. For cross-domain few-shot learning, we suggest representation fusion to unify different abstraction levels of a deep neural network into one representation. We propose Cross-domain Hebbian Ensemble Few-shot learning (CHEF), which achieves representation fusion by an ensemble of Hebbian learners acting on different layers of a deep neural network. Ablation studies show that representation fusion is a decisive factor to boost cross-domain few-shot learning. On the few-shot datasets miniImagenet and tieredImagenet with small domain shifts, CHEF is competitive with state-of-the-art methods. On cross-domain few-shot benchmark challenges with larger domain shifts, CHEF establishes novel state-of-the-art results in all categories. We further apply CHEF on a real-world cross-domain application in drug discovery. We consider a domain shift from bioactive molecules to environmental chemicals and drugs with twelve associated toxicity prediction tasks. On these tasks, that are highly relevant for computational drug discovery, CHEF significantly outperforms all its competitors. Github: https://github.com/ml-jku/chef


title: LASSR: Effective Super-Resolution Method for Plant Disease Diagnosis

abstract: The collection of high-resolution training data is crucial in building robust plant disease diagnosis systems, since such data have a significant impact on diagnostic performance. However, they are very difficult to obtain and are not always available in practice. Deep learning-based techniques, and particularly generative adversarial networks (GANs), can be applied to generate high-quality super-resolution images, but these methods often produce unexpected artifacts that can lower the diagnostic performance. In this paper, we propose a novel artifact-suppression super-resolution method that is specifically designed for diagnosing leaf disease, called Leaf Artifact-Suppression Super Resolution (LASSR). Thanks to its own artifact removal module that detects and suppresses artifacts to a considerable extent, LASSR can generate much more pleasing, high-quality images compared to the state-of-the-art ESRGAN model. Experiments based on a five-class cucumber disease (including healthy) discrimination model show that training with data generated by LASSR significantly boosts the performance on an unseen test dataset by nearly 22% compared with the baseline, and that our approach is more than 2% better than a model trained with images generated by ESRGAN.


title: The geometry of the moduli Space of non-cyclic biquadratic field extensions

abstract: In this paper, we investigate the existence of an elementary abelian closure in characteristic not $2$ for biquadratic extensions. We discover that it exists for any non-cyclic extension. We make use of it to obtain a classification for this class of extensions up to isomorphism via descent. This permits us to describe the geometry of this moduli space in group theoretic terms. We also provide two families of polynomials with two parameters that can describe any quartic extensions.


title: Weather Radar in Nepal: Opportunities and Challenges in Mountainous Region

abstract: Extreme rainfall is one of the major causes of natural hazards (for example flood, landslide, and debris flow) in the central Himalayan region, Nepal. The performance of strategies to manage these risks relies on the accuracy of quantitative rainfall estimates. Rain gauges have traditionally been used to measure the amount of rainfall at a given location. The point measurement often misrepresents the basin estimates, because of limited density and high spatial variability of rainfall fields across the Himalayas. The Department of Hydrology and Meteorology (DHM), Nepal has planned to install a network of three weather radars that cover the entire country. So far, the first weather radar has been installed in 2019 in the western region of the country. Two more radars will be added for the planned radar network in the near future covering the central and eastern regions of the country. Here we introduce the first installed weather radar in Nepal. We highlight both the opportunities and challenges with the radar observation in the mountainous regions. Radar rainfall estimates across the Himalayas are critical to issue severe weather warnings; forecast floods and landslides; and inform decision making in a broad range of sectors, including water and energy, construction, transportation, and agriculture.


title: Experimental Demonstration of Optoelectronic Equalization for Short-reach Transmission with Reservoir Computing

abstract: A receiver with shared complexity between optical and digital domains is experimentally demonstrated. Reservoir computing is used to equalize up to 4 directly-detected optically filtered spectral slices of a 32 GBd OOK signal over up to 80 km of SMF.


title: Transfer Learning and SpecAugment applied to SSVEP Based BCI Classification

abstract: Objective: We used deep convolutional neural networks (DCNNs) to classify electroencephalography (EEG) signals in a steady-state visually evoked potentials (SSVEP) based single-channel brain-computer interface (BCI), which does not require calibration on the user. Methods: EEG signals were converted to spectrograms and served as input to train DCNNs using the transfer learning technique. We also modified and applied a data augmentation method, SpecAugment, generally employed for speech recognition. Furthermore, for comparison purposes, we classified the SSVEP dataset using Support-vector machines (SVMs) and Filter Bank canonical correlation analysis (FBCCA). Results: Excluding the evaluated user's data from the fine-tuning process, we reached 82.2% mean test accuracy and 0.825 mean F1-Score on 35 subjects from an open dataset, using a small data length (0.5 s), only one electrode (Oz) and the DCNN with transfer learning, window slicing (WS) and SpecAugment's time masks. Conclusion: The DCNN results surpassed SVM and FBCCA performances, using a single electrode and a small data length. Transfer learning provided minimal accuracy change, but made training faster. SpecAugment created a small performance improvement and was successfully combined with WS, yielding higher accuracies. Significance: We present a new methodology to solve the problem of SSVEP classification using DCNNs. We also modified a speech recognition data augmentation technique and applied it to the context of BCIs. The presented methodology surpassed performances obtained with FBCCA and SVMs (more traditional SSVEP classification methods) in BCIs with small data lengths and one electrode. This type of BCI can be used to develop small and fast systems.


title: Single-Sideband Time-Modulated Phased Array With 2-bit Phased Shifters

abstract: A novel single-sideband (SSB) time-modulated technique with 2-bit phase shifters is proposed. The timemodulated module is implemented by adding periodic phase modulation to 2-bit phase shifters, which is simpler without performance loss compared to existing SSB time-modulated method. During one modulation period, four phase states (0, {\pi}/2, {\pi}, 3{\pi}/2) of 2-bit phase shifters are switched in sequence. After the modulation, the SSB time modulation is realized and the main power is distributed to the first harmonic component. The feasibility of the proposed method is verified by experiments. The undesired harmonics are efficiently suppressed. Meanwhile, 80{\deg} beam scanning range are realized through the proposed module.


title: A Lean and Highly-automated Model-Based Software Development Process Based on DO-178C/DO-331

abstract: The emergence of a global market for urban air mobility and unmanned aerial systems has attracted many startups across the world. These organizations have little training or experience in the traditional processes used in civil aviation for the development of software and electronic hardware. They are also constrained in the resources they can allocate for dedicated teams of professionals to follow these standardized processes. To fill this gap, this paper presents a custom workflow based on a subset of objectives derived from the foundational standards for safety critical software DO-178C/DO-331. The selection of objectives from the standards is based on the importance, degree of automation, and reusability of specific objectives. This custom workflow is intended to establish a lean and highly automated development life cycle resulting in higher quality software with better maintainability characteristics for research and prototype aircraft. It can also be proposed as means of compliance for software of certain applications such as unmanned aircraft systems, urban air mobility and general aviation. By producing the essential set of development and verification artifacts, the custom workflow also provides a scalable basis for potential future certification in compliance with DO-178C/DO-331. The custom workflow is demonstrated in a case study of an Autopilot Manual Disconnection System.


title: On invariant rank two vector bundles on $\mathbb{P}^2$

abstract: In this paper we characterize the rank two vector bundles on $\mathbb{P}^2$ which are invariant under the actions of the parabolic subgroups $G_p:=\mathrm{Stab}_p(\mathrm{PGL}(3))$ fixing a point in the projective plane, $G_L:=\mathrm{Stab}_L(\mathrm{PGL}(3))$ fixing a line, and when $p\in L$, the Borel subgroup $\mathbf{B} = G_p \cap G_L$ of $\mathrm{PGL}(3)$. Moreover, we prove that the geometrical configuration of the jumping locus induced by the invariance does not, on the other hand, characterize the invariance itself. Indeed, we find infinite families that are almost uniform but not almost homogeneous.


title: Machine Learning of Partial Differential Equations from Noise Data

abstract: Machine learning of partial differential equations from data is a potential breakthrough to solve the lack of physical equations in complex dynamic systems, but because numerical differentiation is ill-posed to noise data, noise has become the biggest obstacle in the application of partial differential equation identification method. To overcome this problem, we propose Frequency Domain Identification method based on Fourier transforms, which effectively eliminates the influence of noise by using the low frequency component of frequency domain data to identify partial differential equations in frequency domain. We also propose a new sparse identification criterion, which can accurately identify the terms in the equation from low signal-to-noise ratio data. Through identifying a variety of canonical equations spanning a number of scientific domains, the proposed method is proved to have high accuracy and robustness for equation structure and parameters identification for low signal-to-noise ratio data. The method provides a promising technique to discover potential partial differential equations from noisy experimental data.


title: Super duality for quantum affine algebras of type $A$

abstract: We introduce a new approach to the study of finite-dimensional representations of the quantum group of the affine Lie superalgebra $\mathrm{L}\mathfrak{gl}_{M|N}=\mathbb{C}[t,t^{-1}]\otimes\mathfrak{gl}_{M|N}$ ($M\neq N$). We explain how the representations of the quantum group of $\mathrm{L}\mathfrak{gl}_{M|N}$ are directly related to those of the quantum affine algebra of type $A$, using an exact monoidal functor called truncation. This can be viewed as an affine analogue of super duality of type $A$.


title: Approximation of Integral Fractional Laplacian and Fractional PDEs via sinc-Basis

abstract: Fueled by many applications in random processes, imaging science, geophysics, etc., fractional Laplacians have recently received significant attention. The key driving force behind the success of this operator is its ability to capture non-local effects while enforcing less smoothness on functions. In this paper, we introduce a spectral method to approximate this operator employing a sinc basis. Using our scheme, the evaluation of the operator and its application onto a vector has complexity of $\mathcal O(N\log(N))$ where $N$ is the number of unknowns. Thus, using iterative methods such as CG, we provide an efficient strategy to solve fractional partial differential equations with exterior Dirichlet conditions on arbitrary Lipschitz domains. Our implementation works in both $2d$ and $3d$. We also recover the FEM rates of convergence on benchmark problems. We further illustrate the efficiency of our approach by applying it to fractional Allen-Cahn and image denoising problems.


title: Piece-wise Matching Layer in Representation Learning for ECG Classification

abstract: This paper proposes piece-wise matching layer as a novel layer in representation learning methods for electrocardiogram (ECG) classification. Despite the remarkable performance of representation learning methods in the analysis of time series, there are still several challenges associated with these methods ranging from the complex structures of methods, the lack of generality of solutions, the need for expert knowledge, and large-scale training datasets. We introduce the piece-wise matching layer that works based on two levels to address some of the aforementioned challenges. At the first level, a set of morphological, statistical, and frequency features and comparative forms of them are computed based on each periodic part and its neighbors. At the second level, these features are modified by predefined transformation functions based on a receptive field scenario. Several scenarios of offline processing, incremental processing, fixed sliding receptive field, and event-based triggering receptive field can be implemented based on the choice of length and mechanism of indicating the receptive field. We propose dynamic time wrapping as a mechanism that indicates a receptive field based on event triggering tactics. To evaluate the performance of this method in time series analysis, we applied the proposed layer in two publicly available datasets of PhysioNet competitions in 2015 and 2017 where the input data is ECG signal. We compared the performance of our method against a variety of known tuned methods from expert knowledge, machine learning, deep learning methods, and the combination of them. The proposed approach improves the state of the art in two known completions 2015 and 2017 around 4% and 7% correspondingly while it does not rely on in advance knowledge of the classes or the possible places of arrhythmia.


title: Low rank compression in the numerical solution of the nonequilibrium Dyson equation

abstract: We propose a method to improve the computational and memory efficiency of numerical solvers for the nonequilibrium Dyson equation in the Keldysh formalism. It is based on the empirical observation that the nonequilibrium Green's functions and self energies arising in many problems of physical interest, discretized as matrices, have low rank off-diagonal blocks, and can therefore be compressed using a hierarchical low rank data structure. We describe an efficient algorithm to build this compressed representation on the fly during the course of time stepping, and use the representation to reduce the cost of computing history integrals, which is the main computational bottleneck. For systems with the hierarchical low rank property, our method reduces the computational complexity of solving the nonequilibrium Dyson equation from cubic to near quadratic, and the memory complexity from quadratic to near linear. We demonstrate the full solver for the Falicov-Kimball model exposed to a rapid ramp and Floquet driving of system parameters, and are able to increase feasible propagation times substantially. We present examples with 262144 time steps, which would require approximately five months of computing time and 2.2 TB of memory using the direct time stepping method, but can be completed in just over a day on a laptop with less than 4 GB of memory using our method. We also confirm the hierarchical low rank property for the driven Hubbard model in the weak coupling regime within the GW approximation, and in the strong coupling regime within dynamical mean-field theory.


title: Transforming Neural Network Visual Representations to Predict Human Judgments of Similarity

abstract: Deep-learning vision models have shown intriguing similarities and differences with respect to human vision. We investigate how to bring machine visual representations into better alignment with human representations. Human representations are often inferred from behavioral evidence such as the selection of an image most similar to a query image. We find that with appropriate linear transformations of deep embeddings, we can improve prediction of human binary choice on a data set of bird images from 72% at baseline to 89%. We hypothesized that deep embeddings have redundant, high (4096) dimensional representations; however, reducing the rank of these representations results in a loss of explanatory power. We hypothesized that the dilation transformation of representations explored in past research is too restrictive, and indeed we found that model explanatory power can be significantly improved with a more expressive linear transform. Most surprising and exciting, we found that, consistent with classic psychological literature, human similarity judgments are asymmetric: the similarity of X to Y is not necessarily equal to the similarity of Y to X, and allowing models to express this asymmetry improves explanatory power.


title: Representations of orthogonal and symplectic Yangians

abstract: Exteded Yangian algebras of orthogonal and symplectic types are defined by the Yang-Baxter RLL relation involving the fundamental R-matrix with $so(n)$ or $sp(2m)$ symmetry. We study representations of highest weight characterized by weight function ratios. We consider the algebra relations for the linear and the quadratic evaluations and the resulting conditions imposed on the representation weights. We present expressions of L-operators constructed on underlying Clifford and Heisenberg algebras and characterize their representations.


title: Observables of spheroidal magnetized Strange Stars

abstract: We study stable spheroidal configurations of magnetized Strange Stars using an axially symmetric metric in spherical coordinates that uses a gamma parameter to link the anisotropy in the Equation of State due to the magnetic field with the deformation of the star. The stars are composed by magnetized Strange Quark Matter described within the framework of the MIT-Bag model. Their masses, radii, eccentricity, redshift and mass quadrupole moment are computed. Results are compared with spherical Strange Stars solutions obtained with TOV equations and observational data of Strange Stars candidates. In the spheroidal model the observables depend directly on the deformation of the stars, and even though it is small, the observables strongly deviate from the corresponding spherical configurations. Thus, the highest values of the mass quadrupole moment correspond to the intermediate mass regime. These differences might allow to discriminate between models with/without magnetic field when compared with observations.


title: Batch-sequential design and heteroskedastic surrogate modeling for delta smelt conservation

abstract: Delta smelt is an endangered fish species in the San Francisco estuary that have shown an overall population decline over the past 30 years. Researchers have developed a stochastic, agent-based simulator to virtualize the system, with the goal of understanding the relative contribution of natural and anthropogenic factors suggested as playing a role in their decline. However, the input configuration space is high-dimensional, running the simulator is time-consuming, and its noisy outputs change nonlinearly in both mean and variance. Getting enough runs to effectively learn input--output dynamics requires both a nimble modeling strategy and parallel supercomputer evaluation. Recent advances in heteroskedastic Gaussian process (HetGP) surrogate modeling helps, but little is known about how to appropriately plan experiments for highly distributed simulator evaluation. We propose a batch sequential design scheme, generalizing one-at-a-time variance-based active learning for HetGP surrogates, as a means of keeping multi-core cluster nodes fully engaged with expensive runs. Our acquisition strategy is carefully engineered to favor selection of replicates which boost statistical and computational efficiencies when training surrogates to isolate signal in high noise regions. Design and modeling performance is illustrated on a range of toy examples before embarking on a large-scale smelt simulation campaign and downstream high-fidelity input sensitivity analysis.


title: Approximation of free convolutions by free infinitely divisible laws

abstract: Based on the~method of subordinating functions we prove bounds for the minimal error of approximations of $n$-fold convolutions of probability measures by free infinitely divisible probability measures.


title: CrimAnalyzer: Understanding Crime Patterns in S\~ao Paulo City

abstract: S\~ao Paulo is the largest city in South America, with high criminality rates. The number and type of crimes varies considerably around the city, assuming different patterns depending on urban and social characteristics. In this scenario, enabling tools to explore particular locations of the city is very important for domain experts to understand how urban features as to mobility, passersby behavior, and urban infrastructures can influence the quantity and type of crimes. In present work, we present CrimAnalyzer, a visualization assisted analytic tool that allows users to analyze crime behavior in specific regions of a city, providing new methodologies to identify local crime hotspots and their corresponding patterns over time. CrimAnalyzer has been developed from the demand of experts in criminology and it deals with three major challenges: i) flexibility to explore local regions and understand their crime patterns, ii) Identification of not only prevalent hotspots in terms of number of crimes but also hotspots where crimes are frequent but not in large amount, and iii) understand the dynamic of crime patterns over time. The effectiveness and usefulness of the proposed system are demonstrated by qualitative/quantitative comparisons as well as case studies involving real data and run by domain experts.


title: A Seamless Phase I/II Platform Design with a Time-To-Event Efficacy Endpoint for Potential COVID-19 Therapies

abstract: In the search for effective treatments for COVID-19, initial emphasis has been on re-purposed treatments. To maximise the chances of finding successful treatments, novel treatments that have been developed for this disease in particular, are needed. In this manuscript we describe and evaluate the statistical design of the AGILE platform, an adaptive randomized seamless Phase I/II trial platform that seeks to quickly establish a safe range of doses and investigates treatments for potential efficacy using a Bayesian sequential trial design. Both single agent and combination treatments are considered. We find that the design can identify potential treatments that are safe and efficacious reliably with small to moderate sample sizes.


title: A full-chain tube-based constitutive model for living linear polymers

abstract: We present a new strategy for introducing population balances into full-chain constitutive models of living polymers with linear chain architectures. We provide equations to describe a range of stress relaxation processes covering both unentangled systems (Rouse-like motion) and well entangled systems (reptation, contour length fluctuations, chain retraction, and constraint release). Special attention is given to the solutions that emerge when the 'breaking time' of the chain becomes fast compared to various stress relaxation processes. In these 'fast breaking' limits, we reproduce previously known results (with some corrections) and also present new results for nonlinear stress relaxation dynamics. Our analysis culminates with a fully developed constitutive model for the fast breaking regime in which stress relaxation is dominated by contour length fluctuations. Linear and nonlinear rheology predictions of the model are presented and discussed.


title: Deep Learning for Recognizing Mobile Targets in Satellite Imagery

abstract: There is an increasing demand for software that automatically detects and classifies mobile targets such as airplanes, cars, and ships in satellite imagery. Applications of such automated target recognition (ATR) software include economic forecasting, traffic planning, maritime law enforcement, and disaster response. This paper describes the extension of a convolutional neural network (CNN) for classification to a sliding window algorithm for detection. It is evaluated on mobile targets of the xView dataset, on which it achieves detection and classification accuracies higher than 95%.


title: Autotuning Search Space for Loop Transformations

abstract: One of the challenges for optimizing compilers is to predict whether applying an optimization will improve its execution speed. Programmers may override the compiler's profitability heuristic using optimization directives such as pragmas in the source code. Machine learning in the form of autotuning can assist users in finding the best optimizations for each platform. In this paper we propose a loop transformation search space that takes the form of a tree, in contrast to previous approaches that usually use vector spaces to represent loop optimization configurations. We implemented a simple autotuner exploring the search space and applied it to a selected set of PolyBench kernels. While the autotuner is capable of representing every possible sequence of loop transformations and their relations, the results motivate the use of better search strategies such as Monte Carlo tree search to find sophisticated loop transformations such as multilevel tiling.


title: Elementary construction of minimal free resolutions of the Specht ideals of shapes $(n-2,2)$ and $(d,d,1)$

abstract: For a partition $\lambda$ of $n \in \mathbb{N}$, let $I^{\rm Sp}_\lambda$ be the ideal of $R=K[x_1,\ldots,x_n]$ generated by all Specht polynomials of shape $\lambda$. We assume that ${\rm char}(K)=0$. Then $R/I^{\rm Sp}_{(n-2,2)}$ is Gorenstein, and $R/I^{\rm Sp}_{(d,d,1)}$ is a Cohen-Macaulay ring with a linear free resolution. In this paper, we construct minimal free resolutions of these rings. Berkesch Zamaere, Griffeth, and Sam had already studied minimal free resolutions of $R/I^{\rm Sp}_{(n-d,d)}$, which are also Cohen-Macaulay, using heighly advanced technique of the representation theory. However we only use the basic theory of Specht modules, and explicitly describe the differential maps.


title: Nonthermal element abundances at astrophysical shocks

abstract: The nonthermal source abundances of elements play a crucial role in the understanding of cosmic ray (CR) phenomena from a few GeV up to several tens of EeV. We present a first systematic approach to describe the change of the abundances from the thermal to the nonthermal state via diffusive shock acceleration by a temporally evolving shock. We consider hereby not only ionization states of elements contained in the ambient gas, which we allow to be time dependent due to shock heating, but also elements condensed on solid, charged dust grains which can be injected into the acceleration process as well. Our generic parametrized model is then applied to the case of particle acceleration by supernova remnants in various ISM phases, for which we use state-of-the-art computation packages to calculate the ionization states of all elements. The resulting predictions for low energy cosmic ray (LECR) source abundances are compared with the data obtained by various experiments. We obtain excellent agreement for shocks in warm ionized ISM environments, which include HII regions, if dust grains are injected into the diffusive shock acceleration process with a much higher efficiency than ions. Less dependence of the fit quality is found on the mass-to-charge ratio of ions. For neutral environments, assuming that there are shocks in the weakly ionized component, and for the hot ionized medium we obtain generally inferior fits, but except for the cold neutral medium we do not exclude them as subdominant sites of Galactic CR production. The key challenge is found to be putting the LECR abundance of pure gas phase elements like neon and the (semi-)volatile elements phosphorus, sulfur and chlorine into the right balance with silicon, calcium and elements of the iron group. Finally, a brief outlook to the potential consequences for the understanding of the CR composition at higher energies is presented.


title: Principles for Designing Computer Music Controllers

abstract: This paper will present observations on the design, artistic, and human factors of creating digital music controllers. Specific projects will be presented, and a set of design principles will be supported from those examples.


title: A Deep Learning Forecaster with Exogenous Variables for Day-Ahead Locational Marginal Price

abstract: Several approaches have been proposed to forecast day-ahead locational marginal price (daLMP) in deregulated energy markets. The rise of deep learning has motivated its use in energy price forecasts but most deep learning approaches fail to accommodate for exogenous variables, which have significant influence in the peaks and valleys of the daLMP. Accurate forecasts of the daLMP valleys are of crucial importance for power generators since one of the most important decisions they face is whether to sell power at a loss to prevent incurring in shutdown and start-up costs, or to bid at production cost and face the risk of shutting down. In this article we propose a deep learning model that incorporates both the history of daLMP and the effect of exogenous variables (e.g., forecasted load, weather data). A numerical study at the PJM independent system operator (ISO) illustrates how the proposed model outperforms traditional time series techniques while supporting risk-based analysis of shutdown decisions.


title: Minimal universal model for chaos in laser with feedback

abstract: We revisit the model of the laser with feedback and the minimal nonlinearity leading to chaos. Although the model has its origin in laser physics, with peculiarities related to the CO2 laser, it belongs to the class of the three dimensional paradigmatic nonlinear oscillator models giving chaos. The proposed model contains three key nonlinearities, two of which are of the type xy, where x and y are the fast and slow variables. The third one is of the type xz^2, where z is an intermediate feedback variable. We analytically demonstrate that it is essential for producing chaos via local or global homoclinic bifurcations. Its electronic implementation in the range of kilo Hertz region confirms its potential in describing phenomena evolving on different time scales.


title: On a question of Teissier

abstract: Using a result of Demailly and Pham on log canonical thresholds, we give an upper bound for polar invariants from a question of Teissier on hypersurface singularities. This provides a weaker alternative upper bound compared to the one conjectured by Teissier.


title: An order on circular permutations

abstract: Motivation coming from the study of affine Weyl groups, a structure of ranked poset is defined on the set of circular permutations in $S_n$ (that is, $n$-cycles). It is isomorphic to the poset of so-called admitted vectors, and to an interval in the affine symmetric group $\tilde S_n$ with the weak order. The poset is a semidistributive lattice, and the rank function, whose range is cubic in $n$, is computed by some special formula involving inversions. We prove also some links with Eulerian numbers, triangulations of an $n$-gon, and Young's lattice.


title: On the Fairness of Causal Algorithmic Recourse

abstract: Algorithmic fairness is typically studied from the perspective of predictions. Instead, here we investigate fairness from the perspective of recourse actions suggested to individuals to remedy an unfavourable classification. We propose two new fairness criteria at the group and individual level, which -- unlike prior work on equalising the average group-wise distance from the decision boundary -- explicitly account for causal relationships between features, thereby capturing downstream effects of recourse actions performed in the physical world. We explore how our criteria relate to others, such as counterfactual fairness, and show that fairness of recourse is complementary to fairness of prediction. We study theoretically and empirically how to enforce fair causal recourse by altering the classifier and perform a case study on the Adult dataset. Finally, we discuss whether fairness violations in the data generating process revealed by our criteria may be better addressed by societal interventions as opposed to constraints on the classifier.


title: Strong anisotropy of electron-phonon interaction in NbP probed by magnetoacoustic quantum oscillations

abstract: In this study, we report on the observation of de Haas-van Alphen-type quantum oscillations (QO) in the ultrasound velocity of NbP as well as `giant QO' in the ultrasound attenuation in pulsed magnetic fields. The difference of the QO amplitude for different acoustic modes reveals a strong anisotropy of the effective deformation potential, which we estimate to be as high as $9\,\mathrm{eV}$ for certain parts of the Fermi surface. Furthermore, the natural filtering of QO frequencies and the tracing of the individual Landau levels to the quantum limit allows for a more detailed investigation of the Fermi surface of NbP as was previously achieved by means of analyzing QO observed in magnetization or electrical resistivity.


title: Impact of Representation Learning in Linear Bandits

abstract: We study how representation learning can improve the efficiency of bandit problems. We study the setting where we play $T$ linear bandits with dimension $d$ concurrently, and these $T$ bandit tasks share a common $k (\ll d)$ dimensional linear representation. For the finite-action setting, we present a new algorithm which achieves $\widetilde{O}(T\sqrt{kN} + \sqrt{dkNT})$ regret, where $N$ is the number of rounds we play for each bandit. When $T$ is sufficiently large, our algorithm significantly outperforms the naive algorithm (playing $T$ bandits independently) that achieves $\widetilde{O}(T\sqrt{d N})$ regret. We also provide an $\Omega(T\sqrt{kN} + \sqrt{dkNT})$ regret lower bound, showing that our algorithm is minimax-optimal up to poly-logarithmic factors. Furthermore, we extend our algorithm to the infinite-action setting and obtain a corresponding regret bound which demonstrates the benefit of representation learning in certain regimes. We also present experiments on synthetic and real-world data to illustrate our theoretical findings and demonstrate the effectiveness of our proposed algorithms.


title: Transport properties of band engineered p-n heterostructures of epitaxial Bi$_2$Se$_3$/(Bi$_{1-x}$Sb$_x$)$_2$(Te$_{1-y}$Se$_y$)$_3$ topological insulators

abstract: The challenge of parasitic bulk doping in Bi-based 3D topological insulator materials is still omnipresent, especially when preparing samples by molecular beam epitaxy (MBE). Here, we present a heterostructure approach for epitaxial BSTS growth. A thin n-type Bi$_2$Se$_3$ (BS) layer is used as an epitaxial and electrostatic seed which drastically improves the crystalline and electronic quality and reproducibility of the sample properties. In heterostructures of BS with p-type (Bi$_{1-x}$Sb$_x$)$_2$(Te$_{1-y}$Se$_y$)$_3$ (BSTS) we demonstrate intrinsic band bending effects to tune the electronic properties solely by adjusting the thickness of the respective layer. The analysis of weak anti-localization features in the magnetoconductance indicates a separation of top and bottom conduction layers with increasing BSTS thickness. By temperature- and gate-dependent transport measurements, we show that the thin BS seed layer can be completely depleted within the heterostructure and demonstrate electrostatic tuning of the bands via a back-gate throughout the whole sample thickness.


title: Generic aspects of the resource theory of quantum coherence

abstract: The class of incoherent operations induces a pre-order on the set of quantum pure states, defined by the possibility of converting one state into the other by transformations within the class. We prove that if two $n$-dimensional pure states are chosen independently according to the natural uniform distribution, then the probability that they are comparable vanishes as $n\rightarrow\infty$. We also study the maximal success probability of incoherent conversions and find an explicit formula for its large-$n$ asymptotic distribution. Our analysis is based on the observation that the extreme values (largest and smallest components) of a random point uniformly sampled from the unit simplex are distributed asymptotically as certain explicit homogeneous Markov chains.

categories: quant-ph cond-mat.stat-mech math-ph math.MP math.PR


title: Multipole decomposition of the general luminosity distance 'Hubble law' -- a new framework for observational cosmology

abstract: We present the luminosity distance series expansion to third order in redshift for a general space-time with no assumption on the metric tensor or the field equations prescribing it. It turns out that the coefficients of this general 'Hubble law' can be expressed in terms of a finite number of physically interpretable multipole coefficients. The multipole terms can be combined into effective direction dependent parameters replacing the Hubble constant, deceleration parameter, curvature parameter, and 'jerk' parameter of the Friedmann-Lema\^{\i}tre-Robertson-Walker (FLRW) class of metrics. Due to the finite number of multipole coefficients, the exact anisotropic Hubble law is given by 9, 25, 61 degrees of freedom in the $\mathcal{O}(z)$, $\mathcal{O}(z^2)$, $\mathcal{O}(z^3)$ vicinity of the observer respectively, where $z\!:=\,$redshift. This makes possible model independent determination of dynamical degrees of freedom of the cosmic neighbourhood of the observer and direct testing of the FLRW ansatz. We argue that the derived multipole representation of the general Hubble law provides a new framework with broad applications in observational cosmology.


title: Two Candidate High-Redshift X-ray Jets Without Coincident Radio Jets

abstract: We report the detection of extended X-ray emission from two high-redshift radio quasars. These quasars, J1405+0415 at $z$=3.208 and J1610+1811 at $z$=3.118, were observed in a Chandra snapshot survey selected from a complete sample of the radio-brightest quasars in the overlap area of the VLA-FIRST radio survey and the Sloan Digital Sky Survey. The extended X-ray emission is located along the line connecting the core to a radio knot or hotspot, favoring the interpretation of X-ray jets. The inferred rest frame jet X-ray luminosities from 2--30 keV would be of order 10$^{45}$ erg~s$^{-1}$ if emitted isotropically and without relativistic beaming. In the scenario of inverse Compton scattering of the cosmic microwave background (CMB), X-ray jets without a coincident radio counterpart may be common, and should be readily detectable to redshifts even beyond 3.2 due to the (1+$z$)$^4$ increase of the CMB energy density compensating for the (1+$z$)$^{-4}$ cosmological diminution of surface brightness. If these can be X-ray confirmed, they would be the second and third examples of quasar X-ray jets without detection of underlying continuous radio jets.


title: Kartta Labs: Collaborative Time Travel

abstract: We introduce the modular and scalable design of Kartta Labs, an open source, open data, and scalable system for virtually reconstructing cities from historical maps and photos. Kartta Labs relies on crowdsourcing and artificial intelligence consisting of two major modules: Maps and 3D models. Each module, in turn, consists of sub-modules that enable the system to reconstruct a city from historical maps and photos. The result is a spatiotemporal reference that can be used to integrate various collected data (curated, sensed, or crowdsourced) for research, education, and entertainment purposes. The system empowers the users to experience collaborative time travel such that they work together to reconstruct the past and experience it on an open source and open data platform.


title: Can Federated Learning Save The Planet?

abstract: Despite impressive results, deep learning-based technologies also raise severe privacy and environmental concerns induced by the training procedure often conducted in data centers. In response, alternatives to centralized training such as Federated Learning (FL) have emerged. Perhaps unexpectedly, FL, in particular, is starting to be deployed at a global scale by companies that must adhere to new legal demands and policies originating from governments and the civil society for privacy protection. However, the potential environmental impact related to FL remains unclear and unexplored. This paper offers the first-ever systematic study of the carbon footprint of FL. First, we propose a rigorous model to quantify the carbon footprint, hence facilitating the investigation of the relationship between FL design and carbon emissions. Then, we compare the carbon footprint of FL to traditional centralized learning. Our findings show FL, despite being slower to converge, can be a greener technology than data center GPUs. Finally, we highlight and connect the reported results to the future challenges and trends in FL to reduce its environmental impact, including algorithms efficiency, hardware capabilities, and stronger industry transparency.


title: Modeling Atmospheric Data and Identifying Dynamics: Temporal Data-Driven Modeling of Air Pollutants

abstract: Atmospheric modeling has recently experienced a surge with the advent of deep learning. Most of these models, however, predict concentrations of pollutants following a data-driven approach in which the physical laws that govern their behaviors and relationships remain hidden. With the aid of real-world air quality data collected hourly in different stations throughout Madrid, we present an empirical approach using data-driven techniques with the following goals: (1) Find parsimonious systems of ordinary differential equations via sparse identification of nonlinear dynamics (SINDy) that model the concentration of pollutants and their changes over time; (2) assess the performance and limitations of our models using stability analysis; (3) reconstruct the time series of chemical pollutants not measured in certain stations using delay coordinate embedding results. Our results show that Akaike's Information Criterion can work well in conjunction with best subset regression as to find an equilibrium between sparsity and goodness of fit. We also find that, due to the complexity of the chemical system under study, identifying the dynamics of this system over longer periods of time require higher levels of data filtering and smoothing. Stability analysis for the reconstructed ordinary differential equations (ODEs) reveals that more than half of the physically relevant critical points are saddle points, suggesting that the system is unstable even under the idealized assumption that all environmental conditions are constant over time.


title: Tensor lattice field theory with applications to the renormalization group and quantum computing

abstract: We discuss the successes and limitations of statistical sampling for a sequence of models studied in the context of lattice QCD and emphasize the need for new methods to deal with finite-density and real-time evolution. We show that these lattice models can be reformulated using tensorial methods where the field integrations in the path-integral formalism are replaced by discrete sums. These formulations involve various types of duality and provide exact coarse-graining formulas which can be combined with truncations to obtain practical implementations of the Wilson renormalization group program. Tensor reformulations are naturally discrete and provide manageable transfer matrices. Combining truncations with the time continuum limit, we derive Hamiltonians suitable to perform quantum simulation experiments, for instance using cold atoms, or to be programmed on existing quantum computers. We review recent progress concerning the tensor field theory treatment of non-compact scalar models, supersymmetric models, economical four-dimensional algorithms, noise-robust enforcement of Gauss's law, symmetry preserving truncations and topological considerations. We discuss connections with other tensor network approaches.

categories: hep-lat cond-mat.stat-mech hep-th quant-ph


title: Structure-preserving algorithms with uniform error bound and long-time energy conservation for highly oscillatory Hamiltonian systems

abstract: Structure-preserving algorithms and algorithms with uniform error bound have constituted two interesting classes of numerical methods. In this paper, we blend these two kinds of methods for solving nonlinear Hamiltonian systems with highly oscillatory solution, and the blended algorithms inherit and respect the advantage of each method. Two kinds of algorithms are presented to preserve the symplecticity and energy of the Hamiltonian systems, respectively. Moreover, the proposed algorithms are shown to have uniform error bound for the highly oscillatory structure. A numerical experiment is carried out to support the theoretical results established in this paper by showing the performance of the blended algorithms.


title: Out-of-equilibrium Frenkel-Kontorova model

abstract: A 1D model of interacting particles moving over a periodic substrate and in a position dependent temperature profile is considered. When the substrate and the temperature profile are spatially asymmetric a center-of-mass velocity develops, corresponding to a directed transport of the chain. This autonomous system can thus transform heath currents into motion. The model parameters can be tuned such that the particles exhibit a crossover from an ordered configuration on the substrate to a disordered one, the maximal motor effect being reached in such a disordered phase. In this case the manybody motor outperforms the single motor system, showing the great importance of collective effects in microscopic thermal devices. Such collective effects represent thus a free resource that can be exploited to enhance the dynamic and thermodynamic performances in microscopic machines.


title: Analytic energy, gradient, and hessian of electrostatic embedding QM/MM based on electrostatic potential fitted atomic charges scaling linearly with the MM subsystem size

abstract: Electrostatic potential fitting method (ESPF) is a powerful way of defining atomic charges derived from quantum density matrices fitted to reproduce a quantum mechanical charge distribution in the presence of an external electrostatic potential. These can be used in the Hamiltonian to define a robust and efficient electrostatic embedding QM/MM method. The original formulation of ESPF QM/MM contained two main approximations, namely, the neglect of grid derivatives and the non-conservation of the total QM charge. Here, we present a new ESPF atomic charge operator which solves these drawbacks at virtually no extra computational cost. The new charge operators employ atom-centered grids and conserve the total charge when traced with the density matrix. We present an efficient and easy-to-implement analytic form for the energy, gradient, and hessian that scale linearly with the MM subsystem size. We show that grid derivatives and charge conservation are fundamental to preserve the translational invariance properties of energies and its derivatives and exact conditions to be satisfied by the atomic charge derivatives. As proof of concept, we compute the transition state that leads to the formation of hydrogen peroxide during cryptochrome's reoxidation reaction. Last, we show that the construction of the full QM/MM hessian scales linearly with the MM subsystem size.


title: The Hubbard model on the Bethe lattice via variational uniform tree states: metal-insulator transition and a Fermi liquid

abstract: We numerically solve the Hubbard model on the Bethe lattice with finite coordination number $z=3$, and determine its zero-temperature phase diagram. For this purpose, we introduce and develop the `variational uniform tree state' (VUTS) algorithm, a tensor network algorithm which generalizes the variational uniform matrix product state algorithm to tree tensor networks. Our results reveal an antiferromagnetic insulating phase and a paramagnetic metallic phase, separated by a first-order doping-driven metal-insulator transition. We show that the metallic state is a Fermi liquid with coherent quasiparticle excitations for all values of the interaction strength $U$, and we obtain the finite quasiparticle weight $Z$ from the single-particle occupation function of a generalized "momentum" variable. We find that $Z$ decreases with increasing $U$, ultimately saturating to a non-zero, doping-dependent value. Our work demonstrates that tensor-network calculations on tree lattices, and the VUTS algorithm in particular, are a platform for obtaining controlled results for phenomena absent in one dimension, such as Fermi liquids, while avoiding computational difficulties associated with tensor networks in two dimensions. We envision that future studies could observe non-Fermi liquids, interaction-driven metal-insulator transitions, and doped spin liquids using this platform.


title: Real-Time Deep Learning Approach to Visual Servo Control and Grasp Detection for Autonomous Robotic Manipulation

abstract: In order to explore robotic grasping in unstructured and dynamic environments, this work addresses the visual perception phase involved in the task. This phase involves the processing of visual data to obtain the location of the object to be grasped, its pose and the points at which the robot`s grippers must make contact to ensure a stable grasp. For this, the Cornell Grasping dataset is used to train a convolutional neural network that, having an image of the robot`s workspace, with a certain object, is able to predict a grasp rectangle that symbolizes the position, orientation and opening of the robot`s grippers before its closing. In addition to this network, which runs in real-time, another one is designed to deal with situations in which the object moves in the environment. Therefore, the second network is trained to perform a visual servo control, ensuring that the object remains in the robot`s field of view. This network predicts the proportional values of the linear and angular velocities that the camera must have so that the object is always in the image processed by the grasp network. The dataset used for training was automatically generated by a Kinova Gen3 manipulator. The robot is also used to evaluate the applicability in real-time and obtain practical results from the designed algorithms. Moreover, the offline results obtained through validation sets are also analyzed and discussed regarding their efficiency and processing speed. The developed controller was able to achieve a millimeter accuracy in the final position considering a target object seen for the first time. To the best of our knowledge, we have not found in the literature other works that achieve such precision with a controller learned from scratch. Thus, this work presents a new system for autonomous robotic manipulation with high processing speed and the ability to generalize to several different objects.


title: Toward Few-step Adversarial Training from a Frequency Perspective

abstract: We investigate adversarial-sample generation methods from a frequency domain perspective and extend standard $l_{\infty}$ Projected Gradient Descent (PGD) to the frequency domain. The resulting method, which we call Spectral Projected Gradient Descent (SPGD), has better success rate compared to PGD during early steps of the method. Adversarially training models using SPGD achieves greater adversarial accuracy compared to PGD when holding the number of attack steps constant. The use of SPGD can, therefore, reduce the overhead of adversarial training when utilizing adversarial generation with a smaller number of steps. However, we also prove that SPGD is equivalent to a variant of the PGD ordinarily used for the $l_{\infty}$ threat model. This PGD variant omits the sign function which is ordinarily applied to the gradient. SPGD can, therefore, be performed without explicitly transforming into the frequency domain. Finally, we visualize the perturbations SPGD generates and find they use both high and low-frequency components, which suggests that removing either high-frequency components or low-frequency components is not an effective defense.


title: On the invariance of certain types of generalized Cohen-Macaulay modules under Foxby equivalence

abstract: Let R be a local ring and C a semidualizing module of R. We investigate the behavior of certain classes of generalized Cohen-Macaulay R-modules under the Foxby equivalence between the Auslander and Bass classes with respect to C. In particular, we show that generalized Cohen-Macaulay R-modules are invariant under this equivalence and if M is a finitely generated R-module in the Auslander class with respect to C such that C\otimes_RM is surjective Buchsbaum, then M is also surjective Buchsbaum.


title: Coupled magnetic and structural phase transitions in the antiferromagnetic polar metal Pb2CoOsO6 under pressure

abstract: Pb$_2$CoOsO$_6$ is a newly synthesized polar metal in which inversion symmetry is broken by the magnetic frustration in an antiferromagnetic ordering of Co and Os sublattices. The coupled magnetic and structural transition occurs at 45 K at ambient pressure. Here we perform transport measurements and first-principles calculations to study the pressure effects on the magnetic/structural coupled transition of Pb$_2$CoOsO$_6$. Experimentally we monitor the resistivity anomaly at $T_N$ under various pressures up to 11 GPa in a cubic anvil cell apparatus. We find that $T_N$ determined from the resistivity anomaly first increases quickly with pressure in a large slope of $dT_N/dP$ = +6.8(8) K/GPa for $P < 4$ GPa, and then increases with a much reduced slope of 1.8(4) K/GPa above 4 GPa. Our first-principles calculations suggest that the observed discontinuity of $dT_N/dP$ around 4 GPa may be attributed to the vanishing of Os magnetic moment under pressure. Pressure substantially reduces the Os moment and completely suppresses it above a critical value, which relieves the magnetic frustration in the antiferromagnetic ordering of Pb$_2$CoOsO$_6$. The Co and Os polar distortions decrease with the increasing pressure and simultaneously vanish at the critical pressure. Therefore above the critical pressure a new centrosymmetric antiferromagnetic state emerges in Pb$_2$CoOsO$_6$, distinct from the one under ambient pressure, thus showing a discontinuity in $dT_N/dP$.


title: Search for dark matter produced in association with a dark Higgs boson decaying into $W^\pm W^\mp$ or $ZZ$ in fully hadronic final states from $\sqrt{s}=13$ TeV $pp$ collisions recorded with the ATLAS detector

abstract: Several extensions of the Standard Model predict the production of dark matter particles at the LHC. An uncharted signature of dark matter particles produced in association with $VV=W^\pm W^\mp$ or $ZZ$ pairs from a decay of a dark Higgs boson $s$ is searched for using 139 fb$^{-1}$ of $pp$ collisions recorded by the ATLAS detector at a center-of-mass energy of 13 TeV. The $s\to V(q\bar q)V(q\bar q)$ decays are reconstructed with a novel technique aimed at resolving the dense topology from boosted $VV$ pairs using jets in the calorimeter and tracking information. Dark Higgs scenarios with $m_s > 160$ GeV are excluded.


title: Controlling the Interaction Between Generation and Inference in Semi-Supervised Variational Autoencoders Using Importance Weighting

abstract: Even though Variational Autoencoders (VAEs) are widely used for semi-supervised learning, the reason why they work remains unclear. In fact, the addition of the unsupervised objective is most often vaguely described as a regularization. The strength of this regularization is controlled by down-weighting the objective on the unlabeled part of the training set. Through an analysis of the objective of semi-supervised VAEs, we observe that they use the posterior of the learned generative model to guide the inference model in learning the partially observed latent variable. We show that given this observation, it is possible to gain finer control on the effect of the unsupervised objective on the training procedure. Using importance weighting, we derive two novel objectives that prioritize either one of the partially observed latent variable, or the unobserved latent variable. Experiments on the IMDB english sentiment analysis dataset and on the AG News topic classification dataset show the improvements brought by our prioritization mechanism and exhibit a behavior that is inline with our description of the inner working of Semi-Supervised VAEs.


title: Connected Chord Diagrams and the Combinatorics of Asymptotic Expansions

abstract: In this article we study an asymptotic expansion for $C_n$, the number of connected chord diagrams on $n$ chords. The expansion is obtained in earlier work by means of alien derivatives applied to the generating series of connected chord diagrams; we seek a combinatorial interpretation. The main outcome presented here is a new combinatorial interpretation for entry https://oeis.org/A088221 of the OEIS. We will show that https://oeis.org/A088221 counts pairs of connected chord diagrams (allowing empty diagrams). This gives a combinatorial interpretation for the part of the closed form of the asymptotic expansion of $C_n$.


title: Transverse Measures and Best Lipschitz and Least Gradient Maps

abstract: We exhibit the duality between best Lipschitz (infinity harmonic) maps and least gradient maps in the case of maps from surfaces to the circle. We show that given a homotopy class of a map from a surface to the circle the infinity harmonic map defines a geodesic lamination on the surface and the dual least gradient map defines a transverse measure on the lamination. This is the initial step towards an analytic approach to Thurston's work on best Lipschitz maps between hyperbolic surfaces and Thurston's asymmetric metric on Teichmueller space.


title: Neural-network-driven proton decay sensitivity in the $p\rightarrow \bar{\nu} K^{+}$ channel using large liquid argon time projection chambers

abstract: We report on an updated sensitivity for proton decay via $p \rightarrow \bar{\nu} K^+ $ at large, dual phase liquid argon time projection chambers (LAr TPCs). Our work builds on a previous study in which several nucleon decay channels have been simulated and analyzed [arXiv:hep-ph/0701101]. At the time several assumptions were needed to be made on the detector and the backgrounds. Since then, the community has made progress in defining these, and the computing power available enables us to fully simulate and reconstruct large samples in order to perform a better estimate of the sensitivity to proton decay. In this work, we examine the benchmark channel $p\rightarrow \bar{\nu} K^{+}$, which was previously found to be one of the cleanest channels. Using an improved neutrino event generator and a fully simulated LAr TPC detector response combined with a dedicated neural network for kaon identification, we demonstrate that a lifetime sensitivity of $ \tau / \text{Br} \left( p \rightarrow \bar{\nu} K^+ \right) > 7 \times 10^{34} \, \text{years}$ at $90 \, \%$ confidence level can be reached at an exposure of $1 \, \text{megaton} \cdot \text{year}$ in quasi-background-free conditions, confirming the superiority of the LAr TPC over other technologies to address the challenging proton decay modes.


title: Sharp invertibility of random Bernoulli matrices

abstract: Let $p \in (0,1/2)$ be fixed, and let $B_n(p)$ be an $n\times n$ random matrix with i.i.d. Bernoulli random variables with mean $p$. We show that for all $t \ge 0$, \[\mathbb{P}[s_n(B_n(p)) \le tn^{-1/2}] \le C_p t + 2n(1-p)^{n} + C_p (1-p-\epsilon_p)^{n},\] where $s_n(B_n(p))$ denotes the least singular value of $B_n(p)$ and $C_p, \epsilon_p > 0$ are constants depending only on $p$. In particular, \[\mathbb{P}[B_{n}(p) \text{ is singular}] = 2n(1-p)^{n} + C_{p}(1-p-\epsilon_p)^{n},\] which confirms a conjecture of Litvak and Tikhomirov. We also confirm a conjecture of Nguyen by showing that if $Q_{n}$ is an $n\times n$ random matrix with independent rows that are uniformly distributed on the central slice of $\{0,1\}^{n}$, then \[\mathbb{P}[Q_{n} \text{ is singular}] = (1/2 + o_n(1))^{n}.\] This provides, for the first time, a sharp determination of the logarithm of the probability of singularity in any natural model of random discrete matrices with dependent entries.


title: Singularity of discrete random matrices

abstract: Let $\xi$ be a non-constant real-valued random variable with finite support, and let $M_{n}(\xi)$ denote an $n\times n$ random matrix with entries that are independent copies of $\xi$. For $\xi$ which is not uniform on its support, we show that \begin{align*} \mathbb{P}[M_{n}(\xi)\text{ is singular}] &= \mathbb{P}[\text{zero row or column}] + (1+o_n(1))\mathbb{P}[\text{two equal (up to sign) rows or columns}], \end{align*} thereby confirming a folklore conjecture. As special cases, we obtain: (1) For $\xi = \text{Bernoulli}(p)$ with fixed $p \in (0,1/2)$, \[\mathbb{P}[M_{n}(\xi)\text{ is singular}] = 2n(1-p)^{n} + (1+o_n(1))n(n-1)(p^2 + (1-p)^2)^{n},\] which determines the singularity probability to two asymptotic terms. Previously, no result of such precision was available in the study of the singularity of random matrices. (2) For $\xi = \text{Bernoulli}(p)$ with fixed $p \in (1/2,1)$, \[\mathbb{P}[M_{n}(\xi)\text{ is singular}] = (1+o_n(1))n(n-1)(p^2 + (1-p)^2)^{n}.\] Previously, only the much weaker upper bound of $(\sqrt{p} + o_n(1))^{n}$ was known due to the work of Bourgain-Vu-Wood. For $\xi$ which is uniform on its support: (1) We show that \begin{align*} \mathbb{P}[M_{n}(\xi)\text{ is singular}] &= (1+o_n(1))^{n}\mathbb{P}[\text{two rows or columns are equal}]. \end{align*} (2) Perhaps more importantly, we provide a sharp analysis of the contribution of the `compressible' part of the unit sphere to the lower tail of the smallest singular value of $M_{n}(\xi)$.


title: Dipolar-stabilized first and second-order antiskyrmions in ferrimagnetic multilayers

abstract: Skyrmions and antiskyrmions are topologically protected spin structures with opposite topological charge. Particularly in coexisting phases, these two types of magnetic quasi-particles may show fascinating physics and potential for spintronic devices. While skyrmions are observed in a wide range of materials, until now antiskyrmions were exclusive to materials with D2d symmetry. In this work, we show first and second-order antiskyrmions stabilized by magnetic dipole-dipole interaction in Fe/Gd-based multilayers. We modify the magnetic properties of the multilayers by Ir insertion layers. Using Lorentz transmission electron microscopy imaging, we observe coexisting antiskyrmions, Bloch skyrmions, and type-2 bubbles and determine the range of material properties and magnetic fields where the different spin objects form and dissipate. We perform micromagnetic simulations to obtain more insight into the studied system and conclude that the reduction of saturation magnetization and uniaxial anisotropy leads to the existence of this zoo of different spin objects and that they are primarily stabilized by dipolar interaction.


title: High-pressure synthesis of Ba$_2$RhO$_4$, a rhodate analogue of the layered perovskite Sr-ruthenate

abstract: A new layered perovskite-type oxide Ba$_2$RhO$_4$ was synthesized by a high-pressure technique with the support of convex-hull calculations. The crystal and electronic structure were studied by both experimental and computational tools. Structural refinements for powder x-ray diffraction data showed that Ba$_2$RhO$_4$ crystallizes in a K$_2$NiF$_4$-type structure, isostructural to Sr$_2$RuO$_4$ and Ba$_2$IrO$_4$. Magnetic, resistivity, and specific heat measurements for polycrystalline samples of Ba$_2$RhO$_4$ indicate that the system can be characterized as a correlated metal. Despite the close similarity to its Sr$_2$RuO$_4$ counterpart in the electronic specific heat coefficient and the Wilson ratio, Ba$_2$RhO$_4$ shows no signature of superconductivity down to 0.16 K. Whereas the Fermi surface topology has reminiscent pieces of Sr$_2$RuO$_4$, an electron-like e$_g$-($d_{x^2-y^2}$) band descends below the Fermi level, making of this compound unique also as a metallic counterpart of the spin-orbit-coupled Mott insulator Ba$_2$IrO$_4$.


title: Rectified phase-matching equation for fiber mode converter gratings using two-mode interference

abstract: The resonance wavelength, where a fiber mode converter grating written using periodic external perturbations achieves phase matching, is both a critical design parameter and a device parameter. However, a method to precisely predict the resonance wavelength for any new fiber and grating writing apparatus has been missing so far. The missing link was the lack of direct experimental methods to estimate the modified intermodal phase after writing with external perturbations. The presented method can make this estimation from a single experiment, over a broad wavelength range, based on a novel mathematical connection between two-mode interference and mode conversion. Using the novel methods, experimentally measured resonance wavelengths for different pitch and irradiation conditions have been predicted within relative errors of 0.004.


title: Accurate Calibration of Agent-based Epidemiological Models with Neural Network Surrogates

abstract: Calibrating complex epidemiological models to observed data is a crucial step to provide both insights into the current disease dynamics, i.e.\ by estimating a reproductive number, as well as to provide reliable forecasts and scenario explorations. Here we present a new approach to calibrate an agent-based model -- EpiCast -- using a large set of simulation ensembles for different major metropolitan areas of the United States. In particular, we propose: a new neural network based surrogate model able to simultaneously emulate all different locations; and a novel posterior estimation that provides not only more accurate posterior estimates of all parameters but enables the joint fitting of global parameters across regions.


title: Guts, volume and Skein Modules of 3-manifolds

abstract: We consider hyperbolic links that admit alternating projections on surfaces in compact, irreducible 3-manifolds. We show that, under some mild hypotheses, the volume of the complement of such a link is bounded below in terms of a Kauffman bracket function defined on link diagrams on the surface. In the case that the 3-manifold is a thickened surface, this Kauffman bracket function leads to a Jones-type polynomial that is an isotopy invariant of links. We show that coefficients of this polynomial provide 2-sided linear bounds on the volume of hyperbolic alternating links in the thickened surface. As a corollary of the proof of this result, we deduce that the twist number of a reduced, twist reduced, checkerboard alternating link projection with disk regions, is an invariant of the link.


title: Double-peak specific heat anomaly and correlations in the Bose-Hubbard model

abstract: Considering the thermodynamics of bosons in a lattice described by the Bose-Hubbard Hamiltonian, we report the occurrence of anomalous double peaks in their specific heat dependence on temperature. This feature, usually associated with a high geometrical frustration, can also be a consequence of a purely energetic competition. By employing self-energy functional calculations combined with finite-temperature perturbation theory, we propose a mechanism based on ground-state degeneracies expressed as residual entropies. A general decomposition of the specific heat regarding all possible transitions between the system's eingenvalues provides an insight into the nature of each maximum. Furthermore, we address how the model parameters modify the structure of these peaks based on its spectral properties and atom-atom correlation function.


title: Non-equilibrium Growth and Twist of Cross-Linked Collagen Fibrils

abstract: The lysyl oxidase (LOX) enzyme that catalyses cross-link formation during the assembly of collagen fibrils in vivo is too large to diffuse within assembled fibrils, and so is incompatible with a fully equilibrium mechanism for fibril formation. We propose that enzymatic cross-links are formed at the fibril surface during the growth of collagen fibrils; as a consequence no significant reorientation of previously cross-linked collagen molecules occurs inside collagen fibrils during fibril growth in vivo. By imposing local equilibrium only at the fibril surface, we develop a coarse-grained quantitative model of in vivo fibril structure that incorporates a double-twist orientation of collagen molecules and a periodic D-band density modulation along the fibril axis. Radial growth is controlled by the density of available collagen molecules around the fibril. In contrast with earlier equilibrium models of fibril structure, we find that all fibrils can exhibit a core-shell structure that is controlled only by the fibril radius. At small radii a core is developed with a linear double-twist structure as a function of radius. Within the core the double-twist structure is largely independent of the D-band. Within the shell at larger radii, the structure approaches a constant twist configuration that is strongly coupled with the D-band. We suggest a stable radius control mechanism that corneal fibrils can exploit near the edge of the linear core regime; while larger tendon fibrils can use a cruder version of growth control that does not select a preferred radius.


title: Unified lower bounds for interactive high-dimensional estimation under information constraints

abstract: We consider distributed parameter estimation using interactive protocols subject to local information constraints such as bandwidth limitations, local differential privacy, and restricted measurements. We provide a unified framework enabling us to derive a variety of (tight) minimax lower bounds for different parametric families of distributions, both continuous and discrete, under any $\ell_p$ loss. Our lower bound framework is versatile and yields "plug-and-play" bounds that are widely applicable to a large range of estimation problems, and, for the prototypical case of the Gaussian family, circumvents limitations of previous techniques. In particular, our approach recovers bounds obtained using data processing inequalities and Cram\'er--Rao bounds, two other alternative approaches for proving lower bounds in our setting of interest. Further, for the families considered, we complement our lower bounds with matching upper bounds.

categories: cs.DS cs.DM cs.IT cs.LG math.IT math.ST stat.TH


title: Optimal Low-Degree Hardness of Maximum Independent Set

abstract: We study the algorithmic task of finding a large independent set in a sparse Erd\H{o}s-R\'{e}nyi random graph with $n$ vertices and average degree $d$. The maximum independent set is known to have size $(2 \log d / d)n$ in the double limit $n \to \infty$ followed by $d \to \infty$, but the best known polynomial-time algorithms can only find an independent set of half-optimal size $(\log d / d)n$. We show that the class of low-degree polynomial algorithms can find independent sets of half-optimal size but no larger, improving upon a result of Gamarnik, Jagannath, and the author. This generalizes earlier work by Rahman and Vir\'ag, which proved the analogous result for the weaker class of local algorithms.


title: Tensor Train Factorization and Completion under Noisy Data with Prior Analysis and Rank Estimation

abstract: Tensor train (TT) decomposition, a powerful tool for analyzing multidimensional data, exhibits superior performance in many machine learning tasks. However, existing methods for TT decomposition either suffer from noise overfitting, or require extensive fine-tuning of the balance between model complexity and representation accuracy. In this paper, a fully Bayesian treatment of TT decomposition is employed to avoid noise overfitting, by endowing it with the ability of automatic rank determination. In particular, theoretical evidence is established for adopting a Gaussian-product-Gamma prior to induce sparsity on the slices of the TT cores, so that the model complexity is automatically determined even under incomplete and noisy observed data. Furthermore, based on the proposed probabilistic model, an efficient learning algorithm is derived under the variational inference framework. Simulation results on synthetic data show the success of the proposed model and algorithm in recovering the ground-truth TT structure from incomplete noisy data. Further experiments on real-world data demonstrate the proposed algorithm performs better in image completion and image classification, compared to other existing TT decomposition algorithms.


title: Electromagnetic Emission Produced by Three-wave Interactions in a Plasma with Continiously Injected Counterstreaming Electron Beams

abstract: Three-wave interactions between Langmuir and electromagnetic waves in plasma with unstable electron flows are believed to be the main cause for type II and III solar radio emissions. The narrow band of type II bursts requires to assume that this radiation is generated in some local regions of shock fronts traveling in the solar corona, where the specific conditions for the enhancement of electromagnetic emissions near the plasma frequency harmonics are created. The reason for such enhancement at the second harmonic may be the formation of counter-streaming electron beams. There are different opinions in literature on whether the second harmonic electromagnetic emission in the presence of an additional beam can be efficient enough to markedly dominate emissions produced by a single beam. In the present paper, we carry out particle-in-cell simulations of the collision of two symmetric electron beams in plasma with open boundary conditions and show that the efficiency of beam-to-radiation power conversion can be significantly increased compared to models with periodic boundary conditions and reach the level of a few percent if three-wave interactions with electromagnetic waves near the second harmonic of the plasma frequency becomes available for the most unstable oblique beam driven modes.


title: Gravitoviscous Protoplanetary Discs with a Dust Component. IV. Disc Outer Edges, Spectral Indices, and Opacity Gaps

abstract: The crucial initial step in planet formation is the agglomeration of micron-sized dust into macroscopic aggregates. This phase is likely to happen very early during the protostellar disc formation, which is characterised by active gas dynamics. We present numerical simulations of protostellar/protoplanetary disc long-term evolution, which includes gas dynamics with self-gravity in the thin-disc limit, and bidisperse dust grain evolution due to coagulation, fragmentation, and drift through the gas. We show that the decrease of the grain size to the disc periphery leads to sharp outer edges in dust millimetre emission, which are explained by a drop in dust opacity coefficient rather than by dust surface density variations. These visible outer edges are at the location where average grain size $\approx \lambda/2\pi$, where $\lambda$ is the observational wavelength, so discs typically look more compact at longer wavelengths if dust size decreases outwards. This allows a simple recipe for reconstructing grain sizes in disc outer regions. Discs may look larger at longer wavelengths if grain size does not reach $\lambda/2\pi$ for some wavelength. Disc visible sizes evolve non-monotonically over the first million years and differ from dust and gas physical sizes by factor of a few. We compare our model with recent observation data on gas and dust disc sizes, far-infrared fluxes and spectral indices of protoplanetary discs in Lupus. We also show that non-monotonic variations of the grain size in radial direction can cause wavelength-dependent opacity gaps, which are not associated with any physical gaps in the dust density distribution.


title: Conditional Power and Friends: The Why and How of (Un)planned, Unblinded Sample Size Recalculations in Confirmatory Trials

abstract: Adapting the final sample size of a trial to the evidence accruing during the trial is a natural way to address planning uncertainty. Designs with adaptive sample size need to account for their optional stopping to guarantee strict type-I error-rate control. A variety of different methods to maintain type-I error-rate control after unplanned changes of the initial sample size have been proposed in the literature. This makes interim analyses for the purpose of sample size recalculation feasible in a regulatory context. Since the sample size is usually determined via an argument based on the power of the trial, an interim analysis raises the question of how the final sample size should be determined conditional on the accrued information. Conditional power is a concept often put forward in this context. Since it depends on the unknown effect size, we take a strict estimation perspective and compare assumed conditional power, observed conditional power, and predictive power with respect to their properties as estimators of the unknown conditional power. We then demonstrate that pre-planning an interim analysis using methodology for unplanned interim analyses is ineffective and naturally leads to the concept of optimal two-stage designs. We conclude that unplanned design adaptations should only be conducted as reaction to trial-external new evidence, operational needs to violate the originally chosen design, or post hoc changes in the objective criterion. Finally, we show that commonly discussed sample size recalculation rules can lead to paradoxical outcomes and propose two alternative ways of reacting to newly emerging trial-external evidence.


title: A bi-directional approach to comparing the modular structure of networks

abstract: Here we propose a new method to compare the modular structure of a pair of node-aligned networks. The majority of current methods, such as normalized mutual information, compare two node partitions derived from a community detection algorithm yet ignore the respective underlying network topologies. Addressing this gap, our method deploys a community detection quality function to assess the fit of each node partition with respect to the other network's connectivity structure. Specifically, for two networks A and B, we project the node partition of B onto the connectivity structure of A. By evaluating the fit of B's partition relative to A's own partition on network A (using a standard quality function), we quantify how well network A describes the modular structure of B. Repeating this in the other direction, we obtain a two-dimensional distance measure, the bi-directional (BiDir) distance. The advantages of our methodology are three-fold. First, it is adaptable to a wide class of community detection algorithms that seek to optimize an objective function. Second, it takes into account the network structure, specifically the strength of the connections within and between communities, and can thus capture differences between networks with similar partitions but where one of them might have a more defined or robust community structure. Third, it can also identify cases in which dissimilar optimal partitions hide the fact that the underlying community structure of both networks is relatively similar. We illustrate our method for a variety of community detection algorithms, including multi-resolution approaches, and a range of both simulated and real world networks.


title: Comparing the luminosity distance for gravitational waves and electromagnetic signals in a simple model of quadratic gravity

abstract: We compute the modified friction coefficient controlling the propagation of tensor metric perturbations in the context of a generalized cosmological scenario based on a theory of gravity with quadratic curvature corrections. In such a context we discuss the differences between gravitational and electromagnetic luminosity distance, as well as the differences with the standard results based on the Einstein equations. We present numerical estimates of the modified luminosity distance on the cosmic redshift scale typical of Supernovae and standard sirens.


title: Paschen curve for a PCB trace in low vacuum

abstract: For one PCB trace region, where bare high voltage trace goes near GND pad, we estimated an electrical breakdown voltage in low vacuum using simplified geometrical model under some assumptions. Experimental measurement of Paschen curve for the PCB board considered was proceed in pressure range from 2x10-3 mbar to 25 mbar.


title: PQFabric: A Permissioned Blockchain Secure from Both Classical and Quantum Attacks

abstract: Hyperledger Fabric is a prominent and flexible solution for building permissioned distributed ledger platforms. Access control and identity management relies on a Membership Service Provider (MSP) whose cryptographic interface only handles standard PKI methods for authentication: RSA and ECDSA classical signatures. Also, MSP-issued credentials may use only one signature scheme, tying the credential-related functions to classical single-signature primitives. RSA and ECDSA are vulnerable to quantum attacks, with an ongoing post-quantum standardization process to identify quantum-safe drop-in replacements. In this paper, we propose a redesign of Fabric's credential-management procedures and related specifications in order to incorporate hybrid digital signatures, protecting against both classical and quantum attacks using one classical and one quantum-safe signature. We create PQFabric, an implementation of Fabric with hybrid signatures that integrates with the Open Quantum Safe (OQS) library. Our implementation offers complete crypto-agility, with the ability to perform live migration to a hybrid quantum-safe blockchain and select any existing OQS signature algorithm for each node. We perform comparative benchmarks of PQFabric with each of the NIST candidates and alternates, revealing that long public keys and signatures lead to an increase in hashing time that is sometimes comparable to the time spent signing or verifying messages itself. This is a new and potentially significant issue in the migration of blockchains to post-quantum signatures.


title: Does my multimodal model learn cross-modal interactions? It's harder to tell than you might think!

abstract: Modeling expressive cross-modal interactions seems crucial in multimodal tasks, such as visual question answering. However, sometimes high-performing black-box algorithms turn out to be mostly exploiting unimodal signals in the data. We propose a new diagnostic tool, empirical multimodally-additive function projection (EMAP), for isolating whether or not cross-modal interactions improve performance for a given model on a given task. This function projection modifies model predictions so that cross-modal interactions are eliminated, isolating the additive, unimodal structure. For seven image+text classification tasks (on each of which we set new state-of-the-art benchmarks), we find that, in many cases, removing cross-modal interactions results in little to no performance degradation. Surprisingly, this holds even when expressive models, with capacity to consider interactions, otherwise outperform less expressive models; thus, performance improvements, even when present, often cannot be attributed to consideration of cross-modal feature interactions. We hence recommend that researchers in multimodal machine learning report the performance not only of unimodal baselines, but also the EMAP of their best-performing model.


title: On Filaments, Prolate Halos and Rotation Curves

abstract: We propose a simple geometrical mechanism for the flattening of galactic rotation curves, the local compression of field lines around their planes induced either by the presence of thin string-like objects at the centers of galaxies or by elongated dark-matter halos, and elaborate on its possible role in Nature. We fit 83 rotation curves from the SPARC database with logarithmic potentials produced by a thin "wire" at the origin and then, after selecting 2 galaxies that yield the most interesting fits, analyze them with an alternative model, deformed versions of two popular models of dark-matter halos. Our conclusion is that the presence of a filament clearly improves the fit quality in a number of cases, while bulged dark matter profiles have a lesser effect. If taken at face value, these results would imply the presence of elongated mass distributions away from the galactic plane in a number of galaxies, and may also have some indirect impact on the controversy between cold dark matter (CDM), self-interacting dark matter (SiDM), and modified Newtonian dynamics (MOND).


title: IMPECCABLE: Integrated Modeling PipelinE for COVID Cure by Assessing Better LEads

abstract: The drug discovery process currently employed in the pharmaceutical industry typically requires about 10 years and $2-3 billion to deliver one new drug. This is both too expensive and too slow, especially in emergencies like the COVID-19 pandemic. In silicomethodologies need to be improved to better select lead compounds that can proceed to later stages of the drug discovery protocol accelerating the entire process. No single methodological approach can achieve the necessary accuracy with required efficiency. Here we describe multiple algorithmic innovations to overcome this fundamental limitation, development and deployment of computational infrastructure at scale integrates multiple artificial intelligence and simulation-based approaches. Three measures of performance are:(i) throughput, the number of ligands per unit time; (ii) scientific performance, the number of effective ligands sampled per unit time and (iii) peak performance, in flop/s. The capabilities outlined here have been used in production for several months as the workhorse of the computational infrastructure to support the capabilities of the US-DOE National Virtual Biotechnology Laboratory in combination with resources from the EU Centre of Excellence in Computational Biomedicine.


title: The norm map on the Bloch group for quadratic extensions

abstract: D. Rudenko proved the homotopy invariance of the truncated polylogarithmic complexes. It follows that on these complexes there is the norm map with good proprieties. We apply his result and get the explicit formula for the norm map in the case of quadratic extensions.


title: High-speed calibration and characterization of superconducting quantum processors without qubit reset

abstract: To Characterize and calibrate quantum processing devices a large amount of measurement data has to be collected. Active qubit reset increases the speed at which data can be gathered but requires additional hardware and/or calibration. The experimental apparatus can, however, be operated at elevated repetition rates without reset. In this case, the outcome of a first measurement serves as the initial state for the next experiment. Rol. $\textit{et al.}$ used this restless operation mode to accelerate the calibration of a single-qubit gate by measuring fixed-length sequences of Clifford gates which compose to $X$ gates [Phys. Rev. Appl. $7$, 041001 (2017)]. However, we find that, when measuring pulse sequences which compose to arbitrary operations, a distortion appears in the measured data. Here, we extend the restless methodology by showing how to efficiently analyze restless measurements and correct distortions to achieve an identical outcome and accuracy as compared to measurements in which the superconducting qubits are reset. This allows us to rapidly characterize and calibrate qubits. We illustrate our data collection and analysis method by measuring a Rabi oscillation at a $250~\rm{kHz}$ repetition rate without any reset, for a qubit with a decay rate of $1/2\pi T_1=3~\rm{kHz}$. We also show that we can measure a single- and a two-qubit average gate fidelity with Randomized Benchmarking 20 and 8 times faster, respectively, than measurements that reset the qubits through $T_1$-decay.


title: Non-orientable link cobordisms and torsion order in Floer homologies

abstract: We use unoriented versions of instanton and knot Floer homology to prove inequalities involving the Euler characteristic and the number of local maxima appearing in unorientable cobordisms, which mirror results of a recent paper by Juhasz, Miller, and Zemke concerning orientable cobordisms. Most of the subtlety in our argument lies in the fact that maps for non-orientable cobordisms require more complicated decorations than their orientable counterparts. We introduce unoriented versions of the band unknotting number and the refined cobordism distance and apply our results to give bounds on these based on the torsion orders of the Floer homologies. Finally, we show that the difference between the unoriented refined cobordism distance of a knot $K$ from the unknot and the non-orientable slice genus of $K$ can be arbitrarily large.


title: Refined pointwise estimates for solutions to the 1D barotropic compressible Navier--Stokes equations: An application to the long-time behavior of a point mass

abstract: We study the long-time behavior of a point mass moving in a one-dimensional viscous compressible fluid. Previously, we showed that the velocity of the point mass $V(t)$ satisfies a decay estimate $V(t)=O(t^{-3/2})$~[K. Koike, J. Differential Equations \textbf{271} (2021) 356--413]. This result was obtained as a corollary to pointwise estimates of solutions to a free boundary problem of barotropic compressible Navier--Stokes equations. In this paper, we give a simple necessary and sufficient condition on the initial data for the decay estimate $V(t)=O(t^{-3/2})$ to be optimal. This is achieved by refining the pointwise estimates previously obtained: we make use of \textit{inter-diffusion waves} that, together with the classical \textit{diffusion waves}, give an improved approximation of the fluid behavior around the point mass; this then leads to a sharper understanding of the long-time behavior of the point mass.


title: Fantastic Features and Where to Find Them: Detecting Cognitive Impairment with a Subsequence Classification Guided Approach

abstract: Despite the widely reported success of embedding-based machine learning methods on natural language processing tasks, the use of more easily interpreted engineered features remains common in fields such as cognitive impairment (CI) detection. Manually engineering features from noisy text is time and resource consuming, and can potentially result in features that do not enhance model performance. To combat this, we describe a new approach to feature engineering that leverages sequential machine learning models and domain knowledge to predict which features help enhance performance. We provide a concrete example of this method on a standard data set of CI speech and demonstrate that CI classification accuracy improves by 2.3% over a strong baseline when using features produced by this method. This demonstration provides an ex-ample of how this method can be used to assist classification in fields where interpretability is important, such as health care.


title: Scenic: A Language for Scenario Specification and Data Generation

abstract: We propose a new probabilistic programming language for the design and analysis of cyber-physical systems, especially those based on machine learning. Specifically, we consider the problems of training a system to be robust to rare events, testing its performance under different conditions, and debugging failures. We show how a probabilistic programming language can help address these problems by specifying distributions encoding interesting types of inputs, then sampling these to generate specialized training and test data. More generally, such languages can be used to write environment models, an essential prerequisite to any formal analysis. In this paper, we focus on systems like autonomous cars and robots, whose environment at any point in time is a 'scene', a configuration of physical objects and agents. We design a domain-specific language, Scenic, for describing scenarios that are distributions over scenes and the behaviors of their agents over time. As a probabilistic programming language, Scenic allows assigning distributions to features of the scene, as well as declaratively imposing hard and soft constraints over the scene. We develop specialized techniques for sampling from the resulting distribution, taking advantage of the structure provided by Scenic's domain-specific syntax. Finally, we apply Scenic in a case study on a convolutional neural network designed to detect cars in road images, improving its performance beyond that achieved by state-of-the-art synthetic data generation methods.


title: Geomechanical simulation of energy storage in salt formations

abstract: A promising option for storing large-scale quantities of green gases (e.g., hydrogen) is in subsurface rock salt caverns. The mechanical performance of salt caverns utilized for long-term subsurface energy storage plays a significant role in long-term stability and serviceability. However, rock salt undergoes non-linear creep deformation due to long-term loading caused by subsurface storage. Salt caverns have complex geometries and the geological domain surrounding salt caverns has a vast amount of material heterogeneity. To safely store gases in caverns, a thorough analysis of the geological domain becomes crucial. To date, few studies have attempted to analyze the influence of geometrical and material heterogeneity on the state of stress in salt caverns subjected to long-term loading. In this work, we present a rigorous and systematic modeling study to quantify the impact of heterogeneity on the deformation of salt caverns and quantify the state of stress around the caverns. A 2D finite element simulator was developed to consistently account for the non-linear creep deformation and also to model tertiary creep. The computational scheme was benchmarked with the already existing experimental study. The impact of cyclic loading on the cavern was studied considering maximum and minimum pressure that depends on lithostatic pressure. The influence of geometric heterogeneity such as irregularly-shaped caverns and material heterogeneity, which involves different elastic and creep properties of the different materials in the geological domain, is rigorously studied and quantified. Moreover, multi-cavern simulations are conducted to investigate the influence of a cavern on the adjacent caverns.