title: Google Summer of Code: Student Motivations and Contributions

abstract: Several open source software (OSS) projects expect to foster newcomers' onboarding and to receive contributions by participating in engagement programs, like Summers of Code. However, there is little empirical evidence showing why students join such programs. In this paper, we study the well-established Google Summer of Code (GSoC), which is a 3-month OSS engagement program that offers stipends and mentors to students willing to contribute to OSS projects. We combined a survey (students and mentors) and interviews (students) to understand what motivates students to enter GSoC. Our results show that students enter GSoC for an enriching experience, not necessarily to become frequent contributors. Our data suggest that, while the stipends are an important motivator, the students participate for work experience and the ability to attach the name of the supporting organization to their resum\'es. We also discuss practical implications for students, mentors, OSS projects, and Summer of Code programs.


title: Optical functions in de Sitter

abstract: This paper addresses pure gauge questions in the study of (asymptotically) de Sitter spacetimes. We construct global solutions to the eikonal equation on de Sitter, whose level sets give rise to double null foliations, and give detailed estimates for the structure coefficients in this gauge. We show two results which are relevant for the foliations used in the stability problem of the expanding region of Schwarzschild de Sitter spacetimes [Schlue (2016)]: (i) Small perturbations of round spheres on the cosmological horizons lead to spheres that pinch off at infinity. (ii) Globally well behaved double null foliations can be constructed from infinity using a choice of spheres related to the level sets of a mass aspect function. While (i) shows that in the above stability problem a final gauge choice is necessary, the proof of (ii) already outlines a strategy for the case of spacetimes with decaying, instead of vanishing, conformal Weyl curvature.


title: Improving Precision through Adjustment for Prognostic Variables in Group Sequential Trial Designs: Impact of Baseline Variables, Short-Term Outcomes, and Treatment Effect Heterogeneity

