title: EPIE Dataset: A Corpus For Possible Idiomatic Expressions

abstract: Idiomatic expressions have always been a bottleneck for language comprehension and natural language understanding, specifically for tasks like Machine Translation(MT). MT systems predominantly produce literal translations of idiomatic expressions as they do not exhibit generic and linguistically deterministic patterns which can be exploited for comprehension of the non-compositional meaning of the expressions. These expressions occur in parallel corpora used for training, but due to the comparatively high occurrences of the constituent words of idiomatic expressions in literal context, the idiomatic meaning gets overpowered by the compositional meaning of the expression. State of the art Metaphor Detection Systems are able to detect non-compositional usage at word level but miss out on idiosyncratic phrasal idiomatic expressions. This creates a dire need for a dataset with a wider coverage and higher occurrence of commonly occurring idiomatic expressions, the spans of which can be used for Metaphor Detection. With this in mind, we present our English Possible Idiomatic Expressions(EPIE) corpus containing 25206 sentences labelled with lexical instances of 717 idiomatic expressions. These spans also cover literal usages for the given set of idiomatic expressions. We also present the utility of our dataset by using it to train a sequence labelling module and testing on three independent datasets with high accuracy, precision and recall scores.


title: Existence of global weak solutions to the compressible Ericksen-Leslie system in dimension one

abstract: We consider the compressible Ericksen-Leslie system of liquid crystal flows in one dimension. A global weak solution is constructed with initial density $\rho_0\geq 0$ and $\rho_0\in L^\gamma$ for $\gamma>1$.


title: Efficient nonparametric statistical inference on population feature importance using Shapley values

abstract: The true population-level importance of a variable in a prediction task provides useful knowledge about the underlying data-generating mechanism and can help in deciding which measurements to collect in subsequent experiments. Valid statistical inference on this importance is a key component in understanding the population of interest. We present a computationally efficient procedure for estimating and obtaining valid statistical inference on the Shapley Population Variable Importance Measure (SPVIM). Although the computational complexity of the true SPVIM scales exponentially with the number of variables, we propose an estimator based on randomly sampling only $\Theta(n)$ feature subsets given $n$ observations. We prove that our estimator converges at an asymptotically optimal rate. Moreover, by deriving the asymptotic distribution of our estimator, we construct valid confidence intervals and hypothesis tests. Our procedure has good finite-sample performance in simulations, and for an in-hospital mortality prediction task produces similar variable importance estimates when different machine learning algorithms are applied.


title: Lattice dynamics localization in low-angle twisted bilayer graphene

abstract: A low twist angle between the two stacked crystal networks in bilayer graphene enables self-organized lattice reconstruction with the formation of a periodic domain. This superlattice modulates the vibrational and electronic structures, imposing new rules for electron-phonon coupling and the eventual observation of strong correlation and superconductivity. Direct optical images of the crystal superlattice in reconstructed twisted bilayer graphene are reported here, generated by the inelastic scattering of light in a nano-Raman spectroscope. The observation of the crystallographic structure with visible light is made possible due to lattice dynamics localization, the images resembling spectral variations caused by the presence of strain solitons and topological points. The results are rationalized by a nearly-free-phonon model and electronic calculations that highlight the relevance of solitons and topological points, particularly pronounced for structures with small twist angles. We anticipate our discovery to play a role in understanding Jahn-Teller effects and electronic Cooper pairing, among many other important phonon-related effects, and it may be useful for characterizing devices in the most prominent platform for the field of twistronics.


title: Can tumor location on pre-treatment MRI predict likelihood of pseudo-progression versus tumor recurrence in Glioblastoma? A feasibility study

abstract: A significant challenge in Glioblastoma (GBM) management is identifying pseudo-progression (PsP), a benign radiation-induced effect, from tumor recurrence, on routine imaging following conventional treatment. Previous studies have linked tumor lobar presence and laterality to GBM outcomes, suggesting that disease etiology and progression in GBM may be impacted by tumor location. Hence, in this feasibility study, we seek to investigate the following question: Can tumor location on treatment-na\"ive MRI provide early cues regarding likelihood of a patient developing pseudo-progression versus tumor recurrence? In this study, 74 pre-treatment Glioblastoma MRI scans with PsP (33) and tumor recurrence (41) were analyzed. First, enhancing lesion on Gd-T1w MRI and peri-lesional hyperintensities on T2w/FLAIR were segmented by experts and then registered to a brain atlas. Using patients from the two phenotypes, we construct two atlases by quantifying frequency of occurrence of enhancing lesion and peri-lesion hyperintensities, by averaging voxel intensities across the population. Analysis of differential involvement was then performed to compute voxel-wise significant differences (p-value<0.05) across the atlases. Statistically significant clusters were finally mapped to a structural atlas to provide anatomic localization of their location. Our results demonstrate that patients with tumor recurrence showed prominence of their initial tumor in the parietal lobe, while patients with PsP showed a multi-focal distribution of the initial tumor in the frontal and temporal lobes, insula, and putamen. These preliminary results suggest that lateralization of pre-treatment lesions towards certain anatomical areas of the brain may allow to provide early cues regarding assessing likelihood of occurrence of pseudo-progression from tumor recurrence on MRI scans.


title: Partial Policy Iteration for L1-Robust Markov Decision Processes

abstract: Robust Markov decision processes (MDPs) allow to compute reliable solutions for dynamic decision problems whose evolution is modeled by rewards and partially-known transition probabilities. Unfortunately, accounting for uncertainty in the transition probabilities significantly increases the computational complexity of solving robust MDPs, which severely limits their scalability. This paper describes new efficient algorithms for solving the common class of robust MDPs with s- and sa-rectangular ambiguity sets defined by weighted $L_1$ norms. We propose partial policy iteration, a new, efficient, flexible, and general policy iteration scheme for robust MDPs. We also propose fast methods for computing the robust Bellman operator in quasi-linear time, nearly matching the linear complexity the non-robust Bellman operator. Our experimental results indicate that the proposed methods are many orders of magnitude faster than the state-of-the-art approach which uses linear programming solvers combined with a robust value iteration.


title: Symmetry Abstractions for Hybrid Systems and their Applications

abstract: A symmetry of a dynamical system is a map that transforms one trajectory to another trajectory. We introduce a new type of abstraction for hybrid automata based on symmetries. The abstraction combines different modes in a concrete automaton A, whose trajectories are related by symmetries, into a single mode in the abstract automaton B. The abstraction sets the guard and reset of an abstract edge to be the union of the symmetry-transformed guards and resets of the concrete edges. We establish the soundness of the abstraction using a forward simulation relation (FSR) and present several examples. Our abstraction results in simpler automata, that are more amenable for formal analysis and design. We illustrate an application of this abstraction in making reachability analysis faster and enabling unbounded time safety verification. We show how a fixed point of the reachable set computation of B can be used to answer reachability queries for A, even if the latter visits an infinite and unbounded sequences of modes. We present our implementation of the abstraction construction, the fixed point check, and the map that transforms abstract reachable sets to concrete ones in a software tool. Finally, we show the advantage of our method over existing ones, and the different aspects of our abstraction, in a sequence of experiments including scenarios with linear and nonlinear agents following waypoints.


title: Convergence of Meta-Learning with Task-Specific Adaptation over Partial Parameters

abstract: Although model-agnostic meta-learning (MAML) is a very successful algorithm in meta-learning practice, it can have high computational cost because it updates all model parameters over both the inner loop of task-specific adaptation and the outer-loop of meta initialization training. A more efficient algorithm ANIL (which refers to almost no inner loop) was proposed recently by Raghu et al. 2019, which adapts only a small subset of parameters in the inner loop and thus has substantially less computational cost than MAML as demonstrated by extensive experiments. However, the theoretical convergence of ANIL has not been studied yet. In this paper, we characterize the convergence rate and the computational complexity for ANIL under two representative inner-loop loss geometries, i.e., strongly-convexity and nonconvexity. Our results show that such a geometric property can significantly affect the overall convergence performance of ANIL. For example, ANIL achieves a faster convergence rate for a strongly-convex inner-loop loss as the number $N$ of inner-loop gradient descent steps increases, but a slower convergence rate for a nonconvex inner-loop loss as $N$ increases. Moreover, our complexity analysis provides a theoretical quantification on the improved efficiency of ANIL over MAML. The experiments on standard few-shot meta-learning benchmarks validate our theoretical findings.


title: $E^3$: Visual Exploration of Spatiotemporal Energy Demand

abstract: Understanding demand-side energy behaviour is critical for making efficiency responses for energy demand management. We worked closely with energy experts and identified the key elements of the energy demand problem including temporal and spatial demand and shifts in spatiotemporal demand. To our knowledge, no previous research has investigated the shifts in spatiotemporal demand. To fill this research gap, we propose a unified visual analytics approach to support exploratory demand analysis; we developed E3, a highly interactive tool that support users in making and verifying hypotheses through human-client-server interactions. A novel potential flow based approach was formalized to model shifts in energy demand and integrated into a server-side engine. Experts then evaluated and affirmed the usefulness of this approach through case studies of real-world electricity data. In the future, we will improve the modelling algorithm, enhance visualisation, and expand the process to support more forms of energy data.


title: Circuits: An abstract viewpoint

abstract: Our primary purpose is to isolate the abstract, mathematical properties of circuits -- both classical Boolean circuits and quantum circuits -- that are essential for their computational interpretation. A secondary purpose is to clarify the similarities and differences between the classical and quantum situations. The general philosophy in this note is to include the mathematically essential aspects of circuits but to omit any of the additional structures that are usually included for convenience. We shall, however, retain the assumption that circuits are finite; this assumption does no harm to the applicability of our approach and is necessary for some of our work.


title: Brane inflation: swampland criteria, TCC, and reheating predictions

