title: Neuromorphic Computing is Turing-Complete

abstract: Neuromorphic computing is a non-von Neumann computing paradigm that performs computation by emulating the human brain. Neuromorphic systems are extremely energy-efficient and known to consume thousands of times less power than CPUs and GPUs. They have the potential to drive critical use cases such as autonomous vehicles, edge computing and internet of things in the future. For this reason, they are sought to be an indispensable part of the future computing landscape. Neuromorphic systems are mainly used for spike-based machine learning applications, although there are some non-machine learning applications in graph theory, differential equations, and spike-based simulations. These applications suggest that neuromorphic computing might be capable of general-purpose computing. However, general-purpose computability of neuromorphic computing has not been established yet. In this work, we prove that neuromorphic computing is Turing-complete and therefore capable of general-purpose computing. Specifically, we present a model of neuromorphic computing, with just two neuron parameters (threshold and leak), and two synaptic parameters (weight and delay). We devise neuromorphic circuits for computing all the {\mu}-recursive functions (i.e., constant, successor and projection functions) and all the {\mu}-recursive operators (i.e., composition, primitive recursion and minimization operators). Given that the {\mu}-recursive functions and operators are precisely the ones that can be computed using a Turing machine, this work establishes the Turing-completeness of neuromorphic computing.


title: Unraveling the Unconventional Order of a High-Mobility Indacenodithiophene-Benzothiadiazole Copolymer

abstract: A new class of donor-acceptor (D-A) copolymers found to produce high charge carrier mobilities competitive with amorphous silicon ($> 1 cm^{2}V^{-1}s^{-1}$) exhibits the puzzling microstructure of substantial local order, however lacking long-range order and crystallinity previously deemed necessary for achieving high mobility. Here, we demonstrate the application of low-dose transmission electron microscopy to image and quantify the nanoscale and mesoscale organization of an archetypal D-A copolymer across areas comparable to electronic devices (~ $9 {\mu}m^{2}$). The local structure is spatially resolved by mapping the backbone (001) spacing reflection, revealing nanocrystallites of aligned polymer chains over nearly the entire film. Analysis of the nanoscale structure of its ordered domains suggests significant short- and medium-range order and preferential grain boundary orientations. Moreover, we provide insights into the rich, interconnected mesoscale organization of this new family of D-A copolymers by analysis of the local orientational spatial autocorrelations.


title: Neutron scattering study of the kagome metal Sc3Mn3Al7Si5

abstract: Sc3Mn3Al7Si5 is a rare example of a correlated metal in which the Mn moments form a kagome lattice. The absence of magnetic ordering to the lowest temperatures suggests that geometrical frustration of magnetic interactions may lead to strong magnetic fluctuations. We have performed inelastic neutron scattering measurements on Sc3Mn3Al7Si5, finding that phonon scattering dominates for energies from ~20 - 50 meV. These results are in good agreement with ab initio calculations of the phonon dispersions and densities of states, and as well reproduce the measured specific heat. A weak magnetic signal was detected at energies less than ~10 meV, present only at the lowest temperatures. The magnetic signal is broad and quasielastic, as expected for metallic paramagnets.


title: Freeze-in sterile neutrino dark matter in a class of U$(1)^\prime$ models with inverse seesaw

abstract: We study a class of general U$(1)^\prime$ models to explain the observed dark matter relic abundance and light neutrino masses. The model contains three right handed neutrinos and three gauge singlet Majorana fermions to generate the light neutrino mass via the inverse seesaw mechanism. We assign one pair of degenerate sterile neutrinos to be the dark matter candidate whose relic density is generated by the freeze-in mechanism. We consider different regimes of the masses of the dark matter particle and the ${Z^\prime}$ gauge boson. The production of the dark matter can occur at different reheating temperatures in various scenarios depending on the masses of the ${Z^\prime}$ boson and the dark matter candidate. We also note that if the mass of the sterile neutrino dark matter is $\gtrsim 1 \rm{MeV}$ and if the $Z^\prime$ is heavier than the dark matter, the decay of the dark matter candidate into positrons can explain the long standing puzzle of the galactic $511\rm{keV}$ line in the Milky Way center observed by the INTEGRAL satellite. We constrain the model parameters from the dark matter analysis, vacuum stability and the collider searches of heavy ${Z^\prime}$ at the LHC. For the case with light $Z^\prime$, we also compare how far the parameter space allowed from dark matter relic density can be probed by the future lifetime frontier experiments SHiP and FASERs in the special case of $U(1)_{B-L}$ model.


title: The Lowest Broad Alpha Cluster Resonances in $^{19}$F

abstract: There is a deep astrophysical interest in the structure of $^{19}$F states close to the alpha decay threshold. The nuclear structure of these states is important for understanding of the development of $\alpha$ clustering in the $^{20}$Ne region. Emergence of clustered states and generally states that favor coupling to reaction channels near the corresponding decay thresholds is currently of special interest in theoretical physics. Excitation function for $^{15}$N($\alpha$,$\alpha$) elastic scattering was measured by the TTIK method. These new data together with old, high energy resolution data, were analyzed using the R matrix approach. $^{19}$F nuclear structure was calculated using configuration interaction methods with the recently developed effective interaction Hamiltonian. The parameters of broad low spin $\ell = 0$ and 1 relative partial wave resonances close to the $\alpha$ decay threshold in $^{19}$F were identified. Detailed theoretical analysis was carried out identifying all states coupled to the $\ell = 0$ and 1 alpha cluster channels. Considering hierarchy of states with different harmonic oscillator shell excitations allows to evaluate coupling to the alpha channels with different number of nodes in the relative wave function and helps to explain the distribution of the clustering strength and emergence of broad clustering resonances. Comparison of clustering in $^{20}$Ne into $^{16}$O+$\alpha$ and consideration of spin-orbit splitting of the $^{15}$N+$\alpha$ channel provides additional evidence. Detailed analysis of new and old experimental data allows to identify a series of $\alpha$ clustering resonances in $^{19}$F and to assess the distribution of the clustering strength which is of importance to questions of astrophysics and for theoretical understanding of many-body physics and emergence of clustering in loosely bound or unstable nuclei.


title: On the maximum helium content of multiple populations in the globular cluster NGC6752

abstract: Multiple populations in globular clusters are usually explained by the formation of stars out of material with a chemical composition that is polluted to different degrees by the ejecta of short-lived, massive stars of various type. Among other things, these polluters differ by the amount of helium they spread in the surrounding medium. In this study we investigate whether the present-day photometric method used to infer the helium content of multiple populations indeed gives the true value or underestimates it by missing very He-rich, but rare stars. We focus on the specific case of NGC6752. We compute atmosphere models and synthetic spectra along isochrones produced for this cluster for a very broad range of He abundances covering the predictions of different pollution scenarios, including the extreme case of the fast-rotating massive star (FRMS) scenario. We calculate synthetic photometry in HST filters best suited to study the helium content. We subsequently build synthetic clusters with various distributions of stars. We finally determine the maximum helium mass fraction of these synthetic clusters using a method similar to that applied to observational data. We build toy models of clusters with various distributions of multiple populations and ensure that we are able to recover the input maximum Y. We then build synthetic clusters with the populations predicted by the FRMS scenario and find that while we slightly underestimate the maximum Y value, we are still able to detect stars much more He-rich than the current observed maximum Y. It is easier to determine the maximum Y on main sequence stars than on red giant branch stars, but qualitatively the results are unaffected by the sample choice. We show that in NGC6752 it is unlikely that stars more He-rich than the current observational limit of about 0.3 are present.


title: Survival probabilities of charmonia as a clue to measure transient magnetic fields

abstract: We investigate time evolution of $S$-wave charmonium populations under a time-dependent homogeneous magnetic field and evaluate survival probabilities of the low-lying charmonia to the goal of estimating the magnetic field strength at heavy-ion collisions. Our approach implements mixing between different spin eigenstates and transitions to radially excited states. We show that the survival probabilities can change even by an extremely short magnetic field. Furthermore, we find that the survival probabilities depend on the initial spin states. We propose the sum of the survival probabilities over spin partners as an observable insensitive to the initial states. We also find that the sum can be approximately given as a function of $\sigma B_0^2$ with a duration time $\sigma$ and the maximum strength of the magnetic field $B_0$.


title: A review of quarkonia under strong magnetic fields

abstract: We review the properties of quarkonia under strong magnetic fields. The main phenomena are (i) mixing between different spin eigenstates, (ii) quark Landau levels and deformation of wave function, (iii) modification of $\bar{Q}Q$ potential, and (iv) the motional Stark effect. For theoretical approaches, we review (i) constituent quark models, (ii) effective Lagrangians, (iii) QCD sum rules, and (iv) holographic approaches.


title: Pure hydrocarbon cycles in TMC-1: Discovery of ethynyl cyclopropenylidene, cyclopentadiene and indene

abstract: We report the detection for the first time in space of three new pure hydrocarbon cycles in TMC-1: c-C3HCCH (ethynyl cyclopropenylidene), c-C5H6 (cyclopentadiene) and c-C9H8 (indene). We derive a column density of 3.1e11 cm-2 for the former cycle and similar values, in the range (1-2)e13 cm-2, for the two latter molecules. This means that cyclopentadiene and indene, in spite of their large size, are exceptionally abundant, only a factor of five less abundant than the ubiquitous cyclic hydrocarbon c-C3H2. The high abundance found for these two hydrocarbon cycles, together with the high abundance previously found for the propargyl radical (CH2CCH) and other hydrocarbons like vinyl and allenyl acetylene (Agundez et al. 2021; Cernicharo et al. 2021a,b), start to allow us to quantify the abundant content of hydrocarbon rings in cold dark clouds and to identify the intermediate species that are probably behind the in situ bottom-up synthesis of aromatic cycles in these environments. While c-C3HCCH is most likely formed through the reaction between the radical CCH and c-C3H2, the high observed abundances of cyclopentadiene and indene are difficult to explain through currently proposed chemical mechanisms. Further studies are needed to identify how are five- and six-membered rings formed under the cold conditions of clouds like TMC-1.


