title: Learning Multi-Agent Coordination through Connectivity-driven Communication

abstract: In artificial multi-agent systems, the ability to learn collaborative policies is predicated upon the agents' communication skills: they must be able to encode the information received from the environment and learn how to share it with other agents as required by the task at hand. We present a deep reinforcement learning approach, Connectivity Driven Communication (CDC), that facilitates the emergence of multi-agent collaborative behaviour only through experience. The agents are modelled as nodes of a weighted graph whose state-dependent edges encode pair-wise messages that can be exchanged. We introduce a graph-dependent attention mechanisms that controls how the agents' incoming messages are weighted. This mechanism takes into full account the current state of the system as represented by the graph, and builds upon a diffusion process that captures how the information flows on the graph. The graph topology is not assumed to be known a priori, but depends dynamically on the agents' observations, and is learnt concurrently with the attention mechanism and policy in an end-to-end fashion. Our empirical results show that CDC is able to learn effective collaborative policies and can over-perform competing learning algorithms on cooperative navigation tasks.


title: Visualizing modular forms

abstract: We examine several currently used techniques for visualizing complex-valued functions applied to modular forms. We plot several examples and study the benefits and limitations of each technique. We then introduce a method of visualization that can take advantage of colormaps in Python's matplotlib library, describe an implementation, and give more examples. Much of this discussion applies to general visualizations of complex-valued functions in the plane.


title: Image-to-Image Translation with Text Guidance

abstract: The goal of this paper is to embed controllable factors, i.e., natural language descriptions, into image-to-image translation with generative adversarial networks, which allows text descriptions to determine the visual attributes of synthetic images. We propose four key components: (1) the implementation of part-of-speech tagging to filter out non-semantic words in the given description, (2) the adoption of an affine combination module to effectively fuse different modality text and image features, (3) a novel refined multi-stage architecture to strengthen the differential ability of discriminators and the rectification ability of generators, and (4) a new structure loss to further improve discriminators to better distinguish real and synthetic images. Extensive experiments on the COCO dataset demonstrate that our method has a superior performance on both visual realism and semantic consistency with given descriptions.


title: Volkov-Pankratov states in topological superconductors

abstract: We study the in-gap states that appear at the boundaries of both 1D and 2D topological superconductors. While the massless Majorana quasiparticles are guaranteed to arise by the bulk-edge correspondence, we find that they could be accompanied by massive Volkov-Pankratov (VP) states which are present only when the interface is sufficiently smooth. These predictions can be tested in an s-wave superconductor with Rashba spin-orbit coupling placed on top of a magnetic domain wall. We calculate the spin-resolved local density of states of the VP states about the band inversion generated by a magnetic domain wall and find that they are oppositely spin-polarized on either side of the topological phase boundary. We also demonstrate that the spatial position, energy-level spacing, and spin polarization of the VP states can be modified by the introduction of in-plane electric fields.


title: A Study of Potential Code Borrowing and License Violations in Java Projects on GitHub

abstract: With an ever-increasing amount of open source software, the popularity of services like GitHub that facilitate code reuse, and common misconceptions about the licensing of open source software, the problem of license violations in the code is getting more and more prominent. In this study, we compile an extensive corpus of popular Java projects from GitHub, search it for code clones and perform an original analysis of possible code borrowing and license violations on the level of code fragments. We chose Java as a language because of its popularity in industry, where the plagiarism problem is especially relevant because of possible legal action. We analyze and discuss distribution of 94 different discovered and manually evaluated licenses in files and projects, differences in the licensing of files, distribution of potential code borrowing between licenses, various types of possible license violations, most violated licenses, etc. Studying possible license violations in specific blocks of code, we have discovered that 29.6% of them might be involved in potential code borrowing and 9.4% of them could potentially violate original licenses.


title: The structure of multigranular rough sets

abstract: We study multigranulation spaces of two equivalences. The lattice-theoretical properties of so-called "optimistic" and "pessimistic" multigranular approximation systems are given. We also consider the ordered sets of rough sets determined by these approximation pairs.


title: Complexity Analysis of Generalized and Fractional Hypertree Decompositions

abstract: Hypertree decompositions (HDs), as well as the more powerful generalized hypertree decompositions (GHDs), and the yet more general fractional hypertree decompositions (FHDs) are hypergraph decomposition methods successfully used for answering conjunctive queries and for solving constraint satisfaction problems. Every hypergraph $H$ has a width relative to each of these methods: its hypertree width $hw(H)$, its generalized hypertree width $ghw(H)$, and its fractional hypertree width $fhw(H)$, respectively. It is known that $hw(H)\leq k$ can be checked in polynomial time for fixed $k$, while checking $ghw(H)\leq k$ is NP-complete for $k \geq 3$. The complexity of checking $fhw(H)\leq k$ for a fixed $k$ has been open for over a decade. We settle this open problem by showing that checking $fhw(H)\leq k$ is NP-complete, even for $k=2$. The same construction allows us to prove also the NP-completeness of checking $ghw(H)\leq k$ for $k=2$. After that, we identify meaningful restrictions which make checking for bounded $ghw$ or $fhw$ tractable or allow for an efficient approximation of the $fhw$.


title: The Multiplayer Colonel Blotto Game

abstract: We initiate the study of the natural multiplayer generalization of the classic continuous Colonel Blotto game. The two-player Blotto game, introduced by Borel as a model of resource competition across $n$ simultaneous fronts, has been studied extensively for a century and seen numerous applications throughout the social sciences. Our work defines the multiplayer Colonel Blotto game and derives Nash equilibria for various settings of $k$ (number of players) and $n$. We also introduce a "Boolean" version of Blotto that becomes interesting in the multiplayer setting. The main technical difficulty of our work, as in the two-player theoretical literature, is the challenge of coupling various marginal distributions into a joint distribution satisfying a strict sum constraint. In contrast to previous works in the continuous setting, we derive our couplings algorithmically in the form of efficient sampling algorithms.


title: Runaway Electron Seed Formation at Reactor-Relevant Temperature

abstract: Systematic variation of the pre-disruption core electron temperature (Te) from 1 to 12 keV using an internal transport barrier scenario reveals a dramatic increase in the production of seed runaway electrons (REs), ultimately accessing near-complete conversion of the pre-disruption current into sub-MeV RE current. Injected Ar pellets are observed to ablate more intensely and promptly as Te rises. At high Te, the observed ablation exceeds predictions from published thermal ablation models. Simultaneously, the thermal quench (TQ) is observed to significantly shorten with increasing Te -- a surprising result. While the reason for the shorter TQ is not yet understood, candidate mechanisms include: insufficiently accurate thermal ablation models, enhanced ablation driven by the seed RE population, or significant parallel heat transport along stochastic fields. Kinetic modeling that self-consistently treats the plasma cooling via radiation, the induced electric field, and the formation of the seed RE is performed. Including the combined effect of the inherent dependence of hot-tail RE seeding on Te together with the shortened TQ, modeling recovers the progression towards near-complete conversion of the pre-disruption current to RE current as Te rises. Measurement of the HXR spectrum during the early current quench (CQ) reveals a trend of decreasing energy with pre-disruption Te. At the very highest Te ($\approx$ 12 keV), $\approx$ 100% conversion of the thermal current to runaway current is found. The energy of this peculiar RE beam is inferred to be sub-MeV as it emits vanishingly few MeV hard X-rays (HXRs). These measurements demonstrate novel TQ dynamics as Te is varied and illustrate the limitations of treating the RE seed formation problem without considering the inter-related dependencies of the pellet ablation, radiative energy loss, and resultant variations of the TQ duration.


title: Leveraging Affect Transfer Learning for Behavior Prediction in an Intelligent Tutoring System

abstract: In this work, we propose a video-based transfer learning approach for predicting problem outcomes of students working with an intelligent tutoring system (ITS). By analyzing a student's face and gestures, our method predicts the outcome of a student answering a problem in an ITS from a video feed. Our work is motivated by the reasoning that the ability to predict such outcomes enables tutoring systems to adjust interventions, such as hints and encouragement, and to ultimately yield improved student learning. We collected a large labeled dataset of student interactions with an intelligent online math tutor consisting of 68 sessions, where 54 individual students solved 2,749 problems. The dataset is public and available at https://www.cs.bu.edu/faculty/betke/research/learning/ . Working with this dataset, our transfer-learning challenge was to design a representation in the source domain of pictures obtained "in the wild" for the task of facial expression analysis, and transferring this learned representation to the task of human behavior prediction in the domain of webcam videos of students in a classroom environment. We developed a novel facial affect representation and a user-personalized training scheme that unlocks the potential of this representation. We designed several variants of a recurrent neural network that models the temporal structure of video sequences of students solving math problems. Our final model, named ATL-BP for Affect Transfer Learning for Behavior Prediction, achieves a relative increase in mean F-score of 50% over the state-of-the-art method on this new dataset.


title: Irradiated ocean planets bridge super-Earth and sub-Neptune populations

abstract: Small planets ($\sim$1--3.9 $\Rearth$) constitute more than half of the inventory of the 4000-plus exoplanets discovered so far. Smaller planets are sufficiently dense to be rocky, but those with radii larger than $\sim$1.6 $\Rearth$ are thought to display in many cases hydrogen/helium gaseous envelopes up to $\sim$30\% of the planetary mass. These low-mass planets are highly irradiated and the question of their origin, evolution, and possible links remains open. Here we show that close-in ocean planets affected by greenhouse effect display hydrospheres in supercritical state, which generate inflated atmospheres without invoking the presence of large hydrogen/helium gaseous envelopes. We present a new set of mass-radius relationships for ocean planets with different compositions and different equilibrium temperatures, which are found to be well adapted to low-density sub-Neptune planets. Our model suggests that super-Earths and water-rich sub-Neptunes could belong to the same family of planets, i.e. hydrogen/helium-free planets, with differences between their interiors simply resulting from the variation in the water content.


title: Linearizable boundary value problems for the elliptic sine-Gordon and the elliptic Ernst equations