abstract: We consider inflation in a five -dimensional space time with the inflaton field confined to live on a brane world. In this scenario, we study different types of potentials for the inflaton, discuss their observational consequences, and compare with data. We find that some class of potentials are in good agreement with observation and that the value of the inflaton field can be sub-Planckian. Moreover, we investigate the swampland criteria in this scenario and determine the consistency of the model with the conjectures. Doing so, we could determine models that simultaneously satisfy both observational data and swampland criteria. More constraints are applied by studying the reheating phase where the acceptable range for the reheating temperature imposes some bounds on the models. As the last step, the result of trans-Planckian censorship conjecture for the model is considered where it is shown the constraint of TCC will be very strong and it could be used to applied limit on the brane tension.


title: Nash Equilibrium Problems of Polynomials

abstract: This paper studies Nash equilibrium problems that are given by polynomial functions. We formulate efficient polynomial optimization problems for computing Nash equilibria. The Lasserre type Moment-SOS relaxations are used to solve them. Under generic assumptions, the method can find a Nash equilibrium if there is one. Moreover, it can find all Nash equilibria if there are finitely many ones of them. The method can also detect nonexistence if there is no Nash equilibrium.


title: The transition matrix between the Specht and $\mathfrak{sl}_3$ web bases is unitriangular with respect to shadow containment

abstract: Webs are planar graphs with boundary that describe morphisms in a diagrammatic representation category for $\mathfrak{sl}_k$. They are studied extensively by knot theorists because braiding maps provide a categorical way to express link diagrams in terms of webs, producing quantum invariants like the well-known Jones polynomial. One important question in representation theory is to identify the relationships between different bases; coefficients in the change-of-basis matrix often describe combinatorial, algebraic, or geometric quantities (like, e.g., Kazhdan-Lusztig polynomials). By "flattening" the braiding maps, webs can also be viewed as the basis elements of a symmetric-group representation. In this paper, we define two new combinatorial structures for webs: band diagrams and their one-dimensional projections, shadows, that measure depths of regions inside the web. As an application, we resolve an open conjecture that the change-of-basis between the so-called Specht basis and web basis of this symmetric-group representation is unitriangular for $\mathfrak{sl}_3$-webs. We do this using band diagrams and shadows to construct a new partial order on webs that is a refinement of the usual partial order. In fact, we prove that for $\mathfrak{sl}_2$-webs, our new partial order coincides with the tableau partial order on webs studied by the authors and others. We also prove that though the new partial order for $\mathfrak{sl}_3$-webs is a refinement of the previously-studied tableau order, the two partial orders do not agree for $\mathfrak{sl}_3$.


title: Convex cones spanned by regular polytopes

abstract: We study three families of polyhedral cones whose sections are regular simplices, cubes, and crosspolytopes. We compute solid angles and conic intrinsic volumes of these cones. We show that several quantities appearing in stochastic geometry can be expressed through these conic intrinsic volumes. A list of such quantities includes internal and external solid angles of regular simplices and crosspolytopes, the probability that a (symmetric) Gaussian random polytope or the Gaussian zonotope contains a given point, the expected number of faces of the intersection of a regular polytope with a random linear subspace passing through its centre, and the expected number of faces of the projection of a regular polytope onto a random linear subspace.


title: Stopper-Controller Games embedded in Single-Player Control Problems

abstract: In 2002, Benjamin Jourdain and Claude Martini discovered that for a class of payoff functions, the pricing problem for American options can be reduced to pricing of European options for an appropriately associated payoff, all within a Black-Scholes framework. This discovery has been investigated in great detail by S\"oren Christensen, Jan Kallsen and Matthias Lenga in a recent work in 2020. In the present work we prove that this phenomenon can be observed in a wider context, and even holds true in a setup of non-linear stochastic processes. We analyse this problem from both probabilistic and analytic viewpoints. In the classical situation, Jourdain and Martini used this method to approximate prices of American put options. The broader applicability now potentially covers non-linear frameworks such as model uncertainty and controller-and-stopper-games.


title: Barium Selective Chemosensing by Diazacrown Ether Naphthalimide Turn-on Fluorophores for Single Ion Barium Tagging

abstract: Single molecule fluorescence detection of barium is investigated for enhancing the sensitivity and robustness of a neutrinoless double beta decay ($0\nu\beta\beta$) search in $^{136}$Xe, the discovery of which would alter our understanding of the nature of neutrinos and the early history of the Universe. A key developmental step is the synthesis of barium-selective chemosensors capable of incorporation into ongoing experiments in high-pressure $^{136}$Xe gas. Here we report turn-on fluorescent naphthalimide chemosensors containing monoaza- and diaza-crown ethers as agents for single Ba$^{2+}$ detection. Monoaza-18-crown-6 ether naphthalimide sensors showed sensitivity primarily to Ba$^{2+}$ and Hg$^{2+}$, whereas two diaza-18-crown-6 ether naphthalimides revealed a desirable selectivity toward Ba$^{2+}$. Solution-phase fluorescence and NMR experiments support a photoinduced electron transfer mechanism enabling turn-on fluorescence sensing in the presence of barium ions. Changes in ion-receptor interactions enable effective selectivity between competitive barium, mercury, and potassium ions, with detailed calculations correctly predicting fluorescence responses. With these molecules, dry-phase single Ba$^{2+}$ ion imaging with turn-on fluorescence is realized using oil-free microscopy techniques. This represents a significant advance toward a practical method of single Ba$^{2+}$ detection within large volumes of $^{136}$Xe, plausibly enabling a background-free technique to search for the hypothetical process of $0\nu\beta\beta$.

categories: physics.ins-det hep-ex nucl-ex physics.chem-ph


title: Block orbital-selective Mott insulators: a spin excitation analysis

abstract: We present a comprehensive study of the spin excitations - as measured by the dynamical spin structure factor $S(q,\omega)$ - of the so-called block-magnetic state of low-dimensional orbital-selective Mott insulators. We realize this state via both a multi-orbital Hubbard model and a generalized Kondo-Heisenberg Hamiltonian. Due to various competing energy scales present in the models, the system develops periodic ferromagnetic islands of various shapes and sizes, which are antiferromagnetically coupled. The 2$\times$2 particular case was already found experimentally in the ladder material BaFe$_2$Se$_3$ that becomes superconducting under pressure. Here we discuss the electronic density as well as Hubbard and Hund coupling dependence of $S(q,\omega)$ using density matrix renormalization group method. Several interesting features were identified: (1) An acoustic (dispersive spin-wave) mode develops. (2) The spin-wave bandwidth establishes a new energy scale that is strongly dependent on the size of the magnetic island and becomes abnormally small for large clusters. (3) Optical (dispersionless spin excitation) modes are present for all block states studied here. In addition, a variety of phenomenological spin Hamiltonians have been investigated but none matches entirely our results that were obtained primarily at intermediate Hubbard $U$ strengths. Our comprehensive analysis provides theoretical guidance and motivation to crystal growers to search for appropriate candidate materials to realize the block states, and to neutron scattering experimentalists to confirm the exotic dynamical magnetic properties unveiled here, with a rich mixture of acoustic and optical features.


title: Testing higher-order quantum interference with many-particle states

abstract: Quantum theory permits interference between indistinguishable paths but, at the same time, restricts its order. Single-particle interference, for instance, is limited to the second order, that is, to pairs of single-particle paths. To date, all experimental efforts to search for higher-order interferences beyond those compatible with quantum mechanics have been based on such single-particle schemes. However, quantum physics is not bounded to single-particle interference. We here experimentally study many-particle higher-order interference using a two-photon-five-slit setup. We observe nonzero two-particle interference up to fourth order, corresponding to the interference of two distinct two-particle paths. We further show that fifth-order interference is restricted to $10^{-3}$ in the intensity-correlation regime and to $10^{-2}$ in the photon-correlation regime, thus providing novel bounds on higher-order quantum interference.


title: Task-agnostic Exploration in Reinforcement Learning

abstract: Efficient exploration is one of the main challenges in reinforcement learning (RL). Most existing sample-efficient algorithms assume the existence of a single reward function during exploration. In many practical scenarios, however, there is not a single underlying reward function to guide the exploration, for instance, when an agent needs to learn many skills simultaneously, or multiple conflicting objectives need to be balanced. To address these challenges, we propose the \textit{task-agnostic RL} framework: In the exploration phase, the agent first collects trajectories by exploring the MDP without the guidance of a reward function. After exploration, it aims at finding near-optimal policies for $N$ tasks, given the collected trajectories augmented with \textit{sampled rewards} for each task. We present an efficient task-agnostic RL algorithm, \textsc{UCBZero}, that finds $\epsilon$-optimal policies for $N$ arbitrary tasks after at most $\tilde O(\log(N)H^5SA/\epsilon^2)$ exploration episodes. We also provide an $\Omega(\log (N)H^2SA/\epsilon^2)$ lower bound, showing that the $\log$ dependency on $N$ is unavoidable. Furthermore, we provide an $N$-independent sample complexity bound of \textsc{UCBZero} in the statistically easier setting when the ground truth reward functions are known.


title: Mechanical Deformation Induced Continuously Variable Emission for Radiative Cooling

abstract: Passive radiative cooling drawing the heat energy of objects to the cold outer space through the atmospheric transparent window (8 um - 13 um) is significant for reducing the energy consumption of buildings. Daytime and nighttime radiative cooling have been extensively investigated in the past. However, radiative cooling which can continuously regulate its cooling temperature, like a valve, according to human need is rarely reported. In this study, we present a concept of reconfigurable photonic structure for the adaptive radiative cooling by continuously varying the emission spectra in the atmospheric window region. This is realized by the deformation of the one-dimensional PDMS grating and the nanoparticles embedded PDMS thin film when subjected to mechanical strain. The proposed structure reaches different stagnation temperatures under certain strains. A dynamic exchange between two different strains results in the fluctuation of the photonic structure's temperature around a set temperature.


title: Walk Message Passing Neural Networks and Second-Order Graph Neural Networks