title: Challenges of Adopting SAFe in the Banking Industry -- A Study Two Years after its Introduction

abstract: The Scaled Agile Framework (SAFe) is a framework for scaling agile methods in large organizations. We have found several experience reports and white papers describing SAFe adoptions in different banks, which indicates that SAFe is being used in the banking industry. However, there is a lack of academic publications on the topic, the banking industry is missing in the scientific reports analyzing SAFe transformations. To fill this gap, we present a study on the main challenges with a SAFe transformation at a large full-service bank. We identify the challenges in the bank under study and compare the findings with experience reports from other banks, as well as with research on SAFe transformations in other domains. Many of the challenges reported in this paper overlap with the generic SAFe challenges, including management and organization, education and training, culture and mindset, requirements engineering, quality assurance, and systems architecture. However, we also report some novel challenges specific to the banking domain, e.g., the risk of jeopardizing customer relations, stability, and trust of external stakeholders. This study validates several SAFe-related challenges reported in previous work in the banking context. It also brings up some novel challenges specific to the banking industry. Therefore, we believe our results are particularly useful to practitioners responsible for SAFe transformations at other banks.


title: Filter Distribution Templates in Convolutional Networks for Image Classification Tasks

abstract: Neural network designers have reached progressive accuracy by increasing models depth, introducing new layer types and discovering new combinations of layers. A common element in many architectures is the distribution of the number of filters in each layer. Neural network models keep a pattern design of increasing filters in deeper layers such as those in LeNet, VGG, ResNet, MobileNet and even in automatic discovered architectures such as NASNet. It remains unknown if this pyramidal distribution of filters is the best for different tasks and constrains. In this work we present a series of modifications in the distribution of filters in four popular neural network models and their effects in accuracy and resource consumption. Results show that by applying this approach, some models improve up to 8.9% in accuracy showing reductions in parameters up to 54%.


title: Thermal diffusivity and its lower bound in orthorhombic SnSe

abstract: The orthorhombic monochalcogenide SnSe has attracted much attention in recent years as a promising high-temperature thermoelectric material. We present a study of its thermal conductivity and specific heat of SnSe between 2~K and 300~K and quantify its anisotropic thermal diffusivity, $D$. For both crystallographic orientations, thermal diffusivity remains above the recently identified Planckian limit ($D > v_s^2 \tau_P$, where $v_s$ is the sound velocity and $\tau_P= \hbar/k_BT$) and its anisotropy in $D$ is set by the anisotropy of $v_s$. Comparison with cubic members of the IV-VI family leads to a consistent picture, where the diffusivity in all members of the family is set by the product of v$_s$, $\tau_P$ and the 'melting' velocity derived from the melting temperature.


title: Overcoming timestep limitations in boosted-frame Particle-In-Cell simulations of plasma-based acceleration

abstract: Explicit electromagnetic Particle-In-Cell (PIC) codes are typically limited by the Courant- Friedrichs-Lewy (CFL) condition, which implies that the timestep multiplied by the speed of light must be smaller than the smallest cell size. In the case of boosted-frame PIC simulations of plasma-based acceleration, this limitation can be a major hinderance as the cells are often very elongated along the longitudinal direction and the timestep is thus limited by the small, transverse cell size. This entails many small-timestep PIC iterations, and can limit the potential speed-up of the boosted-frame technique. Here, by using a CFL-free analytical spectral solver, and by mitigating additional numerical instabilities that arise at large timestep, we show that it is possible to overcome traditional limitations on the timestep and thereby realize the full potential of the boosted-frame technique over a much wider range of parameters.


title: Surface theory of a second-order topological insulator beyond the Dirac approximation

abstract: We study the surface states and chiral hinge states of a 3D second-order topological insulator in the presence of an external magnetic gauge field. Surfaces pierced by flux host Landau levels, while surfaces parallel to the applied field are not significantly affected. The chiral hinge modes mediate spectral flow between neighbouring surfaces. As the magnetic field strength is increased, the surface Landau quantization deviates from that of a massive Dirac cone. Quantitatively, the $n = 0$ Landau level falls inside the surface Dirac gap, and not at the gap edge. The $n \ne 0$ levels exhibit a further, qualitative discrepancy: while the massive Dirac cone is expected to produce pairs of levels ($\pm n$) which are symmetric around zero energy, the $n$ and $-n$ levels become asymmetric in our lattice model -- one of the pair may even be absent from the spectrum, or hybridized with the continuum. In order to resolve the issue, we extend the standard 2D massive Dirac surface theory, by including additional Hamiltonian terms at $\mathcal{O} (k^2)$. While these terms do not break particle-hole symmetry in the absence of magnetic field, they lead to the aforementioned Landau level asymmetry once the magnetic field is applied. We argue that similar $\mathcal{O}(k^2)$ correction terms are generically expected in lattice models containing gapped Dirac fermions, using the BHZ model of a 2D topological insulator as an example.


title: MAGMA: An Optimization Framework for Mapping Multiple DNNs on Multiple Accelerator Cores

abstract: As Deep Learning continues to drive a variety of applications in edge and cloud data centers, there is a growing trend towards building large accelerators with several sub-accelerator cores/chiplets. This work looks at the problem of supporting multi-tenancy on such accelerators. In particular, we focus on the problem of mapping jobs from several DNNs simultaneously on an accelerator. Given the extremely large search space, we formulate the search as an optimization problem and develop an optimization framework called M3E. In addition, we develop a specialized optimization algorithm called MAGMA with custom operators to enable structured sample-efficient exploration. We quantitatively compare MAGMA with several state-of-the-art methods, black-box optimization, and reinforcement learning methods across different accelerator settings (large/small accelerators) and different sub-accelerator configurations (homogeneous/heterogeneous), and observe MAGMA can consistently find better mappings.


title: RELIKE: Reionization Effective Likelihood from Planck 2018 Data

abstract: We release RELIKE (Reionization Effective Likelihood), a fast and accurate effective likelihood code based on the latest Planck 2018 data that allows one constrain any model for reionization between $6 < z < 30$ using five constraints from the CMB reionization principal components (PC). We tested the code on two example models which showed excellent agreement with sampling the exact Planck likelihoods using either a simple Gaussian PC likelihood or its full kernel density estimate. This code enables a fast and consistent means for combining Planck constraints with other reionization data sets, such as kinetic Sunyaev-Zeldovich effects, line-intensity mapping, luminosity function, star formation history, quasar spectra, etc, where the redshift dependence of the ionization history is important. Since the PC technique tests any reionization history in the given range, we also derive model-independent constraints for the total Thomson optical depth $\tau_{\rm PC} = 0.0619^{+0.0056}_{-0.0068}$ and its $15\le z \le 30$ high redshift component $\tau_{\rm PC}(15, 30) < 0.020 $ (95\% C.L.). The upper limits on the high-redshift optical depth is a factor of $\sim3$ larger than those reported in the Planck 2018 cosmological parameter paper using the FlexKnot method and we validate our results with a direct analysis of a two-step model which permits this small high-$z$ component.


title: Safety-Augmented Operation of Mobile Robots Using Variable Structure Control

abstract: The design process and complexity of existing safety controls are heavily determined by the geometrical properties of the environment, which affects the proof of convergence, design scalability, performance robustness, and numerical efficiency of the control. Hence, this paper proposes a variable structure control to isolate the environment's geometrical complexity from the control structure. A super-twisting algorithm is used to achieve accurate trajectory tracking and robust safety control. The safety control is designed solely based on distance measurement. First, a nominal safety model for obstacle avoidance is derived, where realistic system constraints are considered. The nominal model is well-suited for safety control design for obstacle avoidance, geofencing, and border patrol with analytically proven stability results. The safety control utilizes distance measurement to maintain a safe distance by compensating the robot's angular velocity. A supervisory logic is constructed to guarantee the overall stability and safety of the system. Operational safety and precision tracking are proven under parametric uncertainty and environmental uncertainty. The proposed design is modular with minimal tuning parameters, which reduces the computational burden and improves the control scalability. The effectiveness of the proposed method is verified against various case studies.


title: A Novel Transmission Policy for Intelligent Reflecting Surface Assisted Wireless Powered Sensor Networks

abstract: This paper proposes a novel transmission policy for an intelligent reflecting surface (IRS) assisted wireless powered sensor network (WPSN). An IRS is deployed to enhance the performance of wireless energy transfer (WET) and wireless information transfer (WIT) by intelligently adjusting phase shifts of each reflecting elements. To achieve its self-sustainability, the IRS needs to collect energy from energy station to support its control circuit operation. Our proposed policy for the considered WPSN is called IRS assisted harvest-then-transmit time switching, which is able to schedule the transmission time slots by switching between energy collection and energy reflection modes. We study the achievable sum throughput of the proposed transmission policy and investigate a joint design of the transmission time slots, the power allocation, as well as the discrete phase shifts of the WET and WIT. This formulates the problem as a mixed-integer non-linear program, which is NP-hard and non-convex. We first relax it to one with continuous phase shifts, and then propose a two-step approach and decompose the original problem into two sub-problems. We solve the first sub-problem with respect to the phase shifts of the WIT in terms of closed-form expression. For the second sub-problem, we consider a special case without the circuit power of each sensor node, the Lagrange dual method and the KKT conditions are applied to derive the optimal closed-form transmission time slots, power allocation, and phase shift of the WET. Then we generalise the case with the circuit power of each sensor node, which can be solved via employing a semi-definite programming relaxation. The optimal discrete phase shifts can be obtained by quantizing the continuous values. Numerical results demonstrate the effectiveness of the proposed policy and validate the beneficial role of the IRS in comparison to the benchmark schemes.