abstract: By employing a novel generalization of the inverse scattering transform method known as the unified transform or Fokas method, it can be shown that the solution of certain physically significant boundary value problems for the elliptic sine-Gordon equation, as well as for the elliptic version of the Ernst equation, can be expressed in terms of the solution of appropriate $2 \times 2$-matrix Riemann--Hilbert (RH) problems. These RH problems are defined in terms of certain functions, called spectral functions, which involve the given boundary conditions, but also unknown boundary values. For arbitrary boundary conditions, the determination of these unknown boundary values requires the analysis of a nonlinear Fredholm integral equation. However, there exist particular boundary conditions, called linearizable, for which it is possible to bypass this nonlinear step and to characterize the spectral functions directly in terms of the given boundary conditions. Here, we review the implementation of this effective procedure for the following linearizable boundary value problems: (a) the elliptic sine-Gordon equation in a semi-strip with zero Dirichlet boundary values on the unbounded sides and with constant Dirichlet boundary value on the bounded side; (b) the elliptic Ernst equation with boundary conditions corresponding to a uniformly rotating disk of dust; (c) the elliptic Ernst equation with boundary conditions corresponding to a disk rotating uniformly around a central black hole; (d) the elliptic Ernst equation with vanishing Neumann boundary values on a rotating disk.


title: Maximin Fairness with Mixed Divisible and Indivisible Goods

abstract: We study fair resource allocation when the resources contain a mixture of divisible and indivisible goods, focusing on the well-studied fairness notion of maximin share fairness (MMS). With only indivisible goods, a full MMS allocation may not exist, but a constant multiplicative approximate allocation always does. We analyze how the MMS approximation guarantee would be affected when the resources to be allocated also contain divisible goods. In particular, we show that the worst-case MMS approximation guarantee with mixed goods is no worse than that with only indivisible goods. However, there exist problem instances to which adding some divisible resources would strictly decrease the MMS approximation ratio of the instance. On the algorithmic front, we propose a constructive algorithm that will always produce an $\alpha$-MMS allocation for any number of agents, where $\alpha$ takes values between $1/2$ and $1$ and is a monotone increasing function determined by how agents value the divisible goods relative to their MMS values.


title: Sterile neutrino oscillometry with Jinping

abstract: The existence of sterile neutrino is an open question in neutrino physics up to now. The method of neutrino oscillometry provides a powerful tool to test the common 3+1 sterile neutrino hypothesis, i.e. three active flavors and one sterile falvor. There are several antineutrino sources can be used for this method. One of them is the well known isotope chain of $^{144}{\rm Ce}-$$^{144}{\rm Pr}$ with initial activity around 50-100 kCi. It has compact size and might be installed either outside or inside the detector. Another one is the short-lived isotope $\rm^8Li$, which can be produced in nuclear reaction of a proton beam hitting beryllium target. The Lithium source has only the out-of-detector option due to its large size. The proposed Jinping water-based liquid scintillator detector will be used as a detection volume. Above experimental setups will allow us to cover the current best fit values of oscillation parameters with 90% C.L. At the same time, it is sensitive to the region of Neutrino-4 result.


title: On Khovanov Homology and Related Invariants

abstract: This paper begins with a survey of some applications of Khovanov homology to low-dimensional topology, with an eye toward extending these results to $\mathfrak{sl}(n)$ homologies. We extend Levine and Zemke's ribbon concordance obstruction from Khovanov homology to $\mathfrak{sl}(n)$ homology for $n \geq 2$, including the universal $\mathfrak{sl}(2)$ and $\mathfrak{sl}(3)$ homology theories. Inspired by Alishahi and Dowlin's bounds for the unknotting number coming from Khovanov homology and relying on spectral sequence arguments, we produce bounds on the alternation number of a knot. Lee and Bar-Natan spectral sequences also provide lower bounds on Turaev genus.


title: A Reflection Formula for the Gaussian Hypergeometric Function of Matrix Argument

abstract: We obtain a reflection formula for the Gaussian hypergeometric function of real symmetric matrix argument. We also show that this result extends to the Gaussian hypergeometric function defined over the symmetric cones, and even to generalizations of the Gaussian hypergeometric function defined in terms of series of Jack polynomials. Finally, we obtain a quadratic transformation formula for the Gaussian hypergeometric function of Hermitian $2 \times 2$ matrix argument.


title: A Competing Risks Model with Binary Time Varying Covariates for Estimation of Breast Cancer Risks in BRCA1 Families

abstract: Mammographic screening and prophylactic surgery such as risk-reducing salpingo oophorectomy (RRSO) can potentially reduce breast cancer risks among mutation carriers of BRCA families. The evaluation of these interventions is usually complicated by the fact that their effects on breast cancer may change over time and by the presence of competing risks. We introduce a correlated competing risks model to model breast and ovarian cancer risks within BRCA1 families that accounts for time-varying covariates (TVCs). Different parametric forms for these TVCs are proposed for more flexibility and a correlated gamma frailty model is specified to account for the correlated competing events. We also introduced a new ascertainment correction approach that accounts for the selection of families through probands affected with either breast or ovarian cancer, or unaffected. Our simulation studies demonstrate the good performances of our proposed approach in terms of bias and precision of the estimators of model parameters and cause-specific penetrances over different levels of familial correlations. We apply our new approach to 498 BRCA1 mutation carrier families recruited through the Breast Cancer Family Registry. Our results demonstrate the importance of the functional form of the TVC when assessing the role of RRSO on breast cancer. In particular, under the best fitting TVC model, the overall effect of RRSO on breast cancer risk was statistically significant in women with BRCA1.


title: Probabilistic aspects of $\Lambda$-coalescents in equilibrium and in evolution

abstract: We present approximation methods which lead to law of large numbers and fluctuation results for functionals of $\Lambda$-coalescents, both in the dust-free case and in the case with a dust component. Our focus is on the tree length (or total branch length) and the total external branch length, as well as the time to the most recent common ancestor and the size of the last merger. In the second part we discuss evolving coalescents. For certain Beta-coalescents we analyse fluctuations of a class of functionals in appropriate time scales. Finally we review results of Gufler on the representation of evolving $\Lambda$-coalescents in terms of the lookdown space.


title: A Step Forward from High-Entropy Ceramics to Compositionally Complex Ceramics: A New Perspective

abstract: High-entropy ceramics (HECs) have quickly gained attention since 2015. To date, nearly all work has focused on five-component, equimolar compositions. This perspective article briefly reviews different families of HECs and selected properties. Following a couple of our most recent studies, we propose a step forward to expand HECs to Compositionally Complex ceramics (CCCs) to include medium-entropy and non-equimolar compositions. Using defective fluorite and ordered pyrochlore oxides as two primary examples, we further consider the complexities of aliovalent cations and anion vacancies as well as ordered structures with two cation sublattices. Better thermally-insulating yet stiff CCCs have been found in non-equimolar compositions with optimal amounts of oxygen vacancies and in ordered pyrochlores with substantial size disorder. It is demonstrated that medium-entropy ceramics (MECs) can prevail over their high-entropy counterparts. The diversifying classes of CCCs provide even more possibilities than HECs to tailor the composition, defects, disorder/order, and, consequently, various properties.


title: Computing Shapley Values for Mean Width in 3-D

abstract: The Shapley value is a common tool in game theory to evaluate the importance of a player in a cooperative setting. In a geometric context, it provides a way to measure the contribution of a geometric object in a set towards some function on the set. Recently, Cabello and Chan (SoCG 2019) presented algorithms for computing Shapley values for a number of functions for point sets in the plane. More formally, a coalition game consists of a set of players $N$ and a characteristic function $v: 2^N \to \mathbb{R}$ with $v(\emptyset) = 0$. Let $\pi$ be a uniformly random permutation of $N$, and $P_N(\pi, i)$ be the set of players in $N$ that appear before player $i$ in the permutation $\pi$. The Shapley value of the game is defined to be $\phi(i) = \mathbb{E}_\pi[v(P_N(\pi, i) \cup \{i\}) - v(P_N(\pi, i))]$. More intuitively, the Shapley value represents the impact of player $i$'s appearance over all insertion orders. We present an algorithm to compute Shapley values in 3-D, where we treat points as players and use the mean width of the convex hull as the characteristic function. Our algorithm runs in $O(n^3\log^2{n})$ time and $O(n)$ space. Our approach is based on a new data structure for a variant of the dynamic convolution problem $(u, v, p)$, where we want to answer $u\cdot v$ dynamically. Our data structure supports updating $u$ at position $p$, incrementing and decrementing $p$ and rotating $v$ by $1$. We present a data structure that supports $n$ operations in $O(n\log^2{n})$ time and $O(n)$ space. Moreover, the same approach can be used to compute the Shapley values for the mean volume of the convex hull projection onto a uniformly random $(d - 2)$-subspace in $O(n^d\log^2{n})$ time and $O(n)$ space for a point set in $d$-dimensional space ($d \geq 3$).


title: Bounds on direct and indirect effects under treatment/mediator endogeneity and outcome attrition

abstract: Causal mediation analysis aims at disentangling a treatment effect into an indirect mechanism operating through an intermediate outcome or mediator, as well as the direct effect of the treatment on the outcome of interest. However, the evaluation of direct and indirect effects is frequently complicated by non-ignorable selection into the treatment and/or mediator, even after controlling for observables, as well as sample selection/outcome attrition. We propose a method for bounding direct and indirect effects in the presence of such complications using a method that is based on a sequence of linear programming problems. Considering inverse probability weighting by propensity scores, we compute the weights that would yield identification in the absence of complications and perturb them by an entropy parameter reflecting a specific amount of propensity score misspecification to set-identify the effects of interest. We apply our method to data from the National Longitudinal Survey of Youth 1979 to derive bounds on the explained and unexplained components of a gender wage gap decomposition that is likely prone to non-ignorable mediator selection and outcome attrition.


title: Optical Design and Characterization of 40-GHz Detector and Module for the BICEP Array

abstract: Families of cosmic inflation models predict a primordial gravitational-wave background that imprints B-mode polarization pattern in the Cosmic Microwave Background (CMB). High sensitivity instruments with wide frequency coverage and well-controlled systematic errors are needed to constrain the faint B-mode amplitude. We have developed antenna-coupled Transition Edge Sensor (TES) arrays for high-sensitivity polarized CMB observations over a wide range of millimeter-wave bands. BICEP Array, the latest phase of the BICEP/Keck experiment series, is a multi-receiver experiment designed to search for inflationary B-mode polarization to a precision $\sigma$(r) between 0.002 and 0.004 after 3 full years of observations, depending on foreground complexity and the degree of lensing removal. We describe the electromagnetic design and measured performance of BICEP Array low-frequency 40-GHz detector, their packaging in focal plane modules, and optical characterization including efficiency and beam matching between polarization pairs. We summarize the design and simulated optical performance, including an approach to improve the optical efficiency due to mismatch losses. We report the measured beam maps for a new broad-band corrugation design to minimize beam differential ellipticity between polarization pairs caused by interactions with the module housing frame, which helps minimize polarized beam mismatch that converts CMB temperature to polarization ($T \rightarrow P$) anisotropy in CMB maps.