abstract: The expressive power of message passing neural networks (MPNNs) is known to match the expressive power of the 1-dimensional Weisfeiler-Leman graph (1-WL) isomorphism test. To boost the expressive power of MPNNs, a number of graph neural network architectures have recently been proposed based on higher-dimensional Weisfeiler-Leman tests. In this paper we consider the two-dimensional (2-WL) test and introduce a new type of MPNNs, referred to as $\ell$-walk MPNNs, which aggregate features along walks of length $\ell$ between vertices. We show that $2$-walk MPNNs match 2-WL in expressive power. More generally, $\ell$-walk MPNNs, for any $\ell\geq 2$, are shown to match the expressive power of the recently introduced $\ell$-walk refinement procedure (W[$\ell$]). Based on a correspondence between 2-WL and W[$\ell$], we observe that $\ell$-walk MPNNs and $2$-walk MPNNs have the same expressive power, i.e., they can distinguish the same pairs of graphs, but $\ell$-walk MPNNs can possibly distinguish pairs of graphs faster than $2$-walk MPNNs. When it comes to concrete learnable graph neural network (GNN) formalisms that match 2-WL or W[$\ell$] in expressive power, we consider second-order graph neural networks that allow for non-linear layers. In particular, to match W[$\ell$] in expressive power, we allow $\ell-1$ matrix multiplications in each layer. We propose different versions of second-order GNNs depending on the type of features (i.e., coming from a countable set, or coming from an uncountable set) as this affects the number of dimensions needed to represent the features. Our results indicate that increasing non-linearity in layers by means of allowing multiple matrix multiplications does not increase expressive power. At the very best, it results in a faster distinction of input graphs.


title: Logic of Machine Learning

abstract: The main question is: why and how can we ever predict based on a finite sample? The question is not answered by statistical learning theory. Here, I suggest that prediction requires belief in "predictability" of the underlying dependence, and learning involves search for a hypothesis where these beliefs are violated the least given the observations. The measure of these violations ("errors") for given data, hypothesis and particular type of predictability beliefs is formalized as concept of incongruity in modal Logic of Observations and Hypotheses (LOH). I show on examples of many popular textbook learners (from hierarchical clustering to k-NN and SVM) that each of them minimizes its own version of incongruity. In addition, the concept of incongruity is shown to be flexible enough for formalization of some important data analysis problems, not considered as part of ML.


title: On the Inference of Soft Biometrics from Typing Patterns Collected in a Multi-device Environment

abstract: In this paper, we study the inference of gender, major/minor (computer science, non-computer science), typing style, age, and height from the typing patterns collected from 117 individuals in a multi-device environment. The inference of the first three identifiers was considered as classification tasks, while the rest as regression tasks. For classification tasks, we benchmark the performance of six classical machine learning (ML) and four deep learning (DL) classifiers. On the other hand, for regression tasks, we evaluated three ML and four DL-based regressors. The overall experiment consisted of two text-entry (free and fixed) and four device (Desktop, Tablet, Phone, and Combined) configurations. The best arrangements achieved accuracies of 96.15%, 93.02%, and 87.80% for typing style, gender, and major/minor, respectively, and mean absolute errors of 1.77 years and 2.65 inches for age and height, respectively. The results are promising considering the variety of application scenarios that we have listed in this work.


title: Settling and Clustering of Snow Particles in Atmospheric Turbulence

abstract: The effect of turbulence on snow precipitation is not incorporated into present weather forecasting models. Here we show evidence that turbulence is in fact a key influence on both fall speed and spatial distribution of settling snow. We consider three snowfall events under vastly different levels of atmospheric turbulence. We characterize the size and morphology of the snow particles, and we simultaneously image their velocity, acceleration, and relative concentration over vertical planes about 30 m2 in area. We find that turbulence-driven settling enhancement explains otherwise contradictory trends between the particle size and velocity. The estimates of the Stokes number and the correlation between vertical velocity and local concentration indicate that the enhanced settling is rooted in the preferential sweeping mechanism. When the snow vertical velocity is large compared to the characteristic turbulence velocity, the crossing trajectories effect results in strong accelerations. When the conditions of preferential sweeping are met, the concentration field is highly non-uniform and clustering appears over a wide range of scales. These clusters, identified for the first time in a naturally occurring flow, display the signature features seen in canonical settings: power-law size distribution, fractal-like shape, vertical elongation, and large fall speed that increases with the cluster size. These findings demonstrate that the fundamental phenomenology of particle-laden turbulence can be leveraged towards a better predictive understanding of snow precipitation and ground snow accumulation. They also demonstrate how environmental flows can be used to investigate dispersed multiphase flows at Reynolds numbers not accessible in laboratory experiments or numerical simulations.


title: Memory-Efficient Pipeline-Parallel DNN Training

abstract: Many state-of-the-art ML results have been obtained by scaling up the number of parameters in existing models. However, parameters and activations for such large models often do not fit in the memory of a single accelerator device; this means that it is necessary to distribute training of large models over multiple accelerators. In this work, we propose PipeDream-2BW, a system that supports memory-efficient pipeline parallelism. PipeDream-2BW uses a novel pipelining and weight gradient coalescing strategy, combined with the double buffering of weights, to ensure high throughput, low memory footprint, and weight update semantics similar to data parallelism. In addition, PipeDream-2BW automatically partitions the model over the available hardware resources, while respecting hardware constraints such as memory capacities of accelerators and interconnect topologies. PipeDream-2BW can accelerate the training of large GPT and BERT language models by up to 20$\times$ with similar final model accuracy.


title: A generalizable saliency map-based interpretation of model outcome

abstract: One of the significant challenges of deep neural networks is that the complex nature of the network prevents human comprehension of the outcome of the network. Consequently, the applicability of complex machine learning models is limited in the safety-critical domains, which incurs risk to life and property. To fully exploit the capabilities of complex neural networks, we propose a non-intrusive interpretability technique that uses the input and output of the model to generate a saliency map. The method works by empirically optimizing a randomly initialized input mask by localizing and weighing individual pixels according to their sensitivity towards the target class. Our experiments show that the proposed model interpretability approach performs better than the existing saliency map-based approaches methods at localizing the relevant input pixels. Furthermore, to obtain a global perspective on the target-specific explanation, we propose a saliency map reconstruction approach to generate acceptable variations of the salient inputs from the space of input data distribution for which the model outcome remains unaltered. Experiments show that our interpretability method can reconstruct the salient part of the input with a classification accuracy of 89%.


title: Temporal clustering network for self-diagnosing faults from vibration measurements

abstract: There is a need to build intelligence in operating machinery and use data analysis on monitored signals in order to quantify the health of the operating system and self-diagnose any initiations of fault. Built-in control procedures can automatically take corrective actions in order to avoid catastrophic failure when a fault is diagnosed. This paper presents a Temporal Clustering Network (TCN) capability for processing acceleration measurement(s) made on the operating system (i.e. machinery foundation, machinery casing, etc.), or any other type of temporal signals, and determine based on the monitored signal when a fault is at its onset. The new capability uses: one-dimensional convolutional neural networks (1D-CNN) for processing the measurements; unsupervised learning (i.e. no labeled signals from the different operating conditions and no signals at pristine vs. damaged conditions are necessary for training the 1D-CNN); clustering (i.e. grouping signals in different clusters reflective of the operating conditions); and statistical analysis for identifying fault signals that are not members of any of the clusters associated with the pristine operating conditions. A case study demonstrating its operation is included in the paper. Finally topics for further research are identified.


title: Origin-to-destination network flow with path preferences and velocity controls: a mean field game-like approach

abstract: In this paper we consider a mean field approach to modeling the agents flow over a transportation network. In particular, beside a standard framework of mean field games, with controlled dynamics by the agents and costs mass-distribution dependent, we also consider a path preferences dynamics obtained as a generalization of the so-called noisy best response dynamics. Such a preferences dynamics says the agents choose their path having access to global information about the network congestion state and based on the observation of the decision of the agents that have preceded. We prove the existence of a mean field equilibrium obtained as a fixed point of a map over a suitable set of time-varying mass-distributions, defined edge by edge in the network. We also address the case where the admissible set of controls is suitably bounded depending on the mass-distribution on the edge itself.


title: Solving the Order Batching and Sequencing Problem using Deep Reinforcement Learning

abstract: In e-commerce markets, on time delivery is of great importance to customer satisfaction. In this paper, we present a Deep Reinforcement Learning (DRL) approach for deciding how and when orders should be batched and picked in a warehouse to minimize the number of tardy orders. In particular, the technique facilitates making decisions on whether an order should be picked individually (pick-by-order) or picked in a batch with other orders (pick-by-batch), and if so with which other orders. We approach the problem by formulating it as a semi-Markov decision process and develop a vector-based state representation that includes the characteristics of the warehouse system. This allows us to create a deep reinforcement learning solution that learns a strategy by interacting with the environment and solve the problem with a proximal policy optimization algorithm. We evaluate the performance of the proposed DRL approach by comparing it with several batching and sequencing heuristics in different problem settings. The results show that the DRL approach is able to develop a strategy that produces consistent, good solutions and performs better than the proposed heuristics.


title: A technique for natural gauge boson masses

abstract: In this work, a novel mechanism for spontaneous symmetry breaking is presented. This mechanism avoids quadratic divergencies and is thus capable of addressing the hierarchy problem in gauge theories. Using the scale-dependent effective action $\Gamma_{k}$ minimally coupled to a gravitational sector, variational parameter setting is applied. This provides a mass and vacuum expectation value as a function of the constants arising in the low scale expansion of Newtons' and cosmological couplings. A comparison with experimental data, such as the Higgs mass, allows putting restrictions on these constants. With this generic approach one can compare with explicit candidates for an effective field theory of gravity. As an example, we use the asymptotic safety scenario, where we find restrictions on the matter content of the theory.


title: Online Algorithms for Weighted Paging with Predictions