title: Verification and Synthesis of Control Barrier Functions

abstract: Control systems often must satisfy strict safety requirements over an extended operating lifetime. Control Barrier Functions (CBFs) are a promising recent approach to constructing simple and safe control policies. This paper proposes a framework for verifying that a CBF guarantees safety for all time and synthesizing CBFs with verifiable safety in polynomial control systems. Our approach is to show that safety of CBFs is equivalent to the non-existence of solutions to a family of polynomial equations, and then prove that this nonexistence is equivalent to a pair of sum-of-squares constraints via the Positivstellensatz of algebraic geometry. We develop this Positivstellensatz to verify CBFs, as well as generalization to high-degree systems and multiple CBF constraints. We then propose a set of heuristics for CBF synthesis, including a general alternating-descent heuristic, a specialized approach for compact safe regions, and an approach for convex unsafe regions. Our approach is illustrated on two numerical examples.


title: Modeling Managerial Search Behavior based on Simon's Concept of Satisficing

abstract: Computational models of managerial search often build on backward-looking search based on hill-climbing algorithms. Regardless of its prevalence, there is some evidence that this family of algorithms does not universally represent managers' search behavior. Against this background, the paper proposes an alternative algorithm that captures key elements of Simon's concept of satisficing which received considerable support in behavioral experiments. The paper contrasts the satisficing-based algorithm to two variants of hill-climbing search in an agent-based model of a simple decision-making organization. The model builds on the framework of NK fitness landscapes which allows controlling for the complexity of the decision problem to be solved. The results suggest that the model's behavior may remarkably differ depending on whether satisficing or hill-climbing serves as an algorithmic representation for decision-makers' search. Moreover, with the satisficing algorithm, results indicate oscillating aspiration levels, even to the negative, and intense - and potentially destabilizing - search activities when intra-organizational complexity increases. Findings may shed some new light on prior computational models of decision-making in organizations and point to avenues for future research.


title: On best approximations in Banach spaces from the perspective of orthogonality

abstract: We study best approximations in Banach spaces via Birkhoff-James orthogonality of functionals. To exhibit the usefulness of Birkhoff-James orthogonality techniques in the study of best approximation problems, some algorithms and distance formulae are presented. As an application of our study, we obtain some crucial inequalities, which also strengthen the classical H\"{o}lder's inequality. The relevance of the algorithms and the inequalities are discussed through concrete examples.


title: Optimal relaxation of bump-like solutions of the one-dimensional Cahn-Hilliard equation

abstract: In this paper we derive optimal relaxation rates for the Cahn-Hilliard equation on the one-dimensional torus and the line. We consider initial conditions with a finite (but not small) $L^1$-distance to an appropriately defined bump. The result extends the relaxation method developed previously for a single transition layer (the ``kink'') to the case of two transition layers (the ``bump''). As in the previous work, the tools include Nash-type inequalities, duality arguments, and Schauder estimates. For both the kink and the bump, the energy gap is translation invariant and its decay alone cannot specify to which member of the family of minimizers the solution converges. Whereas in the case of the kink, the conserved quantity singles out the longtime limit, in the case of a bump, a new argument is needed. On the torus, we quantify the (initially algebraic and ultimately exponential) convergence to the bump that is the longtime limit; on the line, the bump-like states are merely metastable and we quantify the initial algebraic relaxation behavior.


title: Unlocking capacities of viral genomics for the COVID-19 pandemic response

abstract: More than any other infectious disease epidemic, the COVID-19 pandemic has been characterized by the generation of large volumes of viral genomic data at an incredible pace due to recent advances in high-throughput sequencing technologies, the rapid global spread of SARS-CoV-2, and its persistent threat to public health. However, distinguishing the most epidemiologically relevant information encoded in these vast amounts of data requires substantial effort across the research and public health communities. Studies of SARS-CoV-2 genomes have been critical in tracking the spread of variants and understanding its epidemic dynamics, and may prove crucial for controlling future epidemics and alleviating significant public health burdens. Together, genomic data and bioinformatics methods enable broad-scale investigations of the spread of SARS-CoV-2 at the local, national, and global scales and allow researchers the ability to efficiently track the emergence of novel variants, reconstruct epidemic dynamics, and provide important insights into drug and vaccine development and disease control. Here, we discuss the tremendous opportunities that genomics offers to unlock the effective use of SARS-CoV-2 genomic data for efficient public health surveillance and guiding timely responses to COVID-19.


title: EmergencyNet: Efficient Aerial Image Classification for Drone-Based Emergency Monitoring Using Atrous Convolutional Feature Fusion

abstract: Deep learning-based algorithms can provide state-of-the-art accuracy for remote sensing technologies such as unmanned aerial vehicles (UAVs)/drones, potentially enhancing their remote sensing capabilities for many emergency response and disaster management applications. In particular, UAVs equipped with camera sensors can operating in remote and difficult to access disaster-stricken areas, analyze the image and alert in the presence of various calamities such as collapsed buildings, flood, or fire in order to faster mitigate their effects on the environment and on human population. However, the integration of deep learning introduces heavy computational requirements, preventing the deployment of such deep neural networks in many scenarios that impose low-latency constraints on inference, in order to make mission-critical decisions in real time. To this end, this article focuses on the efficient aerial image classification from on-board a UAV for emergency response/monitoring applications. Specifically, a dedicated Aerial Image Database for Emergency Response applications is introduced and a comparative analysis of existing approaches is performed. Through this analysis a lightweight convolutional neural network architecture is proposed, referred to as EmergencyNet, based on atrous convolutions to process multiresolution features and capable of running efficiently on low-power embedded platforms achieving upto 20x higher performance compared to existing models with minimal memory requirements with less than 1% accuracy drop compared to state-of-the-art models.


title: Interaction-GCN: A Graph Convolutional Network based framework for social interaction recognition in egocentric videos

abstract: In this paper we propose a new framework to categorize social interactions in egocentric videos, we named InteractionGCN. Our method extracts patterns of relational and non-relational cues at the frame level and uses them to build a relational graph from which the interactional context at the frame level is estimated via a Graph Convolutional Network based approach. Then it propagates this context over time, together with first-person motion information, through a Gated Recurrent Unit architecture. Ablation studies and experimental evaluation on two publicly available datasets validate the proposed approach and establish state of the art results.


title: BayesSUR: An R package for high-dimensional multivariate Bayesian variable and covariance selection in linear regression

abstract: In molecular biology, advances in high-throughput technologies have made it possible to study complex multivariate phenotypes and their simultaneous associations with high-dimensional genomic and other omics data, a problem that can be studied with high-dimensional multi-response regression, where the response variables are potentially highly correlated. To this purpose, we recently introduced several multivariate Bayesian variable and covariance selection models, e.g., Bayesian estimation methods for sparse seemingly unrelated regression for variable and covariance selection. Several variable selection priors have been implemented in this context, in particular the hotspot detection prior for latent variable inclusion indicators, which results in sparse variable selection for associations between predictors and multiple phenotypes. We also propose an alternative, which uses a Markov random field (MRF) prior for incorporating prior knowledge about the dependence structure of the inclusion indicators. Inference of Bayesian seemingly unrelated regression (SUR) by Markov chain Monte Carlo methods is made computationally feasible by factorisation of the covariance matrix amongst the response variables. In this paper we present BayesSUR, an R package, which allows the user to easily specify and run a range of different Bayesian SUR models, which have been implemented in C++ for computational efficiency. The R package allows the specification of the models in a modular way, where the user chooses the priors for variable selection and for covariance selection separately. We demonstrate the performance of sparse SUR models with the hotspot prior and spike-and-slab MRF prior on synthetic and real data sets representing eQTL or mQTL studies and in vitro anti-cancer drug screening studies as examples for typical applications.


title: A variant of the discrete Burgers equation derived from the correlated random walk and its ultradiscretization

abstract: In this paper, we show that a variant of the discrete Burgers equation can be obtained through the Cole--Hopf transformation to a generalized discrete diffusion equation corresponding to the correlated random walk, which is also known as a generalization of the well known random walk. By applying the technique called ultradiscretization, we obtain the generalized ultradiscrete diffusion equation, the ultradiscrete Cole--Hopf transformation and a variant of the ultradiscrete Burgers equation. Moreover, we show that the resulting ultradiscrete Burgers equation yields cellular automata which can be interpreted as a traffic flow model.


title: Deformation and orientational order of chiral membranes with free edges

abstract: Motivated by experiments on colloidal membranes composed of chiral rod-like viruses, we use Monte Carlo methods to determine the phase diagram for the liquid crystalline order of the rods and the membrane shape. We generalize the Lebwohl-Lasher model for a nematic with a chiral coupling to a curved surface with edge tension and a resistance to bending, and include an energy cost for tilting of the rods relative to the local membrane normal. The membrane is represented by a triangular mesh of hard beads joined by bonds, where each bead is decorated by a director. The beads can move, the bonds can reconnect and the directors can rotate at each Monte Carlo step. When the cost of tilt is small, the membrane tends to be flat, with the rods only twisting near the edge for low chiral coupling, and remaining parallel to the normal in the interior of the membrane. At high chiral coupling, the rods twist everywhere, forming a cholesteric state. When the cost of tilt is large, the emergence of the cholesteric state at high values of the chiral coupling is accompanied by the bending of the membrane into a saddle shape. Increasing the edge tension tends to flatten the membrane. These results illustrate the geometric frustration arising from the inability of a surface normal to have twist.


title: Channel Models and Coding Solutions for 1S1R Crossbar Resistive Memory with High Line Resistance

