title: Volumetric parcellation of the right ventricle for regional geometric and functional assessment

abstract: 3D echocardiography is an increasingly popular tool for assessing cardiac remodelling in the right ventricle (RV). It allows quantification of the cardiac chambers without any geometric assumptions, which is the main weakness of 2D echocardiography. However, regional quantification of geometry and function is limited by the lower spatial and temporal resolution and the scarcity of identifiable anatomical landmarks. We developed a technique for regionally assessing the 3 relevant RV regions: apical, inlet and outflow. The method's inputs are end-diastolic (ED) and end-systolic (ES) segmented 3D surface models. The method first defines a partition of the ED endocardium using the geodesic distances from each surface point to apex, tricuspid valve and pulmonary valve: the landmarks that define the 3 regions. The ED surface mesh is then tetrahedralised, and the endocardial-defined partition is interpolated in the blood cavity via the Laplace equation. For obtaining an ES partition, the endocardial partition is transported from ED to ES using a commercial image-based tracking, and then interpolated towards the endocardium, similarly to ED, for computing volumes and ejection fraction (EF). We present a full assessment of the method's validity and reproducibility. First, we assess reproducibility under segmentation variability, obtaining intra- and inter- observer errors (4-10% and 10-23% resp.). Finally, we use a synthetic remodelling dataset to identify the situations in which our method is able to correctly determine the region that has remodelled. This dataset is generated by a novel mesh reconstruction method that deforms a reference mesh, locally imposing a given strain, expressed in anatomical coordinates. We show that the parcellation method is adequate for capturing local circumferential and global circumferential and longitudinal RV remodelling.


title: On the hypotheses of Penrose's singularity theorem under disformal transformations

abstract: We analyze how the hypotheses of Penrose's singularity theorem (1965) are modified by the action of disformal transformations (defined in terms of light-like vectors) upon a given space-time metric. In particular, we investigate the transformation of the null energy condition and the existence of closed trapped surfaces in such scenario, in order to derive conditions upon the background metric and the disformal vector that guarantee the validity of Penrose's theorem for disformal metrics. Then, we explain how to apply this technique for static and spherically symmetric space-times in general.


title: Taking snapshots of a quantum thermalization process: emergent classicality in quantum jump trajectories

abstract: We investigate theoretically the emergence of classical statistical physics in a finite quantum system that is either totally isolated or otherwise subjected to a quantum measurement process. We show via a random matrix theory approach to nonintegrable quantum systems that the set of outcomes of the measurement of a macroscopic observable evolve in time like stochastic variables, whose variance satisfies the celebrated Einstein relation for Brownian diffusion. Our results show how to extend the framework of eigenstate thermalization to the prediction of properties of quantum measurements on an otherwise closed quantum system. We show numerically the validity of the random matrix approach in quantum chain models.


title: Asymptotic normality of consecutive patterns in permutations encoded by generating trees with one-dimensional labels

abstract: We consider uniform random permutations drawn from a family enumerated through generating trees. We develop a new general technique to establish a central limit theorem for the number of consecutive occurrences of a fixed pattern in such permutations. We propose a technique to sample uniform permutations in such families as conditioned random colored walks. Building on that, we derive the behavior of the consecutive patterns in random permutations studying properties of the consecutive increments in the corresponding random walks. The method applies to families of permutations with a one-dimensional-labeled generating tree (together with some technical assumptions) and implies local convergence for random permutations in such families. We exhibit ten different families of permutations, most of them being permutation classes, that satisfy our assumptions. To the best of our knowledge, this is the first work where generating trees - which were introduced to enumerate combinatorial objects - have been used to establish probabilistic results.


title: Stress-Induced Dinoflagellate Bioluminescence at the Single Cell Level

abstract: One of the characteristic features of many marine dinoflagellates is their bioluminescence, which lights up nighttime breaking waves or seawater sliced by a ship's prow. While the internal biochemistry of light production by these microorganisms is well established, the manner by which fluid shear or mechanical forces trigger bioluminescence is still poorly understood. We report controlled measurements of the relation between mechanical stress and light production at the single-cell level, using high-speed imaging of micropipette-held cells of the marine dinoflagellate $Pyrocystis~lunula$ subjected to localized fluid flows or direct indentation. We find a viscoelastic response in which light intensity depends on both the amplitude and rate of deformation, consistent with the action of stretch-activated ion channels. A phenomenological model captures the experimental observations.

categories: cond-mat.soft physics.bio-ph physics.flu-dyn q-bio.CB


title: Hamiltonian theory of classical and quantum gauge invariant perturbations in Bianchi I spacetimes

abstract: We derive a Hamiltonian formulation of the theory of gauge invariant, linear perturbations in anisotropic Bianchi I spacetimes, and describe how to quantize this system. The matter content is assumed to be a minimally coupled scalar field with potential $V(\phi)$. We show that a Bianchi I spacetime generically induces both anisotropies and quantum entanglement on cosmological perturbations, and provide the tools to compute the details of these features. We then apply this formalism to a scenario in which the inflationary era is preceded by an anisotropic Bianchi I phase, and discuss the potential imprints in observable quantities. The formalism developed here paves the road to a simultaneous canonical quantization of both the homogeneous degrees of freedom and the perturbations, a task that we develop in a companion paper.


title: STEm-Seg: Spatio-temporal Embeddings for Instance Segmentation in Videos

abstract: Existing methods for instance segmentation in videos typically involve multi-stage pipelines that follow the tracking-by-detection paradigm and model a video clip as a sequence of images. Multiple networks are used to detect objects in individual frames, and then associate these detections over time. Hence, these methods are often non-end-to-end trainable and highly tailored to specific tasks. In this paper, we propose a different approach that is well-suited to a variety of tasks involving instance segmentation in videos. In particular, we model a video clip as a single 3D spatio-temporal volume, and propose a novel approach that segments and tracks instances across space and time in a single stage. Our problem formulation is centered around the idea of spatio-temporal embeddings which are trained to cluster pixels belonging to a specific object instance over an entire video clip. To this end, we introduce (i) novel mixing functions that enhance the feature representation of spatio-temporal embeddings, and (ii) a single-stage, proposal-free network that can reason about temporal context. Our network is trained end-to-end to learn spatio-temporal embeddings as well as parameters required to cluster these embeddings, thus simplifying inference. Our method achieves state-of-the-art results across multiple datasets and tasks. Code and models are available at https://github.com/sabarim/STEm-Seg.


title: One-meson-loop NJL model: Effect of collective and noncollective excitations on the quark condensate at finite temperature

abstract: We explore the effect of including quantum fluctuations in the two flavor Nambu$-$Jona-Lasinio model at finite temperature. This is accomplished, in a symmetry preserving way, by including collective and noncollective modes in the one-meson-loop gap equation which originate from poles and branch cuts in the complex plane, respectively. The inclusion of a boson cutoff, $\Lambda_b$, is necessary to regularize the meson-loop momenta. This new parameter is used to study the influence of going beyond the usual mean field approximation in the quark condensate. As the temperature increases, chiral symmetry tends to get restored, the collective modes melt and only noncollective modes contribute to the quark condensate. With the inclusion of such modes, the quark condensate at finite temperature has a different behavior than at mean field level that will be explored.


title: Super-Earths and sub-Neptunes are Insensitive to Stellar Metallicity

abstract: Kepler planets (including super-Earths and sub-Neptunes, from 1 to 4 Earth radii) are likely formed before the gaseous proto-planetary disks have dissipated, as are the Jovian planets. If the metal content in these disks resembles that in the host stars, one might expect Kepler planets to occur more frequently, and to be more massive, around metal-rich stars. Contrary to these expectations, we find that the radii of Kepler planets (a proxy for mass) are independent of host metallicity. Previous claims that larger planets prefer more metal-rich stars can be adequately explained by the combined facts that more massive stars tend to host bigger planets, and that more massive stars are also more metal rich in the Kepler sample. We interpret this independence as that the mass of a Kepler planet is not determined by the availability of solids, but is instead regulated by an as yet unknown process. Moreover, we find that the occurrence rates of Kepler planets rise only weakly with stellar metallicity, a trend that is further flattened when the influence of close stellar binaries is accounted for. We explain this weak dependence, in contrast to the strong dependence exhibited by Jovian planets, using a phenomenological model, wherein the masses of proto-planetary disks have a much larger spread than the spread in stellar metallicity, and wherein the formation of Jovian planets requires disks that contain some $5$ times more solid than that needed to form Kepler planets. This model predicts that stars more metal-poor than half-solar should rarely host any Kepler planets.


title: Optically rewritable memory in a graphene/ferroelectric-photovoltaic heterostructure

abstract: Achieving optical operation of logic elements, especially those that involve 2D layers, can open the long sought era of optical computing. However, the efficient optical modulation of the electronic properties of 2D materials including memory effect is currently missing. Here we report a fully optical control of the conductivity of a graphene with write/erase option yet under ultralow optical fluence. The competition between light-induced charge generation in ferroelectric-photovoltaic substrate with subsequent relaxation processes provides the selective photocarrier trapping control affecting the doping of 2D overlayer. These findings open the road to photonic control of 2D devices for all -optical modulators and a variety of all-optical logic circuits, memories and field-effect transistors.


title: Neural Fuzzy Extractors: A Secure Way to Use Artificial Neural Networks for Biometric User Authentication

abstract: Powered by new advances in sensor development and artificial intelligence, the decreasing cost of computation, and the pervasiveness of handheld computation devices, biometric user authentication (and identification) is rapidly becoming ubiquitous. Modern approaches to biometric authentication, based on sophisticated machine learning techniques, cannot avoid storing either trained-classifier details or explicit user biometric data, thus exposing users' credentials to falsification. In this paper, we introduce a secure way to handle user-specific information involved with the use of vector-space classifiers or artificial neural networks for biometric authentication. Our proposed architecture, called a Neural Fuzzy Extractor (NFE), allows the coupling of pre-existing classifiers with fuzzy extractors, through a artificial-neural-network-based buffer called an expander, with minimal or no performance degradation. The NFE thus offers all the performance advantages of modern deep-learning-based classifiers, and all the security of standard fuzzy extractors. We demonstrate the NFE retrofit to a classic artificial neural network for a simple scenario of fingerprint-based user authentication.