abstract: In this paper, we initiate the study of the weighted paging problem with predictions. This continues the recent line of work in online algorithms with predictions, particularly that of Lykouris and Vassilvitski (ICML 2018) and Rohatgi (SODA 2020) on unweighted paging with predictions. We show that unlike unweighted paging, neither a fixed lookahead nor knowledge of the next request for every page is sufficient information for an algorithm to overcome existing lower bounds in weighted paging. However, a combination of the two, which we call the strong per request prediction (SPRP) model, suffices to give a 2-competitive algorithm. We also explore the question of gracefully degrading algorithms with increasing prediction error, and give both upper and lower bounds for a set of natural measures of prediction error.


title: On sparse connectivity, adversarial robustness, and a novel model of the artificial neuron

abstract: Deep neural networks have achieved human-level accuracy on almost all perceptual benchmarks. It is interesting that these advances were made using two ideas that are decades old: (a) an artificial neuron based on a linear summator and (b) SGD training. However, there are important metrics beyond accuracy: computational efficiency and stability against adversarial perturbations. In this paper, we propose two closely connected methods to improve these metrics on contour recognition tasks: (a) a novel model of an artificial neuron, a "strong neuron," with low hardware requirements and inherent robustness against adversarial perturbations and (b) a novel constructive training algorithm that generates sparse networks with $O(1)$ connections per neuron. We demonstrate the feasibility of our approach through experiments on SVHN and GTSRB benchmarks. We achieved an impressive 10x-100x reduction in operations count (10x when compared with other sparsification approaches, 100x when compared with dense networks) and a substantial reduction in hardware requirements (8-bit fixed-point math was used) with no reduction in model accuracy. Superior stability against adversarial perturbations (exceeding that of adversarial training) was achieved without any counteradversarial measures, relying on the robustness of strong neurons alone. We also proved that constituent blocks of our strong neuron are the only activation functions with perfect stability against adversarial attacks.


title: A Large-scale Empirical Analysis of Browser Fingerprints Properties for Web Authentication

abstract: Modern browsers give access to several attributes that can be collected to form a browser fingerprint. Although browser fingerprints have primarily been studied as a web tracking tool, they can contribute to improve the current state of web security by augmenting web authentication mechanisms. In this paper, we investigate the adequacy of browser fingerprints for web authentication. We make the link between the digital fingerprints that distinguish browsers, and the biological fingerprints that distinguish Humans, to evaluate browser fingerprints according to properties inspired by biometric authentication factors. These properties include their distinctiveness, their stability through time, their collection time, their size, and the accuracy of a simple verification mechanism. We assess these properties on a large-scale dataset of 4,145,408 fingerprints composed of 216 attributes, and collected from 1,989,365 browsers. We show that, by time-partitioning our dataset, more than 81.3% of our fingerprints are shared by a single browser. Although browser fingerprints are known to evolve, an average of 91% of the attributes of our fingerprints stay identical between two observations, even when separated by nearly 6 months. About their performance, we show that our fingerprints weigh a dozen of kilobytes, and take a few seconds to collect. Finally, by processing a simple verification mechanism, we show that it achieves an equal error rate of 0.61%. We enrich our results with the analysis of the correlation between the attributes, and of their contribution to the evaluated properties. We conclude that our browser fingerprints carry the promise to strengthen web authentication mechanisms.


title: Visual Chirality

abstract: How can we tell whether an image has been mirrored? While we understand the geometry of mirror reflections very well, less has been said about how it affects distributions of imagery at scale, despite widespread use for data augmentation in computer vision. In this paper, we investigate how the statistics of visual data are changed by reflection. We refer to these changes as "visual chirality", after the concept of geometric chirality - the notion of objects that are distinct from their mirror image. Our analysis of visual chirality reveals surprising results, including low-level chiral signals pervading imagery stemming from image processing in cameras, to the ability to discover visual chirality in images of people and faces. Our work has implications for data augmentation, self-supervised learning, and image forensics.


title: Local photodoping in monolayer MoS2

abstract: Inducing electrostatic doping in 2D materials by laser exposure (photodoping effect) is an exciting route to tune optoelectronic phenomena. However, there is a lack of investigation concerning in what respect the action of photodoping in optoelectronic devices is local. Here, we employ scanning photocurrent microscopy (SPCM) techniques to investigate how a permanent photodoping modulates the photocurrent generation in MoS2 transistors locally. We claim that the photodoping fills the electronic states in MoS2 conduction band, preventing the photon-absorption and the photocurrent generation by the MoS2 sheet. Moreover, by comparing the persistent photocurrent (PPC) generation of MoS2 on top of different substrates, we elucidate that the interface between the material used for the gate and the insulator (gate-insulator interface) is essential for the photodoping generation. Our work gives a step forward to the understanding of the photodoping effect in MoS2 transistors and the implementation of such an effect in integrated devices.


title: Remarks on cosmological bulk viscosity in different epochs

abstract: The intention of this paper is mainly two-fold. \textit{First}, we point out a striking numerical agreement between the bulk viscosity in the lepton era calculated by Husdal (2016) and our own calculations of the present-day bulk viscosity when the functional form is $\zeta\,\sim\sqrt{\rho}$. From a phenomenological point of view, we thus seem to have an ansatz for the viscosity which bridges the infancy of the Universe ($\sim 1$ s) with the present. This can also be looked upon as a kind of symmetry between the early-time cosmology and the present-day cosmology: it is quite remarkable that the kinetic theory-based bulk viscosity in the early universe and the experimentally-based bulk viscosity in the present universe can be covered by the same simple analytical formula. \textit{Second}, we consider the Kasner universe as a typical anisotropic model of Bianchi-type I, investigating whether this geometrical model is compatible with constant viscosity coefficients in the fluid. Perhaps surprisingly, the existence of a shear viscosity turns out to be incompatible with the Kasner model. By contrast, a bulk viscosity is non-problematic in the {\it isotropic} version of the model. In the special case of a Zel'dovich (stiff) fluid, the three equal exponents in the Kasner metric are even determined by the bulk viscosity alone, independent of the value of the fluid energy density. We also give a brief comparison with some other recent approaches to viscous cosmology.


title: Lattice symmetry and emergence of antiferromagnetic quantum Hall states

abstract: Strong local interaction in systems with non-trivial topological bands can stabilize quantum states such as magnetic topological insulators. We investigate the influence of the lattice symmetry on the possible emergence of antiferromagnetic quantum Hall states. We consider the spinful Harper-Hofstadter model extended by a next-nearest-neighbor (NNN) hopping which opens a gap at half-filling and allows for the realization of a quantum Hall insulator. The quantum Hall insulator has the Chern number $\mathcal{C}=2$ as both spin components are in the same quantum Hall state. We add to the system a staggered potential $\Delta$ along the $\hat{x}$-direction favoring a normal insulator and the Hubbard interaction $U$ favoring a Mott insulator. The Mott insulator is a N\'eel antiferromagnet for small and a stripe antiferromagnet for large NNN hopping. We investigate the $U$-$\Delta$ phase diagram of the model for both small and large NNN hoppings. We show that while for large NNN hopping there exists a $\mathcal{C}=1$ stripe antiferromagnetic quantum Hall insulator in the phase diagram, there is no equivalent $\mathcal{C}=1$ N\'eel antiferromagnetic quantum Hall insulator at the small NNN hopping. We discuss that a $\mathcal{C}=1$ antiferromagnetic quantum Hall insulator can emerge only if the effect of the spin-flip transformation cannot be compensated by a space group operation. Our findings can be used as a guideline in future investigations searching for antiferromagnetic quantum Hall states.


title: Growth of perturbations to the peaked periodic waves in the Camassa-Holm equation,

abstract: Peaked periodic waves in the Camassa-Holm equation are revisited. Linearized evolution equations are derived for perturbations to the peaked periodic waves and linearized instability is proven both in $H^1$ and $W^{1,\infty}$ norms. Dynamics of perturbations in $H^1$ is related to the existence of two conserved quantities and is bounded in the full nonlinear system due to these conserved quantities. On the other hand, perturbations to the peaked periodic wave grow in $W^{1,\infty}$ norm and may blow up in a finite time in the nonlinear evolution of the Camassa-Holm equation.

categories: math.AP math-ph math.DS math.MP nlin.PS nlin.SI


title: Linear Last-iterate Convergence in Constrained Saddle-point Optimization

abstract: Optimistic Gradient Descent Ascent (OGDA) and Optimistic Multiplicative Weights Update (OMWU) for saddle-point optimization have received growing attention due to their favorable last-iterate convergence. However, their behaviors for simple bilinear games over the probability simplex are still not fully understood - previous analysis lacks explicit convergence rates, only applies to an exponentially small learning rate, or requires additional assumptions such as the uniqueness of the optimal solution. In this work, we significantly expand the understanding of last-iterate convergence for OGDA and OMWU in the constrained setting. Specifically, for OMWU in bilinear games over the simplex, we show that when the equilibrium is unique, linear last-iterate convergence is achieved with a learning rate whose value is set to a universal constant, improving the result of (Daskalakis & Panageas, 2019b) under the same assumption. We then significantly extend the results to more general objectives and feasible sets for the projected OGDA algorithm, by introducing a sufficient condition under which OGDA exhibits concrete last-iterate convergence rates with a constant learning rate whose value only depends on the smoothness of the objective function. We show that bilinear games over any polytope satisfy this condition and OGDA converges exponentially fast even without the unique equilibrium assumption. Our condition also holds for strongly-convex-strongly-concave functions, recovering the result of (Hsieh et al., 2019). Finally, we provide experimental results to further support our theory.


title: Optimal Transport and Risk Aversion in Kyle's Model of Informed Trading

abstract: We establish connections between optimal transport theory and the dynamic version of the Kyle model, including new characterizations of informed trading profits via conjugate duality and Monge-Kantorovich duality. We use these connections to extend the model to multiple assets, general distributions, and risk-averse market makers. With risk-averse market makers, liquidity is lower, assets exhibit short-term reversals, and risk premia depend on market maker inventories, which are mean reverting. We illustrate the model by showing that implied volatilities predict stock returns when there is informed trading in stocks and options and market makers are risk averse.


