title: A prognostic dynamic model applicable to infectious diseases providing easily visualized guides -- A case study of COVID-19 in the UK

abstract: A reasonable prediction of infectious diseases transmission process under different disease control strategies is an important reference point for policy makers. Here we established a dynamic transmission model via Python and realized comprehensive regulation of disease control measures. We classified government interventions into three categories and introduced three parameters as descriptions for the key points in disease control, these being intraregional growth rate, interregional communication rate, and detection rate of infectors. Our simulation predicts the infection by COVID-19 in the UK would be out of control in 73 days without any interventions; at the same time, herd immunity acquisition will begin from the epicentre. After we introduced government interventions, single intervention is effective in disease control but at huge expense while combined interventions would be more efficient, among which, enhancing detection number is crucial in control strategy of COVID-19. In addition, we calculated requirements for the most effective vaccination strategy based on infection number in real situation. Our model was programmed with iterative algorithms, and visualized via cellular automata, it can be applied to similar epidemics in other regions if the basic parameters are inputted, and is able to synthetically mimick the effect of multiple factors in infectious disease control.


title: Light-enhanced gating effect at the interface of oxide heterostructure

abstract: In semiconducting materials, electrostatic gating and light illumination are widely used stimuli to tune the electronic properties of the system. Here, we show a significant enhancement of photoresponse at the conducting interface of LaVO3-SrTiO3 under the simultaneous application of light and negative gate bias voltage, in comparison to their individual application. On the other hand, the LaVO3-SrTiO3 interface remains largely insensitive to light illumination, when a positive gate bias voltage is applied. Our X-ray diffractometer, Raman spectroscopy and photoemission measurements show that unlike the LaAlO3-SrTiO3 interface, migration of oxygen vacancies is not the prime mechanism for the enhanced photoresponse. Rather, we suggest that the photoresponse of our system is intrinsic and this intrinsic mechanism is a complex interplay between band filling, electric field at the interface, strong electron interaction due to mottness of LaVO3 and modification of conducting channel width.


title: Knowledge-Routed Visual Question Reasoning: Challenges for Deep Representation Embedding

abstract: Though beneficial for encouraging the Visual Question Answering (VQA) models to discover the underlying knowledge by exploiting the input-output correlation beyond image and text contexts, the existing knowledge VQA datasets are mostly annotated in a crowdsource way, e.g., collecting questions and external reasons from different users via the internet. In addition to the challenge of knowledge reasoning, how to deal with the annotator bias also remains unsolved, which often leads to superficial over-fitted correlations between questions and answers. To address this issue, we propose a novel dataset named Knowledge-Routed Visual Question Reasoning for VQA model evaluation. Considering that a desirable VQA model should correctly perceive the image context, understand the question, and incorporate its learned knowledge, our proposed dataset aims to cutoff the shortcut learning exploited by the current deep embedding models and push the research boundary of the knowledge-based visual question reasoning. Specifically, we generate the question-answer pair based on both the Visual Genome scene graph and an external knowledge base with controlled programs to disentangle the knowledge from other biases. The programs can select one or two triplets from the scene graph or knowledge base to push multi-step reasoning, avoid answer ambiguity, and balanced the answer distribution. In contrast to the existing VQA datasets, we further imply the following two major constraints on the programs to incorporate knowledge reasoning: i) multiple knowledge triplets can be related to the question, but only one knowledge relates to the image object. This can enforce the VQA model to correctly perceive the image instead of guessing the knowledge based on the given question solely; ii) all questions are based on different knowledge, but the candidate answers are the same for both the training and test sets.


title: MAGIC observations of the nearby short gamma-ray burst GRB 160821B

abstract: The coincident detection of GW170817 in gravitational waves and electromagnetic radiation spanning the radio to MeV gamma-ray bands provided the first direct evidence that short gamma-ray bursts (GRBs) can originate from binary neutron star (BNS) mergers. On the other hand, the properties of short GRBs in high-energy gamma rays are still poorly constrained, with only $\sim$20 events detected in the GeV band, and none in the TeV band. GRB~160821B is one of the nearest short GRBs known at $z=0.162$. Recent analyses of the multiwavelength observational data of its afterglow emission revealed an optical-infrared kilonova component, characteristic of heavy-element nucleosynthesis in a BNS merger. Aiming to better clarify the nature of short GRBs, this burst was automatically followed up with the MAGIC telescopes, starting from 24 seconds after the burst trigger. Evidence of a gamma-ray signal is found above $\sim$0.5 TeV at a significance of $\sim3\,\sigma$ during observations that lasted until 4 hours after the burst. Assuming that the observed excess events correspond to gamma-ray emission from GRB 160821B, in conjunction with data at other wavelengths, we investigate its origin in the framework of GRB afterglow models. The simplest interpretation with one-zone models of synchrotron-self-Compton emission from the external forward shock has difficulty accounting for the putative TeV flux. Alternative scenarios are discussed where the TeV emission can be relatively enhanced. The role of future GeV-TeV observations of short GRBs in advancing our understanding of BNS mergers and related topics is briefly addressed.


title: The Martini Model in Materials Science

abstract: The Martini model, a coarse-grained force field initially developed with biomolecular simulations in mind, has found an increasing number of applications in the field of soft materials science. The model's underlying building block principle does not pose restrictions on its application beyond biomolecular systems. Here we highlight the main applications to date of the Martini model in materials science, and we give a perspective for the future developments in this field, in particular in light of recent developments such as the new version of the model, Martini 3.


title: Efficient Querying for Cooperative Probabilistic Commitments

abstract: Multiagent systems can use commitments as the core of a general coordination infrastructure, supporting both cooperative and non-cooperative interactions. Agents whose objectives are aligned, and where one agent can help another achieve greater reward by sacrificing some of its own reward, should choose a cooperative commitment to maximize their joint reward. We present a solution to the problem of how cooperative agents can efficiently find an (approximately) optimal commitment by querying about carefully-selected commitment choices. We prove structural properties of the agents' values as functions of the parameters of the commitment specification, and develop a greedy method for composing a query with provable approximation bounds, which we empirically show can find nearly optimal commitments in a fraction of the time methods that lack our insights require.


title: Software Quality Assessment for Robot Operating System

abstract: Robot Operating System (ROS) is widely used in academia and industry, and importantly is leveraged in safety-critical robotic systems. The quality of ROS software can affect the safety and security properties of robotics systems; therefore, reliability and quality are imperative to guarantee. Source code static analysis is a key approach to formally perform software verification. We address two concerns in this paper: (1) conducting a systematic literature review study to provide a complete picture of the existing methods that analyze different aspects of ROS software, (2) performing empirical study to evaluate software properties that can influence the functionality of ROS. We leverage PMD1, an off-the-shelf static analysis tool, to conduct our empirical study over a set of ROS repositories implemented using Java. The survey analysis shows a significant shortcoming in the body of research by the lack of tailored analysis mechanisms for assessing ROS2 code and reveals that the majority of research efforts are centered around ROS1. Our empirical study shows that the Java code of ROS2 does not suffer from serious issues and the majority of the detected alerts are code style issues.


title: Bayes DistNet -- A Robust Neural Network for Algorithm Runtime Distribution Predictions

abstract: Randomized algorithms are used in many state-of-the-art solvers for constraint satisfaction problems (CSP) and Boolean satisfiability (SAT) problems. For many of these problems, there is no single solver which will dominate others. Having access to the underlying runtime distributions (RTD) of these solvers can allow for better use of algorithm selection, algorithm portfolios, and restart strategies. Previous state-of-the-art methods directly try to predict a fixed parametric distribution that the input instance follows. In this paper, we extend RTD prediction models into the Bayesian setting for the first time. This new model achieves robust predictive performance in the low observation setting, as well as handling censored observations. This technique also allows for richer representations which cannot be achieved by the classical models which restrict their output representations. Our model outperforms the previous state-of-the-art model in settings in which data is scarce, and can make use of censored data such as lower bound time estimates, where that type of data would otherwise be discarded. It can also quantify its uncertainty in its predictions, allowing for algorithm portfolio models to make better informed decisions about which algorithm to run on a particular instance.


title: Polar Codes for Quantum Reading

abstract: Quantum reading provides a general framework where to formulate the statistical discrimination of quantum channels. Several paths have been taken for such a problem. However, there is much to be done in the avenue of optimizing channel discrimination using classical codes. At least two open questions can be pointed to: how to construct low complexity encoding schemes that are interesting for channel discrimination and, more importantly, how to develop capacity-achieving protocols. The aim of this paper is to present a solution to these questions using polar codes. Firstly, we characterize the rate and reliability of the channels under polar encoding. We also show that the error probability of the scheme proposed decays exponentially with respect to the code length. Lastly, an analysis of the optimal quantum states to be used as probes is given.


title: On Convergence of Gradient Expected Sarsa($\lambda$)

abstract: We study the convergence of $\mathtt{Expected~Sarsa}(\lambda)$ with linear function approximation. We show that applying the off-line estimate (multi-step bootstrapping) to $\mathtt{Expected~Sarsa}(\lambda)$ is unstable for off-policy learning. Furthermore, based on convex-concave saddle-point framework, we propose a convergent $\mathtt{Gradient~Expected~Sarsa}(\lambda)$ ($\mathtt{GES}(\lambda)$) algorithm. The theoretical analysis shows that our $\mathtt{GES}(\lambda)$ converges to the optimal solution at a linear convergence rate, which is comparable to extensive existing state-of-the-art gradient temporal difference learning algorithms. Furthermore, we develop a Lyapunov function technique to investigate how the step-size influences finite-time performance of $\mathtt{GES}(\lambda)$, such technique of Lyapunov function can be potentially generalized to other GTD algorithms. Finally, we conduct experiments to verify the effectiveness of our $\mathtt{GES}(\lambda)$.


title: Contact Lie poset algebras

abstract: We provide a combinatorial recipe for constructing all posets of height at most two for which the corresponding type-A Lie poset algebra is contact. In the case that such posets are connected, a discrete Morse theory argument establishes that the posets' simplicial realizations are contractible. It follows from a cohomological result of Coll and Gerstenhaber on Lie semi-direct products that the corresponding contact Lie algebras are absolutely rigid.