abstract: Crossbar resistive memory with the 1 Selector 1 Resistor (1S1R) structure is attractive for nonvolatile, high-density, and low-latency storage-class memory applications. As technology scales down to the single-nm regime, the increasing resistivity of wordline/bitline becomes a limiting factor to device reliability. This paper presents write/read communication channels while considering the line resistance and device variabilities by statistically relating the degraded write/read margins and the channel parameters. Binary asymmetric channel (BAC) models are proposed for the write/read operations. Simulations based on these models suggest that the bit-error rate of devices are highly non-uniform across the memory array. These models provide quantitative tools for evaluating the trade-offs between memory reliability and design parameters, such as array size, technology nodes, and aspect ratio, and also for designing coding-theoretic solutions that would be most effective for crossbar memory. Method for optimizing the read threshold is proposed to reduce the raw bit-error rate (RBER). We propose two schemes for efficient channel coding based on Bose-Chaudhuri-Hocquenghem (BCH) codes. An interleaved coding scheme is proposed to mitigate the non-uniformity of reliability and a location dependent coding framework is proposed to leverage this non-uniformity. Both of our proposed coding schemes effectively reduce the undetected bit-error rate (UBER).


title: Simplified Kalman filter for online rating: one-fits-all approach

abstract: In this work, we deal with the problem of rating in sports, where the skills of the players/teams are inferred from the observed outcomes of the games. Our focus is on the online rating algorithms which estimate the skills after each new game by exploiting the probabilistic models of the relationship between the skills and the game outcome. We propose a Bayesian approach which may be seen as an approximate Kalman filter and which is generic in the sense that it can be used with any skills-outcome model and can be applied in the individual -- as well as in the group-sports. We show how the well-know algorithms (such as the Elo, the Glicko, and the TrueSkill algorithms) may be seen as instances of the one-fits-all approach we propose. In order to clarify the conditions under which the gains of the Bayesian approach over the simpler solutions can actually materialize, we critically compare the known and the new algorithms by means of numerical examples using the synthetic as well as the empirical data.


title: Cosmology from pregeometry

abstract: We discuss cosmological solutions for a diffeomorphism invariant gauge theory of the non-compact Lorentz group $SO(1,3)$. Besides the gauge bosons our model of pregeometry contains a vector field in the vector representation of $SO(1,3)$ and a scalar singlet. General relativity and variable gravity emerge as effective theories for large distances and times in Planck units. We propose an approximation to the effective action with up to two derivatives. For a suitable range of parameters the universe approaches for large times stable Minkowski space. For late cosmology the model predicts dynamical dark energy and provides for a candidate for dark matter. Early cosmology is characterized by an inflationary epoch. The beginning of the universe in the infinite past is great emptiness, corresponding to an ultraviolet fixed point with the associated quantum scale symmetry. The beginning universe is a vacuum state with vanishing expectation values and finite non-vanishing correlation functions for the fluctuations of all fields. There is no physical big bang singularity.


title: Algorithmic Factors Influencing Bias in Machine Learning

abstract: It is fair to say that many of the prominent examples of bias in Machine Learning (ML) arise from bias that is there in the training data. In fact, some would argue that supervised ML algorithms cannot be biased, they reflect the data on which they are trained. In this paper we demonstrate how ML algorithms can misrepresent the training data through underestimation. We show how irreducible error, regularization and feature and class imbalance can contribute to this underestimation. The paper concludes with a demonstration of how the careful management of synthetic counterfactuals can ameliorate the impact of this underestimation bias.


title: Asymptotic scaling and universality for skew products with factors in SL(2,R)

abstract: We consider skew-product maps over circle rotations $x\mapsto x+\alpha$ (mod 1) with factors that take values in SL(2,R). This includes maps of almost Mathieu type. In numerical experiments, with $\alpha$ the inverse golden mean, Fibonacci iterates of maps from the almost Mathieu family exhibit asymptotic scaling behavior that is reminiscent of critical phase transitions. In a restricted setup that is characterized by a symmetry, we prove that critical behavior indeed occurs and is universal in an open neighborhood of the almost Mathieu family. This behavior is governed by a periodic orbit of a renormalization transformation. An extension of this transformation is shown to have a second periodic orbit as well, and we present some evidence that this orbit attracts supercritical almost Mathieu maps.


title: Reference based multiple imputation -- what is the right variance and how to estimate it

abstract: Reference based multiple imputation methods have become popular for handling missing data in randomised clinical trials. Rubin's variance estimator is well known to be biased compared to the reference based imputation estimator's true repeated sampling variance. Somewhat surprisingly given the increasingly popularity of these methods, there has been relatively little debate in the literature as to whether Rubin's variance estimator or alternative (smaller) variance estimators targeting the repeated sampling variance are more appropriate. We review the arguments made on both sides of this debate, and conclude that the repeated sampling variance is more appropriate. We review different approaches for estimating the frequentist variance, and suggest a recent proposal for combining bootstrapping with multiple imputation as a widely applicable general solution. At the same time, in light of the consequences of reference based assumptions for frequentist variance, we believe further scrutiny of these methods is warranted to determine whether the the strength of their assumptions are generally justifiable.


title: Functional-renormalization-group approach to strongly coupled Bose-Fermi mixtures in two dimensions

abstract: We study theoretically the phase diagram of strongly coupled two-dimensional Bose-Fermi mixtures interacting with attractive short-range potentials as a function of the particle densities. We focus on the limit where the size of the bound state between a boson and a fermion is small compared to the average interboson separation and develop a functional-renormalization-group approach that accounts for the bound-state physics arising from the extended Fr\"{o}hlich Hamiltonian. By including three-body correlations we are able to reproduce the polaron-to-molecule transition in two-dimensional Fermi gases in the extreme limit of vanishing boson density. We predict frequency- and momentum-resolved spectral functions and study the impact of three-body correlations on quasiparticle properties. At finite boson density, we find that when the bound-state energy exceeds the Fermi energy by a critical value, the fermions and bosons can form a fermionic composite with a well-defined Fermi surface. These composites constitute a Fermi sea of dressed Feshbach molecules in the case of ultracold atoms while in the case of atomically thin semiconductors a trion liquid emerges. As the boson density is increased further, the effective energy gap of the composites decreases, leading to a transition into a strongly correlated phase where polarons are hybridized with molecular degrees of freedom. We highlight the universal connection between two-dimensional semiconductors and ultracold atoms and we discuss perspectives for further exploring the rich structure of strongly coupled Bose-Fermi mixtures in these complementary systems.

categories: cond-mat.quant-gas cond-mat.mes-hall cond-mat.stat-mech cond-mat.str-el physics.atom-ph


title: Federated Identity Management (FIdM) Systems Limitation And Solutions

abstract: Efficient identity management system has become one of the fundamental requirements for ensuring safe, secure, and transparent use of identifiable information and attributes. FIdM allows users to distribute their identity information across security domains which increase the portability of their digital identities. However, it also raises new architectural challenges and significant security and privacy issues that need to be mitigated. In this paper, we presented the limitations and risks in Federated Identity Management system and discuss the results and proposed solutions.


title: Pebble transducers with unary output

abstract: Boja\'nczyk recently initiated an intensive study of deterministic pebble transducers, which are two-way automata that can drop marks (named "pebbles") on their input word, and produce an output word. They describe functions from words to words. Two natural restrictions of this definition have been investigated: marble transducers by Dou\'eneau-Tabot et al., and comparison-free pebble transducers (that we rename here "blind transducers") by Nguy\^en et al. Here, we study the decidability of membership problems between the classes of functions computed by pebble, marble and blind transducers that produce a unary output. First, we show that pebble and marble transducers have the same expressive power when the outputs are unary (which is false over non-unary outputs). Then, we characterize 1-pebble transducers with unary output that describe a function computable by a blind transducer, and show that the membership problem is decidable. These results can be interpreted in terms of automated simplification of programs.


title: Bias field correction of MPRAGE by an external reference -- The poor man's MP2RAGE

abstract: Purpose: To implement and evaluate a sequential approach to obtain semi-quantitative T1-weighted MPRAGE images, unbiased by B1 inhomogeneities at 7T. Methods: In the reference gradient echo used for normalization of the MPRAGE image, flip angle (aGE) and acquisition voxel size (Vref) was varied to optimize tissue contrast and acquisition time (Tacq). The finalized protocol was implemented at three different resolutions and the reproducibility was evaluated. Maps of T1 were derived based on the normalized MPRAGE through forward signal modelling. Results: A good compromise between tissue contrast and SNR was reached at aGE=3{\deg}. A reduction of the reference GE Tacq by a factor of 4, at the cost of negligible bias, was obtained by increasing Vref with a factor of 8 relative the MPRAGE resolution. The coefficient-of-variation in segmented WM was 9+/-5% after normalization, compared to 24+/-12% before. The T1 maps showed no obvious bias and had reasonable values with regard to literature, especially after optional B1 correction through separate flip angle mapping. Conclusion: A non-interleaved acquisition for normalization of MPRAGE offers a simple alternative to MP2RAGE to obtain semi-quantitative purely T1-weighted images. These images can be converted to T1 maps analogously to the established MP2RAGE approach. Scan time can be reduced by increasing Vref which has a miniscule effect on image quality.


title: The G\"odel Fibration

abstract: We introduce the notion of a G\"odel fibration, which is a fibration categorically embodying both the logical principle of traditional Skolemization (we can exchange the order of quantifiers paying the price of a functional) and the existence of a prenex normal form presentation for every logical formula. Building up from Hofstra's earlier fibrational characterization of the de Paiva's categorical Dialectica construction, we show that a fibration is an instance of the Dialectica construction if and only if it is a G\"odel fibration. This result establishes an internal presentation of the dialectica construction. Then we provide a deep structural analysis of the Dialectica construction producing a full description of which categorical structure behaves well with respect to this construction, focusing on (weak) finite products and coproducts. We conclude describing the applications we envisage for this generalized fibrational version of the Dialectica construction.


title: Low-energy optical properties of the non-magnetic kagome metal CsV$_3$Sb$_5$