title: Aligning with Heterogeneous Preferences for Kidney Exchange

abstract: AI algorithms increasingly make decisions that impact entire groups of humans. Since humans tend to hold varying and even conflicting preferences, AI algorithms responsible for making decisions on behalf of such groups encounter the problem of preference aggregation: combining inconsistent and sometimes contradictory individual preferences into a representative aggregate. In this paper, we address this problem in a real-world public health context: kidney exchange. The algorithms that allocate kidneys from living donors to patients needing transplants in kidney exchange matching markets should prioritize patients in a way that aligns with the values of the community they serve, but allocation preferences vary widely across individuals. In this paper, we propose, implement and evaluate a methodology for prioritizing patients based on such heterogeneous moral preferences. Instead of selecting a single static set of patient weights, we learn a distribution over preference functions based on human subject responses to allocation dilemmas, then sample from this distribution to dynamically determine patient weights during matching. We find that this methodology increases the average rank of matched patients in the sampled preference ordering, indicating better satisfaction of group preferences. We hope that this work will suggest a roadmap for future automated moral decision making on behalf of heterogeneous groups.


title: A Higher Weight Analogue of Ogg's Theorem on Weierstrass Points

abstract: For a positive integer $N$, we say that $\infty$ is a Weierstrass point on the modular curve $X_0(N)$ if there is a non-zero cusp form of weight $2$ on $\Gamma_0(N)$ which vanishes at $\infty$ to order greater than the genus of $X_0(N)$. If $p$ is a prime with $p \nmid N$, Ogg proved that $\infty $ is not a Weierstrass point on $X_0(pN)$ if the genus of $X_0(N)$ is $0$. We prove a similar result for even weights $k \geq 4$. We also study the space of weight $k$ cusp forms on $\Gamma_0(N)$ vanishing to order greater than the dimension.


title: From Ancient Contemplative Practice to the App Store: Designing a Digital Container for Mindfulness

abstract: Hundreds of popular mobile apps today market their ties to mindfulness. What activities do these apps support and what benefits do they claim? How do mindfulness teachers, as domain experts, view these apps? We first conduct an exploratory review of 370 mindfulness-related apps on Google Play, finding that mindfulness is presented primarily as a tool for relaxation and stress reduction. We then interviewed 15 U.S. mindfulness teachers from the therapeutic, Buddhist, and Yogic traditions about their perspectives on these apps. Teachers expressed concern that apps that introduce mindfulness only as a tool for relaxation neglect its full potential. We draw upon the experiences of these teachers to suggest design implications for linking mindfulness with further contemplative practices like the cultivation of compassion. Our findings speak to the importance of coherence in design: that the metaphors and mechanisms of a technology align with the underlying principles it follows.


title: Revisiting the Sulfur-Water Chemical System in the Middle Atmosphere of Venus

abstract: Sulfur-water chemistry plays an important role in the middle atmosphere of Venus. Ground based observations have found that simultaneously observed SO2 and H2O at ~64 km vary with time and are temporally anti-correlated. To understand these observations, we explore the sulfur-water chemical system using a one-dimensional chemistry-diffusion model. We find that SO2 and H2O mixing ratios above the clouds are highly dependent on mixing ratios of the two species at the middle cloud top (58 km). The behavior of sulfur-water chemical system can be classified into three regimes but there is no abrupt transition among these regimes. In particular, there is no bifurcation behavior as previously claimed. We also find that the SO2 self-shielding effect causes H2O above the clouds to respond to the middle cloud top in a non-monotonic fashion. Through comparison with observations, we find that mixing ratio variations at the middle cloud top can explain the observed variability of SO2 and H2O. The sulfur-water chemistry in the middle atmosphere is responsible for the H2O-SO2 anti-correlation at 64 km. Eddy transport change alone cannot explain the variations of both species. These results imply that variations of species abundance in the middle atmosphere are significantly influenced by the lower atmospheric processes. Continued ground-based measurements of the co-evolution of SO2 and H2O above the clouds and new spacecraft missions will be crucial for uncover the complicated processes underlying the interaction among the lower atmosphere, the clouds and the middle atmosphere of Venus.


title: Concentration Polarization and Metal Dendrite Initiation in Isolated Electrolyte Microchannels

abstract: Lithium metal penetrations through the liquid-electrolyte-wetted porous separator and solid electrolytes are a major safety concern of next-generation rechargeable metal batteries. The penetrations were frequently discovered to occur through only a few isolated channels, as revealed by "black spots" on both sides of the separator or electrolyte, which manifest a highly localized ionic flux or current density. Predictions of the penetration time have been infeasible due to the hidden and unclear dynamics in these penetration channels. Here, using the glass capillary cells, we investigate for the first time the unexpectedly sensitive influence of channel geometry on the concentration polarization and dendrite initiation processes. The characteristic time for the complete depletion of salt concentration on the surface of the advancing electrode, i.e. Sand's time, exhibits a nonlinear dependence on the curvature of the channel walls along the axial direction. While a positively deviated Sand's time scaling exponent can be used to infer a converging penetration area through the electrolyte, a negatively deviated scaling exponent suggests that diffusion limitation can be avoided in expanding channels, such that the fast-advancing tip-growing dendrites will not be initiated. The safety design of rechargeable metal batteries will benefit from considering the true local current densities and the conduction structures.


title: Hybrid Magneto Photonic Material Structure for Plasmon Assisted Magnetic Switching

abstract: We have proposed the use of surface plasmon resonances at the interface of hybrid magneto-photonic heterostructures [Opt. Mat. Exp., 7, 4316 (2017)] for all-optical control of the macroscopic spin orientation in nanostructures in fs time scales. This requires strong spin-photon coupling for the resonant enhancement of opto-magnetic fields, generated through the inverse Faraday effect, in magnetic nanostructures with perpendicular anisotropy. Here we report on the development of nm thick interlayers to control the growth orientation of hcp-Co alloys grown on refractory plasmonic materials to align the magnetic axis out-of-plane, thereby meeting key requirements for the realization of ultrafast magneto-photonic devices.


title: An algorithm for non-parametric estimation in state-space models

abstract: State-space models are ubiquitous in the statistical literature since they provide a flexible and interpretable framework for analyzing many time series. In most practical applications, the state-space model is specified through a parametric model. However, the specification of such a parametric model may require an important modeling effort or may lead to models which are not flexible enough to reproduce all the complexity of the phenomenon of interest. In such situations, an appealing alternative consists in inferring the state-space model directly from the data using a non-parametric framework. The recent developments of powerful simulation techniques have permitted to improve the statistical inference for parametric state-space models. It is proposed to combine two of these techniques, namely the Stochastic Expectation-Maximization (SEM) algorithm and Sequential Monte Carlo (SMC) approaches, for non-parametric estimation in state-space models. The performance of the proposed algorithm is assessed through simulations on toy models and an application to environmental data is discussed.


title: Cross-lingual Retrieval for Iterative Self-Supervised Training

abstract: Recent studies have demonstrated the cross-lingual alignment ability of multilingual pretrained language models. In this work, we found that the cross-lingual alignment can be further improved by training seq2seq models on sentence pairs mined using their own encoder outputs. We utilized these findings to develop a new approach -- cross-lingual retrieval for iterative self-supervised training (CRISS), where mining and training processes are applied iteratively, improving cross-lingual alignment and translation ability at the same time. Using this method, we achieved state-of-the-art unsupervised machine translation results on 9 language directions with an average improvement of 2.4 BLEU, and on the Tatoeba sentence retrieval task in the XTREME benchmark on 16 languages with an average improvement of 21.5% in absolute accuracy. Furthermore, CRISS also brings an additional 1.8 BLEU improvement on average compared to mBART, when finetuned on supervised machine translation downstream tasks.


title: q-Algebraic Equations, their power series solutions, and the asymptotic behavior of their coefficients

abstract: We give a systematic study of q-algebraic equations. The questions of existence, uniqueness and regularity of the solutions are solved in the space of grid-based Hahn series. The regularity is understood in terms of asymptotic behavior of coefficients, and is the main focus of this work. The results and algorithms are illustrated by many examples.


title: Isotropic Quot schemes of orthogonal bundles over a curve

abstract: We study the isotropic Quot schemes $IQ_e (V)$ parameterizing degree $e$ isotropic subsheaves of maximal rank of an orthogonal bundle $V$ over a curve. The scheme $IQ_e (V)$ contains a compactification of the space $IQ^o_e (V)$ of degree $e$ maximal isotropic subbundles, but behaves quite differently from the classical Quot scheme, and the Lagrangian Quot scheme in [6]. We observe that for certain topological types of $V$, the scheme $IQ_e (V)$ is empty for all $e$. In the remaining cases, for infinitely many $e$ there are irreducible components of $IQ_e (V)$ consisting entirely of nonsaturated subsheaves, and so $IQ_e (V)$ is strictly larger than the closure of $IQ^o_e (V)$. As our main result, we prove that for any orthogonal bundle $V$ and for $e \ll 0$, the closure $\overline{IQ^o_e (V)}$ of $IQ^o_e (V)$ is either empty or consists of one or two irreducible connected components, depending on $\deg(V)$ and $e$. In so doing, we also characterize the nonsaturated part of $\overline{IQ^o_e (V)}$ when $V$ has even rank.


title: Quality Management of Machine Learning Systems

abstract: In the past decade, Artificial Intelligence (AI) has become a part of our daily lives due to major advances in Machine Learning (ML) techniques. In spite of an explosive growth in the raw AI technology and in consumer facing applications on the internet, its adoption in business applications has conspicuously lagged behind. For business/mission-critical systems, serious concerns about reliability and maintainability of AI applications remain. Due to the statistical nature of the output, software 'defects' are not well defined. Consequently, many traditional quality management techniques such as program debugging, static code analysis, functional testing, etc. have to be reevaluated. Beyond the correctness of an AI model, many other new quality attributes, such as fairness, robustness, explainability, transparency, etc. become important in delivering an AI system. The purpose of this paper is to present a view of a holistic quality management framework for ML applications based on the current advances and identify new areas of software engineering research to achieve a more trustworthy AI.