title: Shallow NV centers augmented by exploiting n-type diamond

abstract: Creation of nitrogen-vacancy (NV) centers at the nanoscale surface region in diamond, while retaining their excellent spin and optical properties, is essential for applications in quantum technology. Here, we demonstrate the extension of the spin-coherence time ($\it{T}$${_2}$), the stabilization of the charge state, and an improvement of the creation yield of NV centers formed by the ion-implantation technique at a depth of $\sim$15 nm in phosphorus-doped n-type diamond. The longest $\it{T}$${_2}$ of about 580 $\mu$s of a shallow NV center approaches the one in bulk diamond limited by the nuclear spins of natural abundant $^{13}$C. The averaged $\it{T}$${_2}$ in n-type diamond is over 1.7 times longer than that in pure non-doped diamond. Moreover, the stabilization of the charge state and the more than twofold improvement of the creation yield are confirmed. The enhancements for the shallow NV centers in an n-type diamond-semiconductor are significant for future integrated quantum devices.


title: Frozen-in fermionic singlet dark matter in non-standard cosmology with a decaying fluid

abstract: We perform a detailed study of dark matter production via freeze-in under the assumption that some fluid dominates the early Universe before depositing its energy to the plasma causing entropy injection. As a dark matter candidate we consider a fermionic singlet that is produced through its interactions with a scalar particle in the thermal plasma. The fluid alters the expansion rate of the Universe, as well as the scaling of the temperature, which significantly affects the evolution of both the number density and the mean momentum of the dark matter particle. We identify and discuss in detail the effects of the evolution of these quantities by considering several examples representing dark matter production at different stages of expansion and entropy injection. We find that, since the dark matter density is reduced when the entropy injection to the plasma continues after freeze-in, in order to reproduce its observational value an enhanced rate of dark matter production is required relative to standard cosmology. Furthermore, the impact of the assumed non-standard cosmological history on the dark matter mean momentum can result in either a relaxed or a tightened bound on the dark matter mass from large structure formation data.


title: Lectures on dualities ABC in representation theory

abstract: In these lecture notes for a summer mini-course, we provide an exposition on quantum groups and Hecke algebras, including (quasi) R-matrix, canonical basis, and $q$-Schur duality. Then we formulate their counterparts in the setting of $\imath$quantum groups arising from quantum symmetric pairs, including (quasi) K-matrix, $\imath$-canonical basis, and $\imath$Schur duality. As an application, the ($\imath$-)canonical bases are used to formulate Kazhdan-Lusztig theories and character formulas in the BGG categories for Lie (super)algebras of type A-D. Finally, geometric constructions for $q$-Schur and $\imath$Schur dualities are provided.


title: Generalizations of degeneracy second main theorem and Schmidt's subspace theorem

abstract: In this paper, by introducing the notion of "\textit{distributive constant}" of a family of hypersurfaces with respect to a projective variety, we prove a second main theorem in Nevanlinna theory for meromorphic mappings with arbitrary families of hypersurfaces in projective varieties. Our second main theorem generalizes and improves previous results for meromorphic mappings with hypersurfaces, in particular for algebraically degenerate mappings and for the families of hypersurfaces in subgeneral position. The analogous results for the holomorphic curves with finite growth index from a complex disc into a project variety, and for meromorphic mappings on a complete K\"{a}hler manifold are also given. For the last aim, we will prove a Schmidt's subspace theorem for an arbitrary families of homogeneous polynomials, which is the counterpart in Number theory of our second main theorem. Our these results are generalizations and improvements of all previous results.


title: High-Order Approximation Rates for Shallow Neural Networks with Cosine and ReLU$^k$ Activation Functions

abstract: We study the approximation properties of shallow neural networks with an activation function which is a power of the rectified linear unit. Specifically, we consider the dependence of the approximation rate on the dimension and the smoothness in the spectral Barron space of the underlying function $f$ to be approximated. We show that as the smoothness index $s$ of $f$ increases, shallow neural networks with ReLU$^k$ activation function obtain an improved approximation rate up to a best possible rate of $O(n^{-(k+1)}\log(n))$ in $L^2$, independent of the dimension $d$. The significance of this result is that the activation function ReLU$^k$ is fixed independent of the dimension, while for classical methods the degree of polynomial approximation or the smoothness of the wavelets used would have to increase in order to take advantage of the dimension dependent smoothness of $f$. In addition, we derive improved approximation rates for shallow neural networks with cosine activation function on the spectral Barron space. Finally, we prove lower bounds showing that the approximation rates attained are optimal under the given assumptions.


title: D-LEMA: Deep Learning Ensembles from Multiple Annotations -- Application to Skin Lesion Segmentation

abstract: Medical image segmentation annotations suffer from inter- and intra-observer variations even among experts due to intrinsic differences in human annotators and ambiguous boundaries. Leveraging a collection of annotators' opinions for an image is an interesting way of estimating a gold standard. Although training deep models in a supervised setting with a single annotation per image has been extensively studied, generalizing their training to work with datasets containing multiple annotations per image remains a fairly unexplored problem. In this paper, we propose an approach to handle annotators' disagreements when training a deep model. To this end, we propose an ensemble of Bayesian fully convolutional networks (FCNs) for the segmentation task by considering two major factors in the aggregation of multiple ground truth annotations: (1) handling contradictory annotations in the training data originating from inter-annotator disagreements and (2) improving confidence calibration through the fusion of base models' predictions. We demonstrate the superior performance of our approach on the ISIC Archive and explore the generalization performance of our proposed method by cross-dataset evaluation on the PH2 and DermoFit datasets.


title: Reversing the Curse of Densification in mmWave Networks Through Spatial Multiplexing

abstract: The gold standard of a wireless network is that the throughput increases linearly with the density of access points (APs). However, such a linear throughput gain is suspended in the 5G mmWave network mainly due to the short communication distances in mmWave bands and the dense deployments of mmWave APs. As being operated in the interference-limited regime, the aggregate interference resulted from the increasing mmWave APs will gradually become the network performance bottleneck, which leads to the saturation of the throughput. In this paper, we propose to overcome the densification plateau of a mmWave network by employing spatial multiplexing at APs. To study the effect of spatial multiplexing on mmWave networks, we first derive the coverage probability as a function of spatial multiplexing gain. The fixed-rate coding scheme is then used to provide the network throughput. We also introduce the concept of densification gain to capture the improvement in network throughput achieved through the AP densification. Our results indicate that without the spatial multiplexing at APs, the throughput of mmWave networks will reach the plateau when the density of APs becomes sufficiently large. By enabling the spatial multiplexing at APs, however, the mmWave network can continuously harvest the throughput gains as the number of APs grows. By deriving the upper bound for the throughput of mmWave networks, we quantify the potential throughput improvement as the spatial multiplexing gain increases. However, our numerical results show that such a potential throughput gain cannot be explored by the fixed-rate coding scheme. We then demonstrate the necessity for deploying the multi-rate coding scheme in mmWave networks, especially when the spatial multiplexing gain at mmWave APs is large.


title: INSPIRE: Intensity and spatial information-based deformable image registration

abstract: We present INSPIRE, a top-performing general-purpose method for deformable image registration. INSPIRE brings distance measures which combine intensity and spatial information into an elastic B-splines-based transformation model and incorporates an inverse inconsistency penalization supporting symmetric registration performance. We introduce several theoretical and algorithmic solutions which provide high computational efficiency and thereby applicability of the proposed framework in a wide range of real scenarios. We show that INSPIRE delivers highly accurate, as well as stable and robust registration results. We evaluate the method on a 2D dataset created from retinal images, characterized by presence of networks of thin structures. Here INSPIRE exhibits excellent performance, substantially outperforming the widely used reference methods. {We also evaluate INSPIRE on the Fundus Image Registration Dataset (FIRE), which consists of 134 pairs of separately acquired retinal images. INSPIRE exhibits excellent performance on the FIRE dataset, substantially outperforming several domain-specific methods.} We also evaluate the method on four benchmark datasets of 3D magnetic resonance images of brains, for a total of 2088 pairwise registrations. A comparison with 17 other state-of-the-art methods reveals that INSPIRE provides the best overall performance. Code is available at http://github.com/MIDA-group/inspire


title: Interpretation of XENON1T excess with MeV boosted dark matter

abstract: The XENON1T excess of keV electron recoil events may be induced by the scattering of electrons and long-lived particles with MeV mass and high-speed. We consider a tangible model composed of two scalar MeV dark matter (DM) particles $S_A$ and $S_B$ to interpret the XENON1T keV excess via boosted $S_B$. A small mass splitting $m_{S_A}-m_{S_B}>0$ is introduced and the boosted $S_B$ can be produced by the dark annihilation process of $S_A S_A^\dagger \to \phi \to S_B S_B^\dagger$ via a resonant scalar $\phi$. The $S_B-$electron scattering is intermediated by a vector boson $X$. Although the constraints from BBN, CMB and low-energy experiments set the $X-$mediated $S_B-$electron scattering cross section to be $\lesssim 10^{-35} \mathrm{cm}^2$, the MeV scale DM with a resonance enhanced dark annihilation today can still provide enough boosted $S_B$ and induce the XENON1T keV excess. The relic density of $S_B$ is significantly reduced by the $s$-wave process of $S_B S_B^\dagger \to X X$ which is allowed by the constraints from CMB and 21-cm absorption. A very small relic fraction of $S_B$ is compatible with the stringent bounds on un-boosted $S_B$-electron scattering in DM direct detection and the $S_A$-electron scattering is also allowed.


title: Efficacy of Electrically-Polarized 3D Printed Graphene-blended Spacers on the Flux Enhancement and Scaling Resistance of Water Filtration Membranes

