title: Examining mobility data justice during 2017 Hurricane Harvey

abstract: Natural disasters can significantly disrupt human mobility in urban areas. Studies have attempted to understand and quantify such disruptions using crowdsourced mobility data sets. However, limited research has studied the justice issues of mobility data in the context of natural disasters. The lack of research leaves us without an empirical foundation to quantify and control the possible biases in the data. This study, using 2017 Hurricane Harvey as a case study, explores three aspects of mobility data that could potentially cause injustice: representativeness, quality, and precision. We find representativeness being a major factor contributing to mobility data injustice. There is a persistent disparity of representativeness across neighborhoods of different socioeconomic characteristics before, during, and after the hurricane's landfall. Additionally, we observed significant drops of data precision during the hurricane, adding uncertainty to locate people and understand their movements during extreme weather events. The findings highlight the necessity in understanding and controlling the possible bias of mobility data as well as developing practical tools through data justice lenses in collecting and analyzing data during disasters.


title: Some basic results on finite linear recurring sequence subgroups

abstract: An $f$-subgroup is a linear recurring sequence subgroup, a multiplicative subgroup of a field whose elements can be generated (without repetition) by a linear recurrence relation, with characteristic polynomial $f$. It is called non-standard if it can be generated in a non-cyclic way (that is, not in the order $\alpha^i, \alpha^{i+1}, \alpha^{i+2} \ldots$ for a zero $\alpha$ of $f$), and standard otherwise. We will show that a finite $f$-subgroup is necessarily generated by a subset of the zeros of $f$. We use this result to improve on a recent theorem of Brison and Nogueira. A old question by Brison and Nogueira asks if there exist automatically non-standard $f$-subgroups, $f$-subgroups that cannot be generated by a zero of $f$. We answer that question affirmatively by constructing infinitely many examples.


title: Plasma Spray Process Parameters Configuration using Sample-efficient Batch Bayesian Optimization

abstract: Recent work has shown constrained Bayesian optimization to be a powerful technique for the optimization of industrial processes. In complex manufacturing processes, the possibility to run extensive sequences of experiments with the goal of finding good process parameters is severely limited by the time required for quality evaluation of the produced parts. To accelerate the process parameter optimization, we introduce a parallel acquisition procedure tailored on the process characteristics. We further propose an algorithm that adapts to equipment status to improve run-to-run reproducibility. We validate our optimization method numerically and experimentally, and demonstrate that it can efficiently find input parameters that produce the desired outcome and minimize the process cost.


title: If You Must Choose Among Your Children, Pick the Right One

abstract: Given a simplicial complex $K$ and an injective function $f$ from the vertices of $K$ to $\mathbb{R}$, we consider algorithms that extend $f$ to a discrete Morse function on $K$. We show that an algorithm of King, Knudson and Mramor can be described on the directed Hasse diagram of $K$. Our description has a faster runtime for high dimensional data with no increase in space.


title: Risk Bounds and Rademacher Complexity in Batch Reinforcement Learning

abstract: This paper considers batch Reinforcement Learning (RL) with general value function approximation. Our study investigates the minimal assumptions to reliably estimate/minimize Bellman error, and characterizes the generalization performance by (local) Rademacher complexities of general function classes, which makes initial steps in bridging the gap between statistical learning theory and batch RL. Concretely, we view the Bellman error as a surrogate loss for the optimality gap, and prove the followings: (1) In double sampling regime, the excess risk of Empirical Risk Minimizer (ERM) is bounded by the Rademacher complexity of the function class. (2) In the single sampling regime, sample-efficient risk minimization is not possible without further assumptions, regardless of algorithms. However, with completeness assumptions, the excess risk of FQI and a minimax style algorithm can be again bounded by the Rademacher complexity of the corresponding function classes. (3) Fast statistical rates can be achieved by using tools of local Rademacher complexity. Our analysis covers a wide range of function classes, including finite classes, linear spaces, kernel spaces, sparse linear features, etc.


title: Linear and Nonlinear Optical Propagation in 2D Materials

abstract: Recently, a lot of effort has been dedicated to developing next-generation optoelectronic devices based on two-dimensional materials, thanks to their unique optical properties that are significantly different from those of their bulk counterparts. In order to implement high-performance nanoscale optical devices, an in-depth study of how linear and non-linear propagation occurs in two-dimensional materials is required. Here, we focus on the theory behind the propagation of electromagnetic waves in two-dimensional materials as well as emerging applications in the fields of electronics, optics, sensors, which are summarized and discussed in the paper.


title: Supervised Contrastive Replay: Revisiting the Nearest Class Mean Classifier in Online Class-Incremental Continual Learning

abstract: Online class-incremental continual learning (CL) studies the problem of learning new classes continually from an online non-stationary data stream, intending to adapt to new data while mitigating catastrophic forgetting. While memory replay has shown promising results, the recency bias in online learning caused by the commonly used Softmax classifier remains an unsolved challenge. Although the Nearest-Class-Mean (NCM) classifier is significantly undervalued in the CL community, we demonstrate that it is a simple yet effective substitute for the Softmax classifier. It addresses the recency bias and avoids structural changes in the fully-connected layer for new classes. Moreover, we observe considerable and consistent performance gains when replacing the Softmax classifier with the NCM classifier for several state-of-the-art replay methods. To leverage the NCM classifier more effectively, data embeddings belonging to the same class should be clustered and well-separated from those with a different class label. To this end, we contribute Supervised Contrastive Replay (SCR), which explicitly encourages samples from the same class to cluster tightly in embedding space while pushing those of different classes further apart during replay-based training. Overall, we observe that our proposed SCR substantially reduces catastrophic forgetting and outperforms state-of-the-art CL methods by a significant margin on a variety of datasets.


title: Robust and Accurate Object Detection via Adversarial Learning

abstract: Data augmentation has become a de facto component for training high-performance deep image classifiers, but its potential is under-explored for object detection. Noting that most state-of-the-art object detectors benefit from fine-tuning a pre-trained classifier, we first study how the classifiers' gains from various data augmentations transfer to object detection. The results are discouraging; the gains diminish after fine-tuning in terms of either accuracy or robustness. This work instead augments the fine-tuning stage for object detectors by exploring adversarial examples, which can be viewed as a model-dependent data augmentation. Our method dynamically selects the stronger adversarial images sourced from a detector's classification and localization branches and evolves with the detector to ensure the augmentation policy stays current and relevant. This model-dependent augmentation generalizes to different object detectors better than AutoAugment, a model-agnostic augmentation policy searched based on one particular detector. Our approach boosts the performance of state-of-the-art EfficientDets by +1.1 mAP on the COCO object detection benchmark. It also improves the detectors' robustness against natural distortions by +3.8 mAP and against domain shift by +1.3 mAP. Models are available at https://github.com/google/automl/tree/master/efficientdet/Det-AdvProp.md


title: Adversarial Imitation Learning with Trajectorial Augmentation and Correction

abstract: Deep Imitation Learning requires a large number of expert demonstrations, which are not always easy to obtain, especially for complex tasks. A way to overcome this shortage of labels is through data augmentation. However, this cannot be easily applied to control tasks due to the sequential nature of the problem. In this work, we introduce a novel augmentation method which preserves the success of the augmented trajectories. To achieve this, we introduce a semi-supervised correction network that aims to correct distorted expert actions. To adequately test the abilities of the correction network, we develop an adversarial data augmented imitation architecture to train an imitation agent using synthetic experts. Additionally, we introduce a metric to measure diversity in trajectory datasets. Experiments show that our data augmentation strategy can improve accuracy and convergence time of adversarial imitation while preserving the diversity between the generated and real trajectories.


title: On a theorem of Chernoff on rank one Riemannian symmetric spaces

abstract: In 1975, P.R. Chernoff used iterates of the Laplacian on $\mathbb{R}^n$ to prove an $L^2$ version of the Denjoy-Carleman theorem which provides a sufficient condition for a smooth function on $\mathbb{R}^n$ to be quasi-analytic. In this paper, we prove an exact analogue of Chernoff's theorem for all rank one Riemannian symmetric spaces (of noncompact and compact types) using iterates of the associated Laplace-Beltrami operators.


title: Exponential localization of Steklov eigenfunctions on warped product manifolds: the flea on the elephant phenomenon

abstract: This paper is devoted to the analysis of Steklov eigenvalues and Steklov eigenfunctions on a class of warped product Riemannian manifolds $(M,g)$ whose boundary $\partial M$ consists in two distinct connected components $\Gamma_0$ and $\Gamma_1$. First, we show that the Steklov eigenvalues can be divided into two families $(\lambda_m^\pm)_{m \geq 0}$ which satisfy accurate asymptotics as $m \to \infty$. Second, we consider the associated Steklov eigenfunctions which are the harmonic extensions of the boundary Dirichlet to Neumann eigenfunctions. In the case of symmetric warped product, we prove that the Steklov eigenfunctions are exponentially localized on the whole boundary $\partial M$ as $m \to \infty$. Whenever we add an asymmetric perturbation to a symmetric warped product, we observe a flea on the elephant effect. Roughly speaking, we prove that "half" the Steklov eigenfunctions are exponentially localized on one connected component of the boundary, say $\Gamma_0$, and the other half on the other connected component $\Gamma_1$ as $m \to \infty$.


title: A Multigrid Preconditioner for Jacobian-free Newton-Krylov Methods

abstract: In this work, we propose a multigrid preconditioner for Jacobian-free Newton-Krylov (JFNK) methods. Our multigrid method does not require knowledge of the Jacobian at any level of the multigrid hierarchy. As it is common in standard multigrid methods, the proposed method also relies on three building blocks: transfer operators, smoothers, and a coarse level solver. In addition to the restriction and prolongation operator, we also use a projection operator to transfer the current Newton iterate to a coarser level. The three-level Chebyshev semi-iterative method is employed as a smoother, as it has good smoothing properties and does not require the representation of the Jacobian matrix. We replace the direct solver on the coarsest level with a matrix-free Krylov subspace method, thus giving rise to a truly Jacobian-free multigrid preconditioner. We will discuss all building blocks of our multigrid preconditioner in detail and demonstrate the robustness and the efficiency of the proposed method using several numerical examples.