title: Nonlinear dynamic analysis of an epidemiological model for COVID-19 including public behavior and government action

abstract: This paper is concerned with nonlinear modeling and analysis of the COVID-19 pandemic currently ravaging the planet. There are two objectives: to arrive at an appropriate model that captures the collected data faithfully, and to use that as a basis to explore the nonlinear behavior. We use a nonlinear SEIR (Susceptible, Exposed, Infectious & Removed) transmission model with added behavioral and government policy dynamics. We develop a genetic algorithm technique to identify key model parameters employing COVID19 data from South Korea. Stability, bifurcations and dynamic behavior are analyzed. Parametric analysis reveals conditions for sustained epidemic equilibria to occur. This work points to the value of nonlinear dynamic analysis in pandemic modeling and demonstrates the dramatic influence of social and government behavior on disease dynamics.


title: Bayesian Evaluation of User App Choices in the Presence of Risk Communication on Android Devices

abstract: In the age of ubiquitous technologies, security- and privacy-focused choices have turned out to be a significant concern for individuals and organizations. Risks of such pervasive technologies are extensive and often misaligned with user risk perception, thus failing to help users in taking privacy-aware decisions. Researchers usually try to find solutions for coherently extending trust into our often inscrutable electronic networked environment. To enable security- and privacy-focused decision-making, we mainly focused on the realm of the mobile marketplace, examining how risk indicators can help people choose more secure and privacy-preserving apps. We performed a naturalistic experiment with N=60 participants, where we asked them to select applications on Android tablets with accurate real-time marketplace data. We found that, in aggregate, app selections changed to be more risk-averse in the presence of user risk-perception-aligned visual indicators. Our study design and research propose practical and usable interactions that enable more informed, risk-aware comparisons for individuals during app selections. We include an explicit argument for the role of human decision-making during app selection, beyond the current trend of using machine learning to automate privacy preferences after selection during run-time.


title: BoMaNet: Boolean Masking of an Entire Neural Network

abstract: Recent work on stealing machine learning (ML) models from inference engines with physical side-channel attacks warrant an urgent need for effective side-channel defenses. This work proposes the first $\textit{fully-masked}$ neural network inference engine design. Masking uses secure multi-party computation to split the secrets into random shares and to decorrelate the statistical relation of secret-dependent computations to side-channels (e.g., the power draw). In this work, we construct secure hardware primitives to mask $\textit{all}$ the linear and non-linear operations in a neural network. We address the challenge of masking integer addition by converting each addition into a sequence of XOR and AND gates and by augmenting Trichina's secure Boolean masking style. We improve the traditional Trichina's AND gates by adding pipelining elements for better glitch-resistance and we architect the whole design to sustain a throughput of 1 masked addition per cycle. We implement the proposed secure inference engine on a Xilinx Spartan-6 (XC6SLX75) FPGA. The results show that masking incurs an overhead of 3.5\% in latency and 5.9$\times$ in area. Finally, we demonstrate the security of the masked design with 2M traces.


title: Decomposable Families of Itemsets

abstract: The problem of selecting a small, yet high quality subset of patterns from a larger collection of itemsets has recently attracted lot of research. Here we discuss an approach to this problem using the notion of decomposable families of itemsets. Such itemset families define a probabilistic model for the data from which the original collection of itemsets has been derived from. Furthermore, they induce a special tree structure, called a junction tree, familiar from the theory of Markov Random Fields. The method has several advantages. The junction trees provide an intuitive representation of the mining results. From the computational point of view, the model provides leverage for problems that could be intractable using the entire collection of itemsets. We provide an efficient algorithm to build decomposable itemset families, and give an application example with frequency bound querying using the model. Empirical results show that our algorithm yields high quality results.


title: Towards improving discriminative reconstruction via simultaneous dense and sparse coding

abstract: Discriminative features extracted from the sparse coding model have been shown to perform well for classification. Recent deep learning architectures have further improved reconstruction in inverse problems by considering new dense priors learned from data. We propose a novel dense and sparse coding model that integrates both representation capability and discriminative features. The model studies the problem of recovering a dense vector $\mathbf{x}$ and a sparse vector $\mathbf{u}$ given measurements of the form $\mathbf{y} = \mathbf{A}\mathbf{x}+\mathbf{B}\mathbf{u}$. Our first analysis proposes a geometric condition based on the minimal angle between spanning subspaces corresponding to the matrices $\mathbf{A}$ and $\mathbf{B}$ that guarantees unique solution to the model. The second analysis shows that, under mild assumptions, a convex program recovers the dense and sparse components. We validate the effectiveness of the model on simulated data and propose a dense and sparse autoencoder (DenSaE) tailored to learning the dictionaries from the dense and sparse model. We demonstrate that (i) DenSaE denoises natural images better than architectures derived from the sparse coding model ($\mathbf{B}\mathbf{u}$), (ii) in the presence of noise, training the biases in the latter amounts to implicitly learning the $\mathbf{A}\mathbf{x} + \mathbf{B}\mathbf{u}$ model, (iii) $\mathbf{A}$ and $\mathbf{B}$ capture low- and high-frequency contents, respectively, and (iv) compared to the sparse coding model, DenSaE offers a balance between discriminative power and representation.


title: Multipole Graph Neural Operator for Parametric Partial Differential Equations

abstract: One of the main challenges in using deep learning-based methods for simulating physical systems and solving partial differential equations (PDEs) is formulating physics-based data in the desired structure for neural networks. Graph neural networks (GNNs) have gained popularity in this area since graphs offer a natural way of modeling particle interactions and provide a clear way of discretizing the continuum models. However, the graphs constructed for approximating such tasks usually ignore long-range interactions due to unfavorable scaling of the computational complexity with respect to the number of nodes. The errors due to these approximations scale with the discretization of the system, thereby not allowing for generalization under mesh-refinement. Inspired by the classical multipole methods, we propose a novel multi-level graph neural network framework that captures interaction at all ranges with only linear complexity. Our multi-level formulation is equivalent to recursively adding inducing points to the kernel matrix, unifying GNNs with multi-resolution matrix factorization of the kernel. Experiments confirm our multi-graph network learns discretization-invariant solution operators to PDEs and can be evaluated in linear time.


title: Causal inference of brain connectivity from fMRI with $\psi$-Learning Incorporated Linear non-Gaussian Acyclic Model ($\psi$-LiNGAM)

abstract: Functional connectivity (FC) has become a primary means of understanding brain functions by identifying brain network interactions and, ultimately, how those interactions produce cognitions. A popular definition of FC is by statistical associations between measured brain regions. However, this could be problematic since the associations can only provide spatial connections but not causal interactions among regions of interests. Hence, it is necessary to study their causal relationship. Directed acyclic graph (DAG) models have been applied in recent FC studies but often encountered problems such as limited sample sizes and large number of variables (namely high-dimensional problems), which lead to both computational difficulty and convergence issues. As a result, the use of DAG models is problematic, where the identification of DAG models in general is nondeterministic polynomial time hard (NP-hard). To this end, we propose a $\psi$-learning incorporated linear non-Gaussian acyclic model ($\psi$-LiNGAM). We use the association model ($\psi$-learning) to facilitate causal inferences and the model works well especially for high-dimensional cases. Our simulation results demonstrate that the proposed method is more robust and accurate than several existing ones in detecting graph structure and direction. We then applied it to the resting state fMRI (rsfMRI) data obtained from the publicly available Philadelphia Neurodevelopmental Cohort (PNC) to study the cognitive variance, which includes 855 individuals aged 8-22 years. Therein, we have identified three types of hub structure: the in-hub, out-hub and sum-hub, which correspond to the centers of receiving, sending and relaying information, respectively. We also detected 16 most important pairs of causal flows. Several of the results have been verified to be biologically significant.


title: The automation of robust interatomic-force measurements

abstract: Interatomic-force measurements are regularly performed using frequency-modulation atomic force microscopy. This requires conversion of the observed shift in the resonant frequency of a force-sensing cantilever, to the actual force experienced by its tip. Recently, Sader et al. Nature Nanotechnology 13, 1088 (2018) showed that this force conversion can be unreliable and proposed the inflection point test to identify valid and robust force data. Efficient and user-friendly algorithms are required for its routine practical implementation, which currently do not exist. Here, we (1) advance the theoretical framework of the inflection point test, (2) develop the required efficient algorithms for its complete automation, and (3) demonstrate the utility of this automation by studying two experimental datasets. The principal outcome of this report is the development of user-friendly software that integrates this automation with a standard force conversion methodology. This software provides the enabling technology for practitioners to now seamlessly perform robust nanoscale and interatomic-force measurements.


title: Evaluating and Rewarding Teamwork Using Cooperative Game Abstractions

abstract: Can we predict how well a team of individuals will perform together? How should individuals be rewarded for their contributions to the team performance? Cooperative game theory gives us a powerful set of tools for answering these questions: the Characteristic Function (CF) and solution concepts like the Shapley Value (SV). There are two major difficulties in applying these techniques to real world problems: first, the CF is rarely given to us and needs to be learned from data. Second, the SV is combinatorial in nature. We introduce a parametric model called cooperative game abstractions (CGAs) for estimating CFs from data. CGAs are easy to learn, readily interpretable, and crucially allow linear-time computation of the SV. We provide identification results and sample complexity bounds for CGA models as well as error bounds in the estimation of the SV using CGAs. We apply our methods to study teams of artificial RL agents as well as real world teams from professional sports.


title: AdvMind: Inferring Adversary Intent of Black-Box Attacks