abstract: In this research, an electrically-polarized graphene-polylactic acid (E-GRP) spacer is introduced for the first time by a novel fabrication method, which consists of 3D printing followed by electrical polarization under a high voltage electric field (1.5 kV/cm). The fabricated E-GRP was tested in an osmotic-driven process (forward osmosis system) to evaluate its performance in terms of water flux, reverse solute flux, and ion attraction compared to a 3D printed non-polarized graphene-polylactic acid (GRP) spacer and a polylactic acid (PLA) spacer. The use of the developed E-GRP spacer showed > 50% water flux enhancement (32.4 +- 2 LMH) compared to the system employing the GRP (20.5 +- 2.3 LMH) or PLA (20.8 +- 2.1 LMH) spacer. This increased water flux was attributed to the increased osmotic pressure across the membrane due to the ions adsorbed on the polarized (E-GRP) spacer. The E-GRP spacer also retarded the gypsum scaling on the membrane compared to the GRP spacer due to the dispersion effect of electrostatic forces between the gypsum aggregation and negatively charged surfaces. The electric polarization of the E-GRP spacer was shown to be maintained for > 100 h by observing its salt adsorption properties (in a 3 M NaCl solution).


title: SU(N) q-Toda equations from mass deformed ABJM theory

abstract: It is known that the partition functions of the U(N) x U(N+M) ABJM theory satisfy a set of bilinear relations, which, written in the grand partition function, was recently found to be the q-Painleve III_3 equation. In this paper we have suggested a similar bilinear relation holds for the ABJM theory with N=6 preserving mass deformation for an arbitrary complex value of mass parameter, to which we have provided several non-trivial checks by using the exact values of the partition functions for various N,k,M and the mass parameter. For particular choices of the mass parameters labeled by integers $\nu,a$ as $m_1=m_2=-\pi i(\nu-2a)/\nu$, the bilinear relation corresponds to the q-deformation of the affine SU($\nu$) Toda equation in $\tau$-form.


title: The Noncrossing Bond Poset of a Graph

abstract: The partition lattice and noncrossing partition lattice are well studied objects in combinatorics. Given a graph $G$ on vertex set $\{1,2,\dots, n\}$, its bond lattice, $L_G$, is the subposet of the partition lattice formed by restricting to the partitions whose blocks induce connected subgraphs of $G$. In this article, we introduce a natural noncrossing analogue of the bond lattice, the noncrossing bond poset, $NC_G$, obtained by restricting to the noncrossing partitions of $L_G$. Both the noncrossing partition lattice and the bond lattice have many nice combinatorial properties. We show that, for several families of graphs, the noncrossing bond poset also exhibits these properties. We present simple necessary and sufficient conditions on the graph to ensure the noncrossing bond poset is a lattice. Additionally, for several families of graphs, we give combinatorial descriptions of the M\"obius function and characteristic polynomial of the noncrossing bond poset. These descriptions are in terms of a noncrossing analogue of non-broken circuit (NBC) sets of the graphs and can be thought of as a noncrossing version of Whitney's NBC theorem for the chromatic polynomial. We also consider the shellability and supersolvability of the noncrossing bond poset, providing sufficient conditions for both. We end with some open problems.


title: Subgroups of Classical Groups that are Transitive on Subspaces

abstract: For each finite classical group $G$, we classify the subgroups of $G$ which act transitively on a $G$-invariant set of subspaces of the natural module, where the subspaces are either totally isotropic or nondegenerate. Our proof uses the classification of the maximal factorisations of almost simple groups. As a first application of these results we classify all point-transitive subgroups of automorphisms of finite thick generalised quadrangles.


title: A survey of sketches in traffic measurement: Design, Optimization, Application and Implementation

abstract: Network measurement probes the underlying network to support upper-level decisions such as network management, network update, network maintenance, network defense and beyond. Due to the massive, speedy, unpredictable features of network flows, sketches are widely implemented in measurement nodes to approximately record the frequency or estimate the cardinality of flows. At their cores, sketches usually maintain one or multiple counter array(s), and rely on hash functions to select the counter(s) for each flow. Then the space-efficient sketches from the distributed measurement nodes are aggregated to provide statistics of the undergoing flows. Currently, tremendous redesigns and optimizations have been proposed to improve the sketches for better network measurement performance. However, existing reviews or surveys mainly focus on one particular aspect of measurement tasks. Researchers and engineers in the network measurement community desire an all-in-one survey that covers the entire processing pipeline of sketch-based network measurement. To this end, we present the first comprehensive survey of this area. We first introduce the preparation of flows for measurement, then detail the most recent investigations of design, aggregation, decoding, application and implementation of sketches for network measurement. To summarize the existing efforts, we carry out an in-depth study of the existing literature, covering more than 90 sketch designs and optimization strategies. Furthermore, we conduct a comprehensive analysis and qualitative/quantitative comparison of the sketch designs. Finally,we highlight the open issues for future sketch-based network measurement research.


title: In-Situ Studies of Stress Environment in Amorphous Solids Using Negatively Charged Nitrogen Vacancy Centers in Nanodiamond

abstract: Amorphous solids, which show characteristic differences from crystals, are common in daily usage. Glasses, gels, and polymers are familiar examples, and polymers are particularly important in terms of their role in construction and crafting. Previous studies have mainly focused on the bulk properties of polymeric products, and the local properties are less discussed. Here, we designed a distinctive protocol using the negatively charged nitrogen vacancy center in nanodiamond to study properties inside polymeric products in situ. Choosing the curing of poly dimethylsiloxane and the polymerization of cyanoacrylate as subjects of investigation, we measured the time dependence of local pressure and strain in the materials during the chemical processes. From the measurements, we were able to probe the local shear stress inside the two polymeric substances in situ. By regarding the surprisingly large shear stress as the internal tension, we attempted to provide a microscopic explanation for the ultimate tensile strength of a bulk solid. Our current methodology is applicable to any kind of transparent amorphous solids with the stress in the order of MPa and to the study of in situ properties in nanoscale. With better apparatus, we expect the limit can be pushed to sub-MPa scale.


title: Nuclear spin squeezing in Helium-3 by continuous quantum nondemolition measurement

abstract: We propose a technique to control the macroscopic collective nuclear spin of a Helium-3 vapor in the quantum regime using light. The scheme relies on metastability exchange collisions to mediate interactions between optically accessible metastable states and the ground-state nuclear spin, giving rise to an effective nuclear spin-light quantum nondemolition interaction of the Faraday form. Our technique enables measurement-based quantum control of nuclear spins, such as the preparation of spin-squeezed states. This, combined with the day-long coherence time of nuclear spin states in Helium-3, opens the possibility for a number of applications in quantum technology.


title: A new 2D auxetic CN2 nanostructure with high energy density and mechanical strength

abstract: The existence of a new two dimensional CN2 structure was predicted using ab-initio molecular dynamics (AIMD) and density-functional theory calculations. It consists tetragonal and hexagonal rings with C-N and N-N bonds arranged in a buckling plane, isostructural to tetrahex-carbon allotrope. It is thermodynamically and kinetically stable suggested by its phonon spectrum and AIMD. This nanosheet has high concentration of N and contains N-N single bonds with an energy density of 6.3 kJ/g, indicating potential applications as high energy density materials. It possesses exotic mechanical properties with negative Poisson's ratio and an anisotropic Young's modulus. The modulus in the zigzag direction is predicted to be 340 N/m, stiffer than h-BN and penta-CN2 sheets and comparable to graphene. Its ideal strength of 28.8 N/m outperforms that of penta-graphene. The material maintains phonon stability upon the application of uniaxial strain up to 10% (13%) in the zigzag (armchair) direction or biaxial strain up to 5%. It possesses a wide indirect HSE band gap of 4.57 eV which is tunable between 3.37 - 4.57 eV through strain. Double-layer structures are also explored. Such unique properties may have potential applications in high energy density materials, nanomechanics and electronics.


title: Intelligent Reflecting Surface Aided Full-Duplex Communication: Passive Beamforming and Deployment Design

abstract: This paper investigates the passive beamforming and deployment design for an intelligent reflecting surface (IRS) aided full-duplex (FD) wireless system, where an FD access point (AP) communicates with an uplink (UL) user and a downlink (DL) user simultaneously over the same time-frequency dimension with the help of IRS. Under this setup, we consider three deployment cases: 1) two distributed IRSs placed near the UL user and DL user, respectively; 2) one centralized IRS placed near the DL user; 3) one centralized IRS placed near the UL user. In each case, we aim to minimize the weighted sum transmit power consumption of the AP and UL user by jointly optimizing their transmit power and the passive reflection coefficients at the IRS (or IRSs), subject to the UL and DL users' rate constraints and the uni-modulus constraints on the IRS reflection coefficients. First, we analyze the minimum transmit power required in the IRS-aided FD system under each deployment scheme, and compare it with that of the corresponding half-duplex (HD) system. We show that the FD system outperforms its HD counterpart for all IRS deployment schemes, while the distributed deployment further outperforms the other two centralized deployment schemes. Next, we transform the challenging power minimization problem into an equivalent but more tractable form and propose an efficient algorithm to solve it based on the block coordinate descent (BCD) method. Finally, numerical results are presented to validate our analysis as well as the efficacy of the proposed passive beamforming design.


title: Breaking the Expressive Bottlenecks of Graph Neural Networks

abstract: Recently, the Weisfeiler-Lehman (WL) graph isomorphism test was used to measure the expressiveness of graph neural networks (GNNs), showing that the neighborhood aggregation GNNs were at most as powerful as 1-WL test in distinguishing graph structures. There were also improvements proposed in analogy to $k$-WL test ($k>1$). However, the aggregators in these GNNs are far from injective as required by the WL test, and suffer from weak distinguishing strength, making it become expressive bottlenecks. In this paper, we improve the expressiveness by exploring powerful aggregators. We reformulate aggregation with the corresponding aggregation coefficient matrix, and then systematically analyze the requirements of the aggregation coefficient matrix for building more powerful aggregators and even injective aggregators. It can also be viewed as the strategy for preserving the rank of hidden features, and implies that basic aggregators correspond to a special case of low-rank transformations. We also show the necessity of applying nonlinear units ahead of aggregation, which is different from most aggregation-based GNNs. Based on our theoretical analysis, we develop two GNN layers, ExpandingConv and CombConv. Experimental results show that our models significantly boost performance, especially for large and densely connected graphs.


title: New stable crystal structures of C and GeC2 predicted from first-principles calculations