title: Some Hermitian K-groups via geometric topology

abstract: We compute the first two symplectic quadratic K-theory groups of the integers, or equivalently, the first two stable homology groups of the group of symplectic integral matrices preserving the standard quadratic refinement. The main novelty in our calculation lies in its method, which is based on high-dimensional manifold theory.


title: Change actions: from incremental computation to discrete derivatives

abstract: The goal of this thesis is threefold: first, to provide a general semantic setting for reasoning about incremental computation. Second, to establish and clarify the connection between derivatives in the incremental sense and derivatives in the analytic sense, that is to say, to provide a common definition of derivative of which the previous two are particular instances. Third, to give a theoretically sound calculus for this general setting. To this end we define and explore the notions of change actions and differential maps between change actions and show how these notions relate to incremental computation through the concrete example of the semi-naive evaluation of Datalog queries. We also introduce the notion of a change action model as a setting for higher-order differentiation, and exhibit some interesting examples. Finally, we show how Cartesian difference categories, a family of particularly well-behaved change action models, generalise Cartesian differential categories and give rise to a calculus in the spirit of Ehrhard and Regnier's differential lambda-calculus.


title: Shortest path distance approximation using deep learning techniques

abstract: Computing shortest path distances between nodes lies at the heart of many graph algorithms and applications. Traditional exact methods such as breadth-first-search (BFS) do not scale up to contemporary, rapidly evolving today's massive networks. Therefore, it is required to find approximation methods to enable scalable graph processing with a significant speedup. In this paper, we utilize vector embeddings learnt by deep learning techniques to approximate the shortest paths distances in large graphs. We show that a feedforward neural network fed with embeddings can approximate distances with relatively low distortion error. The suggested method is evaluated on the Facebook, BlogCatalog, Youtube and Flickr social networks.


title: Linear and quadratic magnetoresistance in the semimetal SiP2

abstract: Multiple mechanisms for extremely large magnetoresistance (XMR) found in many topologically nontrivial/trivial semimetals have been theoretically proposed, but experimentally it is unclear which mechanism is responsible in a particular sample. In this article, by the combination of band structure calculations, numerical simulations of magnetoresistance (MR), Hall resistivity and de Haas-van Alphen (dHvA) oscillation measurements, we studied the MR anisotropy of SiP$_{2}$ which is verified to be a topologically trivial, incomplete compensation semimetal. It was found that as magnetic field, $H$, is applied along the $a$ axis, the MR exhibits an unsaturated nearly linear $H$ dependence, which was argued to arise from incomplete carriers compensation. For the $H$ $\parallel$ [101] orientation, an unsaturated nearly quadratic $H$ dependence of MR up to 5.88 $\times$ 10$^{4}$$\%$ (at 1.8 K, 31.2 T) and field-induced up-turn behavior in resistivity were observed, which was suggested due to the existence of hole open orbits extending along the $k_{x}$ direction. Good agreement of the experimental results with the simulations based on the calculated Fermi surface (FS) indicates that the topology of FS plays an important role in its MR.


title: Learning to Generate Levels From Nothing

abstract: Machine learning for procedural content generation has recently become an active area of research. Levels vary in both form and function and are mostly unrelated to each other across games. This has made it difficult to assemble suitably large datasets to bring machine learning to level design in the same way as it's been used for image generation. Here we propose Generative Playing Networks which design levels for itself to play. The algorithm is built in two parts; an agent that learns to play game levels, and a generator that learns the distribution of playable levels. As the agent learns and improves its ability, the space of playable levels, as defined by the agent, grows. The generator targets the agent's playability estimates to then update its understanding of what constitutes a playable level. We call this process of learning the distribution of data found through self-discovery with an environment, self-supervised inductive learning. Unlike previous approaches to procedural content generation, Generative Playing Networks are end-to-end differentiable and do not require human-designed examples or domain knowledge. We demonstrate the capability of this framework by training an agent and level generator for a 2D dungeon crawler game.


title: Rainbow reflection and trapping for energy harvesting

abstract: Important distinctions are made between two related wave control mechanisms that act to spatially separate frequency components; these so-called rainbow mechanisms either slow or reverse guided waves propagating along a graded line array. We demonstrate an important nuance distinguishing rainbow reflection from genuine rainbow trapping and show the implications of this distinction for energy harvesting designs. The difference between these related mechanisms is highlighted using a design methodology, applied to flexural waves on mass loaded thin Kirchhoff-Love elastic plates, and emphasised through simulations for energy harvesting in the setting of elasticity, by elastic metasurfaces of graded line arrays of resonant rods atop a beam. The delineation of these two effects, reflection and trapping, allows us to characterise the behaviour of forced line array systems and predict their capabilities for trapping, conversion and focusing of energy.


title: Electric circuit induced by quantum walk

abstract: We consider the Szegedy walk on graphs adding infinite length tails to a finite internal graph. We assume that on these tails, the dynamics is given by the free quantum walk. We set the $\ell^\infty$-category initial state so that the internal graph receives time independent input from the tails, say $\boldsymbol{\alpha}_{in}$, at every time step. We show that the response of the Szegedy walk to the input, which is the output, say $\boldsymbol{\beta}_{out}$, from the internal graph to the tails in the long time limit, is drastically changed depending on the reversibility of the underlying random walk. If the underlying random walk is reversible, we have $\boldsymbol{\beta}_{out}=\mathrm{Sz}(\boldsymbol{m}_{\delta E})\boldsymbol{\alpha}_{in}$, where the unitary matrix $\mathrm{Sz}(\boldsymbol{m}_{\delta E})$ is the reflection matrix to the unit vector $\boldsymbol{m}_{\delta E}$ which is determined by the boundary of the internal graph $\delta E$. Then the global dynamics so that the internal graph is regarded as one vertex recovers the local dynamics of the Szegedy walk in the long time limit. Moreover if the underlying random walk of the Szegedy walk is reversible, then we obtain that the stationary state is expressed by a linear combination of the reversible measure and the electric current on the electric circuit determined by the internal graph and the random walk's reversible measure. On the other hand, if the underlying random walk is not reversible, then the unitary matrix is just a phase flip; that is, $\boldsymbol{\beta}_{out}=-\boldsymbol{\alpha}_{in}$, and the stationary state is similar to the current flow but satisfies a different type of the Kirchhoff laws.


title: HAN-ECG: An Interpretable Atrial Fibrillation Detection Model Using Hierarchical Attention Networks

abstract: Atrial fibrillation (AF) is one of the most prevalent cardiac arrhythmias that affects the lives of more than 3 million people in the U.S. and over 33 million people around the world and is associated with a five-fold increased risk of stroke and mortality. like other problems in healthcare domain, artificial intelligence (AI)-based algorithms have been used to reliably detect AF from patients' physiological signals. The cardiologist level performance in detecting this arrhythmia is often achieved by deep learning-based methods, however, they suffer from the lack of interpretability. In other words, these approaches are unable to explain the reasons behind their decisions. The lack of interpretability is a common challenge toward a wide application of machine learning-based approaches in the healthcare which limits the trust of clinicians in such methods. To address this challenge, we propose HAN-ECG, an interpretable bidirectional-recurrent-neural-network-based approach for the AF detection task. The HAN-ECG employs three attention mechanism levels to provide a multi-resolution analysis of the patterns in ECG leading to AF. The first level, wave level, computes the wave weights, the second level, heartbeat level, calculates the heartbeat weights, and third level, window (i.e., multiple heartbeats) level, produces the window weights in triggering a class of interest. The detected patterns by this hierarchical attention model facilitate the interpretation of the neural network decision process in identifying the patterns in the signal which contributed the most to the final prediction. Experimental results on two AF databases demonstrate that our proposed model performs significantly better than the existing algorithms. Visualization of these attention layers illustrates that our model decides upon the important waves and heartbeats which are clinically meaningful in the detection task.


title: Ligand-protein interactions in lysozyme investigated through a dual-resolution model

abstract: A fully atomistic modelling of biological macromolecules at relevant length- and time-scales is often cumbersome or not even desirable, both in terms of computational effort required and it a posteriori analysis. This difficulty can be overcome with the use of multi-resolution models, in which different regions of the same system are concurrently described at different levels of detail. In enzymes, computationally expensive atomistic detail is crucial in the modelling of the active site in order to capture e.g. the chemically subtle process of ligand binding. In contrast, important yet more collective properties of the remainder of the protein can be reproduced with a coarser description. In the present work, we demonstrate the effectiveness of this approach through the calculation of the binding free energy of hen egg white lysozyme (HEWL) with the inhibitor di-N-acetylchitotriose. Particular attention is posed to the impact of the mapping, i.e. the selection of atomistic and coarse-grained residues, on the binding free energy. It is shown that, in spite of small variations of the binding free energy with respect to the active site resolution, the separate contributions coming from different energetic terms (such as electrostatic and van der Waals interactions) manifest a stronger dependence on the mapping, thus pointing to the existence of an optimal level of intermediate resolution.


title: On sufficient dimension reduction via principal asymmetric least squares

abstract: In this paper, we introduce principal asymmetric least squares (PALS) as a unified framework for linear and nonlinear sufficient dimension reduction. Classical methods such as sliced inverse regression (Li, 1991) and principal support vector machines (Li, Artemiou and Li, 2011) may not perform well in the presence of heteroscedasticity, while our proposal addresses this limitation by synthesizing different expectile levels. Through extensive numerical studies, we demonstrate the superior performance of PALS in terms of both computation time and estimation accuracy. For the asymptotic analysis of PALS for linear sufficient dimension reduction, we develop new tools to compute the derivative of an expectation of a non-Lipschitz function.


title: Impact of collision models on the physical properties and the stability of lattice Boltzmann methods