title: Convex plumbings in closed hyperbolic 4-manifolds

abstract: We show that every plumbing of disc bundles over surfaces whose genera satisfy a simple inequality may be embedded as a convex submanifold in some closed hyperbolic four-manifold. In particular its interior has a geometrically finite hyperbolic structure that covers a closed hyperbolic four-manifold.


title: Active flow control of a pump-induced wall-normal vortex with steady blowing

abstract: The emergence of a submerged vortex upstream of a pump can reduce pump intake efficiency and cause structural damage. In this study, we consider the use of active flow control with steady blowing to increase the pressure distribution within a single-phase pump-induced wall-normal vortex model, which is based on the Burgers vortex with a no-slip boundary condition prescribed along its symmetry plane. The goal of our control is to modify the vortex core velocity profile. These changes are sought to increase the core pressure such that detrimental effects on the pump are alleviated. Three-dimensional direct numerical simulations (DNS) are performed to examine the dynamics of the vortex with the application of axial momentum injection at and around the root of the vortex. We find that the active flow control approach can effectively modify the wall-normal vortical structure and significantly increase the low-core pressure by up to 81% compared to that of the uncontrolled case. The result shows that the control setup is also effective when it is introduced in an off-centered manner. Compared to the unsteady blowing and suction based actuation from our previous work (Liu et al. 2018), the current steady control technique offers an effective and simple flow control setup that can support robust operations of pumps.


title: Collaborative Distillation for Ultra-Resolution Universal Style Transfer

abstract: Universal style transfer methods typically leverage rich representations from deep Convolutional Neural Network (CNN) models (e.g., VGG-19) pre-trained on large collections of images. Despite the effectiveness, its application is heavily constrained by the large model size to handle ultra-resolution images given limited memory. In this work, we present a new knowledge distillation method (named Collaborative Distillation) for encoder-decoder based neural style transfer to reduce the convolutional filters. The main idea is underpinned by a finding that the encoder-decoder pairs construct an exclusive collaborative relationship, which is regarded as a new kind of knowledge for style transfer models. Moreover, to overcome the feature size mismatch when applying collaborative distillation, a linear embedding loss is introduced to drive the student network to learn a linear embedding of the teacher's features. Extensive experiments show the effectiveness of our method when applied to different universal style transfer approaches (WCT and AdaIN), even if the model size is reduced by 15.5 times. Especially, on WCT with the compressed models, we achieve ultra-resolution (over 40 megapixels) universal style transfer on a 12GB GPU for the first time. Further experiments on optimization-based stylization scheme show the generality of our algorithm on different stylization paradigms. Our code and trained models are available at https://github.com/mingsun-tse/collaborative-distillation.


title: A Corpus of Adpositional Supersenses for Mandarin Chinese

abstract: Adpositions are frequent markers of semantic relations, but they are highly ambiguous and vary significantly from language to language. Moreover, there is a dearth of annotated corpora for investigating the cross-linguistic variation of adposition semantics, or for building multilingual disambiguation systems. This paper presents a corpus in which all adpositions have been semantically annotated in Mandarin Chinese; to the best of our knowledge, this is the first Chinese corpus to be broadly annotated with adposition semantics. Our approach adapts a framework that defined a general set of supersenses according to ostensibly language-independent semantic criteria, though its development focused primarily on English prepositions (Schneider et al., 2018). We find that the supersense categories are well-suited to Chinese adpositions despite syntactic differences from English. On a Mandarin translation of The Little Prince, we achieve high inter-annotator agreement and analyze semantic correspondences of adposition tokens in bitext.


title: Non-Gaussian Entanglement Renormalization for Quantum Fields

abstract: In this work, a non-Gaussian cMERA tensor network for interacting quantum field theories (icMERA) is presented. This consists of a continuous tensor network circuit in which the generator of the entanglement renormalization of the wavefunction is nonperturbatively extended with nonquadratic variational terms. The icMERA circuit nonperturbatively implements a set of scale dependent nonlinear transformations on the fields of the theory, which suppose a generalization of the scale dependent linear transformations induced by the Gaussian cMERA circuit. Here we present these transformations for the case of self-interacting scalar and fermionic field theories. Finally, the icMERA tensor network is fully optimized for the $\lambda \phi^4$ theory in $(1+1)$ dimensions. This allows us to evaluate, nonperturbatively, the connected parts of the two- and four-point correlation functions. Our results show that icMERA wavefunctionals encode proper non-Gaussian correlations of the theory, thus providing a new variational tool to study phenomena related with strongly interacting field theories.


title: Kibble-Zurek behavior in disordered Chern insulators

abstract: Even though no local order parameter in the sense of the Landau theory exists for topological quantum phase transitions in Chern insulators, the highly non-local Berry curvature exhibits critical behavior near a quantum critical point. We investigate the critical properties of its real space analog, the local Chern marker, in weakly disordered Chern insulators. Due to disorder, inhomogeneities appear in the spatial distribution of the local Chern marker. Their size exhibits power-law scaling with the critical exponent matching the one extracted from the Berry curvature of a clean system. We drive the system slowly through such a quantum phase transition. The characteristic size of inhomogeneities in the non-equilibrium post-quench state obeys the Kibble-Zurek scaling. In this setting, the local Chern marker thus does behave in a similar way as a local order parameter for a symmetry breaking second order phase transition. The Kibble-Zurek scaling also holds for the inhomogeneities in the spatial distribution of excitations and of the orbital polarization.


title: Synthesize then Compare: Detecting Failures and Anomalies for Semantic Segmentation

abstract: The ability to detect failures and anomalies are fundamental requirements for building reliable systems for computer vision applications, especially safety-critical applications of semantic segmentation, such as autonomous driving and medical image analysis. In this paper, we systematically study failure and anomaly detection for semantic segmentation and propose a unified framework, consisting of two modules, to address these two related problems. The first module is an image synthesis module, which generates a synthesized image from a segmentation layout map, and the second is a comparison module, which computes the difference between the synthesized image and the input image. We validate our framework on three challenging datasets and improve the state-of-the-arts by large margins, \emph{i.e.}, 6% AUPR-Error on Cityscapes, 7% Pearson correlation on pancreatic tumor segmentation in MSD and 20% AUPR on StreetHazards anomaly segmentation.


title: Detecting Pancreatic Ductal Adenocarcinoma in Multi-phase CT Scans via Alignment Ensemble

abstract: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers among the population. Screening for PDACs in dynamic contrast-enhanced CT is beneficial for early diagnosis. In this paper, we investigate the problem of automated detecting PDACs in multi-phase (arterial and venous) CT scans. Multiple phases provide more information than single phase, but they are unaligned and inhomogeneous in texture, making it difficult to combine cross-phase information seamlessly. We study multiple phase alignment strategies, i.e., early alignment (image registration), late alignment (high-level feature registration), and slow alignment (multi-level feature registration), and suggest an ensemble of all these alignments as a promising way to boost the performance of PDAC detection. We provide an extensive empirical evaluation on two PDAC datasets and show that the proposed alignment ensemble significantly outperforms previous state-of-the-art approaches, illustrating the strong potential for clinical use.


title: Alternating odd pretzel knots and chirally cosmetic surgeries

abstract: A pair of surgeries on a knot is chirally cosmetic if they result in homeomorphic manifolds with opposite orientations. Using recent methods of Ichihara, Ito, and Saito, we show that, except for the (2,5) and (2,7)-torus knots, the genus 2 and 3 alternating odd pretzel knots do not admit any chirally cosmetic surgeries. Further, we show that for a fixed genus, at most finitely many alternating odd pretzel knots admit chirally cosmetic surgeries.


title: The left-right supersymmetric option at a high-energy upgrade of the LHC

abstract: We investigate the possibility that a minimal realization of left-right supersymmetry can be reachable at a high-energy upgrade of the LHC, expected to operate at a center-of-mass energy of 27 TeV. This minimal scenario has a relatively light $SU(2)_R$ doubly-charged Higgs boson, which could decay dominantly into tau-lepton pairs. We explore the associated signals comprised of at least three hadronically-decaying taus, or with at least two hadronic taus and one same-sign-same-flavor charged lepton pair. Our analysis shows that the former signature is challenging to use for getting handles on the signal due to the large corresponding background, and that the latter one can lead to a handful of new physics events in an almost background-free environment. We however find that a signal comprised of three hadronically-decaying tau leptons is likely to be observed at a low luminosity of proton-proton collisions at a 27 TeV upgrade of the LHC.


title: NELA-GT-2019: A Large Multi-Labelled News Dataset for The Study of Misinformation in News Articles

abstract: In this paper, we present an updated version of the NELA-GT-2018 dataset (N{\o}rregaard, Horne, and Adal{\i} 2019), entitled NELA-GT-2019. NELA-GT-2019 contains 1.12M news articles from 260 sources collected between January 1st 2019 and December 31st 2019. Just as with NELA-GT-2018, these sources come from a wide range of mainstream news sources and alternative news sources. Included with the dataset are source-level ground truth labels from 7 different assessment sites covering multiple dimensions of veracity. The NELA-GT-2019 dataset can be found at: https://doi.org/10.7910/DVN/O7FWPO


title: Placement Optimization with Deep Reinforcement Learning