abstract: Temperature-dependent reflectivity measurements on the kagome metal CsV$_3$Sb$_5$ in a broad frequency range of $50-20000$ cm$^{-1}$ down to $T$=10 K are reported. The charge-density wave (CDW) formed below $T_{\rm CDW}$ = 94 K manifests itself in a prominent spectral-weight transfer from low to higher energy regions. The CDW gap of 60-75 meV is observed at the lowest temperature and shows significant deviations from an isotropic BCS-type mean-field behavior. Absorption peaks appear at frequencies as low as 200 cm$^{-1}$ and can be identified with interband transitions according to density-functional calculations. The change in the interband absorption compared to KV$_3$Sb$_5$ reflects the inversion of band saddle points between the K and Cs compounds. Additionally, a broader and strongly temperature-dependent absorption feature is observed below 1000 cm$^{-1}$ and assigned to a displaced Drude peak. It reflects localization effects on charge carriers.


title: Measuring dependence between random vectors via optimal transport

abstract: To quantify the dependence between two random vectors of possibly different dimensions, we propose to rely on the properties of the 2-Wasserstein distance. We first propose two coefficients that are based on the Wasserstein distance between the actual distribution and a reference distribution with independent components. The coefficients are normalized to take values between 0 and 1, where 1 represents the maximal amount of dependence possible given the two multivariate margins. We then make a quasi-Gaussian assumption that yields two additional coefficients rooted in the same ideas as the first two. These different coefficients are more amenable for distributional results and admit attractive formulas in terms of the joint covariance or correlation matrix. Furthermore, maximal dependence is proved to occur at the covariance matrix with minimal von Neumann entropy given the covariance matrices of the two multivariate margins. This result also helps us revisit the RV coefficient by proposing a sharper normalisation. The two coefficients based on the quasi-Gaussian approach can be estimated easily via the empirical covariance matrix. The estimators are asymptotically normal and their asymptotic variances are explicit functions of the covariance matrix, which can thus be estimated consistently too. The results extend to the Gaussian copula case, in which case the estimators are rank-based. The results are illustrated through theoretical examples, Monte Carlo simulations, and a case study involving electroencephalography data.


title: Navier-Stokes Flow past a Rigid Body that Moves by Time-Periodic Motion

abstract: We study existence, uniqueness and asymptotic spatial behavior of time-periodic strong solutions to the Navier-Stokes equations in the exterior of a rigid body, $\mathscr B$, moving by time-periodic motion of given period $T$, when the data are sufficiently regular and small. Our contribution improves all previous ones in several directions. For example, we allow both translational, $\bfxi$, and angular, $\bfomega$, velocities of $\mathscr B$ to depend on time, and do not impose any restriction on the period $T$ nor on the averaged velocity, $\bar{\bfxi}$, of $\mathscr B$. If $\bfxi\not\equiv\0$ we assume that $\bfxi$ and $\bfomega$ are both parallel to a constant direction, while no further assumption is needed if $\bfxi\equiv\0$. We also furnish the spatial asymptotic behavior of the velocity field, $\bfu$, associated to such solutions. In particular, if $\mathscr B$ has a net motion characterized by $\bar{\bfxi}\neq\0$, we then show that, at large distances from $\mathscr B$, $\bfu$ manifests a wake-like behavior in the direction $-\bar{\bfxi}$, entirely similar to that of the velocity field of the steady-state flow occurring when $\mathscr B$ moves with velocity $\bar{\bfxi}$.


title: A Temporal Logic for Asynchronous Hyperproperties

abstract: Hyperproperties are properties of computational systems that require more than one trace to evaluate, e.g., many information-flow security and concurrency requirements. Where a trace property defines a set of traces, a hyperproperty defines a set of sets of traces. The temporal logics HyperLTL and HyperCTL* have been proposed to express hyperproperties. However, their semantics are synchronous in the sense that all traces proceed at the same speed and are evaluated at the same position. This precludes the use of these logics to analyze systems whose traces can proceed at different speeds and allow that different traces take stuttering steps independently. To solve this problem in this paper, we propose an asynchronous variant of HyperLTL. On the negative side, we show that the model-checking problem for this variant is undecidable. On the positive side, we identify a decidable fragment which covers a rich set of formulas with practical applications. We also propose two model-checking algorithms that reduce our problem to the HyperLTL model-checking problem in the synchronous semantics.


title: In-plane critical magnetic fields in magic-angle twisted trilayer graphene

abstract: It has recently been shown that superconductivity in magic-angle twisted trilayer graphene survives to in-plane magnetic fields that are well in excess of the Pauli limit, and much stronger than the in-plane critical magnetic fields of magic-angle twisted bilayer graphene. The difference is surprising because twisted bilayers and trilayers both support the magic-angle flat bands thought to be the fountainhead of twisted graphene superconductivity. We show here that the difference in critical magnetic fields can be traced to a $\mathcal{C}_2 \mathcal{M}_{h}$ symmetry in trilayers that survives in-plane magnetic fields, and also relative displacements between top and bottom layers that are not under experimental control at present. An gate electric field breaks the $\mathcal{C}_2 \mathcal{M}_{h}$ symmetry and therefore limits the in-plane critical magnetic field.


title: Modelling the non-local thermodynamic equilibrium spectra of silylene (SiH2)

abstract: This paper sets out a robust methodology for modelling spectra of polyatomic molecules produced in reactive or dissociative environments, with vibrational populations outside local thermal equilibrium (LTE). The methodology is based on accurate, extensive ro-vibrational line lists containing transitions with high vibrational excitations and relies on the detailed ro-vibrational assignments. The developed methodology is applied to model non-LTE IR and visible spectra of silylene (SiH$_2$) produced in a decomposition of disilane (Si$_2$H$_6$), a reaction of technological importance. Two approaches for non-LTE vibrational populations of the product SiH$_2$ are introduced: a simplistic 1D approach based on the Harmonic approximation and a full 3D model incorporating accurate vibrational wavefunctions of SiH$_2$ computed variationally with the TROVE (Theoretical ROVibrational Energy) program. We show how their non-LTE spectral signatures can be used to trace different reaction channels of molecular dissociations.


title: Analysis of Legal Documents via Non-negative Matrix Factorization Methods

abstract: The California Innocence Project (CIP), a clinical law school program aiming to free wrongfully convicted prisoners, evaluates thousands of mails containing new requests for assistance and corresponding case files. Processing and interpreting this large amount of information presents a significant challenge for CIP officials, which can be successfully aided by topic modeling techniques.In this paper, we apply Non-negative Matrix Factorization (NMF) method and implement various offshoots of it to the important and previously unstudied data set compiled by CIP. We identify underlying topics of existing case files and classify request files by crime type and case status (decision type). The results uncover the semantic structure of current case files and can provide CIP officials with a general understanding of newly received case files before further examinations. We also provide an exposition of popular variants of NMF with their experimental results and discuss the benefits and drawbacks of each variant through the real-world application.


title: Reducing Risk and Uncertainty of Deep Neural Networks on Diagnosing COVID-19 Infection

abstract: Effective and reliable screening of patients via Computer-Aided Diagnosis can play a crucial part in the battle against COVID-19. Most of the existing works focus on developing sophisticated methods yielding high detection performance, yet not addressing the issue of predictive uncertainty. In this work, we introduce uncertainty estimation to detect confusing cases for expert referral to address the unreliability of state-of-the-art (SOTA) DNNs on COVID-19 detection. To the best of our knowledge, we are the first to address this issue on the COVID-19 detection problem. In this work, we investigate a number of SOTA uncertainty estimation methods on publicly available COVID dataset and present our experimental findings. In collaboration with medical professionals, we further validate the results to ensure the viability of the best performing method in clinical practice.


title: Measurement-Robust Control Barrier Functions: Certainty in Safety with Uncertainty in State

abstract: The increasing complexity of modern robotic systems and the environments they operate in necessitates the formal consideration of safety in the presence of imperfect measurements. In this paper we propose a rigorous framework for safety-critical control of systems with erroneous state estimates. We develop this framework by leveraging Control Barrier Functions (CBFs) and unifying the method of Backup Sets for synthesizing control invariant sets with robustness requirements -- the end result is the synthesis of Measurement-Robust Control Barrier Functions (MR-CBFs). This provides theoretical guarantees on safe behavior in the presence of imperfect measurements and improved robustness over standard CBF approaches. We demonstrate the efficacy of this framework both in simulation and experimentally on a Segway platform using an onboard stereo-vision camera for state estimation.


title: Wave-fronts of gravitational waves partially trapped in ultra-compact stars

abstract: We dedicate this work to Dr Omer Blaes, professor of physics at UCSB, on the occasion of his sixtieth birthday. We have been collaborating now and then with Dr. Blaes on problems involving oscillations, waves and stability. Happy birthday, Omer. Enjoy the analytic treatment of damping of the gravitational waves trapped inside ultra compact stars and its possible connection to Quantum Gravity in the context of the LIGO-Virgo efforts in accurately measuring ringdowns and echoes.


title: Randomized Histogram Matching: A Simple Augmentation for Unsupervised Domain Adaptation in Overhead Imagery

abstract: Modern deep neural networks (DNNs) are highly accurate on many recognition tasks for overhead (e.g., satellite) imagery. However, visual domain shifts (e.g., statistical changes due to geography, sensor, or atmospheric conditions) remain a challenge, causing the accuracy of DNNs to degrade substantially and unpredictably when testing on new sets of imagery. In this work, we model domain shifts caused by variations in imaging hardware, lighting, and other conditions as non-linear pixel-wise transformations, and we perform a systematic study indicating that modern DNNs can become largely robust to these types of transformations, if provided with appropriate training data augmentation. In general, however, we do not know the transformation between two sets of imagery. To overcome this, we propose a fast real-time unsupervised training augmentation technique, termed randomized histogram matching (RHM). We conduct experiments with two large benchmark datasets for building segmentation and find that despite its simplicity, RHM consistently yields similar or superior performance compared to state-of-the-art unsupervised domain adaptation approaches, while being significantly simpler and more computationally efficient. RHM also offers substantially better performance than other comparably simple approaches that are widely used for overhead imagery.