abstract: Two novel three-dimensional (3D) crystal structures of carbon (C) and germanium carbide (GeC2) were predicted using first-principles density-functional theory (DFT) calculations. These newly discovered 3D carbon allotrope and GeC2 are in the space group of Fmmm (space group number 69). Their crystal structures have unique tetragonal/hexagonal rings formed by either C or Ge/C atoms. Both structures were proven to be thermodynamically stable through the phonon spectrum calculations. The C allotrope is a semiconductor with a wide band gap of 3.65 eV predicted by the hybrid density functional HSE06 method, while GeC2 is metallic. The new C allotrope has a low mass density of 2.84 g/cm3. More importantly, it is proven to be energetically stable with cohesive energy less than -7.5 eV/atom which is lower than many other carbon allotropes implying the possibility to be fabricated in lab. Such a carbon crystal structure with a low mass density and wide band gap once synthesized would have wide applications in gas adsorption sensors and photo-electronic devices.


title: Intersecting longest paths in chordal graphs

abstract: We consider the size of the smallest set of vertices required to intersect every longest path in a chordal graph. Such sets are known as longest path transversals. We show that if $\omega(G)$ is the clique number of a chordal graph $G$, then there is a transversal of order at most $4\lceil\frac{\omega(G)}{5}\rceil$. We also consider the analogous question for longest cycles, and show that if $G$ is a 2-connected chordal graph then there is a transversal intersecting all longest cycles of order at most $2\lceil\frac{\omega(G)}{3}\rceil$.


title: Performance of Large-Format Deformable Mirrors Constructed with TNO Variable Reluctance Actuators

abstract: Advancements in making high-efficiency actuators are an enabling technology for building the next generation of large-format deformable mirrors. The Netherlands Organization for Applied Scientific Research (TNO) has developed a new style of variable-reluctance actuator that requires approximately eighty times less power to operate as compared to the traditional style of voice-coil actuators. We present the performance results from laboratory testing of TNO's 57-actuator large-format deformable mirror from measuring the influence functions, linearity, hysteresis, natural shape flattening, actuator cross-coupling, creep, repeatability, and actuator lifetime. We measure a linearity of 99.4 +- 0.33% and hysteresis of 2.10 +- 0.23% over a stroke of 10 microns, indicating that this technology has strong potential for use in on-sky adaptive secondary mirrors (ASMs). We summarize plans for future lab prototypes and ASMs that will further demonstrate this technology.


title: Entropy of fully packed hard rigid rods on $d$-dimensional hyper-cubic lattices

abstract: We determine the asymptotic behavior of the entropy of full coverings of a $L \times M$ square lattice by rods of size $k\times 1$ and $1\times k$, in the limit of large $k$. We show that full coverage is possible only if at least one of $L$ and $M$ is a multiple of $k$, and that all allowed configurations can be reached from a standard configuration of all rods being parallel, using only basic flip moves that replace a $k \times k$ square of parallel horizontal rods by vertical rods, and vice versa. In the limit of large $k$, we show that the entropy per site $S_2(k)$ tends to $ A k^{-2} \ln k$, with $A=1$. We conjecture, based on a perturbative series expansion, that this large-$k$ behavior of entropy per site is super-universal and continues to hold on all $d$-dimensional hyper-cubic lattices, with $d \geq 2$.


title: Traffic Rate Network Tomography with Higher-Order Cumulants

abstract: Network tomography aims at estimating source-destination traffic rates from link traffic measurements. This inverse problem was formulated by Vardi in 1996 for Poisson traffic over networks operating under deterministic as well as random routing regimes. In this paper we expand Vardi's second-order moment matching rate estimation approach to higher-order cumulant matching with the goal of increasing the column rank of the mapping and consequently improving the rate estimation accuracy. We develop a systematic set of linear cumulant matching equations and express them compactly in terms of the Khatri-Rao product. Both least squares estimation and iterative minimum I-divergence estimation are considered. We develop an upper bound on the mean squared error (MSE) in least squares rate estimation from empirical cumulants. We demonstrate for the NSFnet that supplementing Vardi's approach with third-order empirical cumulant reduces its averaged normalized MSE relative to the theoretical minimum of the second-order moment matching approach by about 12%-18%. This minimum MSE is obtained when Vardi's second-order moment matching approach is based on the theoretical rather than the empirical moments.


title: Incremental Data-driven Optimization of Complex Systems in Nonstationary Environments

abstract: Existing work on data-driven optimization focuses on problems in static environments, but little attention has been paid to problems in dynamic environments. This paper proposes a data-driven optimization algorithm to deal with the challenges presented by the dynamic environments. First, a data stream ensemble learning method is adopted to train the surrogates so that each base learner of the ensemble learns the time-varying objective function in the previous environments. After that, a multi-task evolutionary algorithm is employed to simultaneously optimize the problems in the past environments assisted by the ensemble surrogate. This way, the optimization tasks in the previous environments can be used to accelerate the tracking of the optimum in the current environment. Since the real fitness function is not available for verifying the surrogates in offline data-driven optimization, a support vector domain description that was designed for outlier detection is introduced to select a reliable solution. Empirical results on six dynamic optimization benchmark problems demonstrate the effectiveness of the proposed algorithm compared with four state-of-the-art data-driven optimization algorithms.


title: Risk Assessment, Threat Modeling and Security Testing in SDLC

abstract: The software development process is considered as one of the key guidelines in the creation of said software and this approach is necessary for providing a more efficient yet satisfactory output. Without separation of work into distinct stages, it may lead to many delays and inefficiency of the project process where this disorganization can directly affect the product quality and reliability. Moreover, with this methodology established as the standard for any project, there are bound to be missteps specifically in regard to the involvement of security due to the lack of awareness. Therefore, the aim of this research is to identify and elaborate the findings and understanding of the security integrated into the process of software development as well as the related individual roles in ensuring that this security is maintained. Through thorough analysis and review of literature, an effort has been made through this paper to showcase the correct processes and ways for securing the software development process. At the same time, certain issues that pertain to this subject have been discussed together with proposing appropriate solutions. Furthermore, in depth discussion is carried out regarding methods such as security testing, risk assessment, threat modeling and other techniques that are able to create a more secure environment and systematic approach in a software development process.


title: Driven dynamics of a quantum dot electron spin coupled to bath of higher-spin nuclei

abstract: The interplay of optical driving and hyperfine interaction between an electron confined in a quantum dot and its surrounding nuclear spin environment produces a range of interesting physics such as mode-locking. In this work, we go beyond the ubiquitous spin 1/2 approximation for nuclear spins and present a comprehensive theoretical framework for an optically driven electron spin in a self-assembled quantum dot coupled to a nuclear spin bath of arbitrary spin. Using a dynamical mean-field approach, we compute the nuclear spin polarization distribution with and without the quadrupolar coupling. We find that while hyperfine interactions drive dynamic nuclear polarization and mode-locking, quadrupolar couplings counteract these effects. The tension between these mechanisms is imprinted on the steady-state electron spin evolution, providing a way to measure the importance of quadrupolar interactions in a quantum dot. Our results show that higher-spin effects such as quadrupolar interactions can have a significant impact on the generation of dynamic nuclear polarization and how it influences the electron spin evolution.


title: Trustworthy Preference Completion in Social Choice

abstract: As from time to time it is impractical to ask agents to provide linear orders over all alternatives, for these partial rankings it is necessary to conduct preference completion. Specifically, the personalized preference of each agent over all the alternatives can be estimated with partial rankings from neighboring agents over subsets of alternatives. However, since the agents' rankings are nondeterministic, where they may provide rankings with noise, it is necessary and important to conduct the trustworthy preference completion. Hence, in this paper firstly, a trust-based anchor-kNN algorithm is proposed to find $k$-nearest trustworthy neighbors of the agent with trust-oriented Kendall-Tau distances, which will handle the cases when an agent exhibits irrational behaviors or provides only noisy rankings. Then, for alternative pairs, a bijection can be built from the ranking space to the preference space, and its certainty and conflict can be evaluated based on a well-built statistical measurement Probability-Certainty Density Function. Therefore, a certain common voting rule for the first $k$ trustworthy neighboring agents based on certainty and conflict can be taken to conduct the trustworthy preference completion. The properties of the proposed certainty and conflict have been studied empirically, and the proposed approach has been experimentally validated compared to state-of-arts approaches with several data sets.


title: Two-dimensional Particle-in-Cell simulations of axisymmetric black hole magnetospheres

abstract: We investigate the temporal evolution of an axisymmetric magnetosphere around a rapidly rotating, stellar-mass black hole, applying a two-dimensional particle-in-cell simulation scheme. Adopting a homogeneous pair production, and assuming that the mass accretion rate is much less than the Eddington limit, we find that the black hole's rotational energy is preferentially extracted from the middle latitudes, and that this outward energy flux exhibits an enhancement that lasts approximately 160 dynamical time scales. It is demonstrated that the magnetohydrodynamic approximations cannot be justified in such a magnetically-dominated magnetosphere, because the Ohm's law completely breaks down, and because the charge-separated electron-positron plasmas are highly non-neutral. An implication is given regarding the collimation of relativistic jets.


title: A Feature Weighted Mixed Naive Bayes Model for Monitoring Anomalies in the Fan System of a Thermal Power Plant

abstract: With the increasing intelligence and integration, a great number of two-valued variables (generally stored in the form of 0 or 1 value) often exist in large-scale industrial processes. However, these variables cannot be effectively handled by traditional monitoring methods such as LDA, PCA and PLS. Recently, a mixed hidden naive Bayesian model (MHNBM) is developed for the first time to utilize both two-valued and continuous variables for abnormality monitoring. Although MHNBM is effective, it still has some shortcomings that need to be improved. For MHNBM, the variables with greater correlation to other variables have greater weights, which cannot guarantee greater weights are assigned to the more discriminating variables. In addition, the conditional probability must be computed based on the historical data. When the training data is scarce, the conditional probability between continuous variables tends to be uniformly distributed, which affects the performance of MHNBM. Here a novel feature weighted mixed naive Bayes model (FWMNBM) is developed to overcome the above shortcomings. For FWMNBM, the variables that are more correlated to the class have greater weights, which makes the more discriminating variables contribute more to the model. At the same time, FWMNBM does not have to calculate the conditional probability between variables, thus it is less restricted by the number of training data samples. Compared with MHNBM, FWMNBM has better performance, and its effectiveness is validated through the numerical cases of a simulation example and a practical case of Zhoushan thermal power plant (ZTPP), China.