abstract: Placement Optimization is an important problem in systems and chip design, which consists of mapping the nodes of a graph onto a limited set of resources to optimize for an objective, subject to constraints. In this paper, we start by motivating reinforcement learning as a solution to the placement problem. We then give an overview of what deep reinforcement learning is. We next formulate the placement problem as a reinforcement learning problem and show how this problem can be solved with policy gradient optimization. Finally, we describe lessons we have learned from training deep reinforcement learning policies across a variety of placement optimization problems.


title: Exploring the relation between students' online learning behavior and course performance by including contextual information in data analysis

abstract: This study examines whether including more contextual information in data analysis could improve our ability to identify the relation between students' online learning behavior and overall performance in an introductory physics course. We created four linear regression models correlating students' pass-fail events in a sequence of online learning modules with their normalized total course score. Each model takes into account an additional level of contextual information than the previous one, such as student learning strategy and duration of assessment attempts. Each of the latter three models is also accompanied by a visual representation of students' interaction states on each learning module. We found that the best performing model is the one that includes the most contextual information, including instruction condition, internal condition, and learning strategy. The model shows that while most students failed on the most challenging learning module, those with normal learning behavior are more likely to obtain higher total course scores, whereas students who resorted to guessing on the assessments of subsequent modules tended to receive lower total scores. Our results suggest that considering more contextual information related to each event can be an effective method to improve the quality of learning analytics, leading to more accurate and actionable recommendations for instructors.


title: Potential Neutralizing Antibodies Discovered for Novel Corona Virus Using Machine Learning

abstract: The fast and untraceable virus mutations take lives of thousands of people before the immune system can produce the inhibitory antibody. Recent outbreak of novel coronavirus infected and killed thousands of people in the world. Rapid methods in finding peptides or antibody sequences that can inhibit the viral epitopes of COVID-19 will save the life of thousands. In this paper, we devised a machine learning (ML) model to predict the possible inhibitory synthetic antibodies for Corona virus. We collected 1933 virus-antibody sequences and their clinical patient neutralization response and trained an ML model to predict the antibody response. Using graph featurization with variety of ML methods, we screened thousands of hypothetical antibody sequences and found 8 stable antibodies that potentially inhibit COVID-19. We combined bioinformatics, structural biology, and Molecular Dynamics (MD) simulations to verify the stability of the candidate antibodies that can inhibit the Corona virus.


title: Model independent tests of the Kerr bound with extreme mass ratio inspirals

abstract: An outstanding prediction of general relativity is the fact that the angular momentum $S$ of an isolated black hole with mass $\mu$ is limited by the Kerr bound, $S\leq G\mu^2/c$. Testing this cornerstone is challenging due to the difficulty in modelling spinning compact objects that violate this bound. We argue that precise, model-independent tests can be achieved by measuring gravitational waves from an extreme mass ratio inspiral around a supermassive object, one of the main targets of the future LISA mission. In the extreme mass ratio limit, the dynamics of the small compact object depends only on its multipole moments, which are free parameters. At variance with the comparable-mass case, accurate waveforms are valid also when the spin of the small object greatly exceeds the Kerr bound. By computing the orbital dephasing and the gravitational-wave signal emitted by a spinning point particle in circular, nonprecessing, equatorial motion around a Kerr black hole, we estimate that LISA will be able to measure the spin of the small compact object at the level of $10\%$. Together with mass measurements, this will allow for theory-agnostic, unprecedented constraints on string-theory inspired objects such as "superspinars", almost in their entire parameter space.


title: Causal Simulation Experiments: Lessons from Bias Amplification

abstract: Recent theoretical work in causal inference has explored an important class of variables which, when conditioned on, may further amplify existing unmeasured confounding bias (bias amplification). Despite this theoretical work, existing simulations of bias amplification in clinical settings have suggested bias amplification may not be as important in many practical cases as suggested in the theoretical literature.We resolve this tension by using tools from the semi-parametric regression literature leading to a general characterization in terms of the geometry of OLS estimators which allows us to extend current results to a larger class of DAGs, functional forms, and distributional assumptions. We further use these results to understand the limitations of current simulation approaches and to propose a new framework for performing causal simulation experiments to compare estimators. We then evaluate the challenges and benefits of extending this simulation approach to the context of a real clinical data set with a binary treatment, laying the groundwork for a principled approach to sensitivity analysis for bias amplification in the presence of unmeasured confounding.


title: A Variational Analysis Approach to Solving the Merton Problem

abstract: We address the Merton problem of maximizing the expected utility of terminal wealth using techniques from variational analysis. Under a general continuous semimartingale market model with stochastic parameters, we obtain a characterization of the optimal portfolio for general utility functions in terms of a forward-backward stochastic differential equation (FBSDE) and derive solutions for a number of well-known utility functions. Our results complement a previous studies conducted on optimal strategies in markets driven by Brownian noise with random drift and volatility parameters.


title: Bolometric Bond Albedo and Thermal Inertia Maps of Mimas

abstract: We use data from Cassini's Composite Infrared Spectrometer CIRS to map Mimas' surface temperatures and its thermophysical properties. This provides a dramatic improvement on the work in Howett et al. (2011), where the values were determined at only two regions on Mimas (one inside, and another outside of the anomalous region). We use all spatially-resolved scans made by CIRS' focal plane 3 (FP3, 600 to 1100 cm-1) of Mimas' surface, which are largely daytime observations but do include one nighttime one. The resulting temperature maps confirm the presence and location of Mimas' previously discovered thermally anomalous region. No other thermally anomalous regions were discovered, although we note that the surface coverage is incomplete on Mimas' leading and anti-Saturn hemisphere. The thermal inertia map confirms that the anomalous region has a notably higher thermal inertia than its surroundings: 98 +/- 42 J m-2 K-1 s-1/2 inside of the anomaly, compared to 34 +/- 32 J m-2 K-1 s-1/2 outside. The albedo inside and outside of the anomalous region agrees within their uncertainty: 0.45+/-0.08 inside compared to 0.41 +/- 0.07 outside the anomaly. Interestingly the albedo appears brighter inside the anomaly region, which may not be surprising given this region does appear brighter at some UV wavelengths (0.338 microns, see Schenk et al., 2011). However, this result should be treated with caution because, as previously stated, statistically the albedo of these two regions is the same when their uncertainties are considered. These thermal inertia and albedo values determined here are consistent with those found by Howett et al. (2011), who determined the thermal inertia inside the anomaly to be 66 +/- 23 J m-2 K-1 s-1/2 and less than 16 J m-2 K-1 s-1/2 outside, with albedos that varied from 0.49 to 0.70.


title: Extensions and deformations of algebras with higher derivations

abstract: Higher derivations on an associative algebra generalizes higher order derivatives. We call a tuple consisting of an algebra and a higher derivation on it by an AssHDer pair. We define a cohomology for AssHDer pairs with coefficients in a representation. Next, we study central extensions, abelian extensions of AssHDer pairs and relate them to the second cohomology group with appropriate coefficients. Deformations of AssHDer pairs are also considered which are governed by the cohomology with self coefficient.


title: Measurement of the shape of the $B_s^0\rightarrow D_s^{*-}\mu^{+}\nu_{\mu}$ differential decay rate

abstract: The shape of the $B_s^0\rightarrow D_s^{*-}\mu^+\nu_{\mu}$ differential decay rate is obtained as a function of the hadron recoil using proton-proton collision data at a centre-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 1.7 fb$^{-1}$ collected by the LHCb detector. The $B_s^0\rightarrow D_s^{*-}\mu^+\nu_{\mu}$ decay is reconstructed through the decays $D_s^{*-}\rightarrow D_s^{-}\gamma$ and $D_s^{-}\rightarrow K^-K^+\pi^-$. The differential decay rate is fitted with the Caprini-Lellouch-Neubert (CLN) and Boyd-Grinstein-Lebed (BGL) parametrisations of the form factors, and the relevant quantities for both are extracted.


title: Local and global densities for Weierstrass models of elliptic curves

abstract: We prove local results on the $p$-adic density of elliptic curves over $\mathbb{Q}_p$ with different reduction types, together with global results on densities of elliptic curves over $\mathbb{Q}$ with specified reduction types at one or more (including infinitely many) primes. These global results include: the density of integral Weierstrass equations which are minimal models of semistable elliptic curves over $\mathbb{Q}$ (that is, elliptic curves with square-free conductor) is $1/\zeta(2)\approx60.79\%$, the same as the density of square-free integers; the density of semistable elliptic curves over $\mathbb{Q}$ is $\zeta(10)/\zeta(2)\approx60.85\%$; the density of integral Weierstrass equations which have square-free discriminant is $\prod_p\left(1-\frac{2}{p^2}+\frac{1}{p^3}\right) \approx 42.89\%$, which is the same (except for a different factor at the prime $2$) as the density of monic integral cubic polynomials with square-free discriminant (and agrees with a previous result of Baier and Browning for short Weierstrass equations); and the density of elliptic curves over $\mathbb{Q}$ with square-free minimal discriminant is $\zeta(10)\prod_p\left(1-\frac{2}{p^2}+\frac{1}{p^3}\right)\approx42.93\%$. The local results derive from a detailed analysis of Tate's Algorithm, while the global ones are obtained through the use of the Ekedahl Sieve, as developed by Poonen, Stoll, and Bhargava.


title: Topological pressure for discontinuous semiflows and a variational principle for impulsive dynamical systems

abstract: We introduce four, a priori different, notions of topological pressure for possibly discontinuous semiflows acting on compact metric spaces and observe that they all agree with the classical one when restricted to the continuous setting. Moreover, for a class of \emph{impulsive semiflows}, which are examples of discontinuous systems, we prove a variational principle. As a consequence, we conclude that for this class of systems the four notions coincide and, moreover, they also coincide with the notion of topological pressure introduced in \cite{ACS17}.


title: Convex Hulls of Random Order Types