abstract: The lattice Boltzmann method (LBM) is known to suffer from stability issues when the collision model relies on the BGK approximation, especially in the zero viscosity limit and for non-vanishing Mach numbers. To tackle this problem, two kinds of solutions were proposed in the literature. They consist in changing either the numerical discretization (finite-volume, finite-difference, spectral-element, etc) of the discrete velocity Boltzmann equation (DVBE), or the collision model. In this work, the latter solution is investigated in details. More precisely, we propose a comprehensive comparison of (static relaxation time based) collision models, in terms of stability, and with preliminary results on their accuracy, for the simulation of isothermal high-Reynolds number flows in the (weakly) compressible regime. It starts by investigating the possible impact of collision models on the macroscopic behavior of stream-and-collide based D2Q9-LBMs, which clarifies the exact physical properties of collision models on LBMs. It is followed by extensive linear and numerical stability analyses, supplemented with an accuracy study based on the transport of vortical structures over long distances. In order to draw conclusions as general as possible, the most common moment spaces (raw, central, Hermite, central Hermite and cumulant), as well as regularized approaches, are considered for the comparative studies. LBMs based on dynamic collision mechanisms (entropic collision, subgrid scale models, explicit filtering, etc) are also briefly discussed.


title: Topological magnonics in the two-dimensional van der Waals magnet CrI3

abstract: In this article, we calculate the magnon spectrum of Kitaev-Heisenberg magnets. This model has been recently proposed as a spin Hamiltonian to model $\text{CrI}_3$ and other two-dimensional magnets. It is a minimal spin Hamiltonian that includes a contribution stemming from a Heisenberg, isotropic exchange, and a contribution arising from a Kitaev interaction, anisotropic and frustrated exchange. Our calculations reveal the topological nature of the magnons and a gap that opens at the $K$ and $K^\prime$ points. These topological properties give rise to effects such as thermal Hall effect. In addition to the bulk properties, we calculate the magnon spectrum of nanoribbons illustrating the corresponding edge states.


title: Combined minivalley and layer control in twisted double bilayer graphene

abstract: Control over minivalley polarization and interlayer coupling is demonstrated in double bilayer graphene twisted with an angle of 2.37$^\circ$. This intermediate angle is small enough for the minibands to form and large enough such that the charge carrier gases in the layers can be tuned independently. Using a dual-gated geometry we identify and control all possible combinations of minivalley polarization via the population of the two bilayers. An applied displacement field opens a band gap in either of the two bilayers, allowing us to even obtain full minivalley polarization. In addition, the wavefunctions of the minivalleys are mixed by tuning through a Lifshitz transition, where the Fermi surface topology changes. The high degree of control makes twisted double bilayer graphene a promising platform for valleytronics devices such as valley valves, filters and logic gates.


title: Generating ${\rm N00N}$-states of surface plasmon-polariton pairs with a nanoparticle

abstract: We consider a generation of two-particle quantum states in the process of spontaneous parametric down-conversion of light by a dielectric nanoparticle with $\chi^{(2)}$ response. As a particular example, we study the generation of surface plasmon-polariton pairs with a ${\rm GaAs}$ nanoparticle located at the silver-air interface. We show that for certain excitation geometries, ${\rm N00N}$-states of surface plasmon-polariton pairs could be obtained. The effect can be physically interpreted as a result of quantum interference between pairs of induced sources, each emitting either signal or idler plasmon. We then relate the resulting ${\rm N00N}$-pattern to the general symmetry properties of dyadic Green's function of a dipole emitter exciting surface waves. It renders the considered effect as a general way towards a robust generation of ${\rm N00N}$-states of surface waves using spontaneous parametric down-conversion in $\chi^{(2)}$ nanoparticles.


title: A Matching Mechanism with Anticipatory Tolls for Congestion Pricing

abstract: This paper presents a matching mechanism for assigning drivers to routes where the drivers pay a toll for the marginal delay they impose on other drivers. The simple matching mechanism is derived from the RANKING algorithm for online bipartite matching proposed by Karp et al. [8]. The toll, which is anticipatory in design, is an adaption of one proposed by Dong et al. [5]. Our research proves that the matching mechanism proposed here is pareto user-optimal and can be adapted to give network optimal results for the minimizing total social cost of travel.


title: Structure and dynamics of hydrodynamically interacting finite-size Brownian particles in a spherical cavity: spheres and cylinders

abstract: The structure and dynamics of confined suspensions of particles of arbitrary shape is of interest in multiple disciplines, from biology to engineering. Theoretical studies are often limited by the complexity of long-range particle-particle and particle-wall forces, including many-body fluctuating hydrodynamic interactions. Here, we report a computational study on the diffusion of spherical and cylindrical particles confined in a spherical cavity. We rely on an Immersed-Boundary General geometry Ewald-like method to capture lubrication and long-range hydrodynamics, and include appropriate non-slip conditions at the confining walls. A Chebyshev polynomial approximation is used to satisfy the fluctuation-dissipation theorem for the Brownian suspension. We explore how lubrication, long-range hydrodynamics, particle volume fraction and shape affect the equilibrium structure and the diffusion of the particles. It is found that once the particle volume fraction is greater than $10\%$, the particles start to form layered aggregates that greatly influence particle dynamics. Hydrodynamic interactions strongly influence the particle diffusion by inducing spatially dependent short-time diffusion coefficients, stronger wall effects on the particle diffusion towards the walls, and a sub-diffusive regime --caused by crowding-- in the long-time particle mobility. The level of asymmetry of the cylindrical particles considered here is enough to induce an orientational order in the layered structure, decreasing the diffusion rate and facilitating a transition to the crowded mobility regime at low particle concentrations. Our results offer fundamental insights into the diffusion and distribution of globular and fibrillar proteins inside cells.


title: HypoML: Visual Analysis for Hypothesis-based Evaluation of Machine Learning Models

abstract: In this paper, we present a visual analytics tool for enabling hypothesis-based evaluation of machine learning (ML) models. We describe a novel ML-testing framework that combines the traditional statistical hypothesis testing (commonly used in empirical research) with logical reasoning about the conclusions of multiple hypotheses. The framework defines a controlled configuration for testing a number of hypotheses as to whether and how some extra information about a "concept" or "feature" may benefit or hinder a ML model. Because reasoning multiple hypotheses is not always straightforward, we provide HypoML as a visual analysis tool, with which, the multi-thread testing data is transformed to a visual representation for rapid observation of the conclusions and the logical flow between the testing data and hypotheses.We have applied HypoML to a number of hypothesized concepts, demonstrating the intuitive and explainable nature of the visual analysis.


title: Health Assessment and Prognostics Based on Higher Order Hidden Semi-Markov Models

abstract: This paper presents a new and flexible prognostics framework based on a higher order hidden semi-Markov model (HOHSMM) for systems or components with unobservable health states and complex transition dynamics. The HOHSMM extends the basic hidden Markov model (HMM) by allowing the hidden state to depend on its more distant history and assuming generally distributed state duration. An effective Gibbs sampling algorithm is designed for statistical inference of an HOHSMM. The performance of the proposed HOHSMM sampler is evaluated by conducting a simulation experiment. We further design a decoding algorithm to estimate the hidden health states using the learned model. Remaining useful life (RUL) is predicted using a simulation approach given the decoded hidden states. The practical utility of the proposed prognostics framework is demonstrated by a case study on NASA turbofan engines. The results show that the HOHSMM-based prognostics framework provides good hidden health state assessment and RUL estimation for complex systems.


title: A Second look at Exponential and Cosine Step Sizes: Simplicity, Adaptivity, and Performance

abstract: Stochastic Gradient Descent (SGD) is a popular tool in training large-scale machine learning models. Its performance, however, is highly variable, depending crucially on the choice of the step sizes. Accordingly, a variety of strategies for tuning the step sizes have been proposed, ranging from coordinate-wise approaches (a.k.a. ``adaptive'' step sizes) to sophisticated heuristics to change the step size in each iteration. In this paper, we study two step size schedules whose power has been repeatedly confirmed in practice: the exponential and the cosine step sizes. For the first time, we provide theoretical support for them proving convergence rates for smooth non-convex functions, with and without the Polyak-\L{}ojasiewicz (PL) condition. Moreover, we show the surprising property that these two strategies are \emph{adaptive} to the noise level in the stochastic gradients of PL functions. That is, contrary to polynomial step sizes, they achieve almost optimal performance without needing to know the noise level nor tuning their hyperparameters based on it. Finally, we conduct a fair and comprehensive empirical evaluation of real-world datasets with deep learning architectures. Results show that, even if only requiring at most two hyperparameters to tune, these two strategies best or match the performance of various finely-tuned state-of-the-art strategies.


title: Solving Missing-Annotation Object Detection with Background Recalibration Loss

abstract: This paper focuses on a novel and challenging detection scenario: A majority of true objects/instances is unlabeled in the datasets, so these missing-labeled areas will be regarded as the background during training. Previous art on this problem has proposed to use soft sampling to re-weight the gradients of RoIs based on the overlaps with positive instances, while their method is mainly based on the two-stage detector (i.e. Faster RCNN) which is more robust and friendly for the missing label scenario. In this paper, we introduce a superior solution called Background Recalibration Loss (BRL) that can automatically re-calibrate the loss signals according to the pre-defined IoU threshold and input image. Our design is built on the one-stage detector which is faster and lighter. Inspired by the Focal Loss formulation, we make several significant modifications to fit on the missing-annotation circumstance. We conduct extensive experiments on the curated PASCAL VOC and MS COCO datasets. The results demonstrate that our proposed method outperforms the baseline and other state-of-the-arts by a large margin. Code available: https://github.com/Dwrety/mmdetection-selective-iou.


title: Lithium Enrichment Signatures of Planetary Engulfment Events in Evolved Stars