title: A Study of the Mathematics of Deep Learning

abstract: "Deep Learning"/"Deep Neural Nets" is a technological marvel that is now increasingly deployed at the cutting-edge of artificial intelligence tasks. This dramatic success of deep learning in the last few years has been hinged on an enormous amount of heuristics and it has turned out to be a serious mathematical challenge to be able to rigorously explain them. In this thesis, submitted to the Department of Applied Mathematics and Statistics, Johns Hopkins University we take several steps towards building strong theoretical foundations for these new paradigms of deep-learning. In chapter 2 we show new circuit complexity theorems for deep neural functions and prove classification theorems about these function spaces which in turn lead to exact algorithms for empirical risk minimization for depth 2 ReLU nets. We also motivate a measure of complexity of neural functions to constructively establish the existence of high-complexity neural functions. In chapter 3 we give the first algorithm which can train a ReLU gate in the realizable setting in linear time in an almost distribution free set up. In chapter 4 we give rigorous proofs towards explaining the phenomenon of autoencoders being able to do sparse-coding. In chapter 5 we give the first-of-its-kind proofs of convergence for stochastic and deterministic versions of the widely used adaptive gradient deep-learning algorithms, RMSProp and ADAM. This chapter also includes a detailed empirical study on autoencoders of the hyper-parameter values at which modern algorithms have a significant advantage over classical acceleration based methods. In the last chapter 6 we give new and improved PAC-Bayesian bounds for the risk of stochastic neural nets. This chapter also includes an experimental investigation revealing new geometric properties of the paths in weight space that are traced out by the net during the training.


title: Dynamic Mode Decomposition in Adaptive Mesh Refinement and Coarsening Simulations

abstract: Dynamic Mode Decomposition (DMD) is a powerful data-driven method used to extract spatio-temporal coherent structures that dictate a given dynamical system. The method consists of stacking collected temporal snapshots into a matrix and mapping the nonlinear dynamics using a linear operator. The standard procedure considers that snapshots possess the same dimensionality for all the observable data. However, this often does not occur in numerical simulations with adaptive mesh refinement/coarsening schemes (AMR/C). This paper proposes a strategy to enable DMD to extract features from observations with different mesh topologies and dimensions, such as those found in AMR/C simulations. For this purpose, the adaptive snapshots are projected onto the same reference function space, enabling the use of snapshot-based methods such as DMD. The present strategy is applied to challenging AMR/C simulations: a continuous diffusion-reaction epidemiological model for COVID-19, a density-driven gravity current simulation, and a bubble rising problem. We also evaluate the DMD efficiency to reconstruct the dynamics and some relevant quantities of interest. In particular, for the SEIRD model and the bubble rising problem, we evaluate DMD's ability to extrapolate in time (short-time future estimates).


title: How Sublimation Delays the Onset of Dusty Debris Disk Formation Around White Dwarf Stars

abstract: Although numerous white dwarf stars host dusty debris disks, the temperature distribution of these stars differs significantly from the white dwarf population as a whole. Dusty debris disks exist exclusively around white dwarfs cooler than 27,000 K. This is all the more enigmatic given that the formation processes of dusty debris disks should favor younger, hotter white dwarfs, which likely host more dynamically unstable planetary systems. Here we apply a sophisticated material sublimation model to white dwarf systems to show that these statistics are actually a natural result of the interplay of thermal and tidal forces, and show how they define the circumstellar regions where dusty debris disks can form. We demonstrate that these processes tend to prevent stability against both sublimative destruction and reaccretion into planetesimals for rocky materials until white dwarfs cool to below ~25,000-32,000 K, in agreement with the observed limit of ~27,000 K. For pure water ice, this critical temperature is less than 2,700 K (requiring a cooling age older the universe); this precludes pure water ice-rich debris disks forming through the accepted two-step mechanism. The critical temperature is size-dependent; more massive white dwarfs could potentially host dusty debris disks at warmer temperatures. Our model suggests that the location of the disks within the PG 0010+280, GD 56, GD 362, and PG 1541+651 systems are consistent with a forsterite-dominated olivine composition. We also find that very cool white dwarfs may simultaneously host multiple, independently formed dusty debris disks, consistent with observations of the LSPM J0207+3331 system.


title: Evaluation of statistical approaches for association testing in noisy drug screening data

abstract: dentifying associations among biological variables is a major challenge in modern quantitative biological research, particularly given the systemic and statistical noise endemic to biological systems. Drug sensitivity data has proven to be a particularly challenging field for identifying associations to inform patient treatment. To address this, we introduce two semi-parametric variations on the commonly used concordance index: the robust concordance index and the kernelized concordance index (rCI, kCI), which incorporate measurements about the noise distribution from the data. We demonstrate that common statistical tests applied to the concordance index and its variations fail to control for false positives, and introduce efficient implementations to compute p-values using adaptive permutation testing. We then evaluate the statistical power of these coefficients under simulation and compare with Pearson and Spearman correlation coefficients. Finally, we evaluate the various statistics in matching drugs across pharmacogenomic datasets. We observe that the rCI and kCI are better powered than the concordance index in simulation and show some improvement on real data. Surprisingly, we observe that the Pearson correlation was the most robust to measurement noise among the different metrics.


title: Assessing the Impact of Deep Neural Network-based Image Denoising on Binary Signal Detection Tasks

abstract: A variety of deep neural network (DNN)-based image denoising methods have been proposed for use with medical images. Traditional measures of image quality (IQ) have been employed to optimize and evaluate these methods. However, the objective evaluation of IQ for the DNN-based denoising methods remains largely lacking. In this work, we evaluate the performance of DNN-based denoising methods by use of task-based IQ measures. Specifically, binary signal detection tasks under signal-known-exactly (SKE) with background-known-statistically (BKS) conditions are considered. The performance of the ideal observer (IO) and common linear numerical observers are quantified and detection efficiencies are computed to assess the impact of the denoising operation on task performance. The numerical results indicate that, in the cases considered, the application of a denoising network can result in a loss of task-relevant information in the image. The impact of the depth of the denoising networks on task performance is also assessed. The presented results highlight the need for the objective evaluation of IQ for DNN-based denoising technologies and may suggest future avenues for improving their effectiveness in medical imaging applications.


title: Riemann-Hilbert problem on an elliptic surface and a uniformly stressed inclusion embedded into a half-plane subjected to antiplane strain

abstract: An inverse problem of elasticity of $n$ elastic inclusions embedded into an elastic half-plane is analyzed. The boundary of the half-plane is free of traction. The half-plane and the inclusions are subjected to antiplane shear, and the conditions of ideal contact hold in the interfaces between the inclusions and the half-plane. The shapes of the inclusions are not prescribed and have to be determined by enforcing uniform stresses inside the inclusions. The method of conformal mappings from a slit domain onto the $(n+1)$-connected physical domain is worked out. It is shown that to recover the map and therefore the inclusions shapes, one needs to solve a vector Riemann-Hilbert problem on a genus-$n$ hyperelliptic surface. In a particular case of loading of a single inclusion in a half-plane, the problem is equivalent to two scalar Riemann-Hilbert problems on two slits on an elliptic surface. In addition to three parameters of the model the conformal map possesses a free geometric parameter. Results of numerical tests which show the impact of these parameters on the inclusion shape are presented.


title: Collision Probabilities in the Edgeworth-Kuiper belt

abstract: Here, we present results on the intrinsic collision probabilities, $ P_I$, and range of collision speeds, $V_I$, as a function of the heliocentric distance, $r$, in the trans-Neptunian region. The collision speed is one of the parameters, that serves as a proxy to a collisional outcome e.g., complete disruption and scattering of fragments, or formation of crater, where both processes are directly related to the impact energy. We utilize an improved and de-biased model of the trans-Neptunian object (TNO) region from the "Outer Solar System Origins Survey" (OSSOS). It provides a well-defined orbital distribution model of TNOs, based on multiple opposition observations of more than 1000 bodies. In this work we compute collisional probabilities for the OSSOS models of the main classical, resonant, detached+outer and scattering TNO populations. The intrinsic collision probabilities and collision speeds are computed using the \"{O}pik's approach, as revised and modified by Wetherill for non-circular and inclined orbits. The calculations are carried out for each of the dynamical TNO groups, allowing for inter-population collisions as well as collisions within each TNO population, resulting in 28 combinations in total. Our results indicate that collisions in the trans-Neptunian region are possible over a wide range in ($r, V_I$) phase space. Although collisions are calculated to happen within $r\sim 20 - 200$~AU and $V_I \sim 0.1$~km/s to as high as $V_I\sim9$~km/s, most of the collisions are likely to happen at low relative velocities $V_I<1$~km/s and are dominated by the main classical belt.


title: Pushing it out of the Way: Interactive Visual Navigation

abstract: We have observed significant progress in visual navigation for embodied agents. A common assumption in studying visual navigation is that the environments are static; this is a limiting assumption. Intelligent navigation may involve interacting with the environment beyond just moving forward/backward and turning left/right. Sometimes, the best way to navigate is to push something out of the way. In this paper, we study the problem of interactive navigation where agents learn to change the environment to navigate more efficiently to their goals. To this end, we introduce the Neural Interaction Engine (NIE) to explicitly predict the change in the environment caused by the agent's actions. By modeling the changes while planning, we find that agents exhibit significant improvements in their navigational capabilities. More specifically, we consider two downstream tasks in the physics-enabled, visually rich, AI2-THOR environment: (1) reaching a target while the path to the target is blocked (2) moving an object to a target location by pushing it. For both tasks, agents equipped with an NIE significantly outperform agents without the understanding of the effect of the actions indicating the benefits of our approach.


title: What Did It Look Like: A service for creating website timelapses using the Memento framework