title: Exact weight cancellation in Monte Carlo eigenvalue transport problems

abstract: Random walks are frequently used as a model for very diverse physical phenomena. The Monte Carlo method is a versatile tool for the study of the properties of systems modelled as random walks. Often, each walker is associated with a statistical weight, used in the estimation of observable quantities. Weights are typically assumed to be positive; nonetheless, some applications require the use of positive and negative weights or complex weights, and often pose particular challenges with convergence. In this paper, we examine such a case from the field of nuclear reactor physics, where the negative particle weights prevent the power iteration algorithm from converging on the sought fundamental eigenstate of the Boltzmann transport equation. We demonstrate how the use of weight cancellation allows convergence on the physical eigenstate. To this end, we develop a novel method to perform weight cancellation in an exact manner, in three spatial dimensions. The viability of this algorithm is then demonstrated on a reactor physics problem.


title: Transmit Beamspace DDMA Based Automotive MIMO Radar

abstract: The time division multiple access (TDMA) technique has been applied in automotive multiple-input multiple-output (MIMO) radar. However, it suffers from the transmit energy loss, and as a result the parameter estimation performance degradation when the number of transmit elements increases. To tackle these problem, a transmit beamspace (TB) Doppler division multiple access (DDMA) approach is proposed. First, a phase modulation matrix with empty Doppler spectrum is introduced. By exploiting the empty Doppler spectrum, a test function based on sequential detection is developed to mitigate the Doppler ambiguity in DDMA waveform. Then, a discrete Fourier transform (DFT)-based TB in slow-time is formed.The proposed method can achieve waveform diversity in Doppler domain and generate a TB in slow-time that concentrates the transmitted power in a fixed spatial region to improve the transmit energy distribution for automotive MIMO radar, which is favored by medium/long range radar (MRR/LRR) applications. As compared to the conventional TDMA technique, the proposed TB DDMA approach can fully exploit the transmission capabilities of all transmit elements to ensure that the emitted power is efficiently used and inherits easy implementation. Moreover, the proposed TB DDMA method avoids the trade-off between the active time for each transmit antenna and the frame time. Simulation results verify the effectiveness of the proposed TB DDMA approach for automotive MIMO radar.


title: Modeling stellar Ca II H & K emission variations. I. Effect of inclination on the S-index

abstract: The emission in the near ultraviolet Ca II H & K lines is modulated by stellar magnetic activity. Although this emission, quantified via the S-index, has been serving as a prime proxy of stellar magnetic activity for several decades, many aspects of the complex relation between stellar magnetism and Ca II H & K emission are still unclear. The amount of measured Ca II H & K emission is suspected to be affected not only by the stellar intrinsic properties but also by the inclination angle of the stellar rotation axis. Until now such an inclination effect on S-index has remained largely unexplored. To fill this gap, we develop a physics-based model to calculate S-index, focusing on the Sun. Using the distributions of solar magnetic features derived from observations together with Ca II H & K spectra synthesized in non-local thermodynamic equilibrium, we validate our model by successfully reconstructing the observed variations of solar S-index over four activity cycles. Further, using the distribution of magnetic features over the visible solar disk obtained from surface flux transport simulations, we obtain S-index time series dating back to 1700 and investigate the effect of inclination on S-index variability, both on the magnetic activity cycle and the rotational timescales. We find that when going from an equatorial to a pole-on view, the amplitude of S-index variations decreases weakly on the activity cycle timescale and strongly on the rotational timescale (by about 22% and 81%, respectively, for a cycle of intermediate strength). The absolute value of S-index depends only weakly on the inclination. We provide analytical expressions that model such dependencies.


title: Boosting Binary Masks for Multi-Domain Learning through Affine Transformations

abstract: In this work, we present a new, algorithm for multi-domain learning. Given a pretrained architecture and a set of visual domains received sequentially, the goal of multi-domain learning is to produce a single model performing a task in all the domains together. Recent works showed how we can address this problem by masking the internal weights of a given original conv-net through learned binary variables. In this work, we provide a general formulation of binary mask based models for multi-domain learning by affine transformations of the original network parameters. Our formulation obtains significantly higher levels of adaptation to new domains, achieving performances comparable to domain-specific models while requiring slightly more than 1 bit per network parameter per additional domain. Experiments on two popular benchmarks showcase the power of our approach, achieving performances close to state-of-the-art methods on the Visual Decathlon Challenge.


title: The $C_2$-equivariant $K(1)$-local sphere

abstract: We compute the $RO(C_2)$-graded Green functor $\underline{\pi}_\star L_{KU_{C_2}/(2)}S_{C_2}$.


title: A unified formalism to study $soft$ as well as $hard$ part of the transverse momentum spectra

abstract: Transverse momentum $p_T$ spectra of final state particles produced in high energy heavy-ion collision can be divided into two distinct regions based on the difference in the underlying particle production process. We have provided a unified formalism to explain both low- and high-$p_T$ regime of spectra in a consistent manner. The $p_T$ spectra of final state particles produced at RHIC and LHC energies have been analysed using unified formalism to test its applicability at different energies, and a good agreement with the data is obtained across all energies. Further, the prospect of extracting the elliptic flow coefficient directly from the transverse momentum spectra is explored.


title: On the Haefliger-Thurston conjecture

abstract: The classifying space for the framed Haefliger structures of codimension $q$ and class $C^r$ is $(2q-1)$-connected, for $1\le r\le\infty$. The corollaries deal with the existence of foliations, with the homology and the perfectness of the diffeomorphism groups, with the existence of foliated products, and of foliated bundles.


title: Recurrent Neural Network for End-to-End Modeling of Laminar-Turbulent Transition

abstract: Accurate prediction of laminar-turbulent transition is a critical element of computational fluid dynamics simulations for aerodynamic design across multiple flow regimes. Traditional methods of transition prediction cannot be easily extended to flow configurations where the transition process depends on a large set of parameters. In comparison, neural network methods allow higher dimensional input features to be considered without compromising the efficiency and accuracy of the traditional data driven models. Neural network methods proposed earlier follow a cumbersome methodology of predicting instability growth rates over a broad range of frequencies, which are then processed to obtain the N-factor envelope, and then, the transition location based on the correlating N-factor. This paper presents an end-to-end transition model based on a recurrent neural network, which sequentially processes the mean boundary-layer profiles along the surface of the aerodynamic body to directly predict the N-factor envelope and the transition locations over a two-dimensional airfoil. The proposed transition model has been developed and assessed using a large database of 53 airfoils over a wide range of chord Reynolds numbers and angles of attack. The sequence-to-sequence transduction model proposed herein provides a more direct approach for accurate predictions of the transition location than the earlier neural network methods, which predict the local amplification rate of a single instability mode at a fixed location along the airfoil. The large universe of airfoils encountered in various applications causes additional difficulties. As such, we provide further insights on selecting training datasets from large amounts of available data.


title: Datasets of the solar quiet (Sq) and solar disturbed (SD) variations of the geomagnetic field from the mid latitudinal Magnetic Observatory of Coimbra (Portugal) obtained by different methods

abstract: Here we present datasets of daily variation obtained from the geomagnetic field raw observations recorded at the Coimbra Magnetic Observatory (COI, Portugal) in 2007-2017, covering almost the entire solar cycle 24. Two methods were used to extract daily variability from the raw geomagnetic hourly data. The first method uses the so-called geomagnetically quiet days to calculate S-type variations as daily means resulting in the data sub-set named IQD Sq and SD. The second method uses the principal component analysis (PCA) to decompose the original series into main variability modes. The first three modes produced by PCA and explaining up to 98 per cent of the variability of the raw data are in the data sub-set named PCA modes. Both methods allow to extract regular geomagnetic field variations related to daily variations (S-type variations) in the ionospheric dynamo region and some magnetospheric currents. The COI location near the mean latitude of the ionospheric Sq current vortex focus allows studying its variability using the S-type regular variations of the geomagnetic field measured near the ground. The S-type variations for the X and Y components of the geomagnetic field obtained at the COI observatory can also be re-scaled and used to analyze field variations obtained at other European geomagnetic observatories at close latitudes. The S-type variations for the Z component of the geomagnetic field obtained at the COI observatory can be compared to similar variations observed at more continental regions to study the so-called coastal effect in the geomagnetic field variations. The dataset described in this paper is analyzed in a companion paper Comparison of the solar variations of the geomagnetic field at the Coimbra Magnetic Observatory (COI) obtained by different methods: effect of the solar and geomagnetic activity by A. Morozova and R. Rebbah submitted to ASR.


title: Logarithmic law of large random correlation matrices

abstract: Consider a random vector $\mathbf{y}=\mathbf{\Sigma}^{1/2}\mathbf{x}$, where the $p$ elements of the vector $\mathbf{x}$ are i.i.d. real-valued random variables with zero mean and finite fourth moment, and $\mathbf{\Sigma}^{1/2}$ is a deterministic $p\times p$ matrix such that the spectral norm of the population correlation matrix $\mathbf{R}$ of $\mathbf{y}$ is uniformly bounded. In this paper, we find that the log determinant of the sample correlation matrix $\hat{\mathbf{R}}$ based on a sample of size $n$ from the distribution of $\mathbf{y}$ satisfies a CLT (central limit theorem) for $p/n\to \gamma\in (0, 1]$ and $p\leq n$. Explicit formulas for the asymptotic mean and variance are provided. In case the mean of $\mathbf{y}$ is unknown, we show that after recentering by the empirical mean the obtained CLT holds with a shift in the asymptotic mean. This result is of independent interest in both large dimensional random matrix theory and high-dimensional statistical literature of large sample correlation matrices for non-normal data. At last, the obtained findings are applied for testing of uncorrelatedness of $p$ random variables. Surprisingly, in the null case $\mathbf{R}=\mathbf{I}$, the test statistic becomes completely pivotal and the extensive simulations show that the obtained CLT also holds if the moments of order four do not exist at all, which conjectures a promising and robust test statistic for heavy-tailed high-dimensional data.