abstract: Planetary engulfment events have long been proposed as a lithium (Li) enrichment mechanism contributing to the population of Li-rich giants (A(Li) >= 1.5 dex). Using MESA stellar models and A(Li) abundance measurements obtained by the GALAH survey, we calculate the strength and observability of the surface Li enrichment signature produced by the engulfment of a hot Jupiter (HJ). We consider solar-metallicity stars in the mass range of 1-2 $\mathrm{M_\odot}$ and the Li supplied by a HJ of 1.0 $\mathrm{M_J}$. We explore engulfment events that occur near the main sequence turn-off (MSTO) and out to orbital separations of 0.1 AU. We map our results onto the Hertzsprung-Russell Diagram, revealing the statistical significance and survival time of lithium enrichment. We identify the parameter space of masses and evolutionary phases where the engulfment of a HJ can lead to Li enrichment signatures at a $5\sigma$ confidence level and with meteoritic abundance strengths. The most compelling strengths and survival times of engulfment-derived Li enrichment are found among host stars of 1.4 $\mathrm{M_\odot}$ near the MSTO. Our calculations indicate that planetary engulfment is not a viable enrichment pathway for stars that have evolved beyond the subgiant branch. For these sources, observed Li enhancements are likely to be produced by other mechanisms, such as the Cameron-Fowler process or the accretion of material from an asymptotic giant branch (AGB) companion. Our results do not account for second-order effects, such as extra mixing processes, which can further dilute Li enrichment signatures.


title: Improved Classical and Quantum Algorithms for Subset-Sum

abstract: We present new classical and quantum algorithms for solving random subset-sum instances. First, we improve over the Becker-Coron-Joux algorithm (EUROCRYPT 2011) from $\tilde{\mathcal{O}}(2^{0.291 n})$ downto $\tilde{\mathcal{O}}(2^{0.283 n})$, using more general representations with values in $\{-1,0,1,2\}$. Next, we improve the state of the art of quantum algorithms for this problem in several directions. By combining the Howgrave-Graham-Joux algorithm (EUROCRYPT 2010) and quantum search, we devise an algorithm with asymptotic cost $\tilde{\mathcal{O}}(2^{0.236 n})$, lower than the cost of the quantum walk based on the same classical algorithm proposed by Bernstein, Jeffery, Lange and Meurer (PQCRYPTO 2013). This algorithm has the advantage of using \emph{classical} memory with quantum random access, while the previously known algorithms used the quantum walk framework, and required \emph{quantum} memory with quantum random access. We also propose new quantum walks for subset-sum, performing better than the previous best time complexity of $\tilde{\mathcal{O}}(2^{0.226 n})$ given by Helm and May (TQC 2018). We combine our new techniques to reach a time $\tilde{\mathcal{O}}(2^{0.216 n})$. This time is dependent on a heuristic on quantum walk updates, formalized by Helm and May, that is also required by the previous algorithms. We show how to partially overcome this heuristic, and we obtain an algorithm with quantum time $\tilde{\mathcal{O}}(2^{0.218 n})$ requiring only the standard classical subset-sum heuristics.


title: Classical limit for the varying-mass Schr\"odinger equation with random inhomogeneities

abstract: The varying-mass Schr\"odinger equation (VMSE) has been successfully applied to model electronic properties of semiconductor hetero-structures, for example, quantum dots and quantum wells. In this paper, we consider VMSE with small random heterogeneities, and derive a radiative transfer equation as its asymptotic limit. The main tool is to systematically apply the Wigner transform in the classical regime when the rescaled Planck constant $\epsilon\ll 1$, and expand the Wigner equation to proper orders of $\epsilon$. As a proof of concept, we numerically compute both VMSE and its limiting radiative transfer equation, and show that their solutions agree well in the classical regime.


title: Delay Effects on Amplitude Death, Oscillation Death, and Renewed Limit Cycle Behavior in Cyclically Coupled Oscillators

abstract: The effects of a distributed 'weak generic kernel' delay on cyclically coupled limit cycle and chaotic oscillators are considered. For coupled Van der Pol oscillators (and in fact, other oscillators as well) the delay can produce transitions from amplitude death(AD) or oscillation death (OD) to Hopf bifurcation-induced periodic behavior, with the delayed limit cycle shrinking or growing as the delay is varied towards or away from the bifurcation point respectively. The transition from AD to OD is mediated here via a pitchfork bifurcation, as seen earlier for other couplings as well. Also, the cyclically coupled undelayed van der Pol system here is already in a state of AD/OD, and introducing the delay allows both oscillations and AD/OD as the delay parameter is varied. This is in contrast to other limit cycle systems, where diffusive coupling alone does not result in the onset of AD/OD. For systems where the individual oscillators are chaotic, such as a Sprott oscillator system or a coupled van der Pol-Rayleigh system with parametric forcing, the delay may produce AD/OD (as in the Sprott case), with the AD to OD transition now occurring via a transcritical bifurcation instead. However, this may not be possible, and the delay might just vary the attractor shape. In either of these situations however, increased delay strength tends to cause the system to have simpler behavior, streamlining the shape of the attractor, or shrinking it in cases with oscillations.


title: Model-based wavefront shaping microscopy

abstract: Wavefront shaping is increasingly being used in modern microscopy to obtain distortion-free, high-resolution images deep inside inhomogeneous media. Wavefront shaping methods typically rely on the presence of a 'guidestar' in order to find the optimal wavefront to mitigate the scattering of light. However, this condition cannot be satisfied in most biomedical applications. Here, we introduce a novel, guidestar-free wavefront shaping method in which the optimal wavefront is computed using a digital model of the sample. The refractive index model of the sample, that serves as the input for the computation, is constructed in-situ by the microscope itself. In a proof of principle imaging experiment, we demonstrate a large improvement in the two-photon fluorescence signal through a diffuse medium, outperforming the state-of-the-art wavefront shaping techniques by a factor of 21.


title: Measurements of Nuclear Reactions that Create and Destroy Li and Be during BBN

abstract: I review measurements of the most important reactions involved in the creation and destruction of Li and Be during big bang nucleosynthesis (BBN) as well as their uncertainties and the relative contributions they make to the uncertainty in the primordial $^7$Li abundance ($^7$Li/H). Examining the sensitivity of calculated $^7$Li/H to these reactions as predicted by different BBN codes I find no significant differences. I compare my calculation of primordial $^7$Li/H to some recently published values and conclude that in the absence of a major undetected experimental blunder, nuclear physics uncertainties cannot account for the cosmological Li problem. With an estimated 13% uncertainty in the calculated abundance, the discrepancy with observation amounts to some $4.6\sigma$.


title: A note on the continuity in the Hurst index of the solution of rough differential equations driven by a fractional Brownian motion

abstract: Within the rough path framework we prove the continuity of the solution to random differential equations driven by fractional Brownian motion with respect to the Hurst parameter $H$ when $H \in (1/3, 1/2]$.


title: A Bounded Measure for Estimating the Benefit of Visualization

abstract: Information theory can be used to analyze the cost-benefit of visualization processes. However, the current measure of benefit contains an unbounded term that is neither easy to estimate nor intuitive to interpret. In this work, we propose to revise the existing cost-benefit measure by replacing the unbounded term with a bounded one. We examine a number of bounded measures that include the Jenson-Shannon divergence and a new divergence measure formulated as part of this work. We use visual analysis to support the multi-criteria comparison, narrowing the search down to those options with better mathematical properties. We apply those remaining options to two visualization case studies to instantiate their uses in practical scenarios, while the collected real world data further informs the selection of a bounded measure, which can be used to estimate the benefit of visualization.


title: Stabilizing Differentiable Architecture Search via Perturbation-based Regularization

abstract: Differentiable architecture search (DARTS) is a prevailing NAS solution to identify architectures. Based on the continuous relaxation of the architecture space, DARTS learns a differentiable architecture weight and largely reduces the search cost. However, its stability has been challenged for yielding deteriorating architectures as the search proceeds. We find that the precipitous validation loss landscape, which leads to a dramatic performance drop when distilling the final architecture, is an essential factor that causes instability. Based on this observation, we propose a perturbation-based regularization - SmoothDARTS (SDARTS), to smooth the loss landscape and improve the generalizability of DARTS-based methods. In particular, our new formulations stabilize DARTS-based methods by either random smoothing or adversarial attack. The search trajectory on NAS-Bench-1Shot1 demonstrates the effectiveness of our approach and due to the improved stability, we achieve performance gain across various search spaces on 4 datasets. Furthermore, we mathematically show that SDARTS implicitly regularizes the Hessian norm of the validation loss, which accounts for a smoother loss landscape and improved performance.


title: On the Coefficients of $(\mathbb{Z}/p)^n$-Equivariant Ordinary Cohomology with Coefficients in $\mathbb{Z}/p$

abstract: This note contains a generalization to $p>2$ of the authors' previous calculations of the coefficients of $(\mathbb{Z}/2)^n$-equivariant ordinary cohomology with coefficients in the constant $\mathbb{Z}/2$-Mackey functor. The algberaic results by S.Kriz allow us to calculate the coefficients of the geometric fixed point spectrum $\Phi^{(\mathbb{Z}/p)^n}H\mathbb{Z}/p$, and more generally, the $\mathbb{Z}$-graded coefficients of the localization of $H\mathbb{Z}/p_{(\mathbb{Z}/p)^n}$ by inverting any chosen set of embeddings $S^0\rightarrow S^{\alpha_i}$ where $\alpha_i$ are non-trivial irreducible representations. We also calculate the $RO(G)^+$-graded coefficients of $H\mathbb{Z}/p_{(\mathbb{Z}/p)^n}$, which means the cohomology of a point indexed by an actual (not virtual) representation. (This is the "non-derived" part, which has a nice algebraic description.)


title: Toward room-temperature nanoscale skyrmions in ultrathin films

abstract: Breaking the dilemma between small size and room temperature stability is a necessary prerequisite for skyrmion-based information technology. Here we demonstrate by means of rate theory and an atomistic spin Hamiltonian that the stability of isolated skyrmions in ultrathin ferromagnetic films can be enhanced by the concerted variation of magnetic interactions while keeping the skyrmion size unchanged. We predict film systems where the lifetime of sub-10 nm skyrmions can reach years at ambient conditions. The long lifetime of such small skyrmions is due to exceptionally large Arrhenius pre-exponential and the stabilizing effect of the energy barrier is insignificant at room temperature. A dramatic increase in the pre-exponential is achieved thanks to softening of magnon modes of the skyrmion, thereby increasing the entropy of the skyrmion with respect to the transition state for collapse. Increasing the number of skyrmion deformation modes should be a guiding principle for the realization of nanoscale, room-temperature stable skyrmions.


title: Enabling high repetition rate nonlinear THz science with a kilowatt-class sub-100 fs laser source