abstract: Popular web pages are archived frequently, which makes it difficult to visualize the progression of the site through the years at web archives. The What Did It Look Like (WDILL) Twitter bot shows web page transitions by creating a timelapse of a given website using one archived copy from each calendar year. Originally implemented in 2015, we recently added new features to WDILL, such as date range requests, diversified memento selection, updated visualizations, and sharing visualizations to Instagram. This would allow scholars and the general public to explore the temporal nature of web archives.


title: Weather and Light Level Classification for Autonomous Driving: Dataset, Baseline and Active Learning

abstract: Autonomous driving is rapidly advancing, and Level 2 functions are becoming a standard feature. One of the foremost outstanding hurdles is to obtain robust visual perception in harsh weather and low light conditions where accuracy degradation is severe. It is critical to have a weather classification model to decrease visual perception confidence during these scenarios. Thus, we have built a new dataset for weather (fog, rain, and snow) classification and light level (bright, moderate, and low) classification. Furthermore, we provide street type (asphalt, grass, and cobblestone) classification, leading to 9 labels. Each image has three labels corresponding to weather, light level, and street type. We recorded the data utilizing an industrial front camera of RCCC (red/clear) format with a resolution of $1024\times1084$. We collected 15k video sequences and sampled 60k images. We implement an active learning framework to reduce the dataset's redundancy and find the optimal set of frames for training a model. We distilled the 60k images further to 1.1k images, which will be shared publicly after privacy anonymization. There is no public dataset for weather and light level classification focused on autonomous driving to the best of our knowledge. The baseline ResNet18 network used for weather classification achieves state-of-the-art results in two non-automotive weather classification public datasets but significantly lower accuracy on our proposed dataset, demonstrating it is not saturated and needs further research.


title: Where to Refuel: Modeling On-the-way Choice of Convenience Outlet

abstract: This paper introduces on-the-way choice of retail outlet as a form of convenience shopping. It presents a model of on-the-way choice of retail outlet and applies the model in the context of fuel retailing to explore its implications for segmentation and spatial competition. The model is a latent class random utility choice model. An application to gas station choices observed in a medium-sized Asian city show the model to fit substantially better than existing models. The empirical results indicate consumers may adopt one of two decision strategies. When adopting an immediacy-oriented strategy they behave in accordance with the traditional gravity-based retail models and tend to choose the most spatially convenient outlet. When following a destination-oriented strategy they focus more on maintaining their overall trip efficiency and so will tend to visit outlets located closer to their main destination and are more susceptible to retail agglomeration effects. The paper demonstrates how the model can be used to inform segmentation and local competition analyses that account for variations in these strategies as well as variations in consumer type, origin and time of travel. Simulations of a duopoly setting further demonstrate the implications.


title: Peierls/Su-Schrieffer-Heeger polarons in two dimensions

abstract: Polarons with different types of electron-phonon coupling have fundamentally different properties. When the dominant interaction is between the electron density and lattice displacement, the momentum of the ground state does not change and the polaron gets exponentially heavy at strong coupling. In contrast, one-dimensional Peierls/Su-Schrieffer-Heeger (PSSH) polarons with interaction originating from displacement-modulated hopping feature a shift of the ground-state momentum to finite values and moderate values of effective mass as coupling is increased REF[Phys. Rev. Lett. 105, 266605 (2010)]. Based on Diagrammatic Monte Carlo method, we investigate whether unusual properties of PSSH polarons depend on the type of the displacement-modulated hopping and to what degree they survive in higher dimension. We study two different PSSH models: with bosonic degrees of freedom residing on sites (model A) and bonds (model B) of the two-dimensional square lattice. For model A, we find that in both adiabatic and intermediate regimes, the momentum of the ground state experiences a continuous transition from zero to a finite value as a function of coupling strength. The transition is driven by quadratic instability of the dispersion function, implying that effective mass diverges at the critical point, and then decreases in an anisotropic fashion with increasing coupling. Unexpectedly, for model B, the momentum of the ground state always stays at zero and the effective mass increases monotonously with coupling. The increase is far from exponential and tends to level-off at strong interaction, resulting in relatively light polarons. Having light polarons in the strong coupling regime is crucial for the bi-polaron mechanism of high-temperature superconductivity REF[Phys. Rev. Lett. 121, 247001 (2018)].


title: Deep Neural Network as an alternative to Boosted Decision Trees for PID

abstract: In this paper we recreate, and improve, the binary classification method for particles proposed in Roe et al. (2005) paper "Boosted decision trees as an alternative to artificial neural networks for particle identification". Such particles are tau neutrinos, which we will refer to as background, and electronic neutrinos: the signal we are interested in. In the original paper the preferred algorithm is a Boosted decision tree. This is due to its low effort tuning and good overall performance at the time. Our choice for implementation is a deep neural network, faster and more promising in performance. We will show how, using modern techniques, we are able to improve on the original result, both in accuracy and in training time.


title: Ten-tier and multi-scale supplychain network analysis of medical equipment: Random failure and intelligent attack analysis

abstract: Motivated by the COVID-19 pandemic, this paper explores the supply chain viability of medical equipment, an industry whose supply chain was put under a crucial test during the pandemic. This paper includes an empirical network-level analysis of supplier reachability under Random Failure Experiment (RFE) and Intelligent Attack Experiment (IAE). Specifically, this study investigates the effect of RFA and IAE across multiple tiers and scales. The global supply chain data was mined and analyzed from about 45,000 firms with about 115,000 intertwined relationships spanning across 10 tiers of the backward supply chain of medical equipment. This complex supply chain network was analyzed at four scales, namely: firm, country-industry, industry, and country. A notable contribution of this study is the application of a supply chain tier optimization tool to identify the lowest tier of the supply chain that can provide adequate resolution for the study of the supply chain pattern. We also developed data-driven-tools to identify the thresholds for breakdown and fragmentation of the medical equipment supply chain when faced with random failures or different intelligent attack scenarios. The novel network analysis tools utilized in the study can be applied to the study of supply chain reachability and viability in other industries.


title: Grauert-Riemenschneider multiplier ideal sheaves and the (optimal) Brian\c{c}on-Skoda number

abstract: The goal of this note is to survey some recent results on the Grauert-Riemenschneider multiplier ideal sheaves on any (reduced) complex space of pure dimension. In particular, we obtain the Brian\c{c}on-Skoda number for any Noetherian ring of weakly holomorphic functions with weakly rational singularities (\emph{not} essentially of finite type over $\mathbb{C}$ and Cohen-Macaulay local rings), which will partially answer a question of Huneke.


title: Characterizing representability by principal congruences for finite distributive lattices with a join-irreducible unit element

abstract: For a finite distributive lattice $D$, let us call $Q \subseteq D$ \emph{principal congruence representable}, if there is a finite lattice $L$ such that the congruence lattice of $L$ is isomorphic to $D$ and the principal congruences of $L$ correspond to $Q$ under this isomorphism. We find a necessary condition for representability by principal congruences and prove that for finite distributive lattices with a join-irreducible unit element this condition is also sufficient.


title: Continuous Decoding of Daily-Life Hand Movements from Forearm Muscle Activity for Enhanced Myoelectric Control of Hand Prostheses

abstract: State-of-the-art motorized hand prostheses are endowed with actuators able to provide independent and proportional control of as many as six degrees of freedom (DOFs). The control signals are derived from residual electromyographic (EMG) activity, recorded concurrently from relevant forearm muscles. Nevertheless, the functional mapping between forearm EMG activity and hand kinematics is only known with limited accuracy. Therefore, no robust method exists for the reliable computation of control signals for the independent and proportional actuation of more than two DOFs. A common approach to deal with this limitation is to pre-program the prostheses for the execution of a restricted number of behaviors (e.g., pinching, grasping, and wrist rotation) that are activated by the detection of specific EMG activation patterns. However, this approach severely limits the range of activities users can perform with the prostheses during their daily living. In this work, we introduce a novel method, based on a long short-term memory (LSTM) network, to continuously map forearm EMG activity onto hand kinematics. Critically, unlike previous work, which often focuses on simple and highly controlled motor tasks, we tested our method on a dataset of activities of daily living (ADLs): the KIN-MUS UJI dataset. To the best of our knowledge, ours is the first reported work on the prediction of hand kinematics that uses this challenging dataset. Remarkably, we show that our network is able to generalize to novel untrained ADLs. Our results suggest that the presented method is suitable for the generation of control signals for the independent and proportional actuation of the multiple DOFs of state-of-the-art hand prostheses.


title: The Hidden cost of the Edge: A Performance Comparison of Edge and Cloud Latencies

abstract: Edge computing has emerged as a popular paradigm for running latency-sensitive applications due to its ability to offer lower network latencies to end-users. In this paper, we argue that despite its lower network latency, the resource-constrained nature of the edge can result in higher end-to-end latency, especially at higher utilizations, when compared to cloud data centers. We study this edge performance inversion problem through an analytic comparison of edge and cloud latencies and analyze conditions under which the edge can yield worse performance than the cloud. To verify our analytic results, we conduct a detailed experimental comparison of the edge and the cloud latencies using a realistic application and real cloud workloads. Both our analytical and experimental results show that even at moderate utilizations, the edge queuing delays can offset the benefits of lower network latencies, and even result in performance inversion where running in the cloud would provide superior latencies. We finally discuss practical implications of our results and provide insights into how application designers and service providers should design edge applications and systems to avoid these pitfalls.


title: Depinning of the Charge-Density Waves in Quasi-2D 1T-TaS2 Devices Operating at Room Temperature