abstract: In randomized trials, appropriately adjusting for baseline variables and short-term outcomes can lead to increased precision and reduced sample size. We examine the impact of such adjustment in group sequential designs, i.e., designs with preplanned interim analyses where enrollment may be stopped early for efficacy or futility. We address the following questions: how much precision gain can be obtained by appropriately adjusting for baseline variables and a short-term outcome? How is this precision gain impacted by factors such as the proportion of pipeline participants (those who enrolled but haven't yet had their primary outcomes measured) and treatment effect heterogeneity? What is the resulting impact on power and average sample size in a group sequential design? We derive an asymptotic formula that decomposes the overall precision gain from adjusting for baseline variables and a short-term outcome into contributions from factors mentioned above, for efficient estimators in the model that only assumes randomization and independent censoring. We use our formula to approximate the precision gain from a targeted minimum loss-based estimator applied to data from a completed trial of a new surgical intervention for stroke. Our formula implies that (for an efficient estimator) adjusting for a prognostic baseline variable leads to at least as much asymptotic precision gain as adjusting for an equally prognostic short-term outcome. In many cases, such as our stroke trial application, the former leads to substantially greater precision gains than the latter. In our simulation study, we show how precision gains from adjustment can be converted into sample size reductions (even when there is no treatment effect).


title: Effect of Voltage Source Converters with Electrochemical Storage Systems on Dynamics of Reduced-inertia Bulk Power Grids

abstract: A major concern associated to the massive connection of distributed energy resources is the increasing share of power electronic interfaces resulting in the global inertia reduction of power systems. The recent literature advocated the use of voltage source converter (VSC) interfaced battery energy storage system (BESS) as a potential way to counterbalance this lack of inertia. However, the impact of VSCs on the dynamics of reduced-inertia grids is not well understood especially with respect to large transmission grids interfacing a mix of rotating machines and resources interfaced with power electronics. In this regards, we propose an extension of the IEEE 39-bus test network used to quantify the impact of VSCs on reduced-inertia grids. In this respect, a reduced-inertia 39-bus system is obtained by replacing 4 synchronous generators in the original 10-synchronous machine system, with 4 wind power plants modeled as aggregated type-3 wind turbines. Then, a large-scale BESS is integrated into the reduced-inertia network via a three-level neutral-point clamped (NPC) converter, thereby to be used for studying the impact of VSC on the dynamics of the inertia-reduced power system, as well as for comparing different VSC controls. The proposed models are implemented on a real-time simulator to conduct post-contingency analysis, respectively, for the original power system and the reduced-inertia one, with and without the BESS-VSC.


title: Universality in the emergence of oscillatory instabilities in turbulent flows

abstract: Spontaneous emergence of periodic oscillations due to self-organization is ubiquitous in turbulent flows. The emergence of such oscillatory instabilities in turbulent fluid mechanical systems is often studied in different system-specific frameworks. We uncover the existence of a universal scaling behaviour during self-organization in turbulent flows leading to oscillatory instability. Our experiments show that the spectral amplitude of the dominant mode of oscillations scales inversely with the Hurst exponent of a fluctuating state variable following an inverse power law relation. Interestingly, we observe the same power law behaviour with a constant exponent near -2 across various turbulent systems such as aeroacoustic, thermoacoustic and aeroelastic systems.


title: Lifetimes and Oscillator Strengths for Ultraviolet Transitions in Neutral Chlorine

abstract: We present lifetime measurements using beam-foil techniques for radiative transitions from the 3$p^4$($^1S$)4$s$ $^2S_{1/2}$, 3$p^4$($^3P$)5$s$ $^2P_{1/2,3/2}$, and 3$p^4$($^3P$)3$d$ $^2F_{5/2}$ levels in Cl I and the corresponding results of the oscillator strengths for transitions at 1004.68, 1079.88, 1090.73, and 1094.77 \AA, respectively. We compare our experimental results with available theoretical calculations and astronomical observations in an effort to resolve discrepancies among them.


title: The Role of Embedding Complexity in Domain-invariant Representations

abstract: Unsupervised domain adaptation aims to generalize the hypothesis trained in a source domain to an unlabeled target domain. One popular approach to this problem is to learn domain-invariant embeddings for both domains. In this work, we study, theoretically and empirically, the effect of the embedding complexity on generalization to the target domain. In particular, this complexity affects an upper bound on the target risk; this is reflected in experiments, too. Next, we specify our theoretical framework to multilayer neural networks. As a result, we develop a strategy that mitigates sensitivity to the embedding complexity, and empirically achieves performance on par with or better than the best layer-dependent complexity tradeoff.


title: On Pro-$2$ Identities of $2\times2$ Linear Groups

abstract: Let $\hat{F}$ be a free pro-$p$ non-abelian group, and let $\Delta$ be a commutative Noetherian complete local ring with a maximal ideal $I$ such that $\textrm{char}(\Delta/I)=p>0$. In [Zu], Zubkov showed that when $p\neq2$, the pro-$p$ congruence subgroup $$GL_{2}^{1}(\Delta)=\ker(GL_{2}(\Delta)\overset{\Delta\to\Delta/I}{\longrightarrow}GL_{2}(\Delta/I))$$ admits a pro-$p$ identity, i.e., there exists an element $1\neq w\in\hat{F}$ that vanishes under any continuous homomorphism $\hat{F}\to GL_{2}^{1}(\Delta)$. In this paper we investigate the case $p=2$. The main result is that when $\textrm{char}(\Delta)=2$, the pro-$2$ group $GL_{2}^{1}(\Delta)$ admits a pro-$2$ identity. This result was obtained by the use of trace identities that originate in PI-theory.


title: Phonon-assisted diffusion in bcc phase of titanium and zirconium from first-principles

abstract: Diffusion is the underlying mechanism for many complicated materials phenomena, and understanding it is basic to the discovery of novel materials with desired physical and mechanical properties. Certain groups of solid phases, such as the bcc phase of IIIB and IVB metals and their alloys, which are only stable when they reach high enough temperatures and experience anharmonic vibration entropic effects, exhibit "anomalously fast diffusion". However, the underlying reason for the observed extraordinary fast diffusion is poorly understood and due to the existence of harmonic vibration instabilities in these phases the standard models fail to predict their diffusivity. Here, we indicate that the anharmonic phonon-phonon coupling effects can accurately describe the anomalously large macroscopic diffusion coefficients in the bcc phase of IVB metals, and therefore yield a new understanding of the underlying mechanism for diffusion in these phases. We utilize temperature-dependent phonon analysis by combining ab initio molecular dynamics with lattice dynamics calculations to provide a new approach to use the transition state theory beyond the harmonic approximation. We validate the diffusivity predictions for the bcc phase of titanium and zirconium with available experimental measurements, while we show that predictions based on harmonic transition state theory severely underestimates diffusivity in these phases.


title: Shared E-scooters: Business, Pleasure, or Transit?

abstract: Shared e-scooters have become a familiar sight in many cities around the world. Yet the role they play in the mobility space is still poorly understood. This paper presents a study of the use of Bird e-scooters in the city of Atlanta. Starting with raw data which contains the location of available Birds over time, the study identifies trips and leverages the Google Places API to associate each trip origin and destination with a Point of Interest (POI). The resulting trip data is then used to understand the role of e-scooters in mobility by clustering trips using 10 collections of POIs, including business, food and recreation, parking, transit, health, and residential. The trips between these POI clusters reveal some surprising, albeit sensible, findings about the role of e-scooters in mobility, as well as the time of the day where they are most popular.


title: Construct $\alpha^{\prime}$ corrected or loop corrected solutions without curvature singularities

abstract: For the bosonic gravi-dilaton system, we provide systematical approaches to construct non-perturbative string cosmological solutions without curvature singularities, which can match the perturbative solution to any order in $\alpha^{\prime}$ expansion. When higher order $\alpha^{\prime}$ corrections are calculated, they can be straightforwardly plugged in to generate compatible non-perturbative evolutions without curvature singularities. We also give a (phenomenological) map between $\alpha^{\prime}$ corrected EOM and loop corrected EOM. This map enables us to easily generate a loop corrected solution from an $\alpha^{\prime}$ corrected solution, and vice versa, therefore substantially enlarges the solution space.


title: Convexity of constant mean curvature graphs in $\mathbb{R}^{n+1}$ with planar boundary

abstract: We study the Dirichlet problem for a graph $\Sigma$ in $\mathbb{R}^{n+1}$ with normalized constant mean curvature $H>0$ and planar boundary $\Gamma=\partial \Omega$. Our main result is that the optimal solvability condition, namely that the normalized mean curvature $h$ of $\Gamma$ satisfies $h\geq H$, also suffices when $\Omega$ is strictly convex, to prove the strict convexity of $\Sigma$.


title: Deep Crowd-Flow Prediction in Built Environments

abstract: Predicting the behavior of crowds in complex environments is a key requirement in a multitude of application areas, including crowd and disaster management, architectural design, and urban planning. Given a crowd's immediate state, current approaches simulate crowd movement to arrive at a future state. However, most applications require the ability to predict hundreds of possible simulation outcomes (e.g., under different environment and crowd situations) at real-time rates, for which these approaches are prohibitively expensive. In this paper, we propose an approach to instantly predict the long-term flow of crowds in arbitrarily large, realistic environments. Central to our approach is a novel CAGE representation consisting of Capacity, Agent, Goal, and Environment-oriented information, which efficiently encodes and decodes crowd scenarios into compact, fixed-size representations that are environmentally lossless. We present a framework to facilitate the accurate and efficient prediction of crowd flow in never-before-seen crowd scenarios. We conduct a series of experiments to evaluate the efficacy of our approach and showcase positive results.


title: Towards Efficient Discrete Integration via Adaptive Quantile Queries

abstract: Discrete integration in a high dimensional space of n variables poses fundamental challenges. The WISH algorithm reduces the intractable discrete integration problem into n optimization queries subject to randomized constraints, obtaining a constant approximation guarantee. The optimization queries are expensive, which limits the applicability of WISH. We propose AdaWISH, which is able to obtain the same guarantee but accesses only a small subset of queries of WISH. For example, when the number of function values is bounded by a constant, AdaWISH issues only O(log n) queries. The key idea is to query adaptively, taking advantage of the shape of the weight function being integrated. In general, we prove that AdaWISH has a regret of only O(log n) relative to an idealistic oracle that issues queries at data-dependent optimal points. Experimentally, AdaWISH gives precise estimates for discrete integration problems, of the same quality as that of WISH and better than several competing approaches, on a variety of probabilistic inference benchmarks. At the same time, it saves substantially on the number of optimization queries compared to WISH. On a suite of UAI inference challenge benchmarks, it saves 81.5% of WISH queries while retaining the quality of results.


title: Spectral identities for Schr\"{o}dinger operators

abstract: We obtain a system of identities relating boundary coefficients and spectral data for the one-dimensional Schr\"{o}dinger equation with boundary conditions containing rational Herglotz--Nevanlinna functions of the eigenvalue parameter. These identities can be thought of as a kind of mini version of the Gelfand--Levitan integral equation for boundary coefficients only.


title: First-order transition in a model of prestige bias

abstract: One of the major benefits of belonging to a prestigious group is that it affects the way you are viewed by others. Here I use a simple mathematical model to explore the implications of this "prestige bias" when candidates undergo repeated rounds of evaluation. In the model, candidates who are evaluated most highly are admitted to a "prestige class", and their membership biases future rounds of evaluation in their favor. I use the language of Bayesian inference to describe this bias, and show that it can lead to a runaway effect in which the weight given to the prior expectation associated with a candidate's class becomes stronger with each round. Most dramatically, the strength of the prestige bias after many rounds undergoes a first-order transition as a function of the precision of the examination on which the evaluation is based.


title: Discovering a sparse set of pairwise discriminating features in high dimensional data

abstract: Extracting an understanding of the underlying system from high dimensional data is a growing problem in science. Discovering informative and meaningful features is crucial for clustering, classification, and low dimensional data embedding. Here we propose to construct features based on their ability to discriminate between clusters of the data points. We define a class of problems in which linear separability of clusters is hidden in a low dimensional space. We propose an unsupervised method to identify the subset of features that define a low dimensional subspace in which clustering can be conducted. This is achieved by averaging over discriminators trained on an ensemble of proposed cluster configurations. We then apply our method to single cell RNA-seq data from mouse gastrulation, and identify 27 key transcription factors (out of 409 total), 18 of which are known to define cell states through their expression levels. In this inferred subspace, we find clear signatures of known cell types that eluded classification prior to discovery of the correct low dimensional subspace.


title: An Efficient Beam and Channel Acquisition via Sparsity Map and Joint Angle-Delay Power Profile Estimation for Wideband Massive MIMO Systems

abstract: In this paper, an efficient beam and channel acquisition scheme together with joint angle-delay power profile (JADPP) construction are proposed for single-carrier mm-wave wideband sparse massive multiple-input multiple-output (MIMO) channels when hybrid beamforming architecture is utilized. We consider two different modes of operation, namely slow-time beam acquisition and fast-time instantaneous channel estimation, for training stage of time division duplex based systems. In the first mode, where pre-structured hybrid beams are formed to scan intended angular sectors, the joint angle-delay sparsity map together with power intensities of each user channels are obtained by using a novel constant false alarm rate thresholding algorithm inspired from adaptive radar detection theory. The proposed thresholding algorithm employs a spatio-temporal adaptive matched filter type estimator, taking the strong interference due to simultaneously active multipath components of different user channels into account, in order to estimate JADPP of each user. After applying the proposed thresholding algorithm on the estimated power profile, the angle-delay sparsity map of the massive MIMO channel is constructed, based on which the channel covariance matrices (CCMs) are formed with significantly reduced amount of training snapshots. Then, by using the estimated CCMs, the analog beamformer is reconstructed by means of a virtual sectorization while taking the inter-group and inter-symbol interference into account. Finally, for the second mode of operation, two novel reduced-rank instantaneous channel estimators, operating in a proper beamspace formed by the hybrid structure, are proposed. The proposed beam and channel acquisition techniques attain the channel estimation accuracy of minimum mean square error filter with true knowledge of CCMs.


title: Homomorphisms from Functional Equations: The Goldie Equation, II

abstract: In this sequel to arXiv1407.4089 by the second author, we extend to multi-dimensional (or infinite-dimensional) settings the Goldie equation arising in the general regular variation of `General regular variation, Popa groups and quantifier weakening', J. Math. Anal. Appl. 483 (2020) 123610, 31 pp. (arXiv1901.05996). We extend the work there in two directions. First, we algebraicize the theory, by systematic use of certain groups -- the Popa groups arising in earlier work by Popa, and their relatives the Javor groups. Secondly, we extend from the original context on the real line to multi-dimensional (or infinite-dimensional) settings.


title: The Goldie Equation: III. Homomorphisms from Functional Equations

abstract: This is the second of three sequels to arXiv1407.4089 -- the third of the resulting quartet -- concerning the real-valued continuous solutions of the multivariate Goldie functional equation (GFE) below of Levi-Civita type. Following on from the preceding paper arXiv1910.05816, in which these solutions are described explicitly, here we characterize (GFE) as expressing homomorphy (in all but some exceptional "improper" cases) between multivariate Popa groups, defined and characterized earlier in the sequence. The group operation involves a form of affine addition (with local scalar acceleration) similar to the circle operation of ring theory. We show the affine action in (GFE) may be replaced by general continuous acceleration yielding a functional equation (GGE) which it emerges has the same solution structure as (GFE). The final member of the sequence arXiv2105.07794 considers the richer framework of a Banach algebra which allows vectorial acceleration, giving the closest possible similarity to the circle operation.


title: A Bayesian Statistics Course for Undergraduates: Bayesian Thinking, Computing, and Research

abstract: We propose a semester-long Bayesian statistics course for undergraduate students with calculus and probability background. We cultivate students' Bayesian thinking with Bayesian methods applied to real data problems. We leverage modern Bayesian computing techniques not only for implementing Bayesian methods, but also to deepen students' understanding of the methods. Collaborative case studies further enrich students' learning and provide experience to solve open-ended applied problems. The course has an emphasis on undergraduate research, where accessible academic journal articles are read, discussed, and critiqued in class. With increased confidence and familiarity, students take the challenge of reading, implementing, and sometimes extending methods in journal articles for their course projects.


title: Modeling Nanoribbon Peeling

abstract: The lifting, peeling and exfoliation of physisorbed ribbons (or flakes) of 2D material such as graphene off a solid surface are common and important manoeuvres in nanoscience. The feature that makes this case peculiar is the structural lubricity generally realized by stiff 2D material contacts. We model theoretically the mechanical peeling of a nanoribbon of graphene as realized by the tip-forced lifting of one of its extremes off a flat crystal surface. The evolution of shape, energy, local curvature and body advancement are ideally expected to follow a succession of regimes: (A) initial prying, (B) peeling with stretching but without sliding (stripping), (C) peeling with sliding, (D) liftoff. In the case where in addition the substrate surface corrugation is small or negligible, then (B) disappears, and we find that the (A)-(C) transition becomes universal, analytical and sharp, determined by the interplay between bending rigidity and adsorption energy. This general two-stage peeling transition is identified as a sharp crossover in published data of graphene nanoribbons pulled off an atomic-scale Au(111) substrate.


title: Resolution of a conjecture on majority dynamics: rapid stabilisation in dense random graphs

abstract: We study majority dynamics on the binomial random graph $G(n,p)$ with $p = d/n$ and $d > \lambda n^{1/2}$, for some large $\lambda>0$. In this process, each vertex has a state in $\{-1,+1 \}$ and at each round every vertex adopts the state of the majority of its neighbours, retaining its state in the case of a tie. We show that with high probability the process reaches unanimity in at most four rounds. This confirms a conjecture of Benjamini, Chan, O' Donnel, Tamuz and Tan.


title: Policy Poisoning in Batch Reinforcement Learning and Control

abstract: We study a security threat to batch reinforcement learning and control where the attacker aims to poison the learned policy. The victim is a reinforcement learner / controller which first estimates the dynamics and the rewards from a batch data set, and then solves for the optimal policy with respect to the estimates. The attacker can modify the data set slightly before learning happens, and wants to force the learner into learning a target policy chosen by the attacker. We present a unified framework for solving batch policy poisoning attacks, and instantiate the attack on two standard victims: tabular certainty equivalence learner in reinforcement learning and linear quadratic regulator in control. We show that both instantiation result in a convex optimization problem on which global optimality is guaranteed, and provide analysis on attack feasibility and attack cost. Experiments show the effectiveness of policy poisoning attacks.


title: Cheeger-Gromoll Splitting Theorem for groups

abstract: We study a notion of curvature for finitely generated groups which serves as a role of Ricci curvature for Riemannian manifolds. We prove an analog of Cheeger-Gromoll splitting theorem. As a consequence, we give a geometric characterization of virtually abelian groups. We also explore the relation between this notion of curvature and the growth of groups.


title: Transition between extinction and blow-up in a generalized Fisher-KPP model

abstract: Stationary solutions of the Fisher-KPP equation with general nonlinear diffusion and arbitrary reactional kinetic orders terms are characterized. Such stationary (separatrix-like) solutions disjoint the blow-up solutions from those showing extinction. In addition a criterion for general parameter values is presented, which allows determining the blow-up or vanishing character of the solutions.

categories: math.AP math-ph math.MP nlin.PS physics.bio-ph physics.chem-ph


title: Global boundedness of multilinear Fourier integral operators

abstract: We establish global regularity of multilinear Fourier integral operators that are associated to nonlinear wave equations on product of $L^p$ spaces by proving endpoint boundedness on suitable products spaces containing combinations of the local Hardy space, the local BMO and the $L^2$ spaces.


title: A minimal pair in the generic degrees

abstract: We show that there is a minimal pair in the nonuniform generic degrees, and hence also in the uniform generic degrees. This fact contrasts with Igusa's result that there are no minimal pairs for relative generic computability, and answers a basic structural question mentioned in several papers in the area.


title: A class of optimization problems motivated by rank estimators in robust regression

abstract: A rank estimator in robust regression is a minimizer of a function which depends (in addition to other factors) on the ordering of residuals but not on their values. Here we focus on the optimization aspects of rank estimators. We distinguish two classes of functions: the class with a continuous and convex objective function (CCC), which covers the class of rank estimators known from statistics, and also another class (GEN), which is far more general. We propose efficient algorithms for both classes. For GEN we propose an enumerative algorithm that works in polynomial time as long as the number of regressors is O(1). The proposed algorithm utilizes the special structure of arrangements of hyperplanes that occur in our problem and is superior to other known algorithms in this area. For the continuous and convex case, we propose an unconditionally polynomial algorithm finding the exact minimizer, unlike the heuristic or approximate methods implemented in statistical packages.


title: Generative Image Translation for Data Augmentation in Colorectal Histopathology Images

abstract: We present an image translation approach to generate augmented data for mitigating data imbalances in a dataset of histopathology images of colorectal polyps, adenomatous tumors that can lead to colorectal cancer if left untreated. By applying cycle-consistent generative adversarial networks (CycleGANs) to a source domain of normal colonic mucosa images, we generate synthetic colorectal polyp images that belong to diagnostically less common polyp classes. Generated images maintain the general structure of their source image but exhibit adenomatous features that can be enhanced with our proposed filtration module, called Path-Rank-Filter. We evaluate the quality of generated images through Turing tests with four gastrointestinal pathologists, finding that at least two of the four pathologists could not identify generated images at a statistically significant level. Finally, we demonstrate that using CycleGAN-generated images to augment training data improves the AUC of a convolutional neural network for detecting sessile serrated adenomas by over 10%, suggesting that our approach might warrant further research for other histopathology image classification tasks.


title: Dynamics of many-body photon bound states in chiral waveguide QED

abstract: We theoretically study the few- and many-body dynamics of photons in chiral waveguides. In particular, we examine pulse propagation through a system of $N$ two-level systems chirally coupled to a waveguide. We show that the system supports correlated multi-photon bound states, which have a well-defined photon number $n$ and propagate through the system with a group delay scaling as $1/n^2$. This has the interesting consequence that, during propagation, an incident coherent state pulse breaks up into different bound state components that can become spatially separated at the output in a sufficiently long system. For sufficiently many photons and sufficiently short systems, we show that linear combinations of $n$-body bound states recover the well-known phenomenon of mean-field solitons in self-induced transparency. For longer systems, however, the solitons break apart through quantum correlated dynamics. Our work thus covers the entire spectrum from few-photon quantum propagation, to genuine quantum many-body (atom and photon) phenomena, and ultimately the quantum-to-classical transition. Finally, we demonstrate that the bound states can undergo elastic scattering with additional photons. Together, our results demonstrate that photon bound states are truly distinct physical objects emerging from the most elementary light-matter interaction between photons and two-level emitters. Our work opens the door to studying quantum many-body physics and soliton physics with photons in chiral waveguide QED.


title: Trajectory construction of Dirac evolution

abstract: We extend our programme of representing the quantum state through exact stand-alone trajectory models to the Dirac equation. We show that the free Dirac equation in the angular coordinate representation is a continuity equation for which the real and imaginary parts of the wavefunction, angular versions of Majorana spinors, define conserved densities. We hence deduce an exact formula for the propagation of the Dirac spinor derived from the self-contained first-order dynamics of two sets of trajectories in 3-space together with a mass-dependent evolution operator. The Lorentz covariance of the trajectory equations is established by invoking the relativity of the trajectory label. We show how these results extend to the inclusion of external potentials. We further show that the angular version of the Dirac equation implies continuity equations for currents with non-negative densities, for which the Dirac current defines the mean flow. This provides an alternative trajectory construction of free evolution. Finally, we examine the polar representation of the Dirac equation, which also implies a non-negative conserved density but does not map into a stand-alone trajectory theory. It reveals how the quantum potential is tacit in the Dirac equation.


title: SPARKS II.: Complex organic molecules in accretion shocks around a hot core precursor

abstract: Classical hot cores are rich in molecular emission, and they show a high abundance of complex organic molecules (COMs). The emergence of molecular complexity is poorly constrained in the early evolution of hot cores. Using the Atacama Large Millimeter Array we put observational constraints on the physical location of COMs in a high-mass protostellar envelope associated with the G328.2551-0.5321 clump. The protostar is single down to ~400au scales and we resolve the emission region of COMs. Using thermodynamic equilibrium modelling of the available 7.5 GHz bandwidth around ~345 GHz, we detect emission from 10 COMs, and identify a line of deuterated water (HDO). The most extended emission originates from methanol, methyl formate and formamide. Together with HDO, these molecules are found to be associated with both the accretion shocks and the inner envelope, for which we estimate a moderate temperature of $T_{\rm kin}\sim$110 K. Our findings reveal a significant difference in the distribution of COMs. O-bearing COMs, such as ethanol, acetone, and ethylene glycol are almost exclusively found and show a higher abundance towards the accretion shocks with $T_{\rm kin}\sim$180 K. Whereas N-bearing COMs with a CN group, such as vinyl and ethyl cyanide peak on the central position, thus the protostar and the accretion disk. This is the first observational evidence for a large column density of COMs seen towards accretion shocks at the centrifugal barrier at the inner envelope. Since the molecular composition is dominated by that of the accretion shocks and the radiatively heated hot inner region is very compact, we propose this source to be a precursor to a classical, radiatively heated hot core.


title: Non-renormalization and operator mixing via on-shell methods

abstract: Using on-shell methods, we present a new perturbative non-renormalization theorem for operator mixing in massless four-dimensional quantum field theories. By examining how unitarity cuts of form factors encode anomalous dimensions we show that longer operators are often restricted from renormalizing shorter operators at the first order where there exist Feynman diagrams. The theorem applies quite generally and depends only on the field content of the operators involved. We apply our theorem to operators of dimension five through seven in the Standard Model Effective Field Theory, including examples of nontrivial zeros in the anomalous-dimension matrix at one through four loops. The zeros at two and higher loops go beyond those previously explained using helicity selection rules. We also include explicit sample calculations at two loops.


title: A Large-Deviation Upperbound on Directed Last-Passage Percolation Growth Rate Based on Entropy of Direction Vector

abstract: This short note provides a large-deviation-based upper bound on the growth rate of directed last passage percolation (LPP) using the entropy of the normalized direction vector.


title: Analytical Solution of (2+1) Dimensional Dirac Equation in Time-Dependent Noncommutative Phase-Space

abstract: In this article, we studied the system of (2+1) dimensional Dirac equation in time-dependent noncommutative phase-space. Exactly, we investigated the analytical solution of the corresponding system by the Lewis-Riesenfeld invariant method based on the construction of the Lewis-Riesenfeld invariant. Knowing that we obtained the time-dependent Dirac Hamiltonian of the problem in question from a time-dependent Bopp-Shift translation, then it used to set the Lewis-Riesenfeld invariant operators. Thereafter, the obtained results used to express the eigenfunctions that lead to determining the general solution of the system.


title: On a strengthening of J\'onssonness for $\aleph_\omega$

abstract: We discuss a system of strengthenings of "$\aleph_\omega$ is J\'onsson" indexed by real numbers, and identify a strongest one. We give a proof of a theorem of Silver and show that there is a barrier to weakening its hypothesis.


title: A multi-label, dual-output deep neural network for automated bug triaging

abstract: Bug tracking enables the monitoring and resolution of issues and bugs within organizations. Bug triaging, or assigning bugs to the owner(s) who will resolve them, is a critical component of this process because there are many incorrect assignments that waste developer time and reduce bug resolution throughput. In this work, we explore the use of a novel two-output deep neural network architecture (Dual DNN) for triaging a bug to both an individual team and developer, simultaneously. Dual DNN leverages this simultaneous prediction by exploiting its own guess of the team classes to aid in developer assignment. A multi-label classification approach is used for each of the two outputs to learn from all interim owners, not just the last one who closed the bug. We make use of a heuristic combination of the interim owners (owner-importance-weighted labeling) which is converted into a probability mass function (pmf). We employ a two-stage learning scheme, whereby the team portion of the model is trained first and then held static to train the team--developer and bug--developer relationships. The scheme employed to encode the team--developer relationships is based on an organizational chart (org chart), which renders the model robust to organizational changes as it can adapt to role changes within an organization. There is an observed average lift (with respect to both team and developer assignment) of 13%-points in 11-fold incremental-learning cross-validation (IL-CV) accuracy for Dual DNN utilizing owner-weighted labels compared with the traditional multi-class classification approach. Furthermore, Dual DNN with owner-weighted labels achieves average 11-fold IL-CV accuracies of 76% (team assignment) and 55% (developer assignment), outperforming reference models by 14%- and 25%-points, respectively, on a proprietary dataset with 236,865 entries.


title: Modelling Resistive and Phase Change Memory with Passive Selector Arrays -- A Matlab Tool

abstract: Memristor devices are crucial for developing neuromorphic computers and next-generation memory technologies. In this work, we provide a comprehensive modelling tool for simulating static DC reading operations of memristor crossbar arrays that use passive selectors with matrix algebra in MATLAB. The software tool was parallel coded and optimized to run with personal computers and distributed computer clusters with minimized CPU and memory consumption. Using the tool, we demonstrate the effect of changing the line resistance, array size, voltage selection scheme, selector diode's ideality factor, reverse saturation current, temperature and sense resistance on the electrical behavior and expected sense margin of one-diode-one-resistor crossbar arrays. We then investigate the effect of single and dual side array biasing and grounding on the dissipated current throughout the array cells. The tool we offer to the memristor community and the studies we present enables the design of larger and more practical memristor arrays for application in data storage and neuromorphic computing.


title: A horseshoe with a discontinuous entropy spectrum

abstract: We study the regularity of the entropy spectrum of the Lyapunov exponents for hyperbolic maps on surfaces. It is well-known that the entropy spectrum is a concave upper semi-continuous function which is analytic on the interior of the set Lyapunov exponents. In this paper we construct a family of horseshoes with a discontinuous entropy spectrum at the boundary of the set of Lyapunov exponents.


title: Virtual Optical Pulling Force

abstract: The tremendous progress in light scattering engineering made it feasible to develop optical tweezers allowing capture, hold, and controllable displacement of submicronsize particles and biological structures. However, the momentum conservation law imposes a fundamental restriction on the optical pressure to be repulsive in paraxial fields. Although different approaches to get around this restriction have been proposed, they are rather sophisticated and rely on either wavefront engineering or utilize active media. Herein, we revisit the issue of optical forces by their analytic continuation to the complex frequency plane and considering their behavior in transient. We show that the exponential excitation at the complex frequency offers an intriguing ability to achieve a pulling force for a passive resonator of any shape and composition even in the paraxial approximation, the remarkable effect which is not reduced to the Fourier transform. The approach is linked to the virtual gain effect when an appropriate transient decay of the excitation signal makes it weaker than the outgoing signal that carries away greater energy and momentum flux density. The approach is implemented for the Fabry-Perot cavity and a high refractive index dielectric nanoparticle, a fruitful platform for intracellular spectroscopy and lab-on-a-chip technologies where the proposed technique may found unprecedented capabilities.


title: RGB-Infrared Cross-Modality Person Re-Identification via Joint Pixel and Feature Alignment

abstract: RGB-Infrared (IR) person re-identification is an important and challenging task due to large cross-modality variations between RGB and IR images. Most conventional approaches aim to bridge the cross-modality gap with feature alignment by feature representation learning. Different from existing methods, in this paper, we propose a novel and end-to-end Alignment Generative Adversarial Network (AlignGAN) for the RGB-IR RE-ID task. The proposed model enjoys several merits. First, it can exploit pixel alignment and feature alignment jointly. To the best of our knowledge, this is the first work to model the two alignment strategies jointly for the RGB-IR RE-ID problem. Second, the proposed model consists of a pixel generator, a feature generator, and a joint discriminator. By playing a min-max game among the three components, our model is able to not only alleviate the cross-modality and intra-modality variations but also learn identity-consistent features. Extensive experimental results on two standard benchmarks demonstrate that the proposed model performs favorably against state-of-the-art methods. Especially, on SYSU-MM01 dataset, our model can achieve an absolute gain of 15.4% and 12.9% in terms of Rank-1 and mAP.


title: Empirical and Constrained Empirical Bayes Variance Estimation Under A One Unit Per Stratum Sample Design

abstract: A single primary sampling unit (PSU) per stratum design is a popular design for estimating the parameter of interest. Although, the point estimator of the design is unbiased and efficient, an unbiased variance estimator does not exist. A common practice to solve this is to collapse or combine the two adjacent strata, but the attained estimator of variance is not design-unbiased, and the bias increases as the population means of the collapsed strata become more variant. Therefore, the one PSU per stratum design with collapsed stratum variance estimator might not be a good choice, and some statisticians prefer a design in which two PSUs per stratum are selected. In this paper, we first compare a one PSU per stratum design to a two PSUs per stratum design. Then, we propose an empirical Bayes estimator for the variance of one PSU per stratum design, where it over-shrinks towards the prior mean. To protect against this, we investigate the potential of a constrained empirical Bayes estimator. Through a simulation study, we show that the empirical Bayes and constrained empirical Bayes estimators outperform the classical collapsed one in terms of empirical relative mean squared error.


title: Assessing the potential of Ge/SiGe quantum dots as hosts for singlet-triplet qubits

abstract: We study double quantum dots in a Ge/SiGe heterostructure and test their maturity towards singlet-triplet ($S-T_0$) qubits. We demonstrate a large range of tunability, from two single quantum dots to a double quantum dot. We measure Pauli spin blockade and study the anisotropy of the $g$-factor. We use an adjacent quantum dot for sensing charge transitions in the double quantum dot at interest. In conclusion, Ge/SiGe possesses all ingredients necessary for building a singlet-triplet qubit.


title: Statistical Topology of Bond Networks with Applications to Silica

abstract: Whereas knowledge of a crystalline material's unit cell is fundamental to understanding the material's properties and behavior, there are not obvious analogues to unit cells for disordered materials despite the frequent existence of considerable medium-range order. This article views a material's structure as a collection of local atomic environments that are sampled from some underlying probability distribution of such environments, with the advantage of offering a unified description of both ordered and disordered materials. Crystalline materials can then be regarded as special cases where the underlying probability distribution is highly concentrated around the traditional unit cell. Four descriptors of local atomic environments suitable for disordered bond networks are proposed and applied to molecular dynamics simulations of silica glasses. Each of them reliably distinguishes the structure of glasses produced at different cooling rates, with the $H_1$ barcode and coordination profile providing the best separation.


title: Regularized Sparse Gaussian Processes

abstract: Gaussian processes are a flexible Bayesian nonparametric modelling approach that has been widely applied but poses computational challenges. To address the poor scaling of exact inference methods, approximation methods based on sparse Gaussian processes (SGP) are attractive. An issue faced by SGP, especially in latent variable models, is the inefficient learning of the inducing inputs, which leads to poor model prediction. We propose a regularization approach by balancing the reconstruction performance of data and the approximation performance of the model itself. This regularization improves both inference and prediction performance. We extend this regularization approach into latent variable models with SGPs and show that performing variational inference (VI) on those models is equivalent to performing VI on a related empirical Bayes model.


title: Open problems in the Kontsevich graph construction of Poisson bracket symmetries

abstract: Poisson brackets admit infinitesimal symmetries which are encoded using oriented graphs; this construction is due to Kontsevich (1996). We formulate several open problems about combinatorial and topological properties of the graphs involved, about integrability and analytic properties of such symmetry flows (in particular, for known classes of Poisson brackets), and about cohomological, differential geometric, and quantum aspects of the theory. Keywords: Poisson bracket, Jacobi identity, infinitesimal symmetry, affine manifold, Poisson cohomology, graph complex.

categories: math-ph math.CO math.MP math.QA math.SG nlin.SI


title: A Pooled Quantile Estimator for Parallel Simulations

abstract: Quantile is an important risk measure quantifying the stochastic system random behaviors. This paper studies a pooled quantile estimator, which is the sample quantile of detailed simulation outputs after directly pooling independent sample paths together. We derive the asymptotic representation of the pooled quantile estimator and further prove its normality. By comparing with the classical quantile estimator used in stochastic simulation, both theoretical and empirical studies demonstrate the advantages of the proposal under the context of parallel simulation.


title: Mathematical Modelling and Analysis of Transmission Dynamics of Lassa fever

abstract: In this work, a periodically-forced seasonal non-autonomous system of a non-linear ordinary differential equation is developed that captures the dynamics of Lassa fever transmission and seasonal variation in the birth of mastomys rodents where time was measured in days to capture seasonality. It was shown that the model is epidemiologically meaningful and mathematically well-posed by using the results from the qualitative properties of the solution of the model. It was established that in order to eliminate Lassa fever disease, treatments with Ribavirin must be provided early to reduce mortality and other preventive measures like an educational campaign, community hygiene, Isolation of infected humans, and culling/destruction of rodents must be applied to also reduce the morbidity of the disease. Finally, the obtained results gave a primer framework for planning and designing cost-effective strategies for good interventions in eliminating Lassa fever.


title: Hierarchical Hidden Markov Jump Processes for Cancer Screening Modeling

abstract: Hidden Markov jump processes are an attractive approach for modeling clinical disease progression data because they are explainable and capable of handling both irregularly sampled and noisy data. Most applications in this context consider time-homogeneous models due to their relative computational simplicity. However, the time homogeneous assumption is too strong to accurately model the natural history of many diseases. Moreover, the population at risk is not homogeneous either, since disease exposure and susceptibility can vary considerably. In this paper, we propose a piece-wise stationary transition matrix to explain the heterogeneity in time. We propose a hierarchical structure for the heterogeneity in population, where prior information is considered to deal with unbalanced data. Moreover, an efficient, scalable EM algorithm is proposed for inference. We demonstrate the feasibility and superiority of our model on a cervical cancer screening dataset from the Cancer Registry of Norway. Experiments show that our model outperforms state-of-the-art recurrent neural network models in terms of prediction accuracy and significantly outperforms a standard hidden Markov jump process in generating Kaplan-Meier estimators.


title: Macdonald Polynomials and level two Demazure modules for affine $\mathfrak{sl}_{n+1}$

abstract: We define a family of symmetric polynomials $G_{\nu,\lambda}(z_1,\cdots, z_{n+1},q)$ indexed by a pair of dominant integral weights. The polynomial $G_{\nu,0}(z,q)$ is the specialized Macdonald polynomial and we prove that $G_{0,\lambda}(z,q)$ is the graded character of a level two Demazure module associated to the affine Lie algebra $\widehat{\mathfrak{sl}}_{n+1}$. Under suitable conditions on $(\nu,\lambda)$ (which includes the case when $\nu=0$ or $\lambda=0$) we prove that $G_{\nu,\lambda}(z,q)$ is Schur positive and give explicit formulae for them in terms of Macdonald polynomials.


title: Bimeromorphic automorphisms groups of certain conic bundles

abstract: We study the groups of biholomorphic and bimeromorphic automorphisms of conic bundles over certain compact complex manifolds of algebraic dimension zero.


title: Topological data analysis approaches to uncovering the timing of ring structure onset in filamentous networks

abstract: Improvements in experimental and computational technologies have led to significant increases in data available for analysis. Topological data analysis (TDA) is an emerging area of mathematical research that can identify structures in these data sets. Here we develop a TDA method to detect physical structures in a cell that persist over time. In most cells, protein filaments (actin) interact with motor proteins (myosins) and organize into polymer networks and higher-order structures. An example of these structures are ring channels that maintain constant diameters over time and play key roles in processes such as cell division, development, and wound healing. The interactions of actin with myosin can be challenging to investigate experimentally in living systems, given limitations in filament visualization \textit{in vivo}. We therefore use complex agent-based models that simulate mechanical and chemical interactions of polymer proteins in cells. To understand how filaments organize into structures, we propose a TDA method that assesses effective ring generation in data consisting of simulated actin filament positions through time. We analyze the topological structure of point clouds sampled along these actin filaments and propose an algorithm for connecting significant topological features in time. We introduce visualization tools that allow the detection of dynamic ring structure formation. This method provides a rigorous way to investigate how specific interactions and parameters may impact the timing of filamentous network organization.


title: Nonstationary Multivariate Gaussian Processes for Electronic Health Records

abstract: We propose multivariate nonstationary Gaussian processes for jointly modeling multiple clinical variables, where the key parameters, length-scales, standard deviations and the correlations between the observed output, are all time dependent. We perform posterior inference via Hamiltonian Monte Carlo (HMC). We also provide methods for obtaining computationally efficient gradient-based maximum a posteriori (MAP) estimates. We validate our model on synthetic data as well as on electronic health records (EHR) data from Kaiser Permanente (KP). We show that the proposed model provides better predictive performance over a stationary model as well as uncovers interesting latent correlation processes across vitals which are potentially predictive of patient risk.


title: Implicit competitive regularization in GANs

abstract: To improve the stability of GAN training we need to understand why they can produce realistic samples. Presently, this is attributed to properties of the divergence obtained under an optimal discriminator. This argument has a fundamental flaw: If we do not impose regularity of the discriminator, it can exploit visually imperceptible errors of the generator to always achieve the maximal generator loss. In practice, gradient penalties are used to regularize the discriminator. However, this needs a metric on the space of images that captures visual similarity. Such a metric is not known, which explains the limited success of gradient penalties in stabilizing GANs. We argue that the performance of GANs is instead due to the implicit competitive regularization (ICR) arising from the simultaneous optimization of generator and discriminator. ICR promotes solutions that look real to the discriminator and thus leverages its inductive biases to generate realistic images. We show that opponent-aware modelling of generator and discriminator, as present in competitive gradient descent (CGD), can significantly strengthen ICR and thus stabilize GAN training without explicit regularization. In our experiments, we use an existing implementation of WGAN-GP and show that by training it with CGD we can improve the inception score (IS) on CIFAR10 for a wide range of scenarios, without any hyperparameter tuning. The highest IS is obtained by combining CGD with the WGAN-loss, without any explicit regularization.


title: Precise radial velocities of giant stars XIII. A second Jupiter orbiting in 4:3 resonance in the 7 CMa system

abstract: We report the discovery of a second planet orbiting the K giant star 7 CMa based on 166 high-precision radial velocities obtained with Lick, HARPS, UCLES and SONG. The periodogram analysis reveals two periodic signals of approximately 745 and 980 d, associated to planetary companions. A double-Keplerian orbital fit of the data reveals two Jupiter-like planets with minimum masses $m_b\sin i \sim 1.9 \,\mathrm{M_{J}}$ and $m_c\sin i \sim 0.9 \,\mathrm{M_{J}}$, orbiting at semi-major axes of $a_b \sim 1.75\,\mathrm{au}$ and $a_c \sim 2.15\,\mathrm{au}$, respectively. Given the small orbital separation and the large minimum masses of the planets close encounters may occur within the time baseline of the observations, thus, a more accurate N-body dynamical modeling of the available data is performed. The dynamical best-fit solution leads to collision of the planets and we explore the long-term stable configuration of the system in a Bayesian framework, confirming that 13% of the posterior samples are stable for at least 10 Myr. The result from the stability analysis indicates that the two-planets are trapped in a low-eccentricity 4:3 mean-motion resonance. This is only the third discovered system to be inside a 4:3 resonance, making it very valuable for planet formation and orbital evolution models.


title: On the Long-Range Dependence of Mixed Fractional Poisson Process

abstract: In this paper, we show that the mixed fractional Poisson process (MFPP) exhibits the long-range dependence (LRD) property. It is proved by establishing an asymptotic result for the covariance of inverse mixed stable subordinator. Also, it is shown that the increments of the MFPP, namely, the mixed fractional Poissonian noise (MFPN) has the short-range dependence (SRD) property.


title: Stallings automata for free-times-abelian groups: intersections and index

abstract: We extend the classical Stallings theory (describing subgroups of free groups as automata) to direct products of free and abelian groups: after introducing enriched automata (i.e., automata with extra abelian labels), we obtain an explicit bijection between subgroups and a certain type of such enriched automata, which - as it happens in the free group - is computable in the finitely generated case. This approach provides a neat geometric description of (even non finitely generated) intersections of finitely generated subgroups within this non-Howson family. In particular, we give a geometric solution to the subgroup intersection problem and the finite index problem, providing recursive bases and transversals respectively.


title: Ultra-luminous quasars at redshift $z>4.5$ from SkyMapper

abstract: The most luminous quasars at high redshift harbour the fastest-growing and most massive black holes in the early Universe. They are exceedingly rare and hard to find. Here, we present our search for the most luminous quasars in the redshift range from $z=4.5$ to $5$ using data from SkyMapper, Gaia and WISE. We use colours to select likely high-redshift quasars and reduce the stellar contamination of the candidate set with parallax and proper motion data. In $\sim$12,500~deg$^2$ of Southern sky, we find 92 candidates brighter than $R_p=18.2$. Spectroscopic follow-up has revealed 21 quasars at $z\ge 4$ (16 of which are within $z=[4.5,5]$), as well as several red quasars, BAL quasars and objects with unusual spectra, which we tentatively label OFeLoBALQSOs at redshifts of $z\approx 1$ to $2$. This work lifts the number of known bright $z\ge 4.5$ quasars in the Southern hemisphere from 10 to 26 and brings the total number of quasars known at $R_p<18.2$ and $z\ge 4.5$ to 42.


title: Improving the Sample and Communication Complexity for Decentralized Non-Convex Optimization: A Joint Gradient Estimation and Tracking Approach

abstract: Many modern large-scale machine learning problems benefit from decentralized and stochastic optimization. Recent works have shown that utilizing both decentralized computing and local stochastic gradient estimates can outperform state-of-the-art centralized algorithms, in applications involving highly non-convex problems, such as training deep neural networks. In this work, we propose a decentralized stochastic algorithm to deal with certain smooth non-convex problems where there are $m$ nodes in the system, and each node has a large number of samples (denoted as $n$). Differently from the majority of the existing decentralized learning algorithms for either stochastic or finite-sum problems, our focus is given to both reducing the total communication rounds among the nodes, while accessing the minimum number of local data samples. In particular, we propose an algorithm named D-GET (decentralized gradient estimation and tracking), which jointly performs decentralized gradient estimation (which estimates the local gradient using a subset of local samples) and gradient tracking (which tracks the global full gradient using local estimates). We show that, to achieve certain $\epsilon$ stationary solution of the deterministic finite sum problem, the proposed algorithm achieves an $\mathcal{O}(mn^{1/2}\epsilon^{-1})$ sample complexity and an $\mathcal{O}(\epsilon^{-1})$ communication complexity. These bounds significantly improve upon the best existing bounds of $\mathcal{O}(mn\epsilon^{-1})$ and $\mathcal{O}(\epsilon^{-1})$, respectively. Similarly, for online problems, the proposed method achieves an $\mathcal{O}(m \epsilon^{-3/2})$ sample complexity and an $\mathcal{O}(\epsilon^{-1})$ communication complexity, while the best existing bounds are $\mathcal{O}(m\epsilon^{-2})$ and $\mathcal{O}(\epsilon^{-2})$, respectively.


title: Deep Kernels with Probabilistic Embeddings for Small-Data Learning

abstract: Gaussian Processes (GPs) are known to provide accurate predictions and uncertainty estimates even with small amounts of labeled data by capturing similarity between data points through their kernel function. However traditional GP kernels are not very effective at capturing similarity between high dimensional data points. Neural networks can be used to learn good representations that encode intricate structures in high dimensional data, and can be used as inputs to the GP kernel. However the huge data requirement of neural networks makes this approach ineffective in small data settings. To solves the conflicting problems of representation learning and data efficiency, we propose to learn deep kernels on probabilistic embeddings by using a probabilistic neural network. Our approach maps high-dimensional data to a probability distribution in a low dimensional subspace and then computes a kernel between these distributions to capture similarity. To enable end-to-end learning, we derive a functional gradient descent procedure for training the model. Experiments on a variety of datasets show that our approach outperforms the state-of-the-art in GP kernel learning in both supervised and semi-supervised settings. We also extend our approach to other small-data paradigms such as few-shot classification where it outperforms previous approaches on mini-Imagenet and CUB datasets.


title: Accelerated Structured Alternating Projections for Robust Spectrally Sparse Signal Recovery

abstract: Consider a spectrally sparse signal $\boldsymbol{x}$ that consists of $r$ complex sinusoids with or without damping. We study the robust recovery problem for the spectrally sparse signal under the fully observed setting, which is about recovering $\boldsymbol{x}$ and a sparse corruption vector $\boldsymbol{s}$ from their sum $\boldsymbol{z}=\boldsymbol{x}+\boldsymbol{s}$. In this paper, we exploit the low-rank property of the Hankel matrix formed by $\boldsymbol{x}$, and formulate the problem as the robust recovery of a corrupted low-rank Hankel matrix. We develop a highly efficient non-convex algorithm, coined Accelerated Structured Alternating Projections (ASAP). The high computational efficiency and low space complexity of ASAP are achieved by fast computations involving structured matrices, and a subspace projection method for accelerated low-rank approximation. Theoretical recovery guarantee with a linear convergence rate has been established for ASAP, under some mild assumptions on $\boldsymbol{x}$ and $\boldsymbol{s}$. Empirical performance comparisons on both synthetic and real-world data confirm the advantages of ASAP, in terms of computational efficiency and robustness aspects.


title: Hot thermal universe endowed with massive dark vector fields and the Hubble tension

abstract: The value of Hubble constant inferred from Planck measurements of anisotropies in the cosmic microwave background is at $4.4\sigma$ tension with direct astronomical measurements at low redshifts. Very recently, it has been conjectured that this discrepancy may be reconciled if a small fraction of the dark matter is described by three mutually orthogonal vector fields of the same mass. We study the thermal description of this model and use the observationally-inferred primordial fractions of baryonic mass in $^4{\rm He}$ to constrain its phase space. We show that while the sterile vector fields may help to alleviate a little bit the existing tension in the measurements of the Hubble parameter, they cannot eliminate the discrepancy between low- and high-redshift observations.


title: Machine Learning for Prediction with Missing Dynamics

abstract: This article presents a general framework for recovering missing dynamical systems using available data and machine learning techniques. The proposed framework reformulates the prediction problem as a supervised learning problem to approximate a map that takes the memories of the resolved and identifiable unresolved variables to the missing components in the resolved dynamics. We demonstrate the effectiveness of the proposed framework with a theoretical guarantee of a path-wise convergence of the resolved variables up to finite time and numerical tests on prototypical models in various scientific domains. These include the 57-mode barotropic stress models with multiscale interactions that mimic the blocked and unblocked patterns observed in the atmosphere, the nonlinear Schr\"{o}dinger equation which found many applications in physics such as optics and Bose-Einstein-Condense, the Kuramoto-Sivashinsky equation which spatiotemporal chaotic pattern formation models trapped ion mode in plasma and phase dynamics in reaction-diffusion systems. While many machine learning techniques can be used to validate the proposed framework, we found that recurrent neural networks outperform kernel regression methods in terms of recovering the trajectory of the resolved components and the equilibrium one-point and two-point statistics. This superb performance suggests that recurrent neural networks are an effective tool for recovering the missing dynamics that involves approximation of high-dimensional functions.


title: Constrained Non-Affine Alignment of Embeddings

abstract: Embeddings are one of the fundamental building blocks for data analysis tasks. Embeddings are already essential tools for large language models and image analysis, and their use is being extended to many other research domains. The generation of these distributed representations is often a data- and computation-expensive process; yet the holistic analysis and adjustment of them after they have been created is still a developing area. In this paper, we first propose a very general quantitatively measure for the presence of features in the embedding data based on if it can be learned. We then devise a method to remove or alleviate undesired features in the embedding while retaining the essential structure of the data. We use a Domain Adversarial Network (DAN) to generate a non-affine transformation, but we add constraints to ensure the essential structure of the embedding is preserved. Our empirical results demonstrate that the proposed algorithm significantly outperforms the state-of-art unsupervised algorithm on several data sets, including novel applications from the industry.


title: Global-Local Metamodel Assisted Two-Stage Optimization via Simulation

abstract: To integrate strategic, tactical and operational decisions, the two-stage optimization has been widely used to guide dynamic decision making. In this paper, we study the two-stage stochastic programming for complex systems with unknown response estimated by simulation. We introduce the global-local metamodel assisted two-stage optimization via simulation that can efficiently employ the simulation resource to iteratively solve for the optimal first- and second-stage decisions. Specifically, at each visited first-stage decision, we develop a local metamodel to simultaneously solve a set of scenario-based second-stage optimization problems, which also allows us to estimate the optimality gap. Then, we construct a global metamodel accounting for the errors induced by: (1) using a finite number of scenarios to approximate the expected future cost occurring in the planning horizon, (2) second-stage optimality gap, and (3) finite visited first-stage decisions. Assisted by the global-local metamodel, we propose a new simulation optimization approach that can efficiently and iteratively search for the optimal first- and second-stage decisions. Our framework can guarantee the convergence of optimal solution for the discrete two-stage optimization with unknown objective, and the empirical study indicates that it achieves substantial efficiency and accuracy.


title: Measurement of $\mathcal{R}(D)$ and $\mathcal{R}(D^*)$ with a semileptonic tagging method

abstract: The experimental results on the ratios of branching fractions $\mathcal{R}(D) = {\cal B}(\bar{B} \to D \tau^- \bar{\nu}_{\tau})/{\cal B}(\bar{B} \to D \ell^- \bar{\nu}_{\ell})$ and $\mathcal{R}(D^*) = {\cal B}(\bar{B} \to D^* \tau^- \bar{\nu}_{\tau})/{\cal B}(\bar{B} \to D^* \ell^- \bar{\nu}_{\ell})$, where $\ell$ denotes an electron or a muon, show a long-standing discrepancy with the Standard Model predictions, and might hint to a violation of lepton flavor universality. We report a new simultaneous measurement of $\mathcal{R}(D)$ and $\mathcal{R}(D^*)$, based on a data sample containing $772 \times 10^6$ $B\bar{B}$ events recorded at the $\Upsilon(4S)$ resonance with the Belle detector at the KEKB $e^+ e^-$ collider. In this analysis the tag-side $B$ meson is reconstructed in a semileptonic decay mode and the signal-side $\tau$ is reconstructed in a purely leptonic decay. The measured values are $\mathcal{R}(D)= 0.307 \pm 0.037 \pm 0.016$ and $\mathcal{R}(D^*) = 0.283 \pm 0.018 \pm 0.014$, where the first uncertainties are statistical and the second are systematic. These results are in agreement with the Standard Model predictions within $0.2$, $1.1$ and $0.8$ standard deviations for $\mathcal{R}(D)$, $\mathcal{R}(D^*)$ and their combination, respectively. This work constitutes the most precise measurements of $\mathcal{R}(D)$ and $\mathcal{R}(D^*)$ performed to date as well as the first result for $\mathcal{R}(D)$ based on a semileptonic tagging method.


title: Neural Program Synthesis By Self-Learning

abstract: Neural inductive program synthesis is a task generating instructions that can produce desired outputs from given inputs. In this paper, we focus on the generation of a chunk of assembly code that can be executed to match a state change inside the CPU and RAM. We develop a neural program synthesis algorithm, AutoAssemblet, learned via self-learning reinforcement learning that explores the large code space efficiently. Policy networks and value networks are learned to reduce the breadth and depth of the Monte Carlo Tree Search, resulting in better synthesis performance. We also propose an effective multi-entropy policy sampling technique to alleviate online update correlations. We apply AutoAssemblet to basic programming tasks and show significant higher success rates compared to several competing baselines.


title: Single-photon pulse induced giant response in N$>$100 qubit system

abstract: The temporal dynamics of large quantum systems perturbed weakly by a single excitation can give rise to unique phenomena at the quantum phase boundaries. Here, we develop a time-dependent model to study the temporal dynamics of a single photon interacting with a defect within a large system of interacting spin qubits (N$>$100). Our model predicts a new quantum resource, giant susceptibility, when the system of qubits is engineered to simulate a first-order quantum phase transition (QPT). We show that the absorption of a single photon pulse by an engineered defect in the large qubit system can nucleate a single shot quantum measurement through spin noise read-out. This concept of a single-shot detection event ("click") is different from parameter estimation which requires repeated measurements. The crucial step of amplifying the weak quantum signal occurs by coupling the defect to a system of interacting qubits biased close to a QPT point. The macroscopic change in long-range order during the QPT generates amplified magnetic noise, which can be read out by a classical device. Our work paves the way for studying the temporal dynamics of large quantum systems interacting with a single-photon pulse.


title: A Talk on 3d Matter Coupled to Chern-Simons Field, Spontaneous Breaking of Scale Invariance, and Fermion-Boson Mapping

abstract: The singlet sector of vector, large $N$, 3d field theory corresponds to Vasiliev higher spin theory on $AdS_4$. Will discuss three dimensional $U(N)$ symmetric field theory with fermion and boson matter coupled to a topological Chern-Simons field. In the presence of a marginal deformation will determine the conditions for the existence of a phase with spontaneous breaking of scale invariance. In this phase the ground state contains massive $U(N)$ quanta and a massless $U(N)$ singlet bound state goldstone boson:~the Dilaton. Will show that such a phase appears only in the presence of a marginal deformation. The massless Dilaton appears in the spectrum provided certain relations between coupling constants are satisfied. Will discuss the fermion-boson mapping and show that the conditions for spontaneous breaking of scale invariance in the boson and fermion theories are copies of each other.


title: Electric dipole moments of neutron and heavy quarks in the B-L symmetric MSSM

abstract: The searching for the electric dipole moments (EDMs) of neutron ($d_n$), $b$ quark ($d_b$) and $c$ quark ($d_c$) gives strict upper bounds on these quantities. And recently, new upper bounds on $d_b$, $d_c$ are obtained by the strict limit on $d_n$. The models of new physics (NP) with additional CP-violating (CPV) sources are constrained strictly by these EDMs. In this work, we focus on the CPV effects on these EDMs in the minimal supersymmetric extension (MSSM) of the SM with local $B-L$ gauge symmetry (B-LSSM). The contributions from one-loop and some two-loop diagrams to the quark EDM are given in general form, which can also be used in the calculation of quark EDM in other models of NP. Considering the constrains from updated experimental data, the numerical results show that the two-loop corrections can make important contributions to these EDMs. Compared with the MSSM, the effects of new CPV phases and new parameters in the B-LSSM on these EDMs are also explored.


title: Many-qutrit Mermin inequalities with three measurement settings

abstract: Mermin inequalities are derived for systems of three-state particles (qutrits) employing three local measurement settings. These establish perfect correlations which violate local realistic bounds more strongly than those previously reported with two bases. The quantum eigenvalue of the Mermin operator grows as the dimension of the Hilbert space, $3^N$, rather than $2^N$, as obtained with two measurement bases. The number of distinct GHZ contradictions also increases as $3^N$.


title: Network Modularity Controls the Speed of Information Diffusion

abstract: The rapid diffusion of information and the adoption of social behaviors are of critical importance in situations as diverse as collective actions, pandemic prevention, or advertising and marketing. Although the dynamics of large cascades have been extensively studied in various contexts, few have systematically examined the impact of network topology on the efficiency of information diffusion. Here, by employing the linear threshold model on networks with communities, we demonstrate that a prominent network feature---the modular structure---strongly affects the speed of information diffusion in complex contagion. Our simulations show that there always exists an optimal network modularity for the most efficient spreading process. Beyond this critical value, either a stronger or a weaker modular structure actually hinders the diffusion speed. These results are confirmed by an analytical approximation. We further demonstrate that the optimal modularity varies with both the seed size and the target cascade size, and is ultimately dependent on the network under investigation. We underscore the importance of our findings in applications from marketing to epidemiology, from neuroscience to engineering, where the understanding of the structural design of complex systems focuses on the efficiency of information propagation.


title: Chazy-Type Asymptotics and Hyperbolic Scattering for the $n$-Body Problem

abstract: We study solutions of the Newtonian $n$-body problem which tend to infinity hyperbolically, that is, all mutual distances tend to infinity with nonzero speed as $t \rightarrow +\infty$ or as $t \rightarrow -\infty$. In suitable coordinates, such solutions form the stable or unstable manifolds of normally hyperbolic equilibrium points in a boundary manifold "at infinity". We show that the flow near these manifolds can be analytically linearized and use this to give a new proof of Chazy's classical asymptotic formulas. We also address the scattering problem, namely, for solutions which are hyperbolic in both forward and backward time, how are the limiting equilibrium points related? After proving some basic theorems about this scattering relation, we use perturbations of our manifold at infinity to study scattering "near infinity", that is, when the bodies stay far apart and interact only weakly.


title: Self-supervised Label Augmentation via Input Transformations

abstract: Self-supervised learning, which learns by constructing artificial labels given only the input signals, has recently gained considerable attention for learning representations with unlabeled datasets, i.e., learning without any human-annotated supervision. In this paper, we show that such a technique can be used to significantly improve the model accuracy even under fully-labeled datasets. Our scheme trains the model to learn both original and self-supervised tasks, but is different from conventional multi-task learning frameworks that optimize the summation of their corresponding losses. Our main idea is to learn a single unified task with respect to the joint distribution of the original and self-supervised labels, i.e., we augment original labels via self-supervision of input transformation. This simple, yet effective approach allows to train models easier by relaxing a certain invariant constraint during learning the original and self-supervised tasks simultaneously. It also enables an aggregated inference which combines the predictions from different augmentations to improve the prediction accuracy. Furthermore, we propose a novel knowledge transfer technique, which we refer to as self-distillation, that has the effect of the aggregated inference in a single (faster) inference. We demonstrate the large accuracy improvement and wide applicability of our framework on various fully-supervised settings, e.g., the few-shot and imbalanced classification scenarios.


title: A Movable Valley Switch Driven by Berry Phase in Bilayer Graphene Resonators

abstract: Since its discovery, Berry phase has been demonstrated to play an important role in many quantum systems. In gapped Bernal bilayer graphene, the Berry phase can be continuously tuned from zero to 2pi, which offers a unique opportunity to explore the tunable Berry phase on the physical phenomena. Here, we report experimental observation of Berry phases-induced valley splitting and crossing in moveable bilayer graphene p-n junction resonators. In our experiment, the bilayer graphene resonators are generated by combining the electric field of scanning tunneling microscope tip with the gap of bilayer graphene. A perpendicular magnetic field changes the Berry phase of the confined bound states in the resonators from zero to 2pi continuously and leads to the Berry phase difference for the two inequivalent valleys in the bilayer graphene. As a consequence, we observe giant valley splitting and unusual valley crossing of the lowest bound states. Our results indicate that the bilayer graphene resonators can be used to manipulate the valley degree of freedom in valleytronics.


title: Effects of Depth, Width, and Initialization: A Convergence Analysis of Layer-wise Training for Deep Linear Neural Networks

abstract: Deep neural networks have been used in various machine learning applications and achieved tremendous empirical successes. However, training deep neural networks is a challenging task. Many alternatives have been proposed in place of end-to-end back-propagation. Layer-wise training is one of them, which trains a single layer at a time, rather than trains the whole layers simultaneously. In this paper, we study a layer-wise training using a block coordinate gradient descent (BCGD) for deep linear networks. We establish a general convergence analysis of BCGD and found the optimal learning rate, which results in the fastest decrease in the loss. More importantly, the optimal learning rate can directly be applied in practice, as it does not require any prior knowledge. Thus, tuning the learning rate is not needed at all. Also, we identify the effects of depth, width, and initialization in the training process. We show that when the orthogonal-like initialization is employed, the width of intermediate layers plays no role in gradient-based training, as long as the width is greater than or equal to both the input and output dimensions. We show that under some conditions, the deeper the network is, the faster the convergence is guaranteed. This implies that in an extreme case, the global optimum is achieved after updating each weight matrix only once. Besides, we found that the use of deep networks could drastically accelerate convergence when it is compared to those of a depth 1 network, even when the computational cost is considered. Numerical examples are provided to justify our theoretical findings and demonstrate the performance of layer-wise training by BCGD.


title: Multiband Nature of the Room-Temperature Superconductivity in Compressed LaH$_{10}$

abstract: Recently, the discovery of room-temperature superconductivity (SC) was experimentally realized in the fcc phase of LaH$_{10}$ under megabar pressures. This SC of compressed LaH$_{10}$ has been explained in terms of strong electron-phonon coupling (EPC), but the mechanism of how the large EPC constant and high superconducting transition temperature $T_{\rm c}$ are attained has not yet been clearly identified. Based on the density-functional theory and the Migdal-Eliashberg formalism, we reveal the presence of two nodeless, anisotropic superconducting gaps on the Fermi surface (FS). Here, the small gap is mostly associated with the hybridized states of H $s$ and La $f$ orbitals on the three outer FS sheets, while the large gap arises mainly from the hybridized state of neighboring H $s$ or $p$ orbitals on the one inner FS sheet. Further, we find that the EPC constant of compressed YH$_{10}$ with the same sodalite-like clathrate structure is enhanced due to the two additional FS sheets, leading to a higher $T_{\rm c}$ than LaH$_{10}$. It is thus demonstrated that the multiband pairing of hybridized electronic states is responsible for the large EPC constant and room-temperature SC in compressed hydrides LaH$_{10}$ and YH$_{10}$.


title: Actor Critic with Differentially Private Critic

abstract: Reinforcement learning algorithms are known to be sample inefficient, and often performance on one task can be substantially improved by leveraging information (e.g., via pre-training) on other related tasks. In this work, we propose a technique to achieve such knowledge transfer in cases where agent trajectories contain sensitive or private information, such as in the healthcare domain. Our approach leverages a differentially private policy evaluation algorithm to initialize an actor-critic model and improve the effectiveness of learning in downstream tasks. We empirically show this technique increases sample efficiency in resource-constrained control problems while preserving the privacy of trajectories collected in an upstream task.


title: Interpretable Deep Neural Networks for Facial Expression and Dimensional Emotion Recognition in-the-wild

abstract: In this project, we created a database with two types of annotations used in the emotion recognition domain : Action Units and Valence Arousal to try to achieve better results than with only one model. The originality of the approach is also based on the type of architecture used to perform the prediction of the emotions : a categorical Generative Adversarial Network. This kind of dual network can generate images based on the pictures from the new dataset thanks to its generative network and decide if an image is fake or real thanks to its discriminative network as well as help to predict the annotations for Action Units and Valence Arousal due to its categorical nature. GANs were trained on the Action Units model only, then the Valence Arousal model only and then on both the Action Units model and Valence Arousal model in order to test different parameters and understand their influence. The generative and discriminative aspects of the GANs have performed interesting results.


title: Manifold Embedded Knowledge Transfer for Brain-Computer Interfaces

abstract: Transfer learning makes use of data or knowledge in one problem to help solve a different, yet related, problem. It is particularly useful in brain-computer interfaces (BCIs), for coping with variations among different subjects and/or tasks. This paper considers offline unsupervised cross-subject electroencephalogram (EEG) classification, i.e., we have labeled EEG trials from one or more source subjects, but only unlabeled EEG trials from the target subject. We propose a novel manifold embedded knowledge transfer (MEKT) approach, which first aligns the covariance matrices of the EEG trials in the Riemannian manifold, extracts features in the tangent space, and then performs domain adaptation by minimizing the joint probability distribution shift between the source and the target domains, while preserving their geometric structures. MEKT can cope with one or multiple source domains, and can be computed efficiently. We also propose a domain transferability estimation (DTE) approach to identify the most beneficial source domains, in case there are a large number of source domains. Experiments on four EEG datasets from two different BCI paradigms demonstrated that MEKT outperformed several state-of-the-art transfer learning approaches, and DTE can reduce more than half of the computational cost when the number of source subjects is large, with little sacrifice of classification accuracy.


title: Learning Lenient Parsing & Typing via Indirect Supervision

abstract: Both professional coders and teachers frequently deal with imperfect (fragmentary, incomplete, ill-formed) code. Such fragments are common in STACKOVERFLOW; students also frequently produce ill-formed code, for which instructors, TAs (or students themselves) must find repairs. In either case, the developer experience could be greatly improved if such code could somehow be parsed & typed; this makes such code more amenable to use within IDEs and allows early detection and repair of potential errors. We introduce a lenient parser, which can parse & type fragments, even ones with simple errors. Training a machine learner to leniently parse and type imperfect code requires a large training set including many pairs of imperfect code and its repair; such training sets are limited by human effort and curation. In this paper, we present a novel, indirectly supervised, approach to train a lenient parser, without access to such human-curated training data. We leverage the huge corpus of mostly correct code available on Github, and the massive, efficient learning capacity of Transformer-based NN architectures. Using GitHub data, we first create a large dataset of fragments of code and corresponding tree fragments and type annotations; we then randomly corrupt the input fragments by seeding errors that mimic corruptions found in STACKOVERFLOW and student data. Using this data, we train high-capacity transformer models to overcome both fragmentation and corruption. With this novel approach, we can achieve reasonable performance on parsing & typing STACKOVERFLOW fragments; we also demonstrate that our approach performs well on shorter student error program and achieves best-in-class performance on longer programs that have more than 400 tokens. We also show that by blending Deepfix and our tool, we could achieve 77% accuracy, which outperforms all previously reported student error correction tools.


title: Positivity of the Rational Function of Gillis, Reznick and Zeilberger

abstract: We prove that the power series expansion of the rational function of Gillis, Reznick and Zeilberger (1983) has only nonnegative coefficients.


title: Bayesian Integration of Multi-resolutional Grid Codes for Spatial Cognition

abstract: Fourier-like summation of several grid cell modules with different spatial frequencies in the medial entorhinal cortex (MEC) has long been proposed to form the contours of place firing fields. Recent experiments largely, but not completely, support this theory. Place fields are obviously expanded by inactivation of dorsal MEC, which fits the hypothesis. However, contrary to the prediction, inactivation of ventral MEC is also weakly broaden the spatial place firing patterns. In this study, we derive the model from grid spatial frequencies represented by Gaussian profiles to a 1D place field by Bayesian inference, and further provide completely theoretical explanations for expansion of place fields and predictions for alignments of grid components. To understand the information transform across between neocortex, entorhinal cortex, and hippocampus, we propose spatial memory indexing theory from hippocampal indexing theory to investigate how neural dynamics work in the entorhinal-hippocampal circuit. The inputs of place cells in CA3 are converged from three grid modules with different grid spacings layer II of MEC by Bayesian mechanism. We resort to the robot system to test Fourier hypothesis and spatial memory indexing theory, and validate our proposed entorhinal-hippocampal model. And then we demonstrate its cognitive mapping capability on the KITTI odometry benchmark dataset. Results suggest that our model provides a rational theoretical explanation for the biological experimental results. Results also show that the proposed model is robust for simultaneous localization and mapping (SLAM) in the large-scale environment. Our proposed model theoretically supports for Fourier hypothesis in a general Bayesian mechanism, which may pertain to other neural systems in addition to spatial cognition.


title: Supertranslation Hair of Schwarzschild Black Hole: A Wilson Line Perspective

abstract: We demonstrate within the quantum field theoretical framework that an asymptotic particle falling into the black hole implants soft graviton hair on the horizon, conforming with the classical proposal of Hawking, Perry and Strominger. A key ingredient to this result is the construction of gravitational Wilson line dressings of an infalling scalar field, carrying a definite horizon supertranslation charge. It is shown that a typical Schwarzschild state is degenerate, and can be labeled by different soft supertranslation hairs parametrized for radial trajectories by the mass and energy of the infalling particle and its asymptotic point of contact with the horizon. The supertranslation zero modes are also obtained in terms of zero-frequency graviton operators, and are shown to be the expected canonical partners of the linearized horizon charge that enlarge the horizon Hilbert space.


title: Parameter-robust Uzawa-type iterative methods for double saddle point problems arising in Biot's consolidation and multiple-network poroelasticity models

abstract: This work is concerned with the iterative solution of systems of quasi-static multiple-network poroelasticity (MPET) equations describing flow in elastic porous media that is permeated by single or multiple fluid networks. Here, the focus is on a three-field formulation of the problem in which the displacement field of the elastic matrix and, additionally, one velocity field and one pressure field for each of the $n \ge 1$ fluid networks are the unknown physical quantities. Generalizing Biot's model of consolidation, which is obtained for $n=1$, the MPET equations for $n\ge1$ exhibit a double saddle point structure. The proposed approach is based on a framework of augmenting and splitting this three-by-three block system in such a way that the resulting block Gauss-Seidel preconditioner defines a fully decoupled iterative scheme for the flux-, pressure-, and displacement fields. In this manner, one obtains an augmented Lagrangian Uzawa-type method, the analysis of which is the main contribution of this work. The parameter-robust uniform linear convergence of this fixed-point iteration is proved by showing that its rate of contraction is strictly less than one independent of all physical and discretization parameters. The theoretical results are confirmed by a series of numerical tests that compare the new fully decoupled scheme to the very popular partially decoupled fixed-stress split iterative method, which decouples only flow--the flux and pressure fields remain coupled in this case--from the mechanics problem. We further test the performance of the block triangular preconditioner defining the new scheme when used to accelerate the GMRES algorithm.


title: Study of $\Xi^-$ hypernuclei in the Skyrme-Hartree-Fock approach

abstract: The properties of $\Xi^-$ hypernuclei are studied systematically using a two-dimensional Skyrme-Hartree-Fock approach combined with three different $\Xi N$ Skyrme forces fitted to reproduce the existing data. We explore the impurity effect of a single $\Xi^-$ hyperon on the radii, deformations, and density distributions of the nuclear core and point out qualitative differences between the different forces. We find that the $\Xi^-$ removal energy of $^{\hskip0.10em13}_{\Xi p}$B [$^{12}$C(g.s.)+ $\Xi^-$(1p)] calculated by the SLX3 force is 0.7 MeV, which is in good agreement with a possible value of $0.82\pm0.17\;$MeV from the KEK E176 experiment. The theoretical prediction for this weakly bound state depends strongly on the deformation of the nuclear core, which is analyzed in detail.


title: Parallelized Training of Restricted Boltzmann Machines using Markov-Chain Monte Carlo Methods

abstract: Restricted Boltzmann Machine (RBM) is a generative stochastic neural network that can be applied to collaborative filtering technique used by recommendation systems. Prediction accuracy of the RBM model is usually better than that of other models for recommendation systems. However, training the RBM model involves Markov-Chain Monte Carlo (MCMC) method, which is computationally expensive. In this paper, we have successfully applied distributed parallel training using Horovod framework to improve the training time of the RBM model. Our tests show that the distributed training approach of the RBM model has a good scaling efficiency. We also show that this approach effectively reduces the training time to little over 12 minutes on 64 CPU nodes compared to 5 hours on a single CPU node. This will make RBM models more practically applicable in recommendation systems.


title: A New Local Transformation Module for Few-shot Segmentation

abstract: Few-shot segmentation segments object regions of new classes with a few of manual annotations. Its key step is to establish the transformation module between support images (annotated images) and query images (unlabeled images), so that the segmentation cues of support images can guide the segmentation of query images. The existing methods form transformation model based on global cues, which however ignores the local cues that are verified in this paper to be very important for the transformation. This paper proposes a new transformation module based on local cues, where the relationship of the local features is used for transformation. To enhance the generalization performance of the network, the relationship matrix is calculated in a high-dimensional metric embedding space based on cosine distance. In addition, to handle the challenging mapping problem from the low-level local relationships to high-level semantic cues, we propose to apply generalized inverse matrix of the annotation matrix of support images to transform the relationship matrix linearly, which is non-parametric and class-agnostic. The result by the matrix transformation can be regarded as an attention map with high-level semantic cues, based on which a transformation module can be built simply.The proposed transformation module is a general module that can be used to replace the transformation module in the existing few-shot segmentation frameworks. We verify the effectiveness of the proposed method on Pascal VOC 2012 dataset. The value of mIoU achieves at 57.0% in 1-shot and 60.6% in 5-shot, which outperforms the state-of-the-art method by 1.6% and 3.5%, respectively.


title: Brownian motion of ellipsoidal particles on a granular magnetic bath

abstract: We study the Brownian motion of ellipsoidal particles lying on an agitated granular bath composed of magnetic particles. We quantify the mobility of different floating ellipsoidal particles using the mean square displacement and the mean square angular displacement, and relate the diffusion coefficients to the bath particle motion. In terms of the particle major radius $R$, we find the translational diffusion coefficient scales roughly as $1/R^2$ and the rotational diffusion coefficient scales as roughly $1/R^4$; this is consistent with the assumption that diffusion arises from random kicks of the bath particles underneath the floating particle. By varying the magnetic forcing, the bath particles diffusivity changes by a factor of ten; over this range, the translational and rotational diffusion of the floating particles change by a factor of fifty. However, the ratio of the two diffusion constants for the floating particles is forcing-independent. Unusual aspects of the floating particle motion include non-Gaussian statistics for their displacements.