abstract: We establish the following two main results on order types of points in general position in the plane (realizable simple planar order types, realizable uniform acyclic oriented matroids of rank $3$): (a) The number of extreme points in an $n$-point order type, chosen uniformly at random from all such order types, is on average $4+o(1)$. For labeled order types, this number has average $4- \frac{8}{n^2 - n +2}$ and variance at most $3$. (b) The (labeled) order types read off a set of $n$ points sampled independently from the uniform measure on a convex planar domain, smooth or polygonal, or from a Gaussian distribution are concentrated, i.e. such sampling typically encounters only a vanishingly small fraction of all order types of the given size. Result (a) generalizes to arbitrary dimension $d$ for labeled order types with the average number of extreme points $2d+o(1)$ and constant variance. We also discuss to what extent our methods generalize to the abstract setting of uniform acyclic oriented matroids. Moreover, our methods allow to show the following relative of the Erd\H{o}s-Szekeres theorem: for any fixed $k$, as $n \to \infty$, a proportion $1 - O(1/n)$ of the $n$-point simple order types contain a triangle enclosing a convex $k$-chain over an edge. For the unlabeled case in (a), we prove that for any antipodal, finite subset of the $2$-dimensional sphere, the group of orientation preserving bijections is cyclic, dihedral or one of $A_4$, $S_4$ or $A_5$ (and each case is possible). These are the finite subgroups of $SO(3)$ and our proof follows the lines of their characterization by Felix Klein.


title: Malicious Experts versus the multiplicative weights algorithm in online prediction

abstract: We consider a prediction problem with two experts and a forecaster. We assume that one of the experts is honest and makes correct prediction with probability $\mu$ at each round. The other one is malicious, who knows true outcomes at each round and makes predictions in order to maximize the loss of the forecaster. Assuming the forecaster adopts the classical multiplicative weights algorithm, we find upper and lower bounds for the value function of the malicious expert. Our results imply that the multiplicative weights algorithm cannot resist the corruption of malicious experts. We also show that an adaptive multiplicative weights algorithm is asymptotically optimal for the forecaster, and hence more resistant to the corruption of malicious experts.


title: My SIM is Leaking My Data: Exposing Self-Login Privacy Breaches in Smartphones

abstract: We expose a new security leak for smartphone users, which allows to stole user personal data by accessing the mobile operator user page when auto-login is employed. We show how any "apparently" genuine app can steal these data from some mobile operators, affecting more than 80% of Italian mobile smartphones.


title: Precise control of magnetic fields and optical polarization in a time-orbiting potential trap

abstract: A time orbiting potential trap confines neutral atoms in a rotating magnetic field. The rotation of the field can be useful for precision measurements, since it can average out some systematic effects. However, the field is more difficult to characterize than a static field, and it makes light applied to the atoms have a time-varying optical polarization relative to the quantization axis. These problems can be overcome using stroboscopic techniques, where either a radio-frequency field or a laser is applied in pulses that are synchronized to the rotating field. Using these methods, the magnetic field can be characterized with a precision of 10 mG and light can be applied with a polarization error of $5\times 10^{-5}$.


title: Testing a thermodynamic approach to collective animal behavior in laboratory fish schools

abstract: Collective behaviors displayed by groups of social animals are observed frequently in nature. Understanding and predicting the behavior of complex biological systems is dependent on developing effective descriptions and models. While collective animal systems are characteristically non-equilibrium, we can employ concepts from equilibrium statistical mechanics to motivate the measurement of material-like properties in laboratory animal aggregates. Here, we present results from a new set of experiments that utilize high speed footage of two-dimensional schooling events, particle tracking, and projected static and dynamic light fields to observe and control the behavior of negatively phototaxic fish schools ($\textit{Hemigrammus}$ $\textit{bleheri}$). First, we use static light fields consisting of dark circular regions to produce visual stimuli that confine the schools to a range of areas. We find that the school has a maximum density independent of group size, and that the swim pressure increases linearly with number density, suggesting that unperturbed schools exist on an isotherm. Next, we use dynamic light fields where the radius of the dark region shrinks linearly with time to compress the schools. We find that the effective temperature parameter depends on the compression time and our results are thus consistent with the school having a constant heat flux. These findings provide further evidence for the utility of effective thermodynamic descriptions descriptions of non-equilibrium systems in collective animal behavior.


title: Stochastic Virtual Battery Modeling of Uncertain Electrical Loads Using Variational Autoencoder

abstract: Effective utilization of flexible loads for grid services, while satisfying end-user preferences and constraints, requires an accurate estimation of the aggregated predictive flexibility offered by the electrical loads. Virtual battery (VB) models are often used to quantify the predictive flexibility in thermostatic loads (e.g. residential air-conditioners, electric water-heaters), which model the temporal evolution of a (virtual) energy state via a first order dynamics including self-dissipation rate, and power and energy capacities as parameters. Uncertainties and lack of information regarding end-usage and equipment models render deterministic VB models impractical. In this paper, we introduce the notion of stochastic VB models, and propose a \textit{variational autoencoder}-based deep learning algorithm to identify the probability distribution of the VB model parameters. Using available sensors and meters data, the proposed algorithm generates not only point estimates of the VB parameters, but also confidence intervals around those values. Effectiveness of the proposed frameworks is demonstrated on a collection of electric water-heater loads, whose operation is driven by uncertain water usage profiles.


title: Semi-supervised few-shot learning for medical image segmentation

abstract: Recent years have witnessed the great progress of deep neural networks on semantic segmentation, particularly in medical imaging. Nevertheless, training high-performing models require large amounts of pixel-level ground truth masks, which can be prohibitive to obtain in the medical domain. Furthermore, training such models in a low-data regime highly increases the risk of overfitting. Recent attempts to alleviate the need for large annotated datasets have developed training strategies under the few-shot learning paradigm, which addresses this shortcoming by learning a novel class from only a few labeled examples. In this context, a segmentation model is trained on episodes, which represent different segmentation problems, each of them trained with a very small labeled dataset. In this work, we propose a novel few-shot learning framework for semantic segmentation, where unlabeled images are also made available at each episode. To handle this new learning paradigm, we propose to include surrogate tasks that can leverage very powerful supervisory signals --derived from the data itself-- for semantic feature learning. We show that including unlabeled surrogate tasks in the episodic training leads to more powerful feature representations, which ultimately results in better generability to unseen tasks. We demonstrate the efficiency of our method in the task of skin lesion segmentation in two publicly available datasets. Furthermore, our approach is general and model-agnostic, which can be combined with different deep architectures.


title: Probing Mid-Infrared Phonon Polaritons in the Aqueous Phase

abstract: Phonon polaritons (PhPs), the collective phonon oscillations with hybridized electromagnetic fields, concentrate optical fields in the mid-infrared frequency range that matches the vibrational modes of molecules. The utilization of PhPs holds the promise for chemical sensing tools and polariton-enhanced nanospectroscopy. However, investigations and innovations on PhPs in the aqueous phase remains stagnant, because of the lack of in situ mid-infrared nano-imaging methods in water. Strong infrared absorption from water prohibits optical delivery and detection in the mid-infrared for scattering-type near-field microscopy. Here, we present our solution: the detection of photothermal responses caused by the excitation of PhPs by liquid phase peak force infrared (LiPFIR) microscopy. Characteristic interference fringes of PhPs in 10B isotope-enriched h-BN were measured in the aqueous phase and their dispersion relationship extracted. LiPFIR enables the measurement of mid-infrared PhPs in the fluid phase, opening possibilities, and facilitating the development of mid-IR phonon polaritonics in water.


title: LARgE Survey -- III. Environments of Ultra-Massive Passive Galaxies at Cosmic Noon: BCG progenitors growing through mergers

abstract: We study the environments of a sample of 61 extremely rare z~1.6 Ultra-Massive Passively Evolving Galaxies (UMPEGs: stellar masses M_stars >10^11.5 M_sun) which -- based on clustering analysis presented in Cheema et al. (2020) -- appear to be associated with very massive (M_halo ~ 10^14.1 h^-1 M_sun) dark matter halos that are likely to be the progenitors of z~0 massive (Coma- and Virgo-like) galaxy clusters. We find that UMPEGs on average have fewer than one satellite galaxy with mass ratio M_sat : M_UMPEG >~ 1:5 (i.e., M_sat >~ 10^10.8 M_sun) within 0.5 Mpc; the large mass gap that we observe between the typical UMPEG and its most massive satellite implies that the z~1.6 UMPEGs assembled through major mergers. Using observed satellite counts with merger timescales from the literature, we estimate the growth rate due to mergers with mass ratio of >~ 1:4 to be ~13% Gyr^-1 (with a ~2x systematic uncertainty). This relatively low growth rate is unlikely to significantly affect the shape of the massive end of the stellar mass function, whose evolution must instead be driven by the quenching of new cohorts of ultra-massive star-forming galaxies. However, this growth rate is high enough that, if sustained to z~0, the typical z~1.6 M_UMPEG=10^11.6 M_sun UMPEG can grow into a M_stars~10^12 M_sun brightest cluster galaxy (BCG) of a present-day massive galaxy cluster. Our observations favour a scenario in which our UMPEGs are main-branch progenitors of some of the present-day BCGs that have first assembled through major mergers at high redshifts and grown further through (likely minor) merging at later times.


title: Topological control of the nonlinear-optical response of hybrid quantum systems