abstract: Manipulating the atomic and electronic structure of matter with strong terahertz (THz) fields while probing the response with ultrafast pulses at x-ray free electron lasers (FELs) has offered unique insights into a multitude of physical phenomena in solid state and atomic physics. Recent upgrades of x-ray FEL facilities are pushing to much higher repetition rates, enabling unprecedented signal to noise for pump probe experiments. This requires the development of suitable THz pump sources that are able to deliver intense pulses at compatible repetition rates. Here we present a high power laser-driven THz source based on optical rectification in LiNbO3 using tilted pulse front pumping. Our source is driven by a kilowatt-level Yb:YAG amplifier system operating at 100 kHz repetition rate and employing nonlinear spectral broadening and recompression to achieve sub-100 fs pulses at 1030 nm wavelength. We demonstrate a maximum of 144 mW average THz power (1.44 uJ pulse energy), consisting of single-cycle pulses centered at 0.6 THz with a peak electric field strength exceeding 150 kV/cm. These high field pulses open up a range of possibilities for nonlinear time-resolved experiments with x-ray probing at unprecedented rates.


title: Geom-GCN: Geometric Graph Convolutional Networks

abstract: Message-passing neural networks (MPNNs) have been successfully applied to representation learning on graphs in a variety of real-world applications. However, two fundamental weaknesses of MPNNs' aggregators limit their ability to represent graph-structured data: losing the structural information of nodes in neighborhoods and lacking the ability to capture long-range dependencies in disassortative graphs. Few studies have noticed the weaknesses from different perspectives. From the observations on classical neural network and network geometry, we propose a novel geometric aggregation scheme for graph neural networks to overcome the two weaknesses. The behind basic idea is the aggregation on a graph can benefit from a continuous space underlying the graph. The proposed aggregation scheme is permutation-invariant and consists of three modules, node embedding, structural neighborhood, and bi-level aggregation. We also present an implementation of the scheme in graph convolutional networks, termed Geom-GCN (Geometric Graph Convolutional Networks), to perform transductive learning on graphs. Experimental results show the proposed Geom-GCN achieved state-of-the-art performance on a wide range of open datasets of graphs. Code is available at https://github.com/graphdml-uiuc-jlu/geom-gcn.


title: Graphs with multi-$4$-cycles and the Barnette's conjecture

abstract: Let ${\cal H}$ denote the family of all graphs with multi-$4$-cycles and suppose that $G \in {\cal H}$. Then, $G$ is a bipartite graph with a vertex bipartition $\{V_{\alpha}, V_{\beta}\}$. We prove that for every vertex $v \in V_{\beta}$ and for every $2$-colouring $V_{\alpha} \rightarrow \{1, 2\}$ there exists a $2$-colouring $V_{\beta} \rightarrow \{1, 2\}$ such that every cycle in $G$ is not monochromatic and $b(v) = 1$ ($b(v) = 2$). Let now $G$ be a simple even plane triangulation with a vertex $3$-partition $\{V_{1}, V_{2}, V_{3}\}$. Denote by $B_{i}$, $i = 1, 2, 3$, the set of all vertices in $V_i$ of degree at least $6$ in $G$. Suppose that $G[B_{1}\cup B_{3}]$ ($G[B_{2}\cup B_{3}]$) is a subgraph of $G$ induced by the set $B_{1}\cup B_{3}$ ($B_{2}\cup B_{3}$, respectively). Let $G^{*}$ be the dual graph of $G$ with the following $3$-face-colouring: a face $f$ of $G^{*}$ is coloured with $i$ if and only if the vertex $v = f^{*} \in V_{i}$. We prove that if $H = G[B_{1}\cup B_{3}] \cup G[B_{2}\cup B_{3}] \in {\cal H}$, then, for any edge chosen on a face coloured $3$ and of size at least $6$ in $G^{*}$, there exists a Hamilton cycle of $G^{*}$ which avoids this edge. Moreover, if every component of $H$ is $2$-connected, then there exists a Hamilton cycle of $G^{*}$ such that for every face coloured $3$ it avoids every second edge of this face or it avoids at most two edges of this face.


title: Multiscale Non-stationary Stochastic Bandits

abstract: Classic contextual bandit algorithms for linear models, such as LinUCB, assume that the reward distribution for an arm is modeled by a stationary linear regression. When the linear regression model is non-stationary over time, the regret of LinUCB can scale linearly with time. In this paper, we propose a novel multiscale changepoint detection method for the non-stationary linear bandit problems, called Multiscale-LinUCB, which actively adapts to the changing environment. We also provide theoretical analysis of regret bound for Multiscale-LinUCB algorithm. Experimental results show that our proposed Multiscale-LinUCB algorithm outperforms other state-of-the-art algorithms in non-stationary contextual environments.


title: Running vacuum model in non-flat universe

abstract: We investigate observational constraints on the running vacuum model (RVM) of $\Lambda=3\nu (H^{2}+K/a^2)+c_0$ in the spatially curved universe, where $\nu$ is the model parameter, $K$ corresponds to the spatial curvature constant, and $c_{0}$ is a constant defined by the boundary conditions. As $\dot{\Lambda} \ne 0$, there are energy exchanges between vacuum, matter and radiation in RVM. We study the "geometrical degeneracy" of RVM on the CMB power spectra. By fitting the cosmological data, we find that the values of $\chi^2$ in RVM and $\Lambda$CDM are similar to each other for the non-flat universe. Explicitly, we obtain the constraints of $\nu\leq O(10^{-4})$ (68 $\%$ C.L.) and $|\Omega_K|\leq O(10^{-2})$ (95 $\%$ C.L.) in our study. In addition, we show that the cosmological constraints of $\Sigma m_{\nu}=0.416^{+0.311}_{-0.407}$ (RVM) and $\Sigma m_{\nu}=0.497^{+0.335}_{-0.387}$ ($\Lambda$CDM) at 95$\%$ C.L. for the neutrino mass sum are relaxed in both models in the spatially curved universe.


title: CSM-NN: Current Source Model Based Logic Circuit Simulation -- A Neural Network Approach

abstract: The miniaturization of transistors down to 5nm and beyond, plus the increasing complexity of integrated circuits, significantly aggravate short channel effects, and demand analysis and optimization of more design corners and modes. Simulators need to model output variables related to circuit timing, power, noise, etc., which exhibit nonlinear behavior. The existing simulation and sign-off tools, based on a combination of closed-form expressions and lookup tables are either inaccurate or slow, when dealing with circuits with more than billions of transistors. In this work, we present CSM-NN, a scalable simulation framework with optimized neural network structures and processing algorithms. CSM-NN is aimed at optimizing the simulation time by accounting for the latency of the required memory query and computation, given the underlying CPU and GPU parallel processing capabilities. Experimental results show that CSM-NN reduces the simulation time by up to $6\times$ compared to a state-of-the-art current source model based simulator running on a CPU. This speedup improves by up to $15\times$ when running on a GPU. CSM-NN also provides high accuracy levels, with less than $2\%$ error, compared to HSPICE.


title: NN-PARS: A Parallelized Neural Network Based Circuit Simulation Framework

abstract: The shrinking of transistor geometries as well as the increasing complexity of integrated circuits, significantly aggravate nonlinear design behavior. This demands accurate and fast circuit simulation to meet the design quality and time-to-market constraints. The existing circuit simulators which utilize lookup tables and/or closed-form expressions are either slow or inaccurate in analyzing the nonlinear behavior of designs with billions of transistors. To address these shortcomings, we present NN-PARS, a neural network (NN) based and parallelized circuit simulation framework with optimized event-driven scheduling of simulation tasks to maximize concurrency, according to the underlying GPU parallel processing capabilities. NN-PARS replaces the required memory queries in traditional techniques with parallelized NN-based computation tasks. Experimental results show that compared to a state-of-the-art current-based simulation method, NN-PARS reduces the simulation time by over two orders of magnitude in large circuits. NN-PARS also provides high accuracy levels in signal waveform calculations, with less than $2\%$ error compared to HSPICE.


title: Improving Efficiency in Neural Network Accelerator Using Operands Hamming Distance optimization

abstract: Neural network accelerator is a key enabler for the on-device AI inference, for which energy efficiency is an important metric. The data-path energy, including the computation energy and the data movement energy among the arithmetic units, claims a significant part of the total accelerator energy. By revisiting the basic physics of the arithmetic logic circuits, we show that the data-path energy is highly correlated with the bit flips when streaming the input operands into the arithmetic units, defined as the hamming distance of the input operand matrices. Based on the insight, we propose a post-training optimization algorithm and a hamming-distance-aware training algorithm to co-design and co-optimize the accelerator and the network synergistically. The experimental results based on post-layout simulation with MobileNetV2 demonstrate on average 2.85X data-path energy reduction and up to 8.51X data-path energy reduction for certain layers.


title: Cooperative Observation of Targets moving over a Planar Graph with Prediction of Positions

abstract: Consider a team with two types of agents: targets and observers. Observers are aerial UAVs that observe targets moving on land with their movements restricted to the paths that form a planar graph on the surface. Observers have limited range of vision and targets do not avoid observers. The objective is to maximize the integral of the number of targets observed in the observation interval. Taking advantage of the fact that the future positions of targets in the short term are predictable, we show in this article a modified hill climbing algorithm that surpasses its previous versions in this new setting of the CTO problem.


title: Exploiting the Matching Information in the Support Set for Few Shot Event Classification

abstract: The existing event classification (EC) work primarily focuseson the traditional supervised learning setting in which models are unableto extract event mentions of new/unseen event types. Few-shot learninghas not been investigated in this area although it enables EC models toextend their operation to unobserved event types. To fill in this gap, inthis work, we investigate event classification under the few-shot learningsetting. We propose a novel training method for this problem that exten-sively exploit the support set during the training process of a few-shotlearning model. In particular, in addition to matching the query exam-ple with those in the support set for training, we seek to further matchthe examples within the support set themselves. This method providesmore training signals for the models and can be applied to every metric-learning-based few-shot learning methods. Our extensive experiments ontwo benchmark EC datasets show that the proposed method can improvethe best reported few-shot learning models by up to 10% on accuracyfor event classification


title: Insights from early mathematical models of 2019-nCoV acute respiratory disease (COVID-19) dynamics