title: Theoretical Analyses of Multiobjective Evolutionary Algorithms on Multimodal Objectives

abstract: The theoretical understanding of MOEAs is lagging far behind their success in practice. In particular, previous theory work considers mostly easy problems that are composed of unimodal objectives. As a first step towards a deeper understanding of how evolutionary algorithms solve multimodal multiobjective problems, we propose the OJZJ problem, a bi-objective problem composed of two objectives isomorphic to the classic jump function benchmark. We prove that SEMO with probability one does not compute the full Pareto front, regardless of the runtime. In contrast, for all problem sizes $n$ and all jump sizes ${k \in [4..\frac n2 - 1]}$, the global SEMO (GSEMO) covers the Pareto front in an expected number of $\Theta((n-2k)n^{k})$ iterations. For $k = o(n)$, we also show the tighter bound $\frac 32 e n^{k+1} \pm o(n^{k+1})$, which might be the first runtime bound for an MOEA that is tight apart from lower-order terms. We also combine the GSEMO with two approaches that showed advantages in single-objective multimodal problems. When using the GSEMO with a heavy-tailed mutation operator, the expected runtime improves by a factor of at least $k^{\Omega(k)}$. When adapting the recent stagnation-detection strategy of Rajabi and Witt (2022) to the GSEMO, the expected runtime also improves by a factor of at least $k^{\Omega(k)}$ and surpasses the heavy-tailed GSEMO by a small polynomial factor in $k$. Via an experimental analysis, we show that these asymptotic differences are visible already for small problem sizes: A factor-$5$ speed-up from heavy-tailed mutation and a factor-$10$ speed-up from stagnation detection can be observed already for jump size~$4$ and problem sizes between $10$ and $50$. Overall, our results show that the ideas recently developed to aid single-objective evolutionary algorithms to cope with local optima can be effectively employed also in multiobjective optimization.


title: Construction of local A-packets

abstract: In this paper, we reformulate Moeglin's explicit construction of local A-packets of split odd special orthogonal groups and symplectic groups. By this reformulation together with results of the previous paper with Minguez, we can compute the A-packets explicitly. Also, we give a non-vanishing criterion of our parametrization, and an algorithm to compute certain derivatives. Finally, we prove a formula for the Aubert duals of irreducible representations of Arthur type.


title: Achieving Adversarial Robustness Requires An Active Teacher

abstract: A new understanding of adversarial examples and adversarial robustness is proposed by decoupling the data generator and the label generator (which we call the teacher). In our framework, adversarial robustness is a conditional concept---the student model is not absolutely robust, but robust with respect to the teacher. Based on the new understanding, we claim that adversarial examples exist because the student cannot obtain sufficient information of the teacher from the training data. Various ways of achieving robustness is compared. Theoretical and numerical evidence shows that to efficiently attain robustness, a teacher that actively provides its information to the student may be necessary.


title: Regularity of fractional heat semigroup associated with Schr\"odinger operators

abstract: Let $L=-\Delta+V$ be a Schr\"odinger operator, where the potential $V$ belongs to the reverse H\"older class. By the subordinative formula, we introduce the fractional heat semigroup $\{e^{-t{L}^\alpha}\}_{t>0}, \alpha>0$, associated with ${L}$. By the aid of the fundamental solution of the heat equation: $$\partial_{t}u+L u=\partial_{t}u -\Delta u+Vu=0,$$ we estimate the gradient and the time-fractional derivatives of the fractional heat kernel $K^{L}_{\alpha,t}(\cdot, \cdot)$, respectively. This method is independent of the Fourier transform, and can be applied to the second order differential operators whose heat kernels satisfying Gaussian upper bounds. As an application, we establish a Carleson measure characterization of the Campanato type space $BMO^{\gamma}_{L}(\mathbb{R}^{n})$ via $\{e^{-t{L}^\alpha}\}_{t>0}$.


title: Fractional 0-1 programming and submodularity

abstract: In this note we study multiple-ratio fractional 0--1 programs, a broad class of NP-hard combinatorial optimization problems. In particular, under some relatively mild assumptions we provide a complete characterization of the conditions, which ensure that a single-ratio function is submodular. Then we illustrate our theoretical results with the assortment optimization and facility location problems, and discuss practical situations that guarantee submodularity in the considered application settings. In such cases, near-optimal solutions for multiple-ratio fractional 0--1 programs can be found via simple greedy algorithms.


title: Multi-Domain Multi-Task Rehearsal for Lifelong Learning

abstract: Rehearsal, seeking to remind the model by storing old knowledge in lifelong learning, is one of the most effective ways to mitigate catastrophic forgetting, i.e., biased forgetting of previous knowledge when moving to new tasks. However, the old tasks of the most previous rehearsal-based methods suffer from the unpredictable domain shift when training the new task. This is because these methods always ignore two significant factors. First, the Data Imbalance between the new task and old tasks that makes the domain of old tasks prone to shift. Second, the Task Isolation among all tasks will make the domain shift toward unpredictable directions; To address the unpredictable domain shift, in this paper, we propose Multi-Domain Multi-Task (MDMT) rehearsal to train the old tasks and new task parallelly and equally to break the isolation among tasks. Specifically, a two-level angular margin loss is proposed to encourage the intra-class/task compactness and inter-class/task discrepancy, which keeps the model from domain chaos. In addition, to further address domain shift of the old tasks, we propose an optional episodic distillation loss on the memory to anchor the knowledge for each old task. Experiments on benchmark datasets validate the proposed approach can effectively mitigate the unpredictable domain shift.


title: Accurate Cell Segmentation in Digital Pathology Images via Attention Enforced Networks

abstract: Automatic cell segmentation is an essential step in the pipeline of computer-aided diagnosis (CAD), such as the detection and grading of breast cancer. Accurate segmentation of cells can not only assist the pathologists to make a more precise diagnosis, but also save much time and labor. However, this task suffers from stain variation, cell inhomogeneous intensities, background clutters and cells from different tissues. To address these issues, we propose an Attention Enforced Network (AENet), which is built on spatial attention module and channel attention module, to integrate local features with global dependencies and weight effective channels adaptively. Besides, we introduce a feature fusion branch to bridge high-level and low-level features. Finally, the marker controlled watershed algorithm is applied to post-process the predicted segmentation maps for reducing the fragmented regions. In the test stage, we present an individual color normalization method to deal with the stain variation problem. We evaluate this model on the MoNuSeg dataset. The quantitative comparisons against several prior methods demonstrate the superiority of our approach.


title: Misspecified Beliefs about Time Lags

abstract: We examine the long-term behavior of a Bayesian agent who has a misspecified belief about the time lag between actions and feedback, and learns about the payoff consequences of his actions over time. Misspecified beliefs about time lags result in attribution errors, which have no long-term effect when the agent's action converges, but can lead to arbitrarily large long-term inefficiencies when his action cycles. Our proof uses concentration inequalities to bound the frequency of action switches, which are useful to study learning problems with history dependence. We apply our methods to study a policy choice game between a policy-maker who has a correctly specified belief about the time lag and the public who has a misspecified belief.


title: Effects of Semitransparent Window AspectRatio on Interaction of Collimated Beam withNatural Convection: Part I

abstract: The effects of the semitransparent winodw's aspect ratio on the interaction of the collimated beam with natural convection have been investigated numerically in the present work. The combination of geometrical parameters of the semitransparent window, i.e., height ratio ($h_r$) and window width ratio ($w_r$) and Planck numbers of the medium have been considered. The other parameters, like flow parameter (Ra$=10^5$), fluid parameter (Pr=0.71), thermal parameter (N), Irradiation (G=1000 $W/m^2$), Angle of incidence ($\phi=135^0$) and geometrical parameter of the geometry ($A_r$=1) and the wall conditions have been kept constant. A collimated beam is irradiated with irradiation value (G=1000 $W/m^2$) on the semitransparent window at an azimuthal angle ($\phi) 135^0$. The cavity is convectively heated from the bottom with heat transfer coefficient 50 $W/m^2 K$ and free stream temperature 305 $K$. A semitransparent window is created on the left wall and isothermal conditions (T=296 $K$) is applied on the semitransparent, left and right vertical walls, wherein adiabatic conditions are applied on upper wall of the cavity. The dynamics of two vortices inside the cavity change considerably by combinations these semitransparent window's aspect ratio and Planck number (Pl) of the medium. The left vortex breaks into two parts and remains confined in upper and lower left corners for some combination of aspect ratios and Planck numbers of the medium. The thermal plume flickers depending on the situation of dynamics of two vortices inside the cavity. The localized hating of the fluid happens mostly for large height ratio of semitransparent window. The conduction; radiation and total Nusselt number are also greatly affected by the semitransparent window's aspect ratio and the Planck number of the medium.


title: Boundedness of differential transforms for fractional Poisson type operators generated by parabolic operators

abstract: In this paper we analyze the convergence of the following type of series $$ T_N^\alpha f(x,t)=\sum_{j=N_1}^{N_2} v_j(P_{a_{j+1}}^\alpha f(x,t)-P_{a_j}^\alpha f(x,t)),\quad (x,t)\in \mathbb R^{n+1}, \ N=(N_1, N_2)\in \mathbb Z^2,\ \alpha>0, $$ where $\{P_{\tau}^\alpha \}_{\tau>0}$ is the fractional Poisson-type operators generated by the parabolic operator $L=\partial_t-\Delta$ with $\Delta$ being the classical Laplacian, $\{v_j\}_{j\in \mathbb Z}$ a bounded real sequences and $\{a_j\}_{j\in \mathbb Z}$ an increasing real sequence. Our analysis will consist {of} the boundedness, in $L^p(\mathbb{R}^n)$ and in $BMO(\mathbb{R}^n)$, of the operators $T^{\alpha}_N$ and its maximal operator $ T^*f(x)= \sup_{N\in \mathbb Z^2} |T^{\alpha}_N f(x)|.$ It is also shown that the local size of the maximal differential transform operators is the same with the order of a singular integral for functions $f$ having local support. Moreover, if $\{v_j\}_{j\in \mathbb Z}\in \ell^p(\mathbb Z)$, we get an intermediate size between the local size of singular integrals and Hardy-Littlewood maximal operator.