abstract: Deep neural networks (DNNs) are inherently susceptible to adversarial attacks even under black-box settings, in which the adversary only has query access to the target models. In practice, while it may be possible to effectively detect such attacks (e.g., observing massive similar but non-identical queries), it is often challenging to exactly infer the adversary intent (e.g., the target class of the adversarial example the adversary attempts to craft) especially during early stages of the attacks, which is crucial for performing effective deterrence and remediation of the threats in many scenarios. In this paper, we present AdvMind, a new class of estimation models that infer the adversary intent of black-box adversarial attacks in a robust and prompt manner. Specifically, to achieve robust detection, AdvMind accounts for the adversary adaptiveness such that her attempt to conceal the target will significantly increase the attack cost (e.g., in terms of the number of queries); to achieve prompt detection, AdvMind proactively synthesizes plausible query results to solicit subsequent queries from the adversary that maximally expose her intent. Through extensive empirical evaluation on benchmark datasets and state-of-the-art black-box attacks, we demonstrate that on average AdvMind detects the adversary intent with over 75% accuracy after observing less than 3 query batches and meanwhile increases the cost of adaptive attacks by over 60%. We further discuss the possible synergy between AdvMind and other defense methods against black-box adversarial attacks, pointing to several promising research directions.


title: COLREG-Compliant Collision Avoidance for Unmanned Surface Vehicle using Deep Reinforcement Learning

abstract: Path Following and Collision Avoidance, be it for unmanned surface vessels or other autonomous vehicles, are two fundamental guidance problems in robotics. For many decades, they have been subject to academic study, leading to a vast number of proposed approaches. However, they have mostly been treated as separate problems, and have typically relied on non-linear first-principles models with parameters that can only be determined experimentally. The rise of Deep Reinforcement Learning (DRL) in recent years suggests an alternative approach: end-to-end learning of the optimal guidance policy from scratch by means of a trial-and-error based approach. In this article, we explore the potential of Proximal Policy Optimization (PPO), a DRL algorithm with demonstrated state-of-the-art performance on Continuous Control tasks, when applied to the dual-objective problem of controlling an underactuated Autonomous Surface Vehicle in a COLREGs compliant manner such that it follows an a priori known desired path while avoiding collisions with other vessels along the way. Based on high-fidelity elevation and AIS tracking data from the Trondheim Fjord, an inlet of the Norwegian sea, we evaluate the trained agent's performance in challenging, dynamic real-world scenarios where the ultimate success of the agent rests upon its ability to navigate non-uniform marine terrain while handling challenging, but realistic vessel encounters.


title: Integers representable as differences of linear recurrence sequences

abstract: Let $\{U_n\}_{n \geqslant 0}$ and $\{G_m\}_{m \geqslant 0}$ be two linear recurrence sequences defined over the integers. We establish an asymptotic formula for the number of integers $c$ in the range $[-x, x]$ which can be represented as differences $ U_n - G_m$, when $x$ goes to infinity. In particular, the density of such integers is $0$.


title: iConViz: Interactive Visual Exploration of the Default Contagion Risk of Networked-Guarantee Loans

abstract: Groups of enterprises can serve as guarantees for one another and form complex networks when obtaining loans from commercial banks. During economic slowdowns, corporate default may spread like a virus and lead to large-scale defaults or even systemic financial crises. To help financial regulatory authorities and banks manage the risk associated with networked loans, we identified the default contagion risk, a pivotal issue in developing preventive measures, and established iConVis, an interactive visual analysis tool that facilitates the closed-loop analysis process. A novel financial metric, the contagion effect, was formulated to quantify the infectious consequences of guarantee chains in this type of network. Based on this metric, we designed and implement a series of novel and coordinated views that address the analysis of financial problems. Experts evaluated the system using real-world financial data. The proposed approach grants practitioners the ability to avoid previous ad hoc analysis methodologies and extend coverage of the conventional Capital Accord to the banking industry.


title: Data Driven Control with Learned Dynamics: Model-Based versus Model-Free Approach

abstract: This paper compares two different types of data-driven control methods, representing model-based and model-free approaches. One is a recently proposed method - Deep Koopman Representation for Control (DKRC), which utilizes a deep neural network to map an unknown nonlinear dynamical system to a high-dimensional linear system, which allows for employing state-of-the-art control strategy. The other one is a classic model-free control method based on an actor-critic architecture - Deep Deterministic Policy Gradient (DDPG), which has been proved to be effective in various dynamical systems. The comparison is carried out in OpenAI Gym, which provides multiple control environments for benchmark purposes. Two examples are provided for comparison, i.e., classic Inverted Pendulum and Lunar Lander Continuous Control. From the results of the experiments, we compare these two methods in terms of control strategies and the effectiveness under various initialization conditions. We also examine the learned dynamic model from DKRC with the analytical model derived from the Euler-Lagrange Linearization method, which demonstrates the accuracy in the learned model for unknown dynamics from a data-driven sample-efficient approach.


title: Supersymmetry of $\mathcal{PT}$- symmetric tridiagonal Hamiltonians

abstract: We extend the study of supersymmetric tridiagonal Hamiltonians to the case of non-Hermitian Hamiltonians with real or complex conjugate eigenvalues. We find the relation between matrix elements of the non-Hermitian Hamiltonian $H$ and its supersymmetric partner $H^{+}$ in a given basis. Moreover, the orthogonal polynomials in the eigenstate expansion problem attached to $H^{+}$ can be recovered from those polynomials arising from the same problem for $H$ with the help of kernel polynomials. Besides its generality, the developed formalism in this work is a natural home for using the numerically powerful Gauss quadrature techniques in probing the nature of some physical quantities such as the energy spectrum of $\mathcal{PT}$-symmetric complex potentials. Finally, we solve the shifted $\mathcal{PT}$-symmetric Morse oscillator exactly in the tridiagonal representation.


title: Go with the Flow: Adaptive Control for Neural ODEs

abstract: Despite their elegant formulation and lightweight memory cost, neural ordinary differential equations (NODEs) suffer from known representational limitations. In particular, the single flow learned by NODEs cannot express all homeomorphisms from a given data space to itself, and their static weight parameterization restricts the type of functions they can learn compared to discrete architectures with layer-dependent weights. Here, we describe a new module called neurally controlled ODE (N-CODE) designed to improve the expressivity of NODEs. The parameters of N-CODE modules are dynamic variables governed by a trainable map from initial or current activation state, resulting in forms of open-loop and closed-loop control, respectively. A single module is sufficient for learning a distribution on non-autonomous flows that adaptively drive neural representations. We provide theoretical and empirical evidence that N-CODE circumvents limitations of previous NODEs models and show how increased model expressivity manifests in several supervised and unsupervised learning problems. These favorable empirical results indicate the potential of using data- and activity-dependent plasticity in neural networks across numerous domains.


title: Instructive Review of Novel SFT with Non-interacting consituents "objects", and attempt Generalization to p-adic theory

abstract: We have constructed a new formalism for describing a situation with {\color{red} several (dual) strings} present at a time, a {\color{red} string field theory}, by means of a constituent / a strings from objects picture similar to, but importantly different from the "bits" by Charles Thorn\cite{Thorn}. Our "objects" (essentially the bits) represent rather a making a lattice in the {\color{red} light cone variables} on the string. The remarkable feature and simplicity of our formalism is, that the "objects" do {\color{red} NOT interact}, basically just run or sit trivially fixed. {\color{red} Scattering is a fake} in our formalism. This opens also up for hoping for generalizations inspired by hadrons with their partons all having Bjorken variable $x=0$, and thus infinitely many constituents. The p-adic string is an example.


title: Non-commutative deformations of perverse coherent sheaves and rational curves

abstract: We consider non-commutative deformations of sheaves on algebraic varieties. We develop some tools to determine parameter algebras of versal non-commutative deformations for partial simple collections and the structure sheaves of smooth rational curves. We apply them to universal flopping contractions of length 2 and higher. We confirm Donovan-Wemyss conjecture in the case of deformations of Laufer's flops.


title: Synthesis and characterization of bulk Nd1-xSrxNiO2 and Nd1-xSrxNiO3

abstract: The recent reports of superconductivity in Nd1-xSrxNiO2/SrTiO3 heterostructures have reinvigorated interest in potential superconductivity of low-valence nickelates. Synthesis of Ni1+-containing compounds is notoriously difficult. In the current work, a combined sol-gel combustion and high-pressure annealing technique was employed to prepare polycrystalline perovskite Nd1-xSrxNiO3 (x = 0, 0.1 and 0.2). Metal nitrates and metal acetates were used as starting materials, and the latter were found to be superior to the former in terms of safety and reactivity. The Nd1-xSrxNiO3 compounds were subsequently reduced to Nd1 xSrxNiO2 using calcium hydride in a sealed, evacuated quartz tube. To understand the synthesis pathway, the evolution from NdNiO3 to NdNiO2 was monitored using in-situ synchrotron X ray diffraction during the reduction process. Electrical transport properties were consistent with an insulator-metal transition occurring between x = 0 and 0.1 for Nd1-xSrxNiO3. Superconductivity was not observed in our bulk samples of Nd1-xSrxNiO2. Neutron diffraction experiments at 3 K and 300 K were performed on Nd0.9Sr0.1NiO2, in which no magnetic Bragg reflections were observed, and the results of structural Rietveld refinement are provided.


title: Learning to Learn with Feedback and Local Plasticity

abstract: Interest in biologically inspired alternatives to backpropagation is driven by the desire to both advance connections between deep learning and neuroscience and address backpropagation's shortcomings on tasks such as online, continual learning. However, local synaptic learning rules like those employed by the brain have so far failed to match the performance of backpropagation in deep networks. In this study, we employ meta-learning to discover networks that learn using feedback connections and local, biologically inspired learning rules. Importantly, the feedback connections are not tied to the feedforward weights, avoiding biologically implausible weight transport. Our experiments show that meta-trained networks effectively use feedback connections to perform online credit assignment in multi-layer architectures. Surprisingly, this approach matches or exceeds a state-of-the-art gradient-based online meta-learning algorithm on regression and classification tasks, excelling in particular at continual learning. Analysis of the weight updates employed by these models reveals that they differ qualitatively from gradient descent in a way that reduces interference between updates. Our results suggest the existence of a class of biologically plausible learning mechanisms that not only match gradient descent-based learning, but also overcome its limitations.