abstract: In December 2019, a novel coronavirus (SARS-CoV-2) has been identified to cause acute respiratory disease in humans. An outbreak of this disease has been reported in mainland China with the city of Wuhan as the recognized epicenter. The disease has also been exported to other countries, including the Philippines, but the level of spread is still under control (as of 08 February 2020). To describe and predict the dynamics of the disease, several preliminary mathematical models are formulated by various international study groups. Here, the insights that can be drawn from these models are discussed, especially as inputs for designing strategies to control the epidemics. Proposed model-based strategies on how to prevent the spread of the disease in local setting, such as during large social gatherings, are also presented. The model shows that the exposure time is a significant factor in spreading the disease. With a basic reproduction number equal to 2, and 14-day infectious period, an infected person staying more than 9 hours in the event could infect other people. Assuming the exposure time is 18 hours, the model recommends that attendees of the social gathering should have a protection with more than 70 percent effectiveness.


title: Solution manifold and Its Statistical Applications

abstract: A solution manifold is the collection of points in a $d$-dimensional space satisfying a system of $s$ equations with $s<d$. Solution manifolds occur in several statistical problems including hypothesis testing, curved-exponential families, constrained mixture models, partial identifications, and nonparametric set estimation. We analyze solution manifolds both theoretically and algorithmically. In terms of theory, we derive five useful results: the smoothness theorem, the stability theorem (which implies the consistency of a plug-in estimator), the convergence of a gradient flow, the local center manifold theorem and the convergence of the gradient descent algorithm. To numerically approximate a solution manifold, we propose a Monte Carlo gradient descent algorithm. In the case of likelihood inference, we design a manifold constraint maximization procedure to find the maximum likelihood estimator on the manifold. We also develop a method to approximate a posterior distribution defined on a solution manifold.


title: Breaking limitation of quantum annealer in solving optimization problems under constraints

abstract: Quantum annealing is a generic solver for optimization problems that uses fictitious quantum fluctuation. The most groundbreaking progress in the research field of quantum annealing is its hardware implementation, i.e., the so-called quantum annealer, using artificial spins. However, the connectivity between the artificial spins is sparse and limited on a special network known as the chimera graph. Several embedding techniques have been proposed, but the number of logical spins, which represents the optimization problems to be solved, is drastically reduced. In particular, an optimization problem including fully or even partly connected spins suffers from low embeddable size on the chimera graph. In the present study, we propose an alternative approach to solve a large-scale optimization problem on the chimera graph via a well-known method in statistical mechanics called the Hubbard-Stratonovich transformation or its variants. The proposed method can be used to deal with a fully connected Ising model without embedding on the chimera graph and leads to nontrivial results of the optimization problem. We tested the proposed method with a number of partition problems involving solving linear equations and the traffic flow optimization problem in Sendai and Kyoto cities in Japan.


title: Depth Descent Synchronization in $\mathrm{SO}(D)$

abstract: We give robust recovery results for synchronization on the rotation group, $\mathrm{SO}(D)$. In particular, we consider an adversarial corruption setting, where a limited percentage of the observations are arbitrarily corrupted. We give a novel algorithm that exploits Tukey depth in the tangent space, which exactly recovers the underlying rotations up to an outlier percentage of $1/(D(D-1)+2)$. This corresponds to an outlier fraction of $1/4$ for $\mathrm{SO}(2)$ and $1/8$ for $\mathrm{SO}(3)$. In the case of $D=2$, we demonstrate that a variant of this algorithm converges linearly to the ground truth rotations. We finish by discussing this result in relation to a simpler nonconvex energy minimization framework based on least absolute deviations, which exhibits spurious fixed points.


title: Fog Computing on Constrained Devices: Paving the Way for the Future IoT

abstract: In the long term, the Internet of Things (IoT) is expected to become an integral part of people's daily lives. In light of this technological advancement, an ever-growing number of objects with limited hardware may become connected to the Internet. In this chapter, we explore the importance of these constrained devices as well as how we can use them in conjunction with fog computing to change the future of the IoT. First, we present an overview of the concepts of constrained devices, IoT, and fog and mist computing, and then we present a classification of applications according to the amount of resources they require (e.g., processing power and memory). After that, we tie in these topics with a discussion of what can be expected in a future where constrained devices and fog computing are used to push the IoT to new limits. Lastly, we discuss some challenges and opportunities that these technologies may bring.


title: Edge State Wave Functions from Momentum-Conserving Tunneling Spectroscopy

abstract: We perform momentum-conserving tunneling spectroscopy using a GaAs cleaved-edge overgrowth quantum wire to investigate adjacent quantum Hall edge states. We use the lowest five wire modes with their distinct wave functions to probe each edge state and apply magnetic fields to modify the wave functions and their overlap. This reveals an intricate and rich tunneling conductance fan structure which is succinctly different for each of the wire modes. We self-consistently solve the Poisson-Schr\"odinger equations to simulate the spectroscopy, reproducing the striking fans in great detail, thus confirming the calculations. Further, the model predicts hybridization between wire states and Landau levels, which is also confirmed experimentally. This establishes momentum-conserving tunneling spectroscopy as a powerful technique to probe edge state wave functions.


title: Low-loss two-dimensional plasmon modes in antimonene

abstract: The effects of spin-orbit (SOC) and electron-phonon coupling on the collective excitation of doped monolayer Sb$_2$ are investigated using density functional and many-body perturbation theories. The spin-orbit coupling is exclusively important for the monolayer Sb$_2$ and it leads to the reconstruction of the electronic band structure. In particular, plasmon modes of monolayer Sb$_2$ are quite sensitive to the SOC and are characterized by very low damping rates owing to small electron-phonon scatterings. Our results show plasmons in antimonene are significantly less damped compared to monolayer graphene when plasmon energies are $\hbar \omega> 0.2$ eV due to smaller plasmon-phonon coupling in the former material.


title: High-Field Magnetization and Magnetic Phase Diagram of Metamagnetic Shape Memory Alloys Ni50-xCoxMn31.5Ga18.5 (x = 9 and 9.7)

abstract: Magnetic phase diagrams of the metamagnetic shape memory alloys Ni50-xCoxMn31.5Ga18.5 (x = 9 and 9.7) were produced from high-field magnetization measurements up to 56 T. For both compounds, magnetic field induced martensitic transformations are observed at various temperatures below 300 K. Hysteresis of the field-induced transformation shows unconventional temperature dependence: it decreases with decreasing temperature after showing a peak. Magnetic susceptibility measurement, microscopy, and X-ray diffraction data suggest a model incorporating the magnetic anisotropy and Zeeman energy in two variants, which qualitatively explains the thermal and the magnetic field history dependence of the hysteresis in these alloys.


title: Predictive Power of Nearest Neighbors Algorithm under Random Perturbation

abstract: We consider a data corruption scenario in the classical $k$ Nearest Neighbors ($k$-NN) algorithm, that is, the testing data are randomly perturbed. Under such a scenario, the impact of corruption level on the asymptotic regret is carefully characterized. In particular, our theoretical analysis reveals a phase transition phenomenon that, when the corruption level $\omega$ is below a critical order (i.e., small-$\omega$ regime), the asymptotic regret remains the same; when it is beyond that order (i.e., large-$\omega$ regime), the asymptotic regret deteriorates polynomially. Surprisingly, we obtain a negative result that the classical noise-injection approach will not help improve the testing performance in the beginning stage of the large-$\omega$ regime, even in the level of the multiplicative constant of asymptotic regret. As a technical by-product, we prove that under different model assumptions, the pre-processed 1-NN proposed in \cite{xue2017achieving} will at most achieve a sub-optimal rate when the data dimension $d>4$ even if $k$ is chosen optimally in the pre-processing step.


title: Interactive Multi-User 3D Visual Analytics in Augmented Reality

abstract: This publication reports on a research project in which we set out to explore the advantages and disadvantages augmented reality (AR) technology has for visual data analytics. We developed a prototype of an AR data analytics application, which provides users with an interactive 3D interface, hand gesture-based controls and multi-user support for a shared experience, enabling multiple people to collaboratively visualize, analyze and manipulate data with high dimensional features in 3D space. Our software prototype, called DataCube, runs on the Microsoft HoloLens - one of the first true stand-alone AR headsets, through which users can see computer-generated images overlaid onto real-world objects in the user's physical environment. Using hand gestures, the users can select menu options, control the 3D data visualization with various filtering and visualization functions, and freely arrange the various menus and virtual displays in their environment. The shared multi-user experience allows all participating users to see and interact with the virtual environment, changes one user makes will become visible to the other users instantly. As users engage together they are not restricted from observing the physical world simultaneously and therefore they can also see non-verbal cues such as gesturing or facial reactions of other users in the physical environment. The main objective of this research project was to find out if AR interfaces and collaborative analysis can provide an effective solution for data analysis tasks, and our experience with our prototype system confirms this.


title: Symplectic invariants of semitoric systems and the inverse problem for quantum systems

abstract: Simple semitoric systems were classified about ten years ago in terms of a collection of invariants, essentially given by a convex polygon with some marked points corresponding to focus-focus singularities. Each marked point is endowed with labels which are symplectic invariants of the system. We will review the construction of these invariants, and explain how they have been generalized or applied in different contexts. One of these applications concerns quantum integrable systems and the corresponding inverse problem, which asks how much information of the associated classical system can be found in the spectrum. An approach to this problem has been to try to compute invariants in the spectrum. We will explain how this has been recently achieved for some of the invariants of semitoric systems, and discuss an open question in this direction.


title: Jet splitting measurements in Pb--Pb and pp collisions at $\sqrt{s}_{\mathrm{NN}} =$ 5.02 TeV with ALICE

abstract: Recent ALICE measurements of jet splittings in Pb--Pb and pp collisions at $\sqrt{s_{\mathrm{NN}}}$ = 5.02 TeV are reported. These measurements scan the phase space of jet emissions in search of medium-induced signals which are expected to emerge at different scales. These include effects such as multiple soft-radiation, single hard emissions, and color coherence. The Lund plane diagram is shown, including projections onto distributions of the splitting scale $k_{\mathrm{T}}$ in intervals of the splitting angle $R_{\mathrm{g}}$. Soft Drop grooming is applied to access hard splittings within the jet, enabling measurements of groomed substructure variables. These include the shared momentum fraction $z_{\mathrm{g}}$ between the two hardest subjets and the number of Soft Drop splittings $n_{\mathrm{SD}}$. The results in Pb--Pb collisions are compared to PYTHIA events embedded into a Pb--Pb background to separate out background from in-medium effects. Measurements of $z_{\mathrm{g}}$ and the normalized splitting angle $\theta_{\mathrm{g}}$ will also be shown in pp collisions at $\sqrt{s}$ = 5.02 TeV for different grooming settings.