title: Deep Optimized Priors for 3D Shape Modeling and Reconstruction

abstract: Many learning-based approaches have difficulty scaling to unseen data, as the generality of its learned prior is limited to the scale and variations of the training samples. This holds particularly true with 3D learning tasks, given the sparsity of 3D datasets available. We introduce a new learning framework for 3D modeling and reconstruction that greatly improves the generalization ability of a deep generator. Our approach strives to connect the good ends of both learning-based and optimization-based methods. In particular, unlike the common practice that fixes the pre-trained priors at test time, we propose to further optimize the learned prior and latent code according to the input physical measurements after the training. We show that the proposed strategy effectively breaks the barriers constrained by the pre-trained priors and could lead to high-quality adaptation to unseen data. We realize our framework using the implicit surface representation and validate the efficacy of our approach in a variety of challenging tasks that take highly sparse or collapsed observations as input. Experimental results show that our approach compares favorably with the state-of-the-art methods in terms of both generality and accuracy.


title: Neighbors From Hell: Voltage Attacks Against Deep Learning Accelerators on Multi-Tenant FPGAs

abstract: Field-programmable gate arrays (FPGAs) are becoming widely used accelerators for a myriad of datacenter applications due to their flexibility and energy efficiency. Among these applications, FPGAs have shown promising results in accelerating low-latency real-time deep learning (DL) inference, which is becoming an indispensable component of many end-user applications. With the emerging research direction towards virtualized cloud FPGAs that can be shared by multiple users, the security aspect of FPGA-based DL accelerators requires careful consideration. In this work, we evaluate the security of DL accelerators against voltage-based integrity attacks in a multitenant FPGA scenario. We first demonstrate the feasibility of such attacks on a state-of-the-art Stratix 10 card using different attacker circuits that are logically and physically isolated in a separate attacker role, and cannot be flagged as malicious circuits by conventional bitstream checkers. We show that aggressive clock gating, an effective power-saving technique, can also be a potential security threat in modern FPGAs. Then, we carry out the attack on a DL accelerator running ImageNet classification in the victim role to evaluate the inherent resilience of DL models against timing faults induced by the adversary. We find that even when using the strongest attacker circuit, the prediction accuracy of the DL accelerator is not compromised when running at its safe operating frequency. Furthermore, we can achieve 1.18-1.31x higher inference performance by over-clocking the DL accelerator without affecting its prediction accuracy.


title: Collisions of false-vacuum bubble walls in a quantum spin chain

abstract: We study the real-time dynamics of a small bubble of "false vacuum" in a quantum spin chain near criticality, where the low-energy physics is described by a relativistic (1+1)-dimensional quantum field theory. Such a bubble can be thought of as a confined kink-antikink pair (a meson). We carefully construct bubbles so that particle production does not occur until the walls collide. To achieve this in the presence of strong correlations, we extend a Matrix Product State (MPS) ansatz for quasiparticle wavepackets [Van Damme et al., arXiv:1907.02474 (2019)] to the case of confined, topological quasiparticles. By choosing the wavepacket width and the bubble size appropriately, we avoid strong lattice effects and observe relativistic kink-antikink collisions. We use the MPS quasiparticle ansatz to detect scattering outcomes: In the Ising model, with transverse and longitudinal fields, we do not observe particle production despite nonintegrability (supporting recent observations of nonthermalizing mesonic states). With additional interactions, we see production of confined and unconfined particle pairs. Although we simulated these low-energy, few-particle events with moderate resources, we observe significant growth of entanglement with energy and with the number of collisions, suggesting that increasing either will ultimately exhaust our methods. Quantum devices, in contrast, are not limited by entanglement production, and promise to allow us to go far beyond classical methods. We anticipate that kink-antikink scattering in 1+1 dimensions will be an instructive benchmark problem for relatively near-term quantum devices.


title: Bayesian Neural Ordinary Differential Equations

abstract: Recently, Neural Ordinary Differential Equations has emerged as a powerful framework for modeling physical simulations without explicitly defining the ODEs governing the system, but instead learning them via machine learning. However, the question: "Can Bayesian learning frameworks be integrated with Neural ODE's to robustly quantify the uncertainty in the weights of a Neural ODE?" remains unanswered. In an effort to address this question, we primarily evaluate the following categories of inference methods: (a) The No-U-Turn MCMC sampler (NUTS), (b) Stochastic Gradient Hamiltonian Monte Carlo (SGHMC) and (c) Stochastic Langevin Gradient Descent (SGLD). We demonstrate the successful integration of Neural ODEs with the above Bayesian inference frameworks on classical physical systems, as well as on standard machine learning datasets like MNIST, using GPU acceleration. On the MNIST dataset, we achieve a posterior sample accuracy of 98.5% on the test ensemble of 10,000 images. Subsequently, for the first time, we demonstrate the successful integration of variational inference with normalizing flows and Neural ODEs, leading to a powerful Bayesian Neural ODE object. Finally, considering a predator-prey model and an epidemiological system, we demonstrate the probabilistic identification of model specification in partially-described dynamical systems using universal ordinary differential equations. Together, this gives a scientific machine learning tool for probabilistic estimation of epistemic uncertainties.


title: Deep Portfolio Optimization via Distributional Prediction of Residual Factors

abstract: Recent developments in deep learning techniques have motivated intensive research in machine learning-aided stock trading strategies. However, since the financial market has a highly non-stationary nature hindering the application of typical data-hungry machine learning methods, leveraging financial inductive biases is important to ensure better sample efficiency and robustness. In this study, we propose a novel method of constructing a portfolio based on predicting the distribution of a financial quantity called residual factors, which is known to be generally useful for hedging the risk exposure to common market factors. The key technical ingredients are twofold. First, we introduce a computationally efficient extraction method for the residual information, which can be easily combined with various prediction algorithms. Second, we propose a novel neural network architecture that allows us to incorporate widely acknowledged financial inductive biases such as amplitude invariance and time-scale invariance. We demonstrate the efficacy of our method on U.S. and Japanese stock market data. Through ablation experiments, we also verify that each individual technique contributes to improving the performance of trading strategies. We anticipate our techniques may have wide applications in various financial problems.


title: Upper Atmosphere Smoke Injection from Large Areal Fires

abstract: Large fires can inject smoke into the upper troposphere and lower stratosphere. Detailed fire simulations allow for assessment of how local weather interacts with these fires and affects smoke lofting. In this study, we employ the fire simulation package in the Weather Research and Forecasting model (WRF-Fire), Version 4.0.1, to explore how smoke lofting from a fire burning a homogeneous fuel bed changes with varying local winds, relative humidity, and atmospheric boundary-layer stability for two different-sized areal fires. We find that the presence of moisture has the greatest influence on the results by raising the altitude of lofting, while faster winds speeds dampen lofting and lower the injection height. Stably stratified conditions further inhibit plume propagation compared with neutrally stratified conditions.


title: Complexes, Graphs, Homotopy, Products and Shannon Capacity

abstract: A finite abstract simplicial complex G defines the Barycentric refinement graph phi(G) = (G,{ (a,b), a subset b or b subset a }) and the connection graph psi(G) = (G,{ (a,b), a intersected with b not empty }). We note here that both functors phi and psi from complexes to graphs are invertible on the image (Theorem 1) and that G,phi(G),psi(G) all have the same automorphism group and that the Cartesian product of G corresponding to the Stanley-Reisner product of phi(G) and the strong Shannon product of psi(G), have the product automorphism groups. Second, we see that if G is a Barycentric refinement, then phi(G) and psi(G) are graph homotopic (Theorem 2). Third, if gamma is the geometric realization functor, assigning to a complex or to a graph the geometric realization of its clique complex, then gamma(G) and gamma(phi(G)) and gamma(psi(G)) are all classically homotopic for a Barycentric refined simplicial complex G (Theorem 3). The Barycentric assumption is necessary in Theorem 2 and 3. There is compatibility with Cartesian products of complexes which manifests in the strong graph product of connection graphs: if two graphs A,A' are homotopic and B,B' are homotopic, then A . B is homotopic to A' . B' (Theorem 4) leading to a commutative ring of homotopy classes of graphs. Finally, we note (Theorem 5) that for all simplicial complexes G as well as product G=G_1 x G_2 ... x G_k, the Shannon capacity Theta(psi(G)) of psi(G) is equal to the number m of zero-dimensional sets in G. An explicit Lowasz umbrella in R^m leads to the Lowasz number theta(G) leq m and so Theta(psi(G))=theta(psi(G))=m making Theta compatible with disjoint union addition and strong multiplication.


title: TDAF: Top-Down Attention Framework for Vision Tasks

abstract: Human attention mechanisms often work in a top-down manner, yet it is not well explored in vision research. Here, we propose the Top-Down Attention Framework (TDAF) to capture top-down attentions, which can be easily adopted in most existing models. The designed Recursive Dual-Directional Nested Structure in it forms two sets of orthogonal paths, recursive and structural ones, where bottom-up spatial features and top-down attention features are extracted respectively. Such spatial and attention features are nested deeply, therefore, the proposed framework works in a mixed top-down and bottom-up manner. Empirical evidence shows that our TDAF can capture effective stratified attention information and boost performance. ResNet with TDAF achieves 2.0% improvements on ImageNet. For object detection, the performance is improved by 2.7% AP over FCOS. For pose estimation, TDAF improves the baseline by 1.6%. And for action recognition, the 3D-ResNet adopting TDAF achieves improvements of 1.7% accuracy.


title: Explicit Formulas for the First Form (q,r)-Dowling Numbers and (q,r)-Whitney-Lah Numbers