title: An encoder-decoder deep surrogate for reverse time migration in seismic imaging under uncertainty

abstract: Seismic imaging faces challenges due to the presence of several uncertainty sources. Uncertainties exist in data measurements, source positioning, and subsurface geophysical properties. Reverse time migration (RTM) is a high-resolution depth migration approach useful for extracting information such as reservoir localization and boundaries. RTM, however, is time-consuming and data-intensive as it requires computing twice the wave equation to generate and store an imaging condition. RTM, when embedded in an uncertainty quantification algorithm (like the Monte Carlo method), shows a many-fold increase in its computational complexity due to the high input-output dimensionality. In this work, we propose an encoder-decoder deep learning surrogate model for RTM under uncertainty. Inputs are an ensemble of velocity fields, expressing the uncertainty, and outputs the seismic images. We show by numerical experimentation that the surrogate model can reproduce the seismic images accurately, and, more importantly, the uncertainty propagation from the input velocity fields to the image ensemble.


title: Time-dependent properties of interacting active matter: dynamical behavior of one-dimensional systems of self-propelled particles

abstract: We study an interacting high-density one-dimensional system of self-propelled particles described by the Active Ornstein-Uhlenbeck particle (AOUP) model where, even in the absence of alignment interactions, velocity and energy domains spontaneously form in analogy with those already observed in two dimensions. Such domains are regions where the individual velocities are spatially correlated as a result of the interplay between self-propulsion and interactions. Their typical size is controlled by a characteristic correlation length. In the present work, we focus on a novel and lesser-known aspect of the model, namely its dynamical behavior. To this purpose, we investigate theoretically and numerically the time-dependent velocity autocorrelation and spatio-temporal velocity correlation functions. The study of these correlations provides a measure of the average life-time and, thus, the stability in time of the velocity domains.


title: An Update on the Student Exoplanet Programme

abstract: An update is given on the exoplanet research collaboration between Nielsen (a marketing research company), Brigham Young University, and NZ universities with the National University of Singapore, which has been expanded to include a community college in the US. Key achievements from the past year are outlined, including density estimates for HD 209458 and Kepler 1 from radial velocity and transit fits. A comparison between the WinFitter optimizer and other techniques is outlined, showing that WinFitter estimated statistical errors are essentially in line (bar a scaling proportion) with those estimated via Markov Chain Monte Carlo techniques.


title: Note on Thermalization of Non-resonantly Produced Sterile Neutrinos

abstract: Using an analytic treatment, we discuss the parameter regions of large active-sterile neutrino mixing angles where sterile neutrinos produced in non-resonant flavour oscillations can approach thermalization in several cosmologies. We show that thermalization affects only large active-sterile neutrino mixing already rejected by different limits. Hence, the allowed sterile neutrino parameter regions are unaffected.


title: Isometric Graph Neural Networks

abstract: Many tasks that rely on representations of nodes in graphs would benefit if those representations were faithful to distances between nodes in the graph. Geometric techniques to extract such representations have poor scaling over large graph size, and recent advances in Graph Neural Network (GNN) algorithms have limited ability to reflect graph distance information beyond the first degree neighborhood. To enable this highly desired capability, we propose a technique to learn Isometric Graph Neural Networks (IGNN), which requires changing the input representation space and loss function to enable any GNN algorithm to generate representations that reflect distances between nodes. We experiment with the isometric technique on several GNN architectures for modeling multiple prediction tasks on multiple datasets. In addition to an improvement in AUC-ROC as high as $43\%$ in these experiments, we observe a consistent and substantial improvement as high as 400% in Kendall's Tau (KT), a measure that directly reflects distance information, demonstrating that the learned embeddings do account for graph distances.


title: Stabilization of one-dimensional Townes solitons by spin-orbit coupling in a dual-core system

abstract: It was recently demonstrated that two-dimensional Townes solitons (TSs) in two-component systems with cubic self-focusing, which are normally made unstable by the critical collapse, can be stabilized by linear spin-orbit coupling (SOC), in Bose-Einstein condensates and optics alike. We demonstrate that one-dimensional TSs, realized as optical spatial solitons in a planar dual-core waveguide with dominant quintic self-focusing, may be stabilized by SOC-like terms emulated by obliquity of the coupling between cores of the waveguide. Thus, SOC offers a universal mechanism for the stabilization of the TSs. A combination of systematic numerical considerations and analytical approximations identifies a vast stability area for skew-symmetric solitons in the system's main (semi-infinite) and annex (finite) bandgaps. Tilted ("moving") solitons are unstable, spontaneously evolving into robust breathers. For broad solitons, diffraction, represented by second derivatives in the system, may be neglected, leading to a simplified model with a finite bandgap. It is populated by skew-antisymmetric gap solitons, which are nearly stable close to the gap's bottom.


title: On conformable fractional Legendre polynomials and their convergence properties with applications

abstract: The main objective of this paper is to give a wide study on the conformable fractional Legendre polynomials (CFLPs). This study is assumed to be a generalization and refinement, in an easy way, of the scalar case into the context of the conformable fractional differentiation. We introduce the CFLPs via different generating functions and provide some of their main properties and convergence results. Subsequently, some pure recurrence and differential recurrence relations, Laplace's first integral formula and orthogonal properties are then developed for CFLPs. We append our study with presenting shifted CFLPs and describing applicable scheme using the collocation method to solve some fractional differential equations (FDEs) in the sense of conformable derivative. Some useful examples of FDEs are treated to support our theoretical results and examining their exact solutions. To the best of our knowledge, the obtained results are newly presented and could enrich the fractional theory of special functions.


title: The $L^1$-contraction principle in optimal transport

abstract: In this work, we use the JKO scheme to approximate a general class of diffusion problems generated by Darcy's law. Although the scheme is now classical, if the energy density is spatially inhomogeneous or irregular, many standard methods fail to apply to establish convergence in the continuum limit. To overcome these difficulties, we analyze the scheme through its dual problem and establish a novel $L^1$-contraction principle for the density variable. Notably, the contraction principle relies only on the existence of an optimal transport map and the convexity structure of the energy. As a result, the principle holds in a very general setting, and opens the door to using optimal-transport-based variational schemes to study a larger class of non-linear inhomogeneous parabolic equations.


title: Darcy's Law with a Source term

abstract: We introduce a novel variant of the JKO scheme to approximate Darcy's law with a pressure dependent source term. By introducing a new variable that implicitly controls the source term, our scheme is still able to use the standard Wasserstein-2-metric even though the total mass changes over time. Leveraging the dual formulation of our scheme, we show that the discrete-in-time approximations satisfy many useful properties expected for the continuum solutions, such as a comparison principle and uniform $L^1$-equicontinuity. Many of these properties are new even in the well-understood case where the growth term is absent. Finally, we show that our discrete approximations converge to a solution of the corresponding PDE system, including a tumor growth model with a general nonlinear source term.


title: The Upgrade I of LHCb VELO -- towards an intelligent monitoring platform

abstract: The Large Hadron Collider beauty (LHCb) detector is designed to detect decays of b- and c- hadrons for the study of CP violation and rare decays. At the end of the LHC Run 2, many of the LHCb measurements remained statistically dominated. In order to increase the trigger yield for purely hadronic channels, the hardware trigger will be removed, and the detector will be read out at 40 MHz. This, in combination with the five-fold increase in luminosity, requires radical changes to LHCb's electronics, and, in some cases, the replacement of entire sub-detectors with state-of-the-art detector technologies. The Vertex Locator (VELO) surrounding the interaction region is used to reconstruct the collision points (primary vertices) and decay vertices of long-lived particles (secondary vertices). The upgraded VELO will be composed of 52 modules placed along the beam axis divided into two retractable halves. The modules will each be equipped with 4 silicon hybrid pixel tiles, each read out by 3 VeloPix ASICs. The total output data rate anticipated for the whole detector will be around 1.6 Tbit/s. The highest occupancy ASICs will have pixel hit rates of approximately 900 Mhit/s, with the corresponding output data rate of 15 Gbit/s. The LHCb upgrade detector will be the first detector to read out at the full LHC rate of 40 MHz. The VELO upgrade will utilize the latest detector technologies to read out at this rate while maintaining the required radiation-hard profile and minimizing the detector material.


title: Airfoil Lift Calculation Using Wind Tunnel Wall Pressures

abstract: An experimental method to calculate lift using static pressure ports on the wind tunnel walls and its associated limits has been explored in this paper. While the wall-pressure measurement (WPM) technique for lift calculation has been implemented by other researchers, there is a lack of literature on the sensitivity of the WPM method to test section size, airfoil chord, and model thickness. Chord sensitivity studies showed that the airfoil chord plays an important role in the accuracy of the measurements and needs to be appropriately sized for a given test section dimensions for optimum performance of the WPM method. A chord sensitivity parameter ($CSP$) was formulated and a lower limit ($=0.025$) was established to relate the ideal chord-length to wind tunnel test-section dimensions to ensure best lift measuring capabilities. Finally, a combination of symmetric and cambered airfoils with thicknesses varying from $6\%-21\%$ were tested and successfully validated against reference data for a freestream chord Reynolds number range of 100,000 to 550,000. The WPM method was found to be sensitive to varying surface flow conditions and airfoil thickness and has been shown to be a viable replacement to traditional lift measurement techniques using load balances or surface pressure ports.


title: Is perturbative study of ground-state correlations valid?

abstract: Perturbative approaches have often been used to include the effects of ground-state correlations in extended theories of the random-phase approximation. Validity of such approaches is investigated for a solvable model where comparison with exact solutions can be made. It is pointed out that there is a case where perturbative approaches give good results in spite of the fact that interaction strength is far beyond a perturbative region.