abstract: We map the topological properties of a one dimensional superlattice to the optical properties of an electronic system. We find that the nonlinear-optical response is optimized for electrons that live in the transitional morphology between topologically protected edge states and delocalized eigenstates. This provides a novel means of tuning the nonlinear-optical response of hybrid quantum systems. We show how these characteristics can be used to mimic saturable absorption and illustrate how `quantum cords' can be used to build an efficient all-optical switch.


title: Fully reversible neural networks for large-scale 3D seismic horizon tracking

abstract: Tracking a horizon in seismic images or 3D volumes is an integral part of seismic interpretation. The last few decades saw progress in using neural networks for this task, starting from shallow networks for 1D traces, to deeper convolutional neural networks for large 2D images. Because geological structures are intrinsically 3D, we hope to see improved horizon tracking by training networks on 3D seismic data cubes. While there are some 3D convolutional neural networks for various seismic interpretation tasks, they are restricted to shallow networks or relatively small 3D inputs because of memory limitations. The required memory for the network states and weights increases with network depth. We present a fully reversible network for horizon tracking that has a memory requirement that is independent of network depth. To tackle memory issues regarding the network weights, we use layers that train in a factorized form directly. Therefore, we can maintain a large number of network channels while keeping the number of convolutional kernels low. We use the saved memory to increase the input size of the data by order of magnitude such that the network can better learn from large structures in the data. A field data example verifies the proposed network structure is suitable for seismic horizon tracking.


title: Modeling electronic, mechanical, optical and thermal properties of graphene-like BC$_6$N materials: Role of prominent BN-bonds

abstract: We model monolayer graphene-like materials with BC$_6$N stoichiometry where the bonding between the B and the N atoms plays an important role for their physical and chemical properties. Two types of BC$_6$N are found based on the BN bonds: In the presence of BN bonds, an even number of $\pi$-bonds emerges indicating an aromatic structure and a large direct bandgap appears, while in the absence of BN bonds, an anti-aromatic structure with an odd-number of $\pi$-bonds is found resulting a direct small bandgap. The stress-strain curves shows high elastic moduli and tensile strength of the structures with BN-bonds, compared to structures without BN-bonds. Self-consistent field calculations demonstrate that BC$_6$N with BN-bonds is energetically more stable than structures without BN-bonds due to a strong binding energy between the B and the N atoms, while their phonon dispersion displays that BC$_6$N without BN-bonds has more dynamical stability. Furthermore, all the BC$_6$N structures considered show a large absorption of electromagnetic radiation with polarization parallel to the monolayers in the visible range. Finer detail of the absorption depend on the actual structures of the layers. A higher electronic thermal conductivity and specific heat are seen in BC$_6$N systems caused by hot carrier--assisted charge transport. This opens up a possible optimization for bolometric applications of graphene based material devices.


title: Fan-Crossing Free Graphs and Their Relationship to other Beyond-Planar Graphs

abstract: A graph is \emph{fan-crossing free} if it has a drawing in the plane so that each edge is crossed by independent edges, that is the crossing edges have distinct vertices. On the other hand, it is \emph{fan-crossing} if the crossing edges have a common vertex, that is they form a fan. Both are prominent examples for beyond-planar graphs. Further well-known beyond-planar classes are the $k$-planar, $k$-gap-planar, quasi-planar, and right angle crossing graphs. We use the subdivision, node-to-circle expansion and path-addition operations to distinguish all these graph classes. In particular, we show that the 2-subdivision and the node-to-circle expansion of any graph is fan-crossing free, which does not hold for fan-crossing and $k$-(gap)-planar graphs, respectively. Thereby, we obtain graphs that are fan-crossing free and neither fan-crossing nor $k$-(gap)-planar. Finally, we show that some graphs have a unique fan-crossing free embedding, that there are thinned maximal fan-crossing free graphs, and that the recognition problem for fan-crossing free graphs is NP-complete.


title: Train, Learn, Expand, Repeat

abstract: High-quality labeled data is essential to successfully train supervised machine learning models. Although a large amount of unlabeled data is present in the medical domain, labeling poses a major challenge: medical professionals who can expertly label the data are a scarce and expensive resource. Making matters worse, voxel-wise delineation of data (e.g. for segmentation tasks) is tedious and suffers from high inter-rater variance, thus dramatically limiting available training data. We propose a recursive training strategy to perform the task of semantic segmentation given only very few training samples with pixel-level annotations. We expand on this small training set having cheaper image-level annotations using a recursive training strategy. We apply this technique on the segmentation of intracranial hemorrhage (ICH) in CT (computed tomography) scans of the brain, where typically few annotated data is available.


title: Imaging domain reversal in an ultrathin van der Waals ferromagnet

abstract: The recent isolation of two-dimensional van der Waals magnetic materials has uncovered rich physics that often differs from the magnetic behaviour of their bulk counterparts. However, the microscopic details of fundamental processes such as the initial magnetization or domain reversal, which govern the magnetic hysteresis, remain largely unknown in the ultrathin limit. Here we employ a widefield nitrogen-vacancy (NV) microscope to directly image these processes in few-layer flakes of magnetic semiconductor vanadium triiodide (VI$_3$). We observe complete and abrupt switching of most flakes at fields $H_c\approx0.5-1$ T (at 5 K) independent of thickness down to two atomic layers, with no intermediate partially-reversed state. The coercive field decreases as the temperature approaches the Curie temperature ($T_c\approx50$ K), however, the switching remains abrupt. We then image the initial magnetization process, which reveals thickness-dependent domain wall depinning fields well below $H_c$. These results point to ultrathin VI$_3$ being a nucleation-type hard ferromagnet, where the coercive field is set by the anisotropy-limited domain wall nucleation field. This work illustrates the power of widefield NV microscopy to investigate magnetization processes in van der Waals ferromagnets, which could be used to elucidate the origin of the hard ferromagnetic properties of other materials and explore field- and current-driven domain wall dynamics.


title: Scientific AI in materials science: a path to a sustainable and scalable paradigm

abstract: Recently there has been an ever-increasing trend in the use of machine learning (ML) and artificial intelligence (AI) methods by the materials science, condensed matter physics, and chemistry communities. This perspective article identifies key scientific, technical, and social opportunities that the materials community must prioritize to consistently develop and leverage Scientific AI to provide a credible path towards the advancement of current materials-limited technologies. Here we highlight the intersections of these opportunities with a series of proposed paths forward. The opportunities are roughly sorted from scientific/technical (e.g., development of robust, physically meaningful multiscale material representations) to social (e.g., promoting an AI-ready workforce). The proposed paths forward range from developing new infrastructure and capabilities to deploying them in industry and academia. We provide a brief introduction to AI in materials science and engineering, followed by detailed discussions of each of the opportunities and paths forward.


title: MINT: Deep Network Compression via Mutual Information-based Neuron Trimming

abstract: Most approaches to deep neural network compression via pruning either evaluate a filter's importance using its weights or optimize an alternative objective function with sparsity constraints. While these methods offer a useful way to approximate contributions from similar filters, they often either ignore the dependency between layers or solve a more difficult optimization objective than standard cross-entropy. Our method, Mutual Information-based Neuron Trimming (MINT), approaches deep compression via pruning by enforcing sparsity based on the strength of the relationship between filters of adjacent layers, across every pair of layers. The relationship is calculated using conditional geometric mutual information which evaluates the amount of similar information exchanged between the filters using a graph-based criterion. When pruning a network, we ensure that retained filters contribute the majority of the information towards succeeding layers which ensures high performance. Our novel approach outperforms existing state-of-the-art compression-via-pruning methods on the standard benchmarks for this task: MNIST, CIFAR-10, and ILSVRC2012, across a variety of network architectures. In addition, we discuss our observations of a common denominator between our pruning methodology's response to adversarial attacks and calibration statistics when compared to the original network.


title: Viewport-Aware Deep Reinforcement Learning Approach for 360$^o$ Video Caching

abstract: 360$^o$ video is an essential component of VR/AR/MR systems that provides immersive experience to the users. However, 360$^o$ video is associated with high bandwidth requirements. The required bandwidth can be reduced by exploiting the fact that users are interested in viewing only a part of the video scene and that users request viewports that overlap with each other. Motivated by the findings of recent works where the benefits of caching video tiles at edge servers instead of caching entire 360$^o$ videos were shown, in this paper, we introduce the concept of virtual viewports that have the same number of tiles with the original viewports. The tiles forming these viewports are the most popular ones for each video and are determined by the users' requests. Then, we propose a proactive caching scheme that assumes unknown videos' and viewports' popularity. Our scheme determines which videos to cache as well as which is the optimal virtual viewport per video. Virtual viewports permit to lower the dimensionality of the cache optimization problem. To solve the problem, we first formulate the content placement of 360$^o$ videos in edge cache networks as a Markov Decision Process (MDP), and then we determine the optimal caching placement using the Deep Q-Network (DQN) algorithm. The proposed solution aims at maximizing the overall quality of the 360$^o$ videos delivered to the end-users by caching the most popular 360$^o$ videos at base quality along with a virtual viewport in high quality. We extensively evaluate the performance of the proposed system and compare it with that of known systems such as LFU, LRU, FIFO, over both synthetic and real 360$^o$ video traces. The results reveal the large benefits coming from proactive caching of virtual viewports instead of the original ones in terms of the overall quality of the rendered viewports, the cache hit ratio, and the servicing cost.


title: TILES-2018, a longitudinal physiologic and behavioral data set of hospital workers

abstract: We present a novel longitudinal multimodal corpus of physiological and behavioral data collected from direct clinical providers in a hospital workplace. We designed the study to investigate the use of off-the-shelf wearable and environmental sensors to understand individual-specific constructs such as job performance, interpersonal interaction, and well-being of hospital workers over time in their natural day-to-day job settings. We collected behavioral and physiological data from $n = 212$ participants through Internet-of-Things Bluetooth data hubs, wearable sensors (including a wristband, a biometrics-tracking garment, a smartphone, and an audio-feature recorder), together with a battery of surveys to assess personality traits, behavioral states, job performance, and well-being over time. Besides the default use of the data set, we envision several novel research opportunities and potential applications, including multi-modal and multi-task behavioral modeling, authentication through biometrics, and privacy-aware and privacy-preserving machine learning.


title: CLASH-VLT: a full dynamical reconstruction of the mass profile of Abell S1063 from 1 kpc out to the virial radius

abstract: Context: The shape of the mass density profiles of cosmological halos informs us of the nature of DM and DM-baryons interactions. Previous estimates of the inner slope of the mass density profiles of clusters of galaxies are in opposition to predictions derived from numerical simulations of CDM. Aims: We determine the inner slope of the DM density profile of a massive cluster of galaxies, Abell S1063 (RXC J2248.7$-$4431) at $z=0.35$, with a dynamical analysis based on an extensive spectroscopic campaign carried out with the VIMOS and MUSE spectrographs at the ESO VLT. This new data set provides an unprecedented sample of 1234 spectroscopic members, 104 of which are located in the cluster core (R < 200 kpc), extracted from the MUSE integral field spectroscopy. The latter also allows the stellar velocity dispersion profile of the brightest cluster galaxy (BCG) to be measured out to 40 kpc. Methods: We used an upgraded version of the MAMPOSSt technique to perform a joint maximum likelihood fit to the velocity dispersion profile of the BCG and to the velocity distribution of cluster member galaxies over a radial range from 1 kpc to the virial radius (~ 2.7 Mpc). Results: We find a value of $\gamma_{DM} =0.99 \pm 0.04$ for the inner logarithmic slope of the DM density profile after marginalizing over all the other parameters. The newly determined dynamical mass profile is found to be in excellent agreement with the mass density profiles obtained from the independent X-ray hydrostatic analysis based on deep Chandra data, as well as the strong and weak lensing analyses. Our value of \gamma_{DM} is in very good agreement with predictions from cosmological CDM simulations. We will extend our analysis to more clusters in future works. If confirmed on a larger cluster sample, our result makes this DM model more appealing than alternative models.


title: Visual link retrieval and knowledge discovery in painting datasets

abstract: Visual arts are of inestimable importance for the cultural, historic and economic growth of our society. One of the building blocks of most analysis in visual arts is to find similarity relationships among paintings of different artists and painting schools. To help art historians better understand visual arts, this paper presents a framework for visual link retrieval and knowledge discovery in digital painting datasets. Visual link retrieval is accomplished by using a deep convolutional neural network to perform feature extraction and a fully unsupervised nearest neighbor mechanism to retrieve links among digitized paintings. Historical knowledge discovery is achieved by performing a graph analysis that makes it possible to study influences among artists. An experimental evaluation on a database collecting paintings by very popular artists shows the effectiveness of the method. The unsupervised strategy makes the method interesting especially in cases where metadata are scarce, unavailable or difficult to collect.


title: Integrating State of the Art Compute, Communication, and Autotuning Strategies to Multiply the Performance of the Application Programm CPMD for Ab Initio Molecular Dynamics Simulations

abstract: We present our recent code modernizations of the of the ab initio molecular dynamics program CPMD (www.cpmd.org) with a special focus on the ultra-soft pseudopotential (USPP) code path. Following the internal instrumentation of CPMD, all time critical routines have been revised to maximize the computational throughput and to minimize the communication overhead for optimal performance. Throughout the program missing hybrid MPI+OpenMP parallelization has been added to optimize scaling. For communication intensive routines, as the multiple distributed 3d FFTs of the electronic states and distributed matrix-matrix multiplications related to the $\beta$-projectors of the pseudopotentials, this MPI+OpenMP parallelization now overlaps computation and communication. The necessary partitioning of the workload is optimized by an auto-tuning algorithm. In addition, the largest global MPI_Allreduce operation has been replaced by highly tuned node-local parallelized operations using MPI shared-memory windows to avoid inter-node communication. A batched algorithm for the multiple 3d FFTs improves the throughput of the MPI_Alltoall communication and, thus, the scalability of the implementation, both for USPP and for the frequently used norm-conserving pseudopotential code path. The enhanced performance and scalability is demonstrated on a mid-sized benchmark system of 256 water molecules and further water systems of from 32 up to 2048 molecules.


title: Correlation between the tolerance factor and phase transition in Ln4-xLnxNi3O10 (Ln and Ln' = La, Pr and Nd; x = 0, 1, 2 and 3)

abstract: We have synthesized a series of the Ruddlesden-Popper nickelate solid solution Ln4-xLn'xNi3O10 (Ln and Ln' = La, Pr and Nd; x = 0, 1, 2 and 3) via the citrate precursor method at different reacting atmospheres. Both the electronic-transport and magnetization measurements on these samples show well-defined phase transitions at temperatures between 135 K and 165 K. These transition temperatures, the room-temperature resistivities, as well as the changes of the Pauli-paramagnetic susceptibilities at the respective phase transitions, strongly correlate with the Goldschmidt tolerance factor t, irrespective of the combination of the magnetic rare-earth ions with unmagnetic La3+. We conclude that these changes of the electronic properties are mostly related to the distortion of the NiO6 octahedra at the phase transition which is strongly correlated with the tolerance factor t, but are rather insensitive to the magnetism of the rare-earth ions Ln3+ and Ln'3+.


title: Bi-stability of SUDR+K model of epidemics and test kits applied to COVID-19

abstract: Motivated with various responses of world governments to COVID-19, here we develop a toy model of the dependence epidemics spreading on the availability of tests for disease. Our model, that we call SUDR+K, is based on usual SIR model, but it splits the total fraction of infected individuals into two components: those that are undetected and those that are detected through tests. Moreover, we assume that available tests increase at a constant rate from the beginning of epidemics but are consumed to detect infected individuals. Strikingly we find a bi-stable behavior between a phase with a giant fraction of infected and a phase with a very small fraction. We show that the separation between these two regimes is governed by a match between the rate of testing and a rate of infection spread at given time. We also show that the existence of two phases does not depend on the mathematical choice of the form of the term describing the rate at which undetected individuals are tested and detected. Presented research implies that a vigorous early testing activity, before the epidemics enters into its giant phase, can potentially keep epidemics under control, and that even a very small change in rate of testing can increase or decrease the size of the whole epidemics of various orders of magnitude. For the real application of realistic model to ongoing epidemics, we would gladly collaborate with field epidemiologists in order to develop quantitative models of testing process.

categories: q-bio.PE cond-mat.stat-mech nlin.AO physics.soc-ph


title: A note on large Kakeya sets

abstract: A Kakeya set $\mathcal{K}$ in an affine plane of order $q$ is the point set covered by a set $\mathcal{L}$ of $q+1$ pairwise non-parallel lines. Large Kakeya sets were studied by Dover and Mellinger; in [6] they showed that Kakeya sets with size at least $q^2-3q+9$ contain a large knot (a point of $\mathcal{K}$ lying on many lines of $\mathcal{L}$). In this paper, we improve on this result by showing that Kakeya set of size at least $\approx q^2-q\sqrt{q}+\frac{3}{2}q$ contain a large knot. Furthermore, we obtain a sharp result for planes of square order containing a Baer subplane.


title: Secant indices of projective varieties

abstract: To each subvariety $X$ in projective $n$-space of codimension $m$ we associate an integer sequence of length $m + 1$ from $1$ to the degree of $X$ recording the maximal cardinalities of finite, reduced intersections of $X$ with linear subvarieties. We call this the sequence of secant indices of $X$. Similar numbers have been studied independently with the aim of classifying subvarieties with extremal secant spaces. Our focus in this note is the study of the combinatorial properties that the secant indices satisfy collectively. We show these sequences are strictly increasing for nondegenerate smooth subvarieties, develop a method to compute term-wise lower bounds for the secant indices, and compute these lower bounds for Veronese and Segre varieties. In the case of Veronese varieties, the truth of the Eisenbud-Green-Harris conjecture would imply the lower bounds we find are in fact equal to the secant indices. Along the way we state several relevant questions and additional conjectures which to our knowledge are open.


title: Tumbling with a Limp: Local Asymmetry in Water's Hydrogen Bond Network and its Consequences

abstract: Ab initio molecular dynamics simulations of liquid water under equilibrium ambient conditions, together with a novel energy decomposition analysis, have recently shown that a substantial fraction of water molecules exhibit a significant asymmetry between the strengths of the two donor and/or the two acceptor interactions. We refer to this recently unraveled aspect as the "local asymmetry in the hydrogen bond network". We discuss how this novel aspect was first revealed, and provide metrics that can be consistently employed on simulated water trajectories to quantify this local heterogeneity in the hydrogen bond network and its dynamics. We then discuss the static aspects of the asymmetry, pertaining to the frozen geometry of liquid water at any given instant of time and the distribution of hydrogen bond strengths therein, and also its dynamic characteristics pertaining to how fast this asymmetry decays and the kinds of molecular motions responsible for this decay. Following this we discuss the spectroscopic manifestations of this asymmetry, from ultrafast X-ray absorption spectra to infrared spectroscopy and down to the much slower terahertz regime. Finally, we discuss the implications of these findings in a broad context and its relation to the current notions about the structure and dynamics of liquid water.


title: Fault Handling in Large Water Networks with Online Dictionary Learning

abstract: Fault detection and isolation in water distribution networks is an active topic due to its model's mathematical complexity and increased data availability through sensor placement. Here we simplify the model by offering a data driven alternative that takes the network topology into account when performing sensor placement and then proceeds to build a network model through online dictionary learning based on the incoming sensor data. Online learning is fast and allows tackling large networks as it processes small batches of signals at a time and has the benefit of continuous integration of new data into the existing network model, be it in the beginning for training or in production when new data samples are encountered. The algorithms show good performance when tested on both small and large-scale networks.


title: Time-delay control for stabilization of the Shapovalov mid-size firm model

abstract: Control and stabilization of irregular and unstable behavior of dynamic systems (including chaotic processes) are interdisciplinary problems of interest to a variety of scientific fields and applications. Using the control methods allows improvements in forecasting the dynamics of unstable economic processes and offers opportunities for governments, central banks, and other policy makers to modify the behaviour of the economic system to achieve its best performance. One effective method for control of chaos and computation of unstable periodic orbits (UPOs) is the unstable delay feedback control (UDFC) approach, suggested by K. Pyragas. This paper proposes the application of the Pyragas' method within framework of economic models. We consider this method through the example of the Shapovalov model, by describing the dynamics of a mid-size firm. The results demonstrate that suppressing chaos is capable in the Shapovalov model, using the UDFC method.


title: Self-Supervised Contextual Bandits in Computer Vision

abstract: Contextual bandits are a common problem faced by machine learning practitioners in domains as diverse as hypothesis testing to product recommendations. There have been a lot of approaches in exploiting rich data representations for contextual bandit problems with varying degree of success. Self-supervised learning is a promising approach to find rich data representations without explicit labels. In a typical self-supervised learning scheme, the primary task is defined by the problem objective (e.g. clustering, classification, embedding generation etc.) and the secondary task is defined by the self-supervision objective (e.g. rotation prediction, words in neighborhood, colorization, etc.). In the usual self-supervision, we learn implicit labels from the training data for a secondary task. However, in the contextual bandit setting, we don't have the advantage of getting implicit labels due to lack of data in the initial phase of learning. We provide a novel approach to tackle this issue by combining a contextual bandit objective with a self supervision objective. By augmenting contextual bandit learning with self-supervision we get a better cumulative reward. Our results on eight popular computer vision datasets show substantial gains in cumulative reward. We provide cases where the proposed scheme doesn't perform optimally and give alternative methods for better learning in these cases.


title: Complete flavor decomposition of the spin and momentum fraction of the proton using lattice QCD simulations at physical pion mass

abstract: We evaluate the gluon and quark contributions to the spin of the proton using an ensemble of gauge configuration generated at physical pion mass. We compute all valence and sea quark contributions to high accuracy. We perform a non-perturbative renormalization for both quark and gluon matrix elements. We find that the contribution of the up, down, strange and charm quarks to the proton intrinsic spin is $\frac{1}{2}\sum_{q=u,d,s,c}\Delta\Sigma^{q^+}=0.191(15)$ and to the total spin $\sum_{q=u,d,s,c}J^{q^+}=0.285(45)$. The gluon contribution to the spin is $J^g=0.187(46)$ yielding $J=J^q+J^g=0.473(71)$ confirming the spin sum. The momentum fraction carried by quarks in the proton is found to be $0.618(60)$ and by gluons $0.427(92)$, the sum of which gives $1.045(118)$ confirming the momentum sum rule. All scale and scheme dependent quantities are given in the $\mathrm{ \overline{MS}}$ scheme at 2 GeV.


title: Effect of different Zinc precursors in Structural and Optical properties of ZnO thin films

abstract: Zinc oxide thin films were deposited by ultrasonic spray pyrolysis technique using different Zinc precursors: Zinc acetate, Zinc nitrate and zinc chloride with different molarities. All deposited films were characterized by various techniques such as X-ray diffraction to determine the films structure, the scanning electron microscopy SEM for the morphology of the surfaces and UV visible spectroscopy to determine the optical proprieties.The obtained results indicated that ZnO nano-crystalline films have a hexagonal structure at type wurtzite for all Zinc precursors. Films deposited with zinc acetate are characterized by a smooth surface, dense network and high transparency, while films deposited with zinc chloride and Zinc Nitrate have a better crystallinity and low optical transmittance.


title: Magnetic Fabry-P\'erot interferometer for valley filtering in a honeycomb-dice model

abstract: Here we theoretically investigate the valley-dependent transmission of particles through a combined electric and magnetic barrier in the $\alpha-T_3$ model which interpolates between the honeycomb and the dice lattices. We put forward that the combination of the Fabry-P\'erot interferences and the magnetic field leads to a perfect transmission for one valley and a suppression of the transmission for the other valley. When only one of the barriers (magnetic or electric) is present, no valley polarized current can be produced. By tuning the Fermi energy, this valley-dependent peculiar behavior can be used as valley filtering. Our results show that highly efficient valley filtering with maximum conductivity and polarization can be achieved by controlling the value of the magnetic field and the electric barrier width and height.


title: An Analysis on the Learning Rules of the Skip-Gram Model

abstract: To improve the generalization of the representations for natural language processing tasks, words are commonly represented using vectors, where distances among the vectors are related to the similarity of the words. While word2vec, the state-of-the-art implementation of the skip-gram model, is widely used and improves the performance of many natural language processing tasks, its mechanism is not yet well understood. In this work, we derive the learning rules for the skip-gram model and establish their close relationship to competitive learning. In addition, we provide the global optimal solution constraints for the skip-gram model and validate them by experimental results.


title: Recent results of nucleon structure & matrix element calculations

abstract: A review of recent lattice calculations of nucleon structure and matrix elements of operators in nucleons is presented. It primarily covers developments in the calculation of the matrix elements of the scalar, tensor, pseudo-scalar, axial-vector and vector operators in the ground state of neutrons and protons in the isospin symmetric limit. Some preliminary calculations of the electric dipole moment, the gravitational moments and stress-energy distribution, and the magnetic polarizability are briefly described.


title: Decay Properties of Conventional and Hybrid $B_c$ Mesons

abstract: Spectrum, radial wave functions at origin, decay constants, weak decay widths, life time and branching ratios for radially excited conventional and hybrid $B_c$ mesons are derived within non-relativistic quark model framework employing Schr$\ddot{\textrm{o}}$dinger equation by shooting method. Calculated results are compared with available theoretical results and experimental observations. This work may help in identifying the new discovered $B_c$ meson states at CDF ,LHCb and ATLAS.


title: Homogenization of Layered Media: Intrinsic and Extrinsic Symmetry Breaking

abstract: A general homogenization procedure for periodic electromagnetic structures, when applied to layered media with asymmetric lattice cells, yields an effective tensor with magnetoelectric coupling. Accurate results for transmission and reflection are obtained even in cases where classical effective medium theory breaks down. Magnetoelectric coupling accounts for symmetry breaking in reflection and transmission when a non-symmetric structure is illuminated from two opposite sides.


title: Lattice QCD Inputs for Nuclear Double Beta Decay

abstract: Second order beta-decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double-beta decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double-beta decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double-beta decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double-beta decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn\to pp e^- e^- \bar\nu_e\bar\nu_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double-beta decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double-beta decay processes themselves or using inputs from LQCD. Future prospects for direct calculations in larger nuclei are also discussed.


title: Progress Extrapolating Algorithmic Learning to Arbitrary Sequence Lengths

abstract: Recent neural network models for algorithmic tasks have led to significant improvements in extrapolation to sequences much longer than training, but it remains an outstanding problem that the performance still degrades for very long or adversarial sequences. We present alternative architectures and loss-terms to address these issues, and our testing of these approaches has not detected any remaining extrapolation errors within memory constraints. We focus on linear time algorithmic tasks including copy, parentheses parsing, and binary addition. First, activation binning was used to discretize the trained network in order to avoid computational drift from continuous operations, and a binning-based digital loss term was added to encourage discretizable representations. In addition, a localized differentiable memory (LDM) architecture, in contrast to distributed memory access, addressed remaining extrapolation errors and avoided unbounded growth of internal computational states. Previous work has found that algorithmic extrapolation issues can also be alleviated with approaches relying on program traces, but the current effort does not rely on such traces.


title: Hydrodynamic limit of the Kob-Andersen model

abstract: This paper concerns with the hydrodynamic limit of the Kob-Andersen model, an interacting particle system that has been introduced by physicists in order to explain glassy behavior, and widely studies since. We will see that the density profile evolves in the hydrodynamic limit according to a non-degenerate hydrodynamic equation, and see how the diffusion coefficient decays as density grows.


title: Electrically small matched antennas with time-periodic and space-uniform modulation

abstract: A conceptual method to realize an impedance-matched electrically small antenna is discussed, using a transmission line modulated only in time. The time-periodic space-uniform modulation of a distributed shunt array of varcapacitors excites a mode with very small frequency but still a small (microwave) wavelength: the matching line then retains a practical size. The modulated transmission line brings the problem of impedance matching to a frequency $f_0$ much higher than the signal, where the large antenna reactance can be easily compensated with a much higher bandwidth for the same $Q$ factor: actually the bandwidth $\Delta f= f_0 / Q$ becomes large when $f_0$ is much larger than the signal frequency. A simple way to obtain the space uniform modulation of the transmission line is illustrated, making use of a Wheatstone-bridge arrangement of the modulated varicaps. The signal and modulating lines, running in parallel on a section of the same length, see two very different equivalent shunt capacitances, so that the modulating voltage is almost uniform while the signal voltage may change significantly across the transmission line.


title: Mixing time of the switch Markov chain and stable degree sequences

abstract: The switch chain is a well-studied Markov chain which can be used to sample approximately uniformly from the set $\Omega(\boldsymbol{d})$ of all graphs with a given degree sequence $\boldsymbol{d}$. Polynomial mixing time (rapid mixing) has been established for the switch chain under various conditions on the degree sequences. Amanatidis and Kleer introduced the notion of strongly stable families of degree sequences, and proved that the switch chain is rapidly mixing for any degree sequence from a strongly stable family. Using a different approach, Erd\H{o}s et al. recently extended this result to the (possibly larger) class of P-stable degree sequences, introduced by Jerrum and Sinclair in 1990. We define a new notion of stability for a given degree sequence, namely $k$-\emph{stability}, and prove that if a degree sequence $\boldsymbol{d}$ is 8-stable then the switch chain on $\Omega(\boldsymbol{d})$ is rapidly mixing. We also provide sufficient conditions for P-stability, strong stability and 8-stability. Using these sufficient conditions, we give the first proof of P-stability for various families of heavy-tailed degree sequences, including power-law degree sequences, and show that the switch chain is rapidly mixing for these families. We further extend these notions and results to directed degree sequences.


title: Directed transport of suspended ferromagnetic nanoparticles under both gradient and uniform magnetic fields

abstract: The suspended ferromagnetic particles subjected to the gradient and uniform magnetic fields experience both the translational force generated by the field gradient and the rotational torque generated by the fields strengths. Although the uniform field does not contribute to the force, it nevertheless influences the translational motion of these particles. This occurs because the translational force depends on the direction of the particle magnetization, which in turn depends on the fields strengths. To study this influence, a minimal set of equations describing the coupled translational and rotational motions of nanosized ferromagnetic particles is introduced and solved in the low Reynolds number approximation. Trajectory analysis reveals that, depending on the initial positions of nanoparticles, there exist four regimes of their directed transport. The intervals of initial positions that correspond to different dynamical regimes are determined, their dependence on the uniform magnetic field is established, and strong impact of this field on the directed transport is demonstrated. The ability and efficiency of the uniform magnetic field to control the separation of suspended ferromagnetic nanoparticles is also discussed.


title: Gaze-Sensing LEDs for Head Mounted Displays

abstract: We introduce a new gaze tracker for Head Mounted Displays (HMDs). We modify two off-the-shelf HMDs to be gaze-aware using Light Emitting Diodes (LEDs). Our key contribution is to exploit the sensing capability of LEDs to create low-power gaze tracker for virtual reality (VR) applications. This yields a simple approach using minimal hardware to achieve good accuracy and low latency using light-weight supervised Gaussian Process Regression (GPR) running on a mobile device. With our hardware, we show that Minkowski distance measure based GPR implementation outperforms the commonly used radial basis function-based support vector regression (SVR) without the need to precisely determine free parameters. We show that our gaze estimation method does not require complex dimension reduction techniques, feature extraction, or distortion corrections due to off-axis optical paths. We demonstrate two complete HMD prototypes with a sample eye-tracked application, and report on a series of subjective tests using our prototypes.


title: The Cost of Privacy in Asynchronous Differentially-Private Machine Learning

abstract: We consider training machine learning models using Training data located on multiple private and geographically-scattered servers with different privacy settings. Due to the distributed nature of the data, communicating with all collaborating private data owners simultaneously may prove challenging or altogether impossible. In this paper, we develop differentially-private asynchronous algorithms for collaboratively training machine-learning models on multiple private datasets. The asynchronous nature of the algorithms implies that a central learner interacts with the private data owners one-on-one whenever they are available for communication without needing to aggregate query responses to construct gradients of the entire fitness function. Therefore, the algorithm efficiently scales to many data owners. We define the cost of privacy as the difference between the fitness of a privacy-preserving machine-learning model and the fitness of trained machine-learning model in the absence of privacy concerns. We prove that we can forecast the performance of the proposed privacy-preserving asynchronous algorithms. We demonstrate that the cost of privacy has an upper bound that is inversely proportional to the combined size of the training datasets squared and the sum of the privacy budgets squared. We validate the theoretical results with experiments on financial and medical datasets. The experiments illustrate that collaboration among more than 10 data owners with at least 10,000 records with privacy budgets greater than or equal to 1 results in a superior machine-learning model in comparison to a model trained in isolation on only one of the datasets, illustrating the value of collaboration and the cost of the privacy. The number of the collaborating datasets can be lowered if the privacy budget is higher.


title: On Broadcasting Time in the Model of Travelling Agents

abstract: Consider the following broadcasting process run on a connected graph $G=(V,E)$. Suppose that $k \ge 2$ agents start on vertices selected from $V$ uniformly and independently at random. One of the agents has a message that she wants to communicate to the other agents. All agents perform independent random walks on $G$, with the message being passed when an agent that knows the message meets an agent that does not know the message. The broadcasting time $\xi(G,k)$ is the time it takes to spread the message to all agents. Our ultimate goal is to gain a better understanding of the broadcasting process run on real-world networks of roads of large cities that might shed some light on the behaviour of future autonomous and connected vehicles. Due to the complexity of road networks, such phenomena have to be studied using simulation in practical applications. In this paper, we study the process on the simplest scenario, i.e., the family of complete graphs, as in this case the problem is analytically tractable. We provide tight bounds for $\xi(K_n,k)$ that hold asymptotically almost surely for the whole range of the parameter $k$. These theoretical results reveal interesting relationships and, at the same time, are also helpful to understand and explain the behaviour we observe in more realistic networks.


title: Reconstructing Sinus Anatomy from Endoscopic Video -- Towards a Radiation-free Approach for Quantitative Longitudinal Assessment

abstract: Reconstructing accurate 3D surface models of sinus anatomy directly from an endoscopic video is a promising avenue for cross-sectional and longitudinal analysis to better understand the relationship between sinus anatomy and surgical outcomes. We present a patient-specific, learning-based method for 3D reconstruction of sinus surface anatomy directly and only from endoscopic videos. We demonstrate the effectiveness and accuracy of our method on in and ex vivo data where we compare to sparse reconstructions from Structure from Motion, dense reconstruction from COLMAP, and ground truth anatomy from CT. Our textured reconstructions are watertight and enable measurement of clinically relevant parameters in good agreement with CT. The source code is available at https://github.com/lppllppl920/DenseReconstruction-Pytorch.


title: Area preserving surface diffeomorphisms with polynomial decay of correlations are ubiquitous

abstract: We show that any surface admits an area preserving $C^{1+\beta}$ diffeomorphism with non-zero Lyapunov exponents which is Bernoulli and has polynomial decay of correlations. We establish both upper and lower polynomial bounds on correlations. In addition, we show that this diffeomorphism satisfies the Central Limit Theorem and has the Large Deviation Property. Finally, we show that the diffeomorphism we constructed possesses a unique hyperbolic Bernoulli measure of maximal entropy with respect to which it has exponential decay of correlations.


title: A One Dimensional Elliptic Distributed Optimal Control Problem with Pointwise Derivative Constraints

abstract: We consider a one dimensional elliptic distributed optimal control problem with pointwise constraints on the derivative of the state. By exploiting the variational inequality satisfied by the derivative of the optimal state, we obtain higher regularity for the optimal state under appropriate assumptions on the data. We also solve the optimal control problem as a fourth order variational inequality by a $C^1$ finite element method, and present the error analysis together with numerical results.


title: A Metric Learning Reality Check

abstract: Deep metric learning papers from the past four years have consistently claimed great advances in accuracy, often more than doubling the performance of decade-old methods. In this paper, we take a closer look at the field to see if this is actually true. We find flaws in the experimental methodology of numerous metric learning papers, and show that the actual improvements over time have been marginal at best.


title: Superconducting Dome in Nd$_{1-x}$Sr$_x$NiO$_2$ Infinite Layer Films

abstract: We report the phase diagram of Nd$_{1-x}$Sr$_x$NiO$_2$ infinite layer thin films grown on SrTiO$_3$. A superconducting dome spanning $0.125 < x < 0.25$ is found, remarkably similar to cuprates, albeit over a narrower doping window. However, while cuprate superconductivity is bounded by an insulator for underdoping and a metal for overdoping, here we observe weakly insulating behavior on either side of the dome. Furthermore, the normal state Hall coefficient is always small and proximate to a continuous zero crossing in doping and in temperature, in contrast to the $\sim 1/x$ dependence observed for cuprates. This suggests the presence of both electron- and hole-like bands, consistent with band structure calculations.


title: Coupled Control Systems: Periodic Orbit Generation with Application to Quadrupedal Locomotion

abstract: A robotic system can be viewed as a collection of lower-dimensional systems that are coupled via reaction forces (Lagrange multipliers) enforcing holonomic constraints. Inspired by this viewpoint, this paper presents a novel formulation for nonlinear control systems that are subject to coupling constraints via virtual "coupling" inputs that abstractly play the role of Lagrange multipliers. The main contribution of this paper is a process---mirroring solving for Lagrange multipliers in robotic systems---wherein we isolate subsystems free of coupling constraints that provably encode the full-order dynamics of the coupled control system from which it was derived. This dimension reduction is leveraged in the formulation of a nonlinear optimization problem for the isolated subsystem that yields periodic orbits for the full-order coupled system. We consider the application of these ideas to robotic systems, which can be decomposed into subsystems. Specifically, we view a quadruped as a coupled control system consisting of two bipedal robots, wherein applying the framework developed allows for gaits (periodic orbits) to be generated for the individual biped yielding a gait for the full-order quadruped. This is demonstrated through walking experiments of a quadrupedal robot in simulation and on rough terrains.


title: An Application of Gaussian Process Modeling for High-order Accurate Adaptive Mesh Refinement Prolongation

abstract: We present a new polynomial-free prolongation scheme for Adaptive Mesh Refinement (AMR) simulations of compressible and incompressible computational fluid dynamics. The new method is constructed using a multi-dimensional kernel-based Gaussian Process (GP) prolongation model. The formulation for this scheme was inspired by the GP methods introduced by A. Reyes et al. (A New Class of High-Order Methods for Fluid Dynamics Simulation using Gaussian Process Modeling, Journal of Scientific Computing, 76 (2017), 443-480; A variable high-order shock-capturing finite difference method with GP-WENO, Journal of Computational Physics, 381 (2019), 189-217). In this paper, we extend the previous GP interpolations and reconstructions to a new GP-based AMR prolongation method that delivers a high-order accurate prolongation of data from coarse to fine grids on AMR grid hierarchies. In compressible flow simulations special care is necessary to handle shocks and discontinuities in a stable manner. To meet this, we utilize the shock handling strategy using the GP-based smoothness indicators developed in the previous GP work by A. Reyes et al. We demonstrate the efficacy of the GP-AMR method in a series of testsuite problems using the AMReX library, in which the GP-AMR method has been implemented.

categories: math.NA astro-ph.IM cs.NA physics.flu-dyn