abstract: In this paper, a q-analogue of r-Whitney-Lah numbers, also known as (q,r)-Whitney-Lah number, denoted by $L_{m,r}[n,k]_q$ is defined using the triangular recurrence relation. Several fundamental properties for the q-analogue are established such as vertical and horizontal recurrence relations, horizontal and exponential generating functions. Moreover, an explicit formula for (q,r)-Whitney-Lah number is derived using the concept of q-difference operator, particularly, the q-analogue of Newton's Interpolation Formula (the umbral version of Taylor series). Furthermore, an explicit formula for the first form (q,r)-Dowling numbers is obtained which is expressed in terms of (q,r)-Whitney-Lah numbers and (q,r)-Whitney numbers of the second kind.


title: Sub-Nyquist computational ghost imaging with orthonormalized colored noise pattern

abstract: Computational ghost imaging generally requires a large number of pattern illumination to obtain a high-quality image. The colored noise speckle pattern was recently proposed to substitute the white noise pattern in a variety of noisy environments and gave a significant signal-to-noise ratio enhancement even with a limited number of patterns. We propose and experimentally demonstrate here an orthonormalization approach based on the colored noise patterns to achieve sub-Nyquist computational ghost imaging. We tested the reconstructed image in quality indicators such as the contrast-to-noise ratio, the mean square error, the peak signal to noise ratio, and the correlation coefficient. The results suggest that our method can provide high-quality images while using a sampling ratio an order lower than the conventional methods.


title: Simulating the radiation loss of superconducting submillimeter wave filters and transmission lines using Sonnet EM

abstract: Superconducting resonators and transmission lines are fundamental building blocks of integrated circuits for millimeter-submillimeter astronomy. Accurate simulation of radiation loss from the circuit is crucial for the design of these circuits because radiation loss increases with frequency, and can thereby deteriorate the system performance. Here we show a stratification for a 2.5-dimensional method-of-moment simulator Sonnet EM that enables accurate simulations of the radiative resonant behavior of submillimeter-wave coplanar resonators and straight coplanar waveguides (CPWs). The Sonnet simulation agrees well with the measurement of the transmission through a coplanar resonant filter at 374.6 GHz. Our Sonnet stratification utilizes artificial lossy layers below the lossless substrate to absorb the radiation, and we use co-calibrated internal ports for de-embedding. With this type of stratification, Sonnet can be used to model superconducting millimeter-submillimeter wave circuits even when radiation loss is a potential concern.


title: Few Shot Adaptive Normalization Driven Multi-Speaker Speech Synthesis

abstract: The style of the speech varies from person to person and every person exhibits his or her own style of speaking that is determined by the language, geography, culture and other factors. Style is best captured by prosody of a signal. High quality multi-speaker speech synthesis while considering prosody and in a few shot manner is an area of active research with many real-world applications. While multiple efforts have been made in this direction, it remains an interesting and challenging problem. In this paper, we present a novel few shot multi-speaker speech synthesis approach (FSM-SS) that leverages adaptive normalization architecture with a non-autoregressive multi-head attention model. Given an input text and a reference speech sample of an unseen person, FSM-SS can generate speech in that person's style in a few shot manner. Additionally, we demonstrate how the affine parameters of normalization help in capturing the prosodic features such as energy and fundamental frequency in a disentangled fashion and can be used to generate morphed speech output. We demonstrate the efficacy of our proposed architecture on multi-speaker VCTK and LibriTTS datasets, using multiple quantitative metrics that measure generated speech distortion and MoS, along with speaker embedding analysis of the generated speech vs the actual speech samples.


title: Characterizations of complete stabilizability

abstract: We present several characterizations, via some weak observability inequalities, on the complete stabilizability for a control system $[A,B]$, i.e., $y'(t)=Ay(t)+Bu(t)$, $t\geq 0$, where $A$ generates a $C_0$-semigroup on a Hilbert space $X$ and $B$ is a linear and bounded operator from another Hilbert space $U$ to $X$. We then extend the aforementioned characterizations in two directions: first, the control operator $B$ is unbounded; second, the control system is time-periodic. We also give some sufficient conditions, from the perspective of the spectral projections, to ensure the weak observability inequalities. As applications, we provide several examples, which are not null controllable, but can be verified, via the weak observability inequalities, to be completely stabilizable.


title: A sharp rise in the detection rate of broad absorption line variations in a quasar SDSS J141955.26+522741.1

abstract: We present an analysis of the variability of broad absorption lines (BALs) in a quasar SDSS J141955.26+522741.1 at $z=2.145$ with 72 observations from the Sloan Digital Sky Survey Data Release 16 (SDSS DR16). The strong correlation between the equivalent widths of BAL and the continuum luminosity, reveals that the variation of BAL trough is dominated by the photoionization. The photoionization model predicts that when the time interval $\Delta T$ between two observations is longer than the recombination timescale $t_{\rm rec}$, the BAL variations can be detected. This can be characterized as a "sharp rise" in the detection rate of BAL variation at $\Delta T=t_{\rm rec}$. For the first time, we detect such a "sharp rise" signature in the detection rate of BAL variations. As a result, we propose that the $t_{\rm rec}$ can be obtained from the "sharp rise" of the detection rate of BAL variation. It is worth mentioning that the BAL variations are detected at the time-intervals less than the $t_{\rm rec}$ for half an order of magnitude in two individual troughs. This result indicates that there may be multiple components with different $t_{\rm rec}$ but the same velocity in an individual trough.


title: Critical condition for electrowetting-induced detachment of a droplet from a curved surface

abstract: Based on energy conservation, we derive a critical condition theoretically for electrowettinginduced droplet detachment from a hydrophobic curved surface. Phase diagrams are constructed in terms of droplet volume, viscosity, Ohnesorge number, friction coefficient at contact line, surface curvature, surface wettability and electrowetting number. The deduced critical condition offers a general and quantitative prediction on when the detachment occurs, a criterion enabling us to gain more insights into how to accurately manipulate the electrowetting-induced detachment of an aqueous droplet on a curved surface. The results obtained in this paper also imply that the detachable regimes of the phase diagrams can be enlarged through increasing droplet volume and surface curvature, and reducing liquid viscosity, friction coefficient, Ohnesorge number and wettability of substrates.


title: Recursive computation of the Hawkes cumulants

abstract: We propose a recursive method for the computation of the cumulants of self-exciting point processes of Hawkes type, based on standard combinatorial tools such as Bell polynomials. This closed-form approach is easier to implement on higher-order cumulants in comparison with existing methods based on differential equations, tree enumeration or martingale arguments. The results are corroborated by Monte Carlo simulations, and also apply to the computation of joint cumulants generated by multidimensional self-exciting processes.


title: A Visual Mining Approach to Improved Multiple-Instance Learning

abstract: Multiple-instance learning (MIL) is a paradigm of machine learning that aims to classify a set (bag) of objects (instances), assigning labels only to the bags. This problem is often addressed by selecting an instance to represent each bag, transforming a MIL problem into standard supervised learning. Visualization can be a useful tool to assess learning scenarios by incorporating the users' knowledge into the classification process. Considering that multiple-instance learning is a paradigm that cannot be handled by current visualization techniques, we propose a multiscale tree-based visualization called MILTree to support MIL problems. The first level of the tree represents the bags, and the second level represents the instances belonging to each bag, allowing users to understand the MIL datasets in an intuitive way. In addition, we propose two new instance selection methods for MIL, which help users improve the model even further. Our methods can handle both binary and multiclass scenarios. In our experiments, SVM was used to build the classifiers. With support of the MILTree layout, the initial classification model was updated by changing the training set, which is composed of the prototype instances. Experimental results validate the effectiveness of our approach, showing that visual mining by MILTree can support exploring and improving models in MIL scenarios and that our instance selection methods outperform the currently available alternatives in most cases.


title: Any multi-index sequence has an interpolating measure

abstract: R. P. Boas showed that any single-index sequence $\left\{ \beta_i \right\}_{i=0}^\infty$ of real numbers can be represented as $\beta_i =\int_0^\infty x^i \, d\mu$ ($i=0,1,2,\ldots$), where $\mu$ is a signed measure. As Boas said his observation seemed to be quite unexpected; however, it is even possible to extend the result to any multi-index sequence of real numbers. In addition, we can also prove that any multi-index finite sequence admits a measure of a similar type.


title: AndroEvolve: Automated Update for Android Deprecated-API Usages

abstract: Android operating system (OS) is often updated, where each new version may involve API deprecation. Usages of deprecated APIs in Android apps need to be updated to ensure the apps' compatibility with the old and new versions of Android OS. In this work, we propose AndroEvolve, an automated tool to update usages of deprecated Android APIs, that addresses the limitations of the state-of-the-art tool, CocciEvolve. AndroEvolve utilizes data flow analysis to solve the problem of out-of-method-boundary variables, and variable denormalization to remove the temporary variables introduced by CocciEvolve. We evaluated the accuracy of AndroEvolve using a dataset of 360 target files and 20 deprecated Android APIs, where AndroEvolve is able to produce 319 correct updates, compared to CocciEvolve which only produces 249 correct updates. We also evaluated the readability of AndroEvolve's update results using a manual and an automatic evaluation. Both evaluations demonstrated that the code produced by AndroEvolve has higher readability than CocciEvolve's. A video demonstration of AndroEvolve is available at https://youtu.be/siU0tuMITXI.


title: Predicting Berborite as Potential Deep-Ultraviolet Nonlinear Optical Crystal from First Principles

abstract: Following our ab initio nonlinear optical (NLO) materials design guidelines, in this Letter, we discovered a novel type of structure to realize potential deep-ultraviolet (DUV) NLO performance in the classical beryllium borate system. By densely stacking the NLO-active layered frameworks, the key design scheme for the structural evolution from the (Be2BO3F2) layers in KBe2BO3F2 (KBBF) to the novel (Be2BO5H3) layers in berborite is illustrated. Based on available experimental results and systematical theoretical evaluation from first principles, the NLO properties of berborite are further obtained as comparable as the only pratical DUV NLO crystal KBBF. It is demonstrated that berborite can achieve available DUV phase-matched output with strong NLO effect for the practically important 177.3 nm and 193.7 nm lasers. Once obtained with sizable single crystal, it can be applied as a promising DUV NLO crystal.