abstract: We report on depinning of nearly-commensurate charge-density waves in 1T-TaS2 thin-films at room temperature. A combination of the differential current-voltage measurements with the low-frequency noise spectroscopy provide unambiguous means for detecting the depinning threshold field in quasi-2D materials. The depinning process in 1T-TaS2 is not accompanied by an observable abrupt increase in electric current - in striking contrast to depinning in the conventional charge-density-wave materials with quasi-1D crystal structure. We explained it by the fact that the current density from the charge-density waves in the 1T-TaS2 devices is orders of magnitude smaller than the current density of the free carriers available in the discommensuration network surrounding the commensurate charge-density-wave islands. The depinning fields in 1T-TaS2 thin-film devices are several orders of magnitude larger than those in quasi-1D van der Waals materials. Obtained results are important for the proposed applications of the charge-density-wave devices in electronics.


title: K\'arm\'an vortex street for the generalized surface quasi-geostrophic equation

abstract: We are concerned with the existence of periodic travelling-wave solutions for the generalized surface quasi-geostrophic (gSQG) equation(including incompressible Euler equation), known as von K\'arm\'an vortex street. These solutions are of $C^1$ type, and are obtained by studying a semilinear problem on an infinite strip whose width equals to the period. By a variational characterization of solutions, we also show the relationship between vortex size, travelling speed and street structure. In particular, the vortices with positive and negative intensity have equal or unequal scaling size in our construction, which constitutes the regularization for K\'arm\'an point vortex street.


title: A positivity-preserving, energy stable scheme for a Ternary Cahn-Hilliard system with the singular interfacial parameters

abstract: In this paper, we construct and analyze a uniquely solvable, positivity preserving and unconditionally energy stable finite-difference scheme for the periodic three-component Macromolecular Microsphere Composite (MMC) hydrogels system, a ternary Cahn-Hilliard system with a Flory-Huggins-deGennes free energy potential. The proposed scheme is based on a convex-concave decomposition of the given energy functional with two variables, and the centered difference method is adopted in space. We provide a theoretical justification that this numerical scheme has a pair of unique solutions, such that the positivity is always preserved for all the singular terms, i.e., not only two phase variables are always between $0$ and $1$, but also the sum of two phase variables is between $0$ and $1$, at a point-wise level. In addition, we use the local Newton approximation and multigrid method to solve this nonlinear numerical scheme, and various numerical results are presented, including the numerical convergence test, positivity-preserving property test, energy dissipation and mass conservation properties.


title: Loss-Based Variational Bayes Prediction

abstract: We propose a new approach to Bayesian prediction that caters for models with a large number of parameters and is robust to model misspecification. Given a class of high-dimensional (but parametric) predictive models, this new approach constructs a posterior predictive using a variational approximation to a generalized posterior that is directly focused on predictive accuracy. The theoretical behavior of the new prediction approach is analyzed and a form of optimality demonstrated. Applications to both simulated and empirical data using high-dimensional Bayesian neural network and autoregressive mixture models demonstrate that the approach provides more accurate results than various alternatives, including misspecified likelihood-based predictions.


title: Traversable Wormholes, Regular Black Holes, and Black-Bounces

abstract: Various spacetime candidates for traversable wormholes, regular black holes, and `black-bounces' are presented and thoroughly explored in the context of the gravitational theory of general relativity. All candidate spacetimes belong to the mathematically simple class of spherically symmetric geometries; the majority are static, with a single dynamical (time-dependent) geometry explored. To the extent possible, the candidates are presented through the use of a global coordinate patch -- some of the prior literature (especially concerning traversable wormholes) has often proposed coordinate systems for desirable solutions to the Einstein equations requiring a multi-patch atlas. The most interesting cases include the so-called `exponential metric' -- well-favoured by proponents of alternative theories of gravity but which actually has a standard classical interpretation, and the `black-bounce' to traversable wormhole case -- where a metric is explored which represents either a traversable wormhole or a regular black hole, depending on the value of the newly introduced scalar parameter $a$. This notion of `black-bounce' is defined as the case where the spherical boundary of a regular black hole forces one to travel towards a one-way traversable `bounce' into a future reincarnation of our own universe. The metric of interest is then explored further in the context of a time-dependent spacetime, where the line element is rephrased with a Vaidya-like time-dependence imposed on the mass of the object, and in terms of outgoing\-/ingoing Eddington-Finkelstein coordinates. Analysing these candidate spacetimes extends the pre-existing discussion concerning the viability of non-singular black hole solutions in the context of general relativity, as well as contributing to the dialogue on whether an arbitrarily advanced civilization would be able to construct a traversable wormhole.


title: Machine Learning Techniques for Software Quality Assurance: A Survey

abstract: Over the last years, machine learning techniques have been applied to more and more application domains, including software engineering and, especially, software quality assurance. Important application domains have been, e.g., software defect prediction or test case selection and prioritization. The ability to predict which components in a large software system are most likely to contain the largest numbers of faults in the next release helps to better manage projects, including early estimation of possible release delays, and affordably guide corrective actions to improve the quality of the software. However, developing robust fault prediction models is a challenging task and many techniques have been proposed in the literature. Closely related to estimating defect-prone parts of a software system is the question of how to select and prioritize test cases, and indeed test case prioritization has been extensively researched as a means for reducing the time taken to discover regressions in software. In this survey, we discuss various approaches in both fault prediction and test case prioritization, also explaining how in recent studies deep learning algorithms for fault prediction help to bridge the gap between programs' semantics and fault prediction features. We also review recently proposed machine learning methods for test case prioritization (TCP), and their ability to reduce the cost of regression testing without negatively affecting fault detection capabilities.


title: Chern conjecture on minimal hypersurfaces

abstract: In this paper, we study $n$-dimensional complete minimal hypersurfaces in a unit sphere. We prove that an $n$-dimensional complete minimal hypersurface with constant scalar curvature in a unit sphere with $f_3$ constant is isometric to the totally geodesic sphere or the Clifford torus if $S\leq 1.8252 n-0.712898$, where $S$ denotes the squared norm of the second fundamental form of this hypersurface.


title: On a class of $k$-entanglement witnesses

abstract: Recently, Yang at al. showed that each 2-positive map acting from $\mathcal{M}_3(\mathbb{C})$ into itself is decomposable. It is equivalent to the statement that each PPT state on $\mathbb{C}^3\otimes\mathbb{C}^3$ has Schmidt number at most 2. It is a generalization of Perez-Horodecki criterion which states that each PPT state on $\mathbb{C}^2\otimes\mathbb{C}^2$ or $\mathbb{C}^2\otimes\mathbb{C}^3$ has Schmidt rank 1 i.e. is separable. Natural question arises whether the result of Yang at al. stays true for PPT states on $\mathbb{C}^3\otimes\mathbb{C}^4$. This question can be considered also in higher dimensions. We construct a positive maps which is suspected for being a counterexample. More generally, we provide a class of positive maps between matrix algebras whose $k$-positivity properties can be easily controlled.


title: The second gap on complete self-shrinkers

abstract: In this paper, we study complete self-shrinkers in Euclidean space and prove that an $n$-dimensional complete self-shrinker in Euclidean space $\mathbb{R}^{n+1}$ is isometric to either $\mathbb{R}^{n}$, $S^{n}(\sqrt{n})$, or $S^k (\sqrt{k})\times\mathbb{R}^{n-k}$, $1\leq k\leq n-1$, if the squared norm $S$ of the second fundamental form, $f_3$ are constant and $S$ satisfies $S<1.83379$. We should remark that the condition of polynomial volume growth is not assumed.


title: WGCN: Graph Convolutional Networks with Weighted Structural Features

abstract: Graph structural information such as topologies or connectivities provides valuable guidance for graph convolutional networks (GCNs) to learn nodes' representations. Existing GCN models that capture nodes' structural information weight in- and out-neighbors equally or differentiate in- and out-neighbors globally without considering nodes' local topologies. We observe that in- and out-neighbors contribute differently for nodes with different local topologies. To explore the directional structural information for different nodes, we propose a GCN model with weighted structural features, named WGCN. WGCN first captures nodes' structural fingerprints via a direction and degree aware Random Walk with Restart algorithm, where the walk is guided by both edge direction and nodes' in- and out-degrees. Then, the interactions between nodes' structural fingerprints are used as the weighted node structural features. To further capture nodes' high-order dependencies and graph geometry, WGCN embeds graphs into a latent space to obtain nodes' latent neighbors and geometrical relationships. Based on nodes' geometrical relationships in the latent space, WGCN differentiates latent, in-, and out-neighbors with an attention-based geometrical aggregation. Experiments on transductive node classification tasks show that WGCN outperforms the baseline models consistently by up to 17.07% in terms of accuracy on five benchmark datasets.


title: Two network Kuramoto-Sakaguchi model under tempered stable L\'evy noise

abstract: We examine a model of two interacting populations of phase oscillators labelled `Blue' and `Red'. To this we apply tempered stable L\'{e}vy noise, a generalisation of Gaussian noise where the heaviness of the tails parametrised by a power law exponent $\alpha$ can be controlled by a tempering parameter $\lambda$. This system models competitive dynamics, where each population seeks both internal phase synchronisation and a phase advantage with respect to the other population, subject to exogenous stochastic shocks. We study the system from an analytic and numerical point of view to understand how the phase lag values and the shape of the noise distribution can lead to steady or noisy behaviour. Comparing the analytic and numerical studies shows that the bulk behaviour of the system can be effectively described by dynamics in the presence of tilted ratchet potentials. Generally, changes in $\alpha$ away from the Gaussian noise limit, $1< \alpha < 2$, disrupts the locking between Blue and Red, while increasing $\lambda$ acts to restore it. However we observe that with further decreases of $\alpha$ to small values, $\alpha\ll 1$, with $\lambda\neq 0$, locking between Blue and Red may be restored. This is seen analytically in a restoration of metastability through the ratchet mechanism, and numerically in transitions between periodic and noisy regions in a fitness landscape using a measure of noise. This non-monotonic transition back to an ordered regime is surprising for a linear variation of a parameter such as the power law exponent and provides a novel mechanism for guiding the collective behaviour of such a complex competitive dynamical system.