title: Efficient Adaptive Experimental Design for Average Treatment Effect Estimation

abstract: The goal of many scientific experiments including A/B testing is to estimate the average treatment effect (ATE), which is defined as the difference between the expected outcomes of two or more treatments. In this paper, we consider a situation where an experimenter can assign a treatment to research subjects sequentially. In adaptive experimental design, the experimenter is allowed to change the probability of assigning a treatment using past observations for estimating the ATE efficiently. However, with this approach, it is difficult to apply a standard statistical method to construct an estimator because the observations are not independent and identically distributed. We thus propose an algorithm for efficient experiments with estimators constructed from dependent samples. We also introduce a sequential testing framework using the proposed estimator. To justify our proposed approach, we provide finite and infinite sample analyses. Finally, we experimentally show that the proposed algorithm exhibits preferable performance.


title: Optimal Epoch Stochastic Gradient Descent Ascent Methods for Min-Max Optimization

abstract: Epoch gradient descent method (a.k.a. Epoch-GD) proposed by Hazan and Kale (2011) was deemed a breakthrough for stochastic strongly convex minimization, which achieves the optimal convergence rate of $O(1/T)$ with $T$ iterative updates for the {\it objective gap}. However, its extension to solving stochastic min-max problems with strong convexity and strong concavity still remains open, and it is still unclear whether a fast rate of $O(1/T)$ for the {\it duality gap} is achievable for stochastic min-max optimization under strong convexity and strong concavity. Although some recent studies have proposed stochastic algorithms with fast convergence rates for min-max problems, they require additional assumptions about the problem, e.g., smoothness, bi-linear structure, etc. In this paper, we bridge this gap by providing a sharp analysis of epoch-wise stochastic gradient descent ascent method (referred to as Epoch-GDA) for solving strongly convex strongly concave (SCSC) min-max problems, without imposing any additional assumption about smoothness or the function's structure. To the best of our knowledge, our result is the first one that shows Epoch-GDA can achieve the optimal rate of $O(1/T)$ for the duality gap of general SCSC min-max problems. We emphasize that such generalization of Epoch-GD for strongly convex minimization problems to Epoch-GDA for SCSC min-max problems is non-trivial and requires novel technical analysis. Moreover, we notice that the key lemma can also be used for proving the convergence of Epoch-GDA for weakly-convex strongly-concave min-max problems, leading to a nearly optimal complexity without resorting to smoothness or other structural conditions.


title: One Planet or Two Planets? The Ultra-sensitive Extreme-magnification Microlensing Event KMT-2019-BLG-1953

abstract: We present the analysis of a very high-magnification ($A\sim 900$) microlensing event KMT-2019-BLG-1953. A single-lens single-source (1L1S) model appears to approximately delineate the observed light curve, but the residuals from the model exhibit small but obvious deviations in the peak region. A binary lens (2L1S) model with a mass ratio $q\sim 2\times 10^{-3}$ improves the fits by $\Delta\chi^2=181.8$, indicating that the lens possesses a planetary companion. From additional modeling by introducing an extra planetary lens component (3L1S model) and an extra source companion (2L2S model), it is found that the residuals from the 2L1S model further diminish, but claiming these interpretations is difficult due to the weak signals with $\Delta\chi^2=16.0$ and $13.5$ for the 3L1S and 2L2L models, respectively. From a Bayesian analysis, we estimate that the host of the planets has a mass of $M_{\rm host}=0.31^{+0.37}_{-0.17}~M_\odot$ and that the planetary system is located at a distance of $D_{\rm L}=7.04^{+1.10}_{-1.33}~{\rm kpc}$ toward the Galactic center. The mass of the securely detected planet is $M_{\rm p}=0.64^{+0.76}_{-0.35}~M_{\rm J}$. The signal of the potential second planet could have been confirmed if the peak of the light curve had been more densely observed by followup observations, and thus the event illustrates the need for intensive followup observations for very high-magnification events even in the current generation of high-cadence surveys.


title: Probing the top-Higgs boson FCNC couplings via the $h\to \gamma\gamma$ channel at the HE-LHC and FCC-hh

abstract: We investigate the sensitivity of future searches for the top-Higgs boson Flavour Changing Neutral Current (FCNC) couplings $tqh$~($q= u, c$) at the proposed High Energy Large Hadron Collider~(HE-LHC) and Future Circular Collider in hadron-hadron mode (FCC-hh). We perform a full simulation for two processes in the $h\to \gamma\gamma$ decay channel (where $h$ is the discovered Higgs state): single top quark FCNC production in association with the Higgs boson (plus a jet) and top quark pair production with FCNC decays $t\to qh$. All the relevant backgrounds are considered in a cut based analysis to obtain the limits on the Branching Ratios (BRs) of $t\to uh$ and $t\to ch$. It is shown that, at the HE-LHC with an integrated luminosity of 15 ab$^{-1}$ and at the FCC-hh with an integrated luminosity of 30 ab$^{-1}$, {{the BR($t\to uh$) (BR($t\to ch$)) can be probed, respectively, to $7.0~(8.5)\times 10^{-5}$ and $2.3~(3.0) \times 10^{-5}$ at the 95\% Confidence Level (CL) (assuming a 10\% systematic uncertainty on the background), which is almost two orders of magnitude better than the current 13 TeV LHC experimental results.


title: Polar state of 0.67BiFeO3-0.33BaTiO3 near the morphotropic phase boundary

abstract: The symmetric studies on the structure-property relationship of the unpoled and poled states of 0.67BiFeO3-0.33BaTiO3 (0.67BF-0.33BT) were conducted to understand the origin of the morphotropic phase boundary (MPB) in BF-BT. A typical relaxor-type dielectric anomaly was observed (Tf, ~627 K). The remnant polarization (Pr) and maximum value of electro-strain (Sm) increase clearly during heating (Pr, ~40 uC/cm2; Sm, 0.191 % under 40 kV/cm at 453 K). The first-cycle electro-strain loops indicate the difference in the polar state between 0.67BF-0.33BT and 0.94BiNaTiO3-0.06BaTiO3. Both the unpoled and poled samples have the similar frequency dispersion behaviors. Even in the poled samples, the transition between the ergodic relaxor state and ferroelectric-like state does not involve a clear dielectric anomaly. Analyses based on the Rietveld refinement of XRD patterns, bright-field images and selected-area electron diffractions (SAED) demonstrated that the formation of the long-range ferroelectric domains was difficult under the poling field.


title: Wide field-of-view study of the Eagle Nebula with the Fourier transform imaging spectrograph SITELLE at CFHT

abstract: We present the very first wide-field, 11 by 11 arcmin, optical spectral mapping of M16, one of the most famous star-forming regions in the Galaxy. The data were acquired with the new imaging Fourier transform spectrograph SITELLE mounted on the Canada-France-Hawaii Telescope (CFHT). We obtained three spectral cubes with R=10'000 (SN1 filter), 1500 (SN2 filter) and 600 (SN3 filter), centered on the Pillars of Creation and the HH216 flow, covering the main optical nebular emission lines: [O II]3726,29 (SN1), Hb, [O III]4959,5007 (SN2), [N II]6548,84, Ha, and [S II]6717,31 (SN3). We validate the performance, calibration, and data reduction of SITELLE, and analyze the structures in the large field-of-view in terms of their kinematics and nebular emission. We compared the SITELLE data to MUSE integral field observations and other spectroscopic and narrow-band imaging data to validate the performance of SITELLE. We computed gas-phase metallicities via the strong-line method, performed a pixel-by-pixel fit to the main emission lines to derive kinematics of the ionized gas, computed the mass-loss rate of the Eastern pillar (the Spire), and combined the SITELLE data with near-infrared narrow-band imaging to characterize the HH216 flow. The comparison with previously published fluxes demonstrates very good agreement. We disentangle the dependence of the gas-phase metallicities (derived via abundance-tracing line ratios) on the degree of ionization and obtain metallicities that are in excellent agreement with the literature. We confirm the bipolar structure of HH216, find evidence for episodic accretion from the source of the flow, and identify its likely driving source. We compute the mass-loss rate of the Spire pillar on the East side of the H II region and find excellent agreement with the correlation between the mass-loss rate and the ionizing photon flux from the nearby cluster NGC6611.


title: Self-supervised learning for audio-visual speaker diarization

abstract: Speaker diarization, which is to find the speech segments of specific speakers, has been widely used in human-centered applications such as video conferences or human-computer interaction systems. In this paper, we propose a self-supervised audio-video synchronization learning method to address the problem of speaker diarization without massive labeling effort. We improve the previous approaches by introducing two new loss functions: the dynamic triplet loss and the multinomial loss. We test them on a real-world human-computer interaction system and the results show our best model yields a remarkable gain of +8%F1-scoresas well as diarization error rate reduction. Finally, we introduce a new large scale audio-video corpus designed to fill the vacancy of audio-video datasets in Chinese.


title: Test of CP invariance in vector-boson fusion production of the Higgs boson in the $H\rightarrow\tau\tau$ channel in proton$-$proton collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

abstract: A test of CP invariance in Higgs boson production via vector-boson fusion is performed in the $H\rightarrow\tau\tau$ decay channel. This test uses the Optimal Observable method and is carried out using 36.1 $\mathrm{fb}^{-1}$ of $\sqrt{s}$ = 13 TeV proton$-$proton collision data collected by the ATLAS experiment at the LHC. Contributions from CP-violating interactions between the Higgs boson and electroweak gauge bosons are described by an effective field theory, in which the parameter $\tilde{d}$ governs the strength of CP violation. No sign of CP violation is observed in the distributions of the Optimal Observable, and $\tilde{d}$ is constrained to the interval $[-0.090, 0.035]$ at the 68% confidence level (CL), compared to an expected interval of $\tilde{d} \in [-0.035,0.033]$ based upon the Standard Model prediction. No constraints can be set on $\tilde{d}$ at 95% CL, while an expected 95% CL interval of $\tilde{d} \in [-0.21,0.15]$ for the Standard Model hypothesis was expected.