title: OCTA-500: A Retinal Dataset for Optical Coherence Tomography Angiography Study

abstract: Optical coherence tomography angiography (OCTA) is a novel imaging modality that has been widely utilized in ophthalmology and neuroscience studies to observe retinal vessels and microvascular systems. However, publicly available OCTA datasets remain scarce. In this paper, we introduce the largest and most comprehensive OCTA dataset dubbed OCTA-500, which contains OCTA imaging under two fields of view (FOVs) from 500 subjects. The dataset provides rich images and annotations including two modalities (OCT/OCTA volumes), six types of projections, four types of text labels (age / gender / eye / disease) and seven types of segmentation labels (large vessel/capillary/artery/vein/2D FAZ/3D FAZ/retinal layers). Then, we propose a multi-object segmentation task called CAVF, which integrates capillary segmentation, artery segmentation, vein segmentation, and FAZ segmentation under a unified framework. In addition, we optimize the 3D-to-2D image projection network (IPN) to IPN-V2 to serve as one of the segmentation baselines. Experimental results demonstrate that IPN-V2 achieves an ~10% mIoU improvement over IPN on CAVF task. Finally, we further study the impact of several dataset characteristics: the training set size, the model input (OCT/OCTA, 3D volume/2D projection), the baseline networks, and the diseases. The dataset and code are publicly available at: https://ieee-dataport.org/open-access/octa-500.


title: Learning Hybrid Representations for Automatic 3D Vessel Centerline Extraction

abstract: Automatic blood vessel extraction from 3D medical images is crucial for vascular disease diagnoses. Existing methods based on convolutional neural networks (CNNs) may suffer from discontinuities of extracted vessels when segmenting such thin tubular structures from 3D images. We argue that preserving the continuity of extracted vessels requires to take into account the global geometry. However, 3D convolutions are computationally inefficient, which prohibits the 3D CNNs from sufficiently large receptive fields to capture the global cues in the entire image. In this work, we propose a hybrid representation learning approach to address this challenge. The main idea is to use CNNs to learn local appearances of vessels in image crops while using another point-cloud network to learn the global geometry of vessels in the entire image. In inference, the proposed approach extracts local segments of vessels using CNNs, classifies each segment based on global geometry using the point-cloud network, and finally connects all the segments that belong to the same vessel using the shortest-path algorithm. This combination results in an efficient, fully-automatic and template-free approach to centerline extraction from 3D images. We validate the proposed approach on CTA datasets and demonstrate its superior performance compared to both traditional and CNN-based baselines.


title: Screw dislocations in the X-cube fracton model

abstract: The X-cube model, a prototypical gapped fracton model, has been shown to have a foliation structure. That is, inside the 3+1D model, there are hidden layers of 2+1D gapped topological states. A screw dislocation in a 3+1D lattice can often reveal nontrivial features associated with a layered structure. In this paper, we study the X-cube model on lattices with screw dislocations. In particular, we find that a screw dislocation results in a finite change in the logarithm of the ground state degeneracy of the model. Part of the change can be traced back to the effect of screw dislocations in a simple stack of 2+1D topological states, hence corroborating the foliation structure in the model. The other part of the change comes from the induced motion of fractons or sub-dimensional excitations along the dislocation, a feature absent in the stack of 2+1D layers.


title: Interpretation of the apparent activation energy of glass transition

abstract: The temperature dependence of the viscosity of glass is a major concern in the field of glass research. Strong deviations from the Arrhenius law make the interpretation of the activation energy difficult. In the present study, a reasonable interpretation of the apparent activation energy is demonstrated along similar lines as those adopted in solid-state physics and chemistry. In contrast to the widely held view that phase transition in glass occurs at the reference temperature $T_{0}$ according to the Vogel--Fulcher--Tammann formula, in the present work the transition observed at the glass-transition temperature $T_{g}$ is regarded as a phase transition from the liquid to solid phases. A distinct feature of glass is that the energy barrier significantly changes in the transition range with width $\Delta T_{g}$. This change in the energy barrier alters the manner in which the apparent activation energy constitutes the Arrhenius form. Analysis of available experimental data showed that the actual energy barrier is significantly smaller than the apparent activation energy, and importantly, the values obtained were in the reasonable range of energy expected for chemical bonds. The overestimation of the apparent activation energy depends on the ratio $T_{g}/\Delta T_{g}$, which explains the existence of two types of glasses strong and fragile glasses. The fragility can be re-interpreted as an indication of the degree of increase in the energy barrier when approaching $T_{g}$ from high temperatures. Since no divergence in viscosity was observed below $T_{g}$, it is unlikely that a transition occurs at $T_{0}$.


title: Optimal Voltage and Current Control of an HVDC System to Improve Real-Time Frequency Regulation

abstract: High-voltage direct-current (HVDC) systems for constant or intermittent power delivery have recently been developed further to support grid frequency regulation (GFR). This paper proposes a new control strategy for a line-commutated converter-based (LCC) HVDC system, wherein the DC-link voltage and current are optimally regulated to improve real-time GFR in both rectifier- and inverter-side AC networks. A dynamic model of an LCC HVDC system is developed using the DC voltage and current as input variables, and is integrated with feedback loops for inertia emulation and droop control. A linear quadratic Gaussian (LQG) controller is also designed for optimal secondary frequency control, while mitigating conflict between the droop controllers of the HVDC converters. An eigenvalue analysis is then conducted, focusing on the effects of model parameters and controller gains on the proposed strategy. Simulation case studies are also performed using the Jeju-Haenam HVDC system as a test bed. The results of the case study confirm that the proposed strategy enables the HVDC system to improve GFR, in coordination with generators in both-side grids, by exploiting the fast dynamics of HVDC converters. The proposed strategy is also effective under various conditions for the LQG weighting coefficients, inertia emulation, and droop control.


title: MADLens, a python package for fast and differentiable non-Gaussian lensing simulations

abstract: We present MADLens a python package for producing non-Gaussian lensing convergence maps at arbitrary source redshifts with unprecedented precision. MADLens is designed to achieve high accuracy while keeping computational costs as low as possible. A MADLens simulation with only $256^3$ particles produces convergence maps whose power agree with theoretical lensing power spectra up to $L{=}10000$ within the accuracy limits of HaloFit. This is made possible by a combination of a highly parallelizable particle-mesh algorithm, a sub-evolution scheme in the lensing projection, and a machine-learning inspired sharpening step. Further, MADLens is fully differentiable with respect to the initial conditions of the underlying particle-mesh simulations and a number of cosmological parameters. These properties allow MADLens to be used as a forward model in Bayesian inference algorithms that require optimization or derivative-aided sampling. Another use case for MADLens is the production of large, high resolution simulation sets as they are required for training novel deep-learning-based lensing analysis tools. We make the MADLens package publicly available under a Creative Commons License (https://github.com/VMBoehm/MADLens).


title: Multi-SpectroGAN: High-Diversity and High-Fidelity Spectrogram Generation with Adversarial Style Combination for Speech Synthesis

abstract: While generative adversarial networks (GANs) based neural text-to-speech (TTS) systems have shown significant improvement in neural speech synthesis, there is no TTS system to learn to synthesize speech from text sequences with only adversarial feedback. Because adversarial feedback alone is not sufficient to train the generator, current models still require the reconstruction loss compared with the ground-truth and the generated mel-spectrogram directly. In this paper, we present Multi-SpectroGAN (MSG), which can train the multi-speaker model with only the adversarial feedback by conditioning a self-supervised hidden representation of the generator to a conditional discriminator. This leads to better guidance for generator training. Moreover, we also propose adversarial style combination (ASC) for better generalization in the unseen speaking style and transcript, which can learn latent representations of the combined style embedding from multiple mel-spectrograms. Trained with ASC and feature matching, the MSG synthesizes a high-diversity mel-spectrogram by controlling and mixing the individual speaking styles (e.g., duration, pitch, and energy). The result shows that the MSG synthesizes a high-fidelity mel-spectrogram, which has almost the same naturalness MOS score as the ground-truth mel-spectrogram.


title: Developing an Analytical Model of Frequency and Voltage Variations for Dynamic Reconfiguration

abstract: This paper develops a new analytical model to estimate real-time variations in grid frequency and voltages resulting from dynamic network reconfiguration (DNR). In the proposed model, switching operations are considered as discrete variations in an admittance matrix, leading to step variations in node injection currents. The network model with discrete admittance variations is then integrated with dynamic models of synchronous generators and voltage-dependent loads, enabling analysis of the dynamic grid operations initiated by the DNR. Case studies are performed to validate the proposed model via comparison with a conventional model and a comprehensive MATLAB/SIMULINK model.


title: Variational State and Parameter Estimation

abstract: This paper considers the problem of computing Bayesian estimates of both states and model parameters for nonlinear state-space models. Generally, this problem does not have a tractable solution and approximations must be utilised. In this work, a variational approach is used to provide an assumed density which approximates the desired, intractable, distribution. The approach is deterministic and results in an optimisation problem of a standard form. Due to the parametrisation of the assumed density selected first- and second-order derivatives are readily available which allows for efficient solutions. The proposed method is compared against state-of-the-art Hamiltonian Monte Carlo in two numerical examples.


title: Strichartz estimates in Wiener amalgam spaces and applications to nonlinear wave equations

abstract: In this paper we obtain some new Strichartz estimates for the wave propagator $e^{it\sqrt{-\Delta}}$ in the context of Wiener amalgam spaces. While it is well understood for the Schr\"odinger case, nothing is known about the wave propagator. This is because there is no such thing as an explicit formula for the integral kernel of the propagator unlike the Schr\"odinger case. To overcome this lack, we instead approach the kernel by rephrasing it as an oscillatory integral involving Bessel functions and then by carefully making use of cancellation in such integrals based on the asymptotic expansion of Bessel functions. Our approach can be applied to the Schr\"odinger case as well. We also obtain some corresponding retarded estimates to give applications to nonlinear wave equations where Wiener amalgam spaces as solution spaces can lead to a finer analysis of the local and global behavior of the solution.