title: Measure Theoretic Weighted Model Integration

abstract: Weighted model counting (WMC) is a popular framework to perform probabilistic inference with discrete random variables. Recently, WMC has been extended to weighted model integration (WMI) in order to additionally handle continuous variables. At their core, WMI problems consist of computing integrals and sums over weighted logical formulas. From a theoretical standpoint, WMI has been formulated by patching the sum over weighted formulas, which is already present in WMC, with Riemann integration. A more principled approach to integration, which is rooted in measure theory, is Lebesgue integration. Lebesgue integration allows one to treat discrete and continuous variables on equal footing in a principled fashion. We propose a theoretically sound measure theoretic formulation of weighted model integration, which naturally reduces to weighted model counting in the absence of continuous variables. Instead of regarding weighted model integration as an extension of weighted model counting, WMC emerges as a special case of WMI in our formulation.


title: Near Real-time Learning and Extraction of Attack Models from Intrusion Alerts

abstract: Critical and sophisticated cyberattacks often take multitudes of reconnaissance, exploitations, and obfuscation techniques to penetrate through well protected enterprise networks. The discovery and detection of attacks, though needing continuous efforts, is no longer sufficient. Security Operation Center (SOC) analysts are overwhelmed by the significant volume of intrusion alerts without being able to extract actionable intelligence. Recognizing this challenge, this paper describes the advances and findings through deploying ASSERT to process intrusion alerts from OmniSOC in collaboration with the Center for Applied Cybersecurity Research (CACR) at Indiana University. ASSERT utilizes information theoretic unsupervised learning to extract and update `attack models' in near real-time without expert knowledge. It consumes streaming intrusion alerts and generates a small number of statistical models for SOC analysts to comprehend ongoing and emerging attacks in a timely manner. This paper presents the architecture and key processes of ASSERT and discusses a few real-world attack models to highlight the use-cases that benefit SOC operations. The research team is developing a light-weight containerized ASSERT that will be shared through a public repository to help the community combat the overwhelming intrusion alerts.


title: Matching tetrads in f(T) gravity

abstract: The procedure underlying the matching of 1-form (tetrad) fields in theories possessing absolute parallelism -- f(T) gravity being within this category -- is addressed and exemplified. We show that the remnant symmetries of the intervening spaces play a central role in the process, because the knowledge of the remnant group of local Lorentz transformations enables one to perform rotations and/or boosts in order to $\mathcal{C}^{1}$-match the corresponding tetrads on the junction surface. This automatically ensures the continuity of the Weitzenb\"{o}ck scalar there, even though this proves to be just a necessary condition in order to obtain a global parallelization of the spacetime.


title: Cosmological implications of the anisotropy of ten galaxy cluster scaling relations

abstract: The hypothesis that the late Universe is isotropic and homogeneous is adopted by most cosmological studies. The expansion rate $H_0$ is thought to be spatially constant, while bulk flows are often presumed to be negligible compared to the Hubble expansion, even at local scales. Their effects on the redshift-distance conversion are hence usually ignored. Any deviation from this consensus can strongly bias the results of such studies and thus the importance of testing these assumptions cannot be understated. Scaling relations of galaxy clusters can be effectively used for that. In previous works, we observed strong anisotropies in cluster scaling relations, whose origins remain ambiguous. By measuring many different cluster properties, several scaling relations with different sensitivities can be built. Nearly independent tests of cosmic isotropy and bulk flows are then feasible. We make use of up to 570 clusters with measured properties at X-ray, microwave, and infrared wavelengths, to construct 10 different cluster scaling relations (five of them presented for the first time) and test the isotropy of the local Universe. Through rigorous tests, we ensure that our analysis is not prone to generally known systematic biases and X-ray absorption issues. By combining all available information, we detect an apparent $9\%$ spatial variation in the local $H_0$ between $(l,b)\sim ({280^{\circ}}^{+35^{\circ}}_{-35^{\circ}},{-15^{\circ}}^{+20^{\circ}}_{-20^{\circ}})$ and the rest of the sky. The observed anisotropy has a nearly dipole form. Using Monte Carlo simulations, we assess the statistical significance of the anisotropy to be $>5\sigma$. This result could also be attributed to a $\sim 900$ km/s bulk flow which seems to extend out to at least $\sim 500$ Mpc. These two effects are indistinguishable until more high$-z$ clusters are observed by future all-sky surveys, such as eROSITA.


title: StyleLess layer: Improving robustness for real-world driving

abstract: Deep Neural Networks (DNNs) are a critical component for self-driving vehicles. They achieve impressive performance by reaping information from high amounts of labeled data. Yet, the full complexity of the real world cannot be encapsulated in the training data, no matter how big the dataset, and DNNs can hardly generalize to unseen conditions. Robustness to various image corruptions, caused by changing weather conditions or sensor degradation and aging, is crucial for safety when such vehicles are deployed in the real world. We address this problem through a novel type of layer, dubbed StyleLess, which enables DNNs to learn robust and informative features that can cope with varying external conditions. We propose multiple variations of this layer that can be integrated in most of the architectures and trained jointly with the main task. We validate our contribution on typical autonomous-driving tasks (detection, semantic segmentation), showing that in most cases, this approach improves predictive performance on unseen conditions (fog, rain), while preserving performance on seen conditions and objects.


title: About exchanging expectation and supremum for conditional Wasserstein GANs

abstract: In cases where a Wasserstein GAN depends on a condition the latter is usually handled via an expectation within the loss function. Depending on the way this is motivated, the discriminator is either required to be Lipschitz-1 in both or in only one of its arguments. For the weaker requirement to become usable one needs to exchange a supremum and an expectation. This is a mathematically perilous operation, which is, so far, only partially justified in the literature. This short mathematical note intends to fill this gap and provides the mathematical rationale for discriminators that are only partially Lipschitz-1 for cases where this approach is more appropriate or successful.


title: On the molecular correlations that result in field-dependent conductivities in electrolyte solutions

abstract: Employing recent advances in response theory and nonequilibrium ensemble reweighting, we study the dynamic and static correlations that give rise to an electric field-dependent ionic conductivity in electrolyte solutions. We consider solutions modeled with both implicit and explicit solvents, with different dielectric properties, and at multiple concentrations. Implicit solvent models at low concentrations and small dielectric constants exhibit strongly field-dependent conductivities. We compared these results to the Onsager-Wilson theory of the Wien effect, which provides a qualitatively consistent prediction at low concentrations and high static dielectric constants, but is inconsistent away from these regimes. The origin of the discrepancy is found to be increased ion correlations under these conditions. Explicit solvent effects act to suppress nonlinear responses, yielding a weakly field-dependent conductivity over the range of physically realizable field strengths. By decomposing the relevant time correlation functions, we find that the insensitivity of the conductivity to the field results from the persistent frictional forces on the ions from the solvent. Our findings illustrate the utility of nonequilibrium response theory in rationalizing nonlinear transport behavior.

categories: cond-mat.soft cond-mat.mes-hall cond-mat.stat-mech physics.chem-ph


title: Kleiss-Kuijf Relations from Momentum Amplituhedron Geometry

abstract: In recent years, it has been understood that color-ordered scattering amplitudes can be encoded as logarithmic differential forms on positive geometries. In particular, amplitudes in maximally supersymmetric Yang-Mills theory in spinor helicity space are governed by the momentum amplituhedron. Due to the group-theoretic structure underlying color decompositions, color-ordered amplitudes enjoy various identities which relate different orderings. In this paper, we show how the Kleiss-Kuijf relations arise from the geometry of the momentum amplituhedron. We also show how similar relations can be realised for the kinematic associahedron, which is the positive geometry of bi-adjoint scalar cubic theory.


title: Regularization by Denoising Sub-sampled Newton Method for Spectral CT Multi-Material Decomposition

abstract: Spectral Computed Tomography (CT) is an emerging technology that enables to estimate the concentration of basis materials within a scanned object by exploiting different photon energy spectra. In this work, we aim at efficiently solving a model-based maximum-a-posterior problem to reconstruct multi-materials images with application to spectral CT. In particular, we propose to solve a regularized optimization problem based on a plug-in image-denoising function using a randomized second order method. By approximating the Newton step using a sketching of the Hessian of the likelihood function, it is possible to reduce the complexity while retaining the complex prior structure given by the data-driven regularizer. We exploit a non-uniform block sub-sampling of the Hessian with inexact but efficient Conjugate gradient updates that require only Jacobian-vector products for denoising term. Finally, we show numerical and experimental results for spectral CT materials decomposition.


title: Deep Underground Neutrino Experiment (DUNE) Near Detector Conceptual Design Report

abstract: This report describes the conceptual design of the DUNE near detector


title: Stable moduli spaces of hermitian forms

abstract: We prove that Grothendieck-Witt spaces of Poincar\'e categories are, in many cases, group completions of certain moduli spaces of hermitian forms. This, in particular, identifies Karoubi's classical hermitian and quadratic K-groups with the genuine Grothendieck-Witt groups from our joint work with Calm\`es, Dotto, Harpaz, Land, Moi, Nardin and Nikolaus, and thereby completes our solution of several conjectures in hermitian K-theory. The method of proof is abstracted from work of Galatius and Randal-Williams on cobordism categories of manifolds using the identification of the Grothendieck-Witt space of a Poincar\'e category as the homotopy type of the associated cobordism category.


title: Existence of weak solutions to the two-dimensional incompressible Euler equations in the presence of sources and sinks

abstract: A classical model for sources and sinks in a two-dimensional perfect incompressible fluid occupying a bounded domain dates back to Yudovich in 1966. In this model, on the one hand, the normal component of the fluid velocity is prescribed on the boundary and is nonzero on an open subset of the boundary, corresponding either to sources (where the flow is incoming) or to sinks (where the flow is outgoing). On the other hand the vorticity of the fluid which is entering into the domain from the sources is prescribed. In this paper we investigate the existence of weak solutions to this system by relying on \textit{a priori} bounds of the vorticity, which satisfies a transport equation associated with the fluid velocity vector field. Our results cover the case where the vorticity has a $L^p$ integrability in space, with $p $ in $[1,+\infty]$, and prove the existence of solutions obtained by compactness methods from viscous approximations. More precisely we prove the existence of solutions which satisfy the vorticity equation in the distributional sense in the case where $p >\frac43$, in the renormalized sense in the case where $p >1$, and in a symmetrized sense in the case where $p =1$.


title: Target formation on the circle by monotone system design

abstract: Positivity and Perron-Frobenius theory provide an elegant framework for the convergence analysis of linear consensus algorithms. Here we consider a generalization of these ideas to the analysis of nonlinear consensus algorithms on the circle and establish tools for the design of consensus protocols that monotonically converge to target formations on the circle.


title: Fixed Point Theorems in M-distance Spaces

abstract: We prove fixed point theorems in a space with a distance function that takes values in a partially ordered monoid. On the one hand, such an approach allows one to generalize some fixed point theorems in a broad class of spaces, including metric and uniform spaces. On the other hand, compared to the so-called cone metric spaces and $K$-metric spaces, we do not require that the distance function range has a linear structure. We also consider several applications of the obtained fixed point theorems. In particular, we consider the questions of the existence of solutions of the Fredholm integral equation in $L$-spaces.


title: An Image is Worth 16x16 Words, What is a Video Worth?

abstract: Leading methods in the domain of action recognition try to distill information from both the spatial and temporal dimensions of an input video. Methods that reach State of the Art (SotA) accuracy, usually make use of 3D convolution layers as a way to abstract the temporal information from video frames. The use of such convolutions requires sampling short clips from the input video, where each clip is a collection of closely sampled frames. Since each short clip covers a small fraction of an input video, multiple clips are sampled at inference in order to cover the whole temporal length of the video. This leads to increased computational load and is impractical for real-world applications. We address the computational bottleneck by significantly reducing the number of frames required for inference. Our approach relies on a temporal transformer that applies global attention over video frames, and thus better exploits the salient information in each frame. Therefore our approach is very input efficient, and can achieve SotA results (on Kinetics dataset) with a fraction of the data (frames per video), computation and latency. Specifically on Kinetics-400, we reach $80.5$ top-1 accuracy with $\times 30$ less frames per video, and $\times 40$ faster inference than the current leading method. Code is available at: https://github.com/Alibaba-MIIL/STAM


title: The phase diagram of a frustrated Heisenberg model: from disorder to order and back again

abstract: We study the effects of bond and site disorder in the classical $J_{1}$-$J_{2}$ Heisenberg model on a square lattice in the order-by-disorder frustrated regime $2J_{2}>\left|J_{1}\right|$. Combining symmetry arguments, numerical energy minimization and large scale Monte Carlo simulations, we establish that the finite temperature Ising-like transition of the clean system is destroyed in the presence of any finite concentration of impurities. We explain this finding via a random-field mechanism which generically emerges in systems where disorder locally breaks the same real-space symmetry spontaneously globally broken by the associated order parameter. We also determine that the phase replacing the clean one is a paramagnet polarized in the nematic glass order with non-trivial magnetic response. This is because disorder also induces non-collinear spin-vortex-crystal order and produces a conjugated transverse dipolar random field. As a result of these many competing effects, the associated magnetic susceptibilities are non-monotonic functions of the temperature. As a further application of our methods, we show the generation of random axes in other frustrated magnets with broken SU(2) symmetry. We also discuss the generality of our findings and their relevance to experiments.

categories: cond-mat.str-el cond-mat.dis-nn cond-mat.stat-mech cond-mat.supr-con


title: Disentanglement-based Cross-Domain Feature Augmentation for Effective Unsupervised Domain Adaptive Person Re-identification

abstract: Unsupervised domain adaptive (UDA) person re-identification (ReID) aims to transfer the knowledge from the labeled source domain to the unlabeled target domain for person matching. One challenge is how to generate target domain samples with reliable labels for training. To address this problem, we propose a Disentanglement-based Cross-Domain Feature Augmentation (DCDFA) strategy, where the augmented features characterize well the target and source domain data distributions while inheriting reliable identity labels. Particularly, we disentangle each sample feature into a robust domain-invariant/shared feature and a domain-specific feature, and perform cross-domain feature recomposition to enhance the diversity of samples used in the training, with the constraints of cross-domain ReID loss and domain classification loss. Each recomposed feature, obtained based on the domain-invariant feature (which enables a reliable inheritance of identity) and an enhancement from a domain specific feature (which enables the approximation of real distributions), is thus an "ideal" augmentation. Extensive experimental results demonstrate the effectiveness of our method, which achieves the state-of-the-art performance.


title: A Characterization of Non-Iterative Normal Modal Logics

abstract: Non-iterative normal modal logics are defined by axioms of modal degree 1. In this paper we use calculations with normal forms to determine the set of all possible non-iterative normal modal logics, unimodal propositional extensions of K. To characterize them, we construct the complete set (lattice) of such logics and we provide the generic specification of their axioms and their semantics.


title: Is there an infant mortality in bacteria?

abstract: This manuscript proposes a significant step in our long-run investigation of infant mortality across species. Since 2016 (Berrut et al. 2016) a succession of studies (Bois et al. 2019) have traced infant mortality from organisms of high complexity (e.g. mammals) down to unicellular organisms. Infant mortality may be considered as a filtering process through which organisms with potentially lethal congenital defects are eliminated. Such defects may have many causes but here we focus particularly on mishaps resulting from non-optimal conditions in the production of proteins, enzymes and other crucial macromolecules. The statistical signature of infant mortality consists in a falling age-specific death rate. The question we address here is whether infant mortality episodes take place in bacteria in the minutes precededing or following cell division. It will be shown that while experiments carried out in the 20th century tried but failed to detect such an effect (mostly because of limited sample size), more recent observations provided consistent evidence of a sizeable mortality, with a rate of the order of 0.7 per 1,000 per hour, in the exponential growth phase of E. coli. A further crucial test will be to measure the age-specific, post-division death rate. An experiment is outlined for that purpose. It is based on the selection of stained cells through flow cytometry and the derivation of their ages at death from their sizes. If an infant mortality effect can be identified in E. coli it can be conjectured that a similar effect also exists in other unicellular organisms, both prokaryote and eukaryote.


title: Osteostatin potentiates the bioactivity of mesoporous glass scaffolds containing Zn2+ ions in human mesenchymal stem cell cultures

abstract: There is an urgent need of biosynthetic bone grafts with enhanced osteogenic capacity. In this study, we describe the design of hierarchical meso-macroporous 3D-scaffolds based on mesoporous bioactive glasses (MBGs), enriched with the peptide osteostatin and Zn2+ ions, and their osteogenic effect on human mesenchymal stem cells (hMSCs) as a preclinical strategy in bone regeneration. By using additive fabrication techniques, scaffolds exhibiting hierarchical porosity: mesopores , macropores and big channels, were prepared. These MBG scaffolds with or without osteostatin were evaluated in cell cultures of hMSCs. Zinc promoted hMSCs colonization (both the surface and inside) of MBG scaffolds. Moreover, Zn2+ ions and osteostatin together, but not independently, in the scaffolds were found to induce the osteoblast differentiation genes runt related transcription factor-2 (RUNX2) and alkaline phosphatase (ALP) in hMSCs after 7 d of culture in the absence of an osteogenic differentiation-promoting medium. These results add credence to the combined use of zinc and osteostatin as an effective strategy for bone regeneration applications.


title: The Resh Programming Language for Multirobot Orchestration

abstract: This paper describes Resh, a new, statically typed, interpreted programming language and associated runtime for orchestrating multirobot systems. The main features of Resh are: (1) It offloads much of the tedious work of programming such systems away from the programmer and into the language runtime; (2) It is based on a small set of temporal and locational operators; and (3) It is not restricted to specific robot types or tasks. The Resh runtime consists of three engines that collaborate to run a Resh program using the available robots in their current environment. This paper describes both Resh and its runtime and gives examples of its use.


title: ScanGAN360: A Generative Model of Realistic Scanpaths for 360$^{\circ}$ Images

abstract: Understanding and modeling the dynamics of human gaze behavior in 360$^\circ$ environments is a key challenge in computer vision and virtual reality. Generative adversarial approaches could alleviate this challenge by generating a large number of possible scanpaths for unseen images. Existing methods for scanpath generation, however, do not adequately predict realistic scanpaths for 360$^\circ$ images. We present ScanGAN360, a new generative adversarial approach to address this challenging problem. Our network generator is tailored to the specifics of 360$^\circ$ images representing immersive environments. Specifically, we accomplish this by leveraging the use of a spherical adaptation of dynamic-time warping as a loss function and proposing a novel parameterization of 360$^\circ$ scanpaths. The quality of our scanpaths outperforms competing approaches by a large margin and is almost on par with the human baseline. ScanGAN360 thus allows fast simulation of large numbers of virtual observers, whose behavior mimics real users, enabling a better understanding of gaze behavior and novel applications in virtual scene design.


title: Maximal Intervals of Decrease and Inflection Points for Node Reliability

abstract: The \textit{node reliability} of a graph $G$ is the probability that at least one node is operational and that the operational nodes can all communicate in the subgraph that they induce, given that the edges are perfectly reliable but each node operates independently with probability $p\in[0,1]$. We show that unlike many other notions of graph reliability, the number of maximal intervals of decrease in $[0,1]$ is unbounded, and that there can be arbitrarily many inflection points in the interval as well.


title: Computational Mechanism for the Effect of Psychosis Community Treatment: A Conceptual Review from Neurobiology to Social Interaction

abstract: The computational underpinnings of positive psychotic symptoms have recently received significant attention. Candidate mechanisms include some combination of maladaptive priors and reduced updating of these priors during perception. A potential benefit of models with such mechanisms is their ability to link multiple levels of explanation. This is key to improving how we understand the experience of psychosis. Moreover, it points us towards more comprehensive avenues for therapeutic research by providing a putative mechanism that could allow for the generation of new treatments from first principles. In order to demonstrate this, our conceptual paper will discuss the application of the insights from previous computational models to an important and complex set of evidence-based clinical interventions with strong social elements, such as coordinated specialty care clinics in early psychosis and assertive community treatment. These interventions may include but also go beyond psychopharmacology, providing, we argue, structure and predictability for patients experiencing psychosis. We develop the argument that this structure and predictability directly counteract the relatively low precision afforded to sensory information in psychosis, while also providing the patient more access to external cognitive resources in the form of providers and the structure of the programs themselves. We discuss how computational models explain the resulting reduction in symptoms, as well as the predictions these models make about potential responses of patients to modifications or to different variations of these interventions. We also link, via the framework of computational models, the experiences of patients and response to interventions to putative neurobiology.


title: An Irreversible Synthetic Route to an Ultra-Strong Two-Dimensional Polymer

abstract: Polymers that extend covalently in two dimensions have attracted recent attention as a means of combining the mechanical strength and in-plane energy conduction of conventional two-dimensional (2D) materials with the low densities, synthetic processability, and organic composition of their one-dimensional counterparts. Efforts to date have proven successful in forms that do not allow full realization of these properties, such as polymerization at flat interfaces or fixation of monomers in immobilized lattices. A frequently employed synthetic approach is to introduce microscopic reversibility, at the cost of bond stability, to achieve 2D crystals after extensive error correction. Herein we demonstrate a synthetic route to 2D irreversible polycondensation directly in the solution phase, resulting in covalently bonded 2D polymer platelets that are chemically stable and highly processable. Further fabrication offers highly oriented, free-standing films which exhibit exceptional 2D elastic modulus and yield strength at 50.9 +- 15.0 GPa and 0.976 +- 0.113 GPa, respectively. Platelet alignment is evidenced by polarized photoluminescence centered at 580 and 680 nm from different dipole transitions. This new synthetic route provides opportunities for 2D polymers in applications ranging from composite structures to molecular sieving membranes.


title: Gamma-convergent projection-free finite element methods for nematic liquid crystals: The Ericksen model

abstract: The Ericksen model for nematic liquid crystals couples a director field with a scalar degree of orientation variable, and allows the formation of various defects with finite energy. We propose a simple but novel finite element approximation of the problem that can be implemented easily within standard finite element packages. Our scheme is projection-free and thus circumvents the use of weakly acute meshes, which are quite restrictive in 3D but are required by recent algorithms for convergence. We prove stability and $\Gamma$-convergence properties of the new method in the presence of defects. We also design an effective nested gradient flow algorithm for computing minimizers that controls the violation of the unit-length constraint of the director. We present several simulations in 2D and 3D that document the performance of the proposed scheme and its ability to capture quite intriguing defects.


title: Casimir Interaction Between a Plane and a Sphere: Correction to the Proximity-Force Approximation at Intermediate Temperatures

abstract: We consider the Casimir interaction energy between a plane and a sphere of radius $R$ at finite temperature $T$ as a function of the distance of closest approach $L$. Typical experimental conditions are such that the thermal wavelength $\lambda_T=\hbar c/k_\mathrm{B}T$ satisfies the condition $L\ll \lambda_T\ll R$. We derive the leading correction to the proximity-force approximation valid for such intermediate temperatures by developing the scattering formula in the plane-wave basis. Our analytical result captures the joint effect of the spherical geometry and temperature and is written as a sum of temperature-dependent logarithmic terms. Surprisingly, two of the logarithmic terms arise from the Matsubara zero-frequency contribution.


title: Light-cone sum rules for proton decay

abstract: We estimate the form factors that parametrise the hadronic matrix elements of proton-to-pion transitions with the help of light-cone sum rules. These form factors are relevant for semi-leptonic proton decay channels induced by baryon-number violating dimension-six operators, as typically studied in the context of grand unified theories. We calculate the form factors in a kinematical regime where the momentum transfer from the proton to the pion is space-like and extrapolate our final results to the regime that is relevant for proton decay. In this way, we obtain estimates for the form factors that show agreement with the state-of-the-art calculations in lattice QCD, if systematic uncertainties are taken into account. Our work is a first step towards calculating more involved proton decay channels where lattice QCD results are not available at present.


title: Multinomial Logit Contextual Bandits: Provable Optimality and Practicality

abstract: We consider a sequential assortment selection problem where the user choice is given by a multinomial logit (MNL) choice model whose parameters are unknown. In each period, the learning agent observes a $d$-dimensional contextual information about the user and the $N$ available items, and offers an assortment of size $K$ to the user, and observes the bandit feedback of the item chosen from the assortment. We propose upper confidence bound based algorithms for this MNL contextual bandit. The first algorithm is a simple and practical method which achieves an $\tilde{\mathcal{O}}(d\sqrt{T})$ regret over $T$ rounds. Next, we propose a second algorithm which achieves a $\tilde{\mathcal{O}}(\sqrt{dT})$ regret. This matches the lower bound for the MNL bandit problem, up to logarithmic terms, and improves on the best known result by a $\sqrt{d}$ factor. To establish this sharper regret bound, we present a non-asymptotic confidence bound for the maximum likelihood estimator of the MNL model that may be of independent interest as its own theoretical contribution. We then revisit the simpler, significantly more practical, first algorithm and show that a simple variant of the algorithm achieves the optimal regret for a broad class of important applications.


title: Pattern Formation in Quantum Ferrofluids: from Supersolids to Superglasses

abstract: Pattern formation is a ubiquitous phenomenon observed in nonlinear and out-of-equilibrium systems. In equilibrium, quantum ferrofluids formed from ultracold atoms were recently shown to spontaneously develop coherent density patterns, manifesting a supersolid. We theoretically investigate the phase diagram of such quantum ferrofluids in oblate trap geometries and find an even wider range of exotic states of matter. Two-dimensional supersolid crystals formed from individual ferrofluid quantum droplets dominate the phase diagram at low densities. For higher densities we find honeycomb and labyrinthine states, as well as a pumpkin phase. We discuss scaling relations which allow us to find these phases for a wide variety of trap geometries, interaction strengths, and atom numbers. Our study illuminates the origin of the various possible patterns of quantum ferrofluids and shows that their occurrence is generic of strongly dipolar interacting systems stabilized by beyond mean-field effects.


title: Infinite Stable Graphs With Large Chromatic Number II

abstract: We prove a version of the strong Taylor's conjecture for stable graphs: if $G$ is a stable graph whose chromatic number is strictly greater than $\beth_2(\aleph_0)$ then $G$ contains all finite subgraphs of Sh$_n(\omega)$ and thus has elementary extensions of unbounded chromatic number. This completes the picture from our previous work. The main new model theoretic ingredient is a generalization of the classical construction of Ehrenfeucht-Mostowski models to an infinitary setting, giving a new characterization of stability.


title: QUOTAS: A new research platform for the data-driven investigation of black holes

abstract: We present QUOTAS, a novel research platform for the data-driven investigation of super-massive black hole (SMBH) populations. While SMBH data sets -- observations and simulations -- have grown rapidly in complexity and abundance, our computational environments and analysis tools have not matured commensurately to exhaust opportunities for discovery. Motivated to explore BH host galaxy and the parent dark matter halo connection, in this pilot version of QUOTAS, we assemble and co-locate the high-redshift, luminous quasar population at $z \geq 3$ alongside simulated data of the same epochs. Leveraging machine learning algorithms (ML) we expand simulation volumes that successfully replicate halo populations beyond the training set. Training ML on the Illustris-TNG300 simulation that includes baryonic physics, we populate the larger LEGACY Expanse dark matter-only box with quasars. Our first science results comparing observational and ML simulated quasars at $z \sim 3$, reveal that while the recovered Black Hole Mass Functions and clustering are in good agreement, simulated SMBHs fail to accrete, shine and grow at high enough rates to match observed quasars. We conclude that sub-grid models of mass accretion and SMBH feedback implemented in Illustris-TNG300 do not reproduce their observed mass growth. QUOTAS, demonstrates the power of ML, both for analyzing large complex datasets, and offering a unique opportunity to interrogate our theoretical model assumptions. We deploy ML again to derive and devise an optimal survey strategy for bringing the undetected lower luminosity quasar population into view. QUOTAS, and all related materials are publicly available at the Google Kaggle platform.

categories: astro-ph.CO astro-ph.GA astro-ph.HE astro-ph.IM


title: Unmanned Aerial Vehicle Visual Detection and Tracking using Deep Neural Networks: A Performance Benchmark

abstract: Unmanned Aerial Vehicles (UAV) can pose a major risk for aviation safety, due to both negligent and malicious use. For this reason, the automated detection and tracking of UAV is a fundamental task in aerial security systems. Common technologies for UAV detection include visible-band and thermal infrared imaging, radio frequency and radar. Recent advances in deep neural networks (DNNs) for image-based object detection open the possibility to use visual information for this detection and tracking task. Furthermore, these detection architectures can be implemented as backbones for visual tracking systems, thereby enabling persistent tracking of UAV incursions. To date, no comprehensive performance benchmark exists that applies DNNs to visible-band imagery for UAV detection and tracking. To this end, three datasets with varied environmental conditions for UAV detection and tracking, comprising a total of 241 videos (331,486 images), are assessed using four detection architectures and three tracking frameworks. The best performing detector architecture obtains an mAP of 98.6% and the best performing tracking framework obtains a MOTA of 96.3%. Cross-modality evaluation is carried out between visible and infrared spectrums, achieving a maximal 82.8% mAP on visible images when training in the infrared modality. These results provide the first public multi-approach benchmark for state-of-the-art deep learning-based methods and give insight into which detection and tracking architectures are effective in the UAV domain.


title: Transport efficiency of biofunctionalized magnetic particles tailored by surfactant concentration

abstract: Controlled transport of surface functionalized magnetic beads in a liquid medium is a central requirement for the handling of captured biomolecular targets in microfluidic lab-on-chip biosensors. Here, the influence of the physiological liquid medium on the transport characteristics of functionalized magnetic particles and on the functionality of the coupled protein is studied. These aspects are theoretically modeled and experimentally investigated for prototype superparamagnetic beads, surface functionalized with green fluorescent protein immersed in buffer solution with different concentrations of a surfactant. The model reports on the tunability of the steady-state particle substrate separation distance to prevent their surface sticking via the choice of surfactant concentration. Experimental and theoretical average velocities are discussed for a ratchet like particle motion induced by a dynamic external field superposed on a static locally varying magnetic field landscape. The developed model and experiment may serve as a basis for quantitative forecasts on the functionality of magnetic particle transport based lab-on-chip devices.


title: Near-optimal approximation methods for elliptic PDEs with lognormal coefficients

abstract: This paper studies numerical methods for the approximation of elliptic PDEs with lognormal coefficients of the form $-{\rm div}(a\nabla u)=f$ where $a=\exp(b)$ and $b$ is a Gaussian random field. The approximant of the solution $u$ is an $n$-term polynomial expansion in the scalar Gaussian random variables that parametrize $b$. We present a general convergence analysis of weighted least-squares approximants for smooth and arbitrarily rough random field, using a suitable random design, for which we prove optimality in the following sense: their convergence rate matches exactly or closely the rate that has been established in \cite{BCDM} for best $n$-term approximation by Hermite polynomials, under the same minimial assumptions on the Gaussian random field. This is in contrast with the current state of the art results for the stochastic Galerkin method that suffers the lack of coercivity due to the lognormal nature of the diffusion field. Numerical tests with $b$ as the Brownian bridge confirm our theoretical findings.


title: Some Fock spaces with depth two action

abstract: The subject of this paper are operators represented on a Fock space which act only on the two leading components of the tensor. We unify the constructions from arXiv:math/0702158, arXiv:0709.4334, arXiv:0812.0895, and arXiv:1003.2998, and extend a number of results from these papers to our more general setting. The results include the quadratic relation satisfied by (the kernel of) the free cumulant generating function, the resolvent form of the generating function for the Wick polynomials, and classification results for the case when the vacuum state on the operator algebra is tracial. We handle the generating functions in infinitely many variables by considering their matrix-valued versions.


title: CUDA Tutorial -- Cryptanalysis of Classical Ciphers Using Modern GPUs and CUDA

abstract: CUDA (formerly an abbreviation of Compute Unified Device Architecture) is a parallel computing platform and API model created by Nvidia allowing software developers to use a CUDA-enabled graphics processing unit (GPU) for general purpose processing. This 90-pages tutorial introduces the CUDA concepts in an easy-to-grasp and interactive way with ready-to-run code samples tested on Windows and Linux. Starting from scratch, a complete stand-alone GPU tool is implemented which automatically performs a ciphertext-only attack on ciphertexts encrypted by monoalphabetic substitution and columnar transposition. Throughout this process, you will learn how to architect the tool, what optimizations could significantly accelerate the routines, why the choice of an adequate metaheuristic is critical, and how to draw sketches to enlighten the design process. This tutorial will be incorporated in the CrypTool book as chapter 13.


title: Double nanowires for hybrid quantum devices

abstract: Parallel one-dimensional semiconductor channels connected by a superconducting strip constitute the core platform in several recent quantum device proposals that rely e.g. on Andreev processes or topological effects. In order to realize these proposals, the actual material systems must have high crystalline purity and the coupling between the different elements should be controllable in terms of their interfaces and geometry. We present a strategy for synthesizing double InAs nanowires by the vapor-liquid-solid mechanism using III-V molecular beam epitaxy. A superconducting layer is deposited onto nanowires without breaking vacuum, ensuring pristine interfaces between the superconductor and the two semiconductor nanowires. The method allows for a high yield of merged as well as separate parallel nanowires, with full or half-shell superconductor coatings. We demonstrate their utility in complex quantum devices by electron transport measurements.


title: Recovering Quantitative Models of Human Information Processing with Differentiable Architecture Search

abstract: The integration of behavioral phenomena into mechanistic models of cognitive function is a fundamental staple of cognitive science. Yet, researchers are beginning to accumulate increasing amounts of data without having the temporal or monetary resources to integrate these data into scientific theories. We seek to overcome these limitations by incorporating existing machine learning techniques into an open-source pipeline for the automated construction of quantitative models. This pipeline leverages the use of neural architecture search to automate the discovery of interpretable model architectures, and automatic differentiation to automate the fitting of model parameters to data. We evaluate the utility of these methods based on their ability to recover quantitative models of human information processing from synthetic data. We find that these methods are capable of recovering basic quantitative motifs from models of psychophysics, learning and decision making. We also highlight weaknesses of this framework and discuss future directions for their mitigation.


title: Reachability and Matching in Single Crossing Minor Free Graphs

abstract: We show that for each single crossing graph $H$, a polynomially bounded weight function for all $H$-minor free graphs $G$ can be constructed in Logspace such that it gives nonzero weights to all the cycles in $G$. This class of graphs subsumes almost all classes of graphs for which such a weight function is known to be constructed in Logspace. As a consequence, we obtain that for the class of $H$-minor free graphs where $H$ is a single crossing graph, reachability can be solved in UL, and bipartite maximum matching can be solved in SPL, which are small subclasses of the parallel complexity class NC. In the restrictive case of bipartite graphs, our maximum matching result improves upon the recent result of Eppstein and Vazirani, where they show an NC bound for constructing perfect matching in general single crossing minor free graphs.


title: SMILE: Self-Distilled MIxup for Efficient Transfer LEarning

abstract: To improve the performance of deep learning, mixup has been proposed to force the neural networks favoring simple linear behaviors in-between training samples. Performing mixup for transfer learning with pre-trained models however is not that simple, a high capacity pre-trained model with a large fully-connected (FC) layer could easily overfit to the target dataset even with samples-to-labels mixed up. In this work, we propose SMILE - Self-Distilled Mixup for EffIcient Transfer LEarning. With mixed images as inputs, SMILE regularizes the outputs of CNN feature extractors to learn from the mixed feature vectors of inputs (sample-to-feature mixup), in addition to the mixed labels. Specifically, SMILE incorporates a mean teacher, inherited from the pre-trained model, to provide the feature vectors of input samples in a self-distilling fashion, and mixes up the feature vectors accordingly via a novel triplet regularizer. The triple regularizer balances the mixup effects in both feature and label spaces while bounding the linearity in-between samples for pre-training tasks. Extensive experiments have been done to verify the performance improvement made by SMILE, in comparisons with a wide spectrum of transfer learning algorithms, including fine-tuning, L2-SP, DELTA, and RIFLE, even with mixup strategies combined. Ablation studies show that the vanilla sample-to-label mixup strategies could marginally increase the linearity in-between training samples but lack of generalizability, while SMILE significantly improve the mixup effects in both label and feature spaces with both training and testing datasets. The empirical observations backup our design intuition and purposes.


title: Visual Grounding Strategies for Text-Only Natural Language Processing

abstract: Visual grounding is a promising path toward more robust and accurate Natural Language Processing (NLP) models. Many multimodal extensions of BERT (e.g., VideoBERT, LXMERT, VL-BERT) allow a joint modeling of texts and images that lead to state-of-the-art results on multimodal tasks such as Visual Question Answering. Here, we leverage multimodal modeling for purely textual tasks (language modeling and classification) with the expectation that the multimodal pretraining provides a grounding that can improve text processing accuracy. We propose possible strategies in this respect. A first type of strategy, referred to as {\it transferred grounding} consists in applying multimodal models to text-only tasks using a placeholder to replace image input. The second one, which we call {\it associative grounding}, harnesses image retrieval to match texts with related images during both pretraining and text-only downstream tasks. We draw further distinctions into both strategies and then compare them according to their impact on language modeling and commonsense-related downstream tasks, showing improvement over text-only baselines.


title: Generalized Gibbs Ensemble of 2D CFTs with U(1) Charge from the AGT Correspondence

abstract: The Generalized Gibbs Ensemble (GGE) is relevant to understand the thermalization of quantum systems with an infinite set of conserved charges. In this work, we analyze the GGE partition function of 2D Conformal Field Theories (CFTs) with a U(1) charge and quantum Benjamin-Ono$_{2}$ (qBO$_{2}$) hierarchy charges. We use the Alday-Gaiotto-Tachikawa (AGT) correspondence to express the thermal trace in terms of the Alba-Fateev-Litvinov-Tarnopolskiy (AFLT) basis of descendants, which diagonalizes all charges. We analyze the GGE partition function in the thermodynamic semiclassical limit, including the first order quantum correction. We find that the equality between GGE averages and primary eigenvalues of the qBO$_{2}$ charges is attainable in the strict large $c$ limit and potentially violated at the subleading $1/c$ order. We also obtain the finite $c$ partition function when only the first non-trivial charge is turned on, expressed in terms of partial theta functions. Our results should be relevant to the eigenstate thermalization hypothesis for charged CFTs, Warped CFTs and effective field theory descriptions of condensed matter systems.


title: Preserve, Promote, or Attack? GNN Explanation via Topology Perturbation

abstract: Prior works on formalizing explanations of a graph neural network (GNN) focus on a single use case - to preserve the prediction results through identifying important edges and nodes. In this paper, we develop a multi-purpose interpretation framework by acquiring a mask that indicates topology perturbations of the input graphs. We pack the framework into an interactive visualization system (GNNViz) which can fulfill multiple purposes: Preserve,Promote, or Attack GNN's predictions. We illustrate our approach's novelty and effectiveness with three case studies: First, GNNViz can assist non expert users to easily explore the relationship between graph topology and GNN's decision (Preserve), or to manipulate the prediction (Promote or Attack) for an image classification task on MS-COCO; Second, on the Pokec social network dataset, our framework can uncover unfairness and demographic biases; Lastly, it compares with state-of-the-art GNN explainer baseline on a synthetic dataset.


title: Explicit asymptotic secret key rate of continuous-variable quantum key distribution with an arbitrary modulation

abstract: We establish an analytical lower bound on the asymptotic secret key rate of continuous-variable quantum key distribution with an arbitrary modulation of coherent states. Previously, such bounds were only available for protocols with a Gaussian modulation, and numerical bounds existed in the case of simple phase-shift-keying modulations. The latter bounds were obtained as a solution of convex optimization problems and our new analytical bound matches the results of Ghorai et al. (2019), up to numerical precision. The more relevant case of quadrature amplitude modulation (QAM) could not be analyzed with the previous techniques, due to their large number of coherent states. Our bound shows that relatively small constellation sizes, with say 64 states, are essentially sufficient to obtain a performance close to a true Gaussian modulation and are therefore an attractive solution for large-scale deployment of continuous-variable quantum key distribution. We also derive similar bounds when the modulation consists of arbitrary states, not necessarily pure.


title: A comparative analysis of local network similarity measurements: application to author citation networks

abstract: Understanding the evolution of paper and author citations is of paramount importance for the design of research policies and evaluation criteria that can promote and accelerate scientific discoveries. Recently many studies on the evolution of science have been conducted in the context of the emergent science of science field. While many studies have probed the link problem in citation networks, only a few works have analyzed the temporal nature of link prediction in author citation networks. In this study we compared the performance of 10 well-known local network similarity measurements to predict future links in author citations networks. Differently from traditional link prediction methods, the temporal nature of the predict links is relevant for our approach. Our analysis revealed interesting results. The Jaccard coefficient was found to be among the most relevant measurements. The preferential attachment measurement, conversely, displayed the worst performance. We also found that the extension of local measurements to their weighted version do not significantly improved the performance of predicting citations. Finally, we also found that a neural network approach summarizing the information from all 10 considered similarity measurements was not able to provide the highest prediction performance.


title: Structure Prediction of Epitaxial Inorganic Interfaces by Lattice and Surface Matching with Ogre

abstract: We present a new version of the Ogre open source Python package with the capability to perform structure prediction of epitaxial inorganic interfaces by lattice and surface matching. In the lattice matching step a scan over combinations of substrate and film Miller indices is performed to identify the domain-matched interfaces with the lowest mismatch. Subsequently, surface matching is conducted by Bayesian optimization to find the optimal interfacial distance and in-plane registry between the substrate and film. For the objective function, a geometric score function is proposed, based on the overlap and empty space between atomic spheres at the interface. The score function reproduces the results of density functional theory (DFT) at a fraction of the computational cost. The optimized interfaces are pre-ranked using a score function based on the similarity of the atomic environment at the interface to the bulk environment. Final ranking of the top candidate structures is performed with DFT. Ogre streamlines DFT calculations of interface energies and electronic properties by automating the construction of interface models. The application of Ogre is demonstrated for two interfaces of interest for quantum computing and spintronics, Al/InAs and Fe/InSb.


title: Composite Test inclusive of Benfords Law, Noise reduction and 0-1 Test for effective detection of Chaos in Rotor-Stator Rub

abstract: Segregating noise from chaos in dynamic systems has been one of the challenging work for the researchers across the globe due to their seemingly similar statistical properties. Even the most used tools such 0-1 test and Lyapunov exponents fail to distinguish chaos when signal is mixed with noise. This paper addresses the issue of segregating the dynamics in a rotor-stator rub system when the vibrations are subjected to different levels of noise. First, the limitation of 0-1 test in segregating chaos from signal mixed with noise has been established. Second, the underexplored Benfords Law and its application to the vibratory dynamical rotor-stator rub system has been introduced for the first time. Using the Benfords Law Compliance Test (BLCT), successful segregation of not only noise from chaos but also very low Signal to Noise Ratio (SNR) signals which are mainly stochastic has been achieved. The Euclidean Distance concept has been used to explore the scale-invariant probability distribution of systems that comply with Benfords Law to separate chaos from noise. Moreover, for moderate bands of noise in signals, we have shown that the Schreibers Nonlinear Noise Reduction technique works effectively in reducing the noise without damaging the dynamic properties of the system. Combining these individual layers (0-1 Test, BLCT and Noise reduction) on a rotor system, a Decision Tree based method to effectively segregate noise from chaos and identify the correct dynamics of any system with time series data set has been proposed.


title: Approximate GCD in Lagrange bases

abstract: For a pair of polynomials with real or complex coefficients, given in any particular basis, the problem of finding their GCD is known to be ill-posed. An answer is still desired for many applications, however. Hence, looking for a GCD of so-called approximate polynomials where this term explicitly denotes small uncertainties in the coefficients has received significant attention in the field of hybrid symbolic-numeric computation. In this paper we give an algorithm, based on one of Victor Ya. Pan, to find an approximate GCD for a pair of approximate polynomials given in a Lagrange basis. More precisely, we suppose that these polynomials are given by their approximate values at distinct known points. We first find each of their roots by using a Lagrange basis companion matrix for each polynomial, cluster the roots of each polynomial to identify multiple roots, and then "marry" the two polynomials to find their GCD. At no point do we change to the monomial basis, thus preserving the good conditioning properties of the original Lagrange basis. We discuss advantages and drawbacks of this method. The computational cost is dominated by the rootfinding step; unless special-purpose eigenvalue algorithms are used, the cost is cubic in the degrees of the polynomials. In principle, this cost could be reduced but we do not do so here.


title: Don't Disturb Me: Challenges of Interacting with SoftwareBots on Open Source Software Projects

abstract: Software bots are used to streamline tasks in Open Source Software (OSS) projects' pull requests, saving development cost, time, and effort. However, their presence can be disruptive to the community. We identified several challenges caused by bots in pull request interactions by interviewing 21 practitioners, including project maintainers, contributors, and bot developers. In particular, our findings indicate noise as a recurrent and central problem. Noise affects both human communication and development workflow by overwhelming and distracting developers. Our main contribution is a theory of how human developers perceive annoying bot behaviors as noise on social coding platforms. This contribution may help practitioners understand the effects of adopting a bot, and researchers and tool designers may leverage our results to better support human-bot interaction on social coding platforms.


title: Comment on Buss et al., Science 2021: An alternative, empirically-supported adjustment for sero-reversion yields a 10 percentage point lower estimate of the cumulative incidence of SARS-CoV-2 in Manaus by October 2020

abstract: The estimate of the cumulative incidence of SARS-CoV-2 of 76% in Manaus by October 2020 by Buss et al. relies on the assumption of an exponentially-declining probability of sero-reversion over time. We present an alternative, empirically-supported approach that is based on the observed dynamics of antibody titers in sero-positive cases. Through this approach we revise the cumulative incidence estimate to 66% (63.3% - 68.5%) by October 2020. This estimate has implications for the proximity to herd immunity and future estimates of fitness advantages of virus variants, such as the P.1 variant of concern. This methodology is also relevant for any other sero-survey analysis that has to adjust for sero-reversion.


title: Estimation of Closest In-Path Vehicle (CIPV) by Low-Channel LiDAR and Camera Sensor Fusion for Autonomous Vehicle

abstract: In autonomous driving, using a variety of sensors to recognize preceding vehicles in middle and long distance is helpful for improving driving performance and developing various functions. However, if only LiDAR or camera is used in the recognition stage, it is difficult to obtain necessary data due to the limitations of each sensor. In this paper, we proposed a method of converting the tracking data of vision into bird's eye view (BEV) coordinates using an equation that projects LiDAR points onto an image, and a method of fusion between LiDAR and vision tracked data. Thus, the newly proposed method was effective through the results of detecting closest in-path vehicle (CIPV) in various situations. In addition, even when experimenting with the EuroNCAP autonomous emergency braking (AEB) test protocol using the result of fusion, AEB performance is improved through improved cognitive performance than when using only LiDAR. In experimental results, the performance of the proposed method was proved through actual vehicle tests in various scenarios. Consequently, it is convincing that the newly proposed sensor fusion method significantly improves the ACC function in autonomous maneuvering. We expect that this improvement in perception performance will contribute to improving the overall stability of ACC.


title: Observations of 12.2 GHz methanol masers towards northern high-mass protostellar objects

abstract: Context. Class II methanol masers at 6.7 and 12.2 GHz occur close to high-mass young stellar objects (HMYSOs). When they are observed simultaneously, such studies may contribute to refining the characterisation of local physical conditions. Aims. We aim to search for the 12.2 GHz methanol emission in 6.7 GHz methanol masers that might have gone undetected in previous surveys of northern sky HMYSOs, mainly due to their variability. Contemporaneous observations of both transitions are used to refine the flux density ratio and examine the physical parameters. Methods. We observed a sample of 153 sites of 6.7 GHz methanol maser emission in the 12.2 GHz methanol line with the Torun 32 m radio telescope, using the newly built X-band receiver. Results. The 12.2 GHz methanol maser emission was detected in 36 HMYSOs, with 4 of them detected for the first time. The 6.7 GHz to 12.2 GHz flux density ratio for spectral features of the contemporaneously observed sources has a median value of 5.1, which is in agreement with earlier reports. The ratio differs significantly among the sources and for the periodic source G107.298+5.639 specifically, the ratio is weakly recurrent from cycle to cycle, but it generally reaches a minimum around the flare peak. This is consistent with the stochastic maser process, where small variations in the physical parameters along the maser path can significantly affect the ratio. A comparison of our data with historical results (from about ten years ago) implies significant (> 50%) variability for about 47% and 14% at 12.2 GHz and 6.7 GHz, respectively. This difference can be explained via the standard model of methanol masers.


title: Developing Apps for Researching the COVID-19 Pandemic with the TrackYourHealth Platform

abstract: Through lockdowns and other severe changes to daily life, almost everyone is affected by the COVID-19 pandemic. Scientists and medical doctors are - among others - mainly interested in researching, monitoring, and improving physical and mental health of the general population. Mobile health apps (mHealth), and apps conducting ecological momentary assessments (EMA) respectively, can help in this context. However, developing such mobile applications poses many challenges like costly software development efforts, strict privacy rules, compliance with ethical guidelines, local laws, and regulations. In this paper, we present TrackYourHealth (TYH), a highly configurable, generic, and modular mobile data collection and EMA platform, which enabled us to develop and release two mobile multi-platform applications related to COVID-19 in just a few weeks. We present TYH and highlight specific challenges researchers and developers of similar apps may also face, especially when developing apps related to the medical field.


title: Nonlinear Estimation for Position-Aided Inertial Navigation Systems

abstract: In this work we solve the position-aided 3D navigation problem using a nonlinear estimation scheme. More precisely, we propose a nonlinear observer to estimate the full state of the vehicle (position, velocity, orientation and gyro bias) from IMU and position measurements. The proposed observer does not introduce additional auxiliary states and is shown to guarantee semi-global exponential stability without any assumption on the acceleration of the vehicle. The performance of the observer is shown, through simulation, to overcome the state-of-the-art approach that assumes negligible accelerations.


title: Shadoks Approach to Low-Makespan Coordinated Motion Planning

abstract: This paper describes the heuristics used by the Shadoks team for the CG:SHOP 2021 challenge. This year's problem is to coordinate the motion of multiple robots in order to reach their targets without collisions and minimizing the makespan. It is a classical multi agent path finding problem with the specificity that the instances are highly dense in an unbounded grid. Using the heuristics outlined in this paper, our team won first place with the best solution to 202 out of 203 instances and optimal solutions to at least 105 of them. The main ingredients include several different strategies to compute initial solutions coupled with a heuristic called Conflict Optimizer to reduce the makespan of existing solutions.


title: The Conformal Anomaly Action to Fourth Order (4T) in $d=4$ in Momentum Space

abstract: We elaborate on the structure of the conformal anomaly effective action up to 4-th order, in an expansion in the gravitational fluctuations $(h)$ of the background metric, in the flat spacetime limit. For this purpose we discuss the renormalization of 4-point functions containing insertions of stress-energy tensors (4T), in conformal field theories in four spacetime dimensions with the goal of identifying the structure of the anomaly action. We focus on a separation of the correlator into its transverse/traceless and longitudinal components, applied to the trace and conservation Ward identities (WI) in momentum space. These are sufficient to identify, from their hierarchical structure, the anomaly contribution, without the need to proceed with a complete determination of all of its independent form factors. Renormalization induces sequential bilinear graviton-scalar mixings on single, double and multiple trace terms, corresponding to $R\square^{-1}$ interactions of the scalar curvature, with intermediate virtual massless exchanges. These dilaton-like terms couple to the conformal anomaly, as for the chiral anomalous WIs. We show that at 4T level a new traceless component appears after renormalization. We comment on future extensions of this result to more general backgrounds, with possible applications to non local cosmologies.


title: Pfaffian paired states for half-integer fractional quantum Hall effect

abstract: In this review the physics of Pfaffian paired states, in the context of fractional quantum Hall effect, is discussed using field-theoretical approaches. The Pfaffian states are prime examples of topological ($p$-wave) Cooper pairing and are characterized by non-Abelian statistics of their quasiparticles. Here we focus on conditions for their realization and competition among them at half-integer filling factors. Using the Dirac composite fermion description, in the presence of a mass term, we study the influence of Landau level mixing in selecting a particular Pfaffian state. While Pfaffian and anti-Pfaffian are selected when Landau level mixing is not strong, and can be taken into account perturbatively, the PH Pfaffian state requires non-perturbative inclusion of at least two Landau levels. Our findings, for small Landau level mixing, are in accordance with numerical investigations in the literature, and call for a non-perturbative approach in the search for PH Pfaffian correlations. We demonstrated that a method based on the Chern-Simons field-theoretical approach can be used to generate characteristic interaction pseudo-potentials for Pfaffian paired states.


title: Interaction energy between two separated charged spheres surrounded inside and outside by electrolyte

abstract: By using the recently generalized version of Newton Shell Theorem analytical equations are derived to calculate the electric interaction energy between two separated charged spheres surrounded outside and inside by electrolyte. This electric interaction energy is calculated as a function of the electrolyte ion concentration, temperature, distance between the spheres and size of the spheres. At the same distance between the spheres the absolute value of the interaction energy decreases with increasing electrolyte ion concentration and increases with increasing temperature. At zero electrolyte ion concentration the derived analytical equation transforms into the Coulomb equation. Finally, the analytical equation is generalized to calculate the electric interaction energy of N separated charged spheres surrounded by electrolyte.


title: Chiral magnetic domain walls under transverse fields: a semi-analytical model

abstract: An analytical model for the domain wall structure in ultrathin films with perpendicular easy axis and interfacial Dzyaloshinskii-Moriya interaction, submitted to an arbitrary in-plane magnetic field, is presented. Its solution is simplified to the numerical minimization of an analytic function of just one variable. The model predictions are compared to numerical micromagnetic simulations, using parameters of existing samples, revealing a very good agreement. Remaining differences are analyzed, and partly corrected. Differences with the predictions of the simplest model, usually found in the literature, in which only the domain wall moment's in-plane orientation can vary, are exemplified. The model allows accurate computations, as a function of in-plane field module and orientation, of the domain wall tension and width, quantities controlling the creep motion of domain walls in such films.


title: The uniqueness of elementary embeddings

abstract: Much of the theory of large cardinals beyond a measurable cardinal concerns the structure of elementary embeddings of the universe of sets into inner models. This paper seeks to answer the question of whether the inner model uniquely determines the elementary embedding.


title: Efficient quantum circuit simulation using a multi-qubit Bloch vector representation of density matrices

abstract: In the Bloch sphere picture, one finds the coefficients for expanding a single-qubit density operator in terms of the identity and Pauli matrices. A generalization to $n$ qubits via tensor products represents a density operator by a real vector of length $4^n$, conceptually similar to a statevector. Here, we study this approach for the purpose of quantum circuit simulation, including noise processes. The tensor structure leads to computationally efficient algorithms for applying circuit gates and performing few-qubit quantum operations. In view of variational circuit optimization, we study ``backpropagation'' through a quantum circuit and gradient computation based on this representation, and generalize our analysis to the Lindblad equation for modeling the (non-unitary) time evolution of a density operator.


title: Design of active network filters as hysteretic sensors

abstract: This work aims to propose and design a class of networks of coupled linear and nonlinear oscillators, in which short bursts of exogenous excitation result in sustained endogenous network activity that returns to a quiescent state only after a characteristic time and along a different path than when originally excited. The desired hysteretic behavior is obtained through the coupling of self-excited oscillations with purposely designed rate laws for slowly-varying nodal parameters, governed only by local interactions in the network. The proposed architecture and the sought dynamics take inspiration from complex biological systems that combine endogenous energy sources with a paradigm for distributed sensing and information processing. In this paper, the network design problem considers arbitrary topologies and investigates the dependence of the desired response on model parameters, as well as on the placement of a single nonlinear node in an otherwise linear network. Perturbation analysis in various asymptotic parameter limits is used to define the proposed internal dynamics. Parameter continuation techniques validate the asymptotic results numerically and demonstrate their robustness over finite ranges of parameter values. Both approaches suggest a nontrivial dependence of the optimal distribution of nonlinearity on the network topology.


title: Deep learning-based photoacoustic imaging of vascular network through thick porous media

abstract: Photoacoustic imaging (PAI) is a promising approach to realize in vivo transcranial cerebral vascular imaging. However, the strong attenuation and distortion of the photoacoustic wave caused by the thick porous skull greatly affect the imaging quality. In this study, we designed a convolutional neural network (CNN) with a U-Net architecture to extract the effective photoacoustic information hidden in the speckle patterns; obtained vascular network images datasets under porous media through simulation and experiment, and trained the network weights respectively. The results show that the proposed neural network can learn the mapping relationship between the speckle pattern and the target; extract the photoacoustic signals of some vessels submerged in noise to reconstruct high-quality images with a sharp outline of the vessel and clean background. Compared with the traditional photoacoustic reconstruction method, the deep learning-based reconstruction algorithms can achieve better performance, the mean absolute error (MAE), Structural SIMilarity (SSIM) and peak signal-to-noise ratio (PSNR) of reconstructed images have been greatly improved. In conclusion, the proposed neural network can effectively extract valid information from highly blurred speckle patterns for the rapid reconstruction of target images, which offers promising applications in transcranial photoacoustic imaging.


title: Putting a price on tenure

abstract: Government employees in Brazil are granted tenure after three years on the job. Firing a tenured government employee is all but impossible, so tenure is a big employee benefit. But exactly how big is it? In other words: how much money is tenure worth to a government employee in Brazil? No one has ever attempted to answer that question. I do that in this paper. I use a modified version of the Sharpe ratio to estimate what the risk-adjusted salaries of government workers should be. The difference between actual salary and risk-adjusted salary gives us an estimate of how much tenure is worth to each employee. I find that in the 2005-2019 period the monthly value of tenure was 3980 reais to the median federal government employee, 1971 reais to the median state government employee, and 500 reais to the median municipal government employee.