title: Continuous spatial field confocal thermometry using lanthanide doped tellurite glass

abstract: Quantifying temperature variations at the micron scale can provide new opportunities in optical sensing. In this paper, we present a novel approach using the temperature-dependent variations in fluorescence of rare-earth doped tellurite glass to provide a micron-scale image of temperature variations over a 200 micrometre field of view. We demonstrate the system by monitoring the evaporation of a water droplet and report a net temperature change of 7.04 K with a sensitivity of at least 0.12 K. These results establish the practicality of this confocal-based approach to provide high-resolution marker-free optical temperature sensing.


title: Orbits around a Kerr black hole and its shadow

abstract: Since the full General Theory of Relativity has been unveiled to the scientific community in 1915, many solutions to the vacuum Einstein field equations have been found and studied. This paper aims at documenting exhaustively the derivation of the shape of the patch of the sky that is left completely black by a spinning black hole described by the Kerr solutions in Boyer-Lindquist coordinates. This dark zone in the observer's sky is called the black hole shadow. Conserved quantities which allow for the analysis of particle orbits are first introduced, with the help of which the trajectories of photons are uniquely described by two impact parameters (Specific angular momentum in the azimuthal direction L and dimensionless Carter's constant Q). We then derive the conditions on those parameters required for a photon to be captured by the black hole. These conditions are then translated into an equation for the black hole shadow. We conclude the paper by drawing out the black hole shadows for an equatorial observer for two separate cases.


title: Modulation of heterogeneous surface charge and flow pattern in electrically gated converging-diverging nanochannel

abstract: The present study aims at utilizing field effect phenomenon to induce heterogeneous surface charge and consequently changing the fluid flow in a solid state nanochannel with converging-diverging periodicity. It is shown that the proposed geometry causes non-uniform radial field adjacent to channel walls which is stronger around the diverging section and weaker next to the converging part of the wall. The later generates heterogeneous surface charge at channel walls depending on the applied gate potential i.e. applying low gate potential enables effective modulation of surface charge with the same polarity of the intrinsic charge at channel walls, while moderate gate potential causes charge inversion in diverging sections of the channel and generates reverse flow and thus results in fluid flow circulation. The potential application of flow circulation for trapping and rejection of particles is also demonstrated.

categories: physics.app-ph physics.bio-ph physics.chem-ph physics.flu-dyn


title: Bioinspired Nanocomposites: 2D Materials Within a 3D Lattice

abstract: Advanced composites are used in a variety of industrial applications and therefore attract much scientific interest. Here we describe the formation of novel carbon-based nanocomposites via incorporation of graphene oxide into the crystal lattice of single crystals of calcite. Incorporation of a 2D organic material into single-crystal lattices has never before been reported. To characterize the resulting nanocomposites we employed high-resolution synchrotron powder X-ray diffraction, electron microscopy, transmission electron microscopy, fluorescence microscopy, and nanoindentation tests. Detailed analysis revealed a layered distribution of graphene oxide sheets incorporated within the calcite host. Moreover, the optical and mechanical properties of the calcite host were altered when a carbon-based nanomaterial was introduced into its lattice. Compared to pure calcite, the composite GO-calcite crystals exhibited lower elastic modulus and higher hardness. The results of this study show that the incorporation of a 2D material within a 3D crystal lattice is not only feasible but also can lead to the formation of hybrid crystals exhibiting new properties.


title: Synchronization properties and reservoir computing capability of hexagonal spintronic oscillator arrays

abstract: The influence of array geometry on synchronization properties of a 2-D oscillator array is investigated based on a comparison between a rectangular and a hexagonal grid. The Kuramoto model is solved for a nearest-neighbor case with periodic boundary conditions and for a small-scale, realistic coupling case with 1/r^3 decay characteristic of spintronic oscillators. In both cases, it is found that the hexagonal grid choice leads to lower synchronization threshold and higher emission power than its rectangular counterpart, which results from increased connectivity, as well as, in the realistic-coupling case, from decreased contributions of the array edges. Additionally, a more general spin-torque oscillator model including both amplitude and phase as degrees of freedom is employed for reservoir computing simulations, showing that by using hexagonal grid one can increase the short-term memory capacity but not the parity-check capacity of the system.


title: Power Law Plateau Inflation Potential In The RS $II$ Braneworld Evading Swampland Conjecture

abstract: In the recent time, inflationary cosmology is facing an existential crisis due to the proposed Swampland criterion which aims to evade any (meta-)stable de Sitter construction within the String landscape. It is been realised that a single field slow roll inflation is inconsistent with the Swampland criterion unless the inflationary model in realised in some non standard scenario such as Warm inflation or the Braneworld scenario. In \cite{Dimopoulos}, Dimopoulos and Owen introduced a new class of model of inflation dubbed as the power law plateau inflation in the standard cold inflationary scenario. But to realise this model in the standard scenario consistent with observation, they had to introduce a phase of thermal inflation. In this paper we have analysed this model in the braneworld scenario to show that for some choice of the parameters defining the model class, one can have an observationally consistent power law plateau without any phase of thermal inflation. We have also shown that, for the correct choice of model parameters, one can easily satisfy the Swampland criterion. Besides, for a particular choice of the potential one can also satisfy the recently proposed Trans-Planckian Censorship Conjecture.


title: An optimal scheduling architecture for accelerating batch algorithms on Neural Network processor architectures

abstract: In neural network topologies, algorithms are running on batches of data tensors. The batches of data are typically scheduled onto the computing cores which execute in parallel. For the algorithms running on batches of data, an optimal batch scheduling architecture is very much needed by suitably utilizing hardware resources - thereby resulting in significant reduction training and inference time. In this paper, we propose to accelerate the batch algorithms for neural networks through a scheduling architecture enabling optimal compute power utilization. The proposed optimal scheduling architecture can be built into HW or can be implemented in SW alone which can be leveraged for accelerating batch algorithms. The results demonstrate that the proposed architecture speeds up the batch algorithms compared to the previous solutions. The proposed idea applies to any HPC architecture meant for neural networks.


title: Artificial-Noise-Aided Secure MIMO Wireless Communications via Intelligent Reflecting Surface

abstract: This paper considers a MIMO secure wireless communication system aided by the physical layer security technique of sending artificial noise (AN). To further enhance the system security performance, the advanced intelligent reflecting surface (IRS) is invoked in the AN-aided communication system, where the base station (BS), legitimate information receiver (IR) and eavesdropper (Eve) are equipped with multiple antennas. With the aim for maximizing the secrecy rate (SR), the transmit precoding (TPC) matrix at the BS, covariance matrix of AN and phase shifts at the IRS are jointly optimized subject to constrains of transmit power limit and unit modulus of IRS phase shifts. Then, the secrecy rate maximization (SRM) problem is formulated, which is a non-convex problem with multiple coupled variables. To tackle it, we propose to utilize the block coordinate descent (BCD) algorithm to alternately update the TPC matrix, AN covariance matrix, and phase shifts while keeping SR non-decreasing. Specifically, the optimal TPC matrix and AN covariance matrix are derived by Lagrangian multiplier method, and the optimal phase shifts are obtained by Majorization-Minimization (MM) algorithm. Since all variables can be calculated in closed form, the proposed algorithm is very efficient. We also extend the SRM problem to the more general multiple-IRs scenario and propose a BCD algorithm to solve it. Finally, simulation results validate the effectiveness of system security enhancement via an IRS.


title: Biological Random Walks: integrating heterogeneous data in disease gene prioritization

abstract: This work proposes a unified framework to leverage biological information in network propagation-based gene prioritization algorithms. Preliminary results on breast cancer data show significant improvements over state-of-the-art baselines, such as the prioritization of genes that are not identified as potential candidates by interactome-based algorithms, but that appear to be involved in/or potentially related to breast cancer, according to a functional analysis based on recent literature.


title: Acoustic Scene Classification Using Bilinear Pooling on Time-liked and Frequency-liked Convolution Neural Network

abstract: The current methodology in tackling Acoustic Scene Classification (ASC) task can be described in two steps, preprocessing of the audio waveform into log-mel spectrogram and then using it as the input representation for Convolutional Neural Network (CNN). This paradigm shift occurs after DCASE 2016 where this framework model achieves the state-of-the-art result in ASC tasks on the (ESC-50) dataset and achieved an accuracy of 64.5%, which constitute to 20.5% improvement over the baseline model, and DCASE 2016 dataset with an accuracy of 90.0% (development) and 86.2% (evaluation), which constitute a 6.4% and 9% improvements with respect to the baseline system. In this paper, we explored the use of harmonic and percussive source separation (HPSS) to split the audio into harmonic audio and percussive audio, which has received popularity in the field of music information retrieval (MIR). Although works have been done in using HPSS as input representation for CNN model in ASC task, this paper further investigate the possibility on leveraging the separated harmonic component and percussive component by curating 2 CNNs which tries to understand harmonic audio and percussive audio in their natural form, one specialized in extracting deep features in time biased domain and another specialized in extracting deep features in frequency biased domain, respectively. The deep features extracted from these 2 CNNs will then be combined using bilinear pooling. Hence, presenting a two-stream time and frequency CNN architecture approach in classifying acoustic scene. The model is being evaluated on DCASE 2019 sub task 1a dataset and scored an average of 65% on development dataset, Kaggle Leadership Private and Public board.


title: Learning Zero-Sum Simultaneous-Move Markov Games Using Function Approximation and Correlated Equilibrium

abstract: We develop provably efficient reinforcement learning algorithms for two-player zero-sum finite-horizon Markov games with simultaneous moves. To incorporate function approximation, we consider a family of Markov games where the reward function and transition kernel possess a linear structure. Both the offline and online settings of the problems are considered. In the offline setting, we control both players and aim to find the Nash Equilibrium by minimizing the duality gap. In the online setting, we control a single player playing against an arbitrary opponent and aim to minimize the regret. For both settings, we propose an optimistic variant of the least-squares minimax value iteration algorithm. We show that our algorithm is computationally efficient and provably achieves an $\tilde O(\sqrt{d^3 H^3 T} )$ upper bound on the duality gap and regret, where $d$ is the linear dimension, $H$ the horizon and $T$ the total number of timesteps. Our results do not require additional assumptions on the sampling model. Our setting requires overcoming several new challenges that are absent in Markov decision processes or turn-based Markov games. In particular, to achieve optimism with simultaneous moves, we construct both upper and lower confidence bounds of the value function, and then compute the optimistic policy by solving a general-sum matrix game with these bounds as the payoff matrices. As finding the Nash Equilibrium of a general-sum game is computationally hard, our algorithm instead solves for a Coarse Correlated Equilibrium (CCE), which can be obtained efficiently. To our best knowledge, such a CCE-based scheme for optimism has not appeared in the literature and might be of interest in its own right.


title: An efficient GFET structure

abstract: A graphene field effect transistor, where the active area is made of monolayer large-area graphene, is simulated including a full 2D Poisson equation and a drift-diffusion model with mobilities deduced by a direct numerical solution of the semiclassical Boltzmann equations for charge transport by a suitable discontinuous Galerkin approach. The critical issue in a graphene field effect transistor is the difficulty of fixing the off state which requires an accurate calibration of the gate voltages. In the present paper we propose and simulate a graphene field effect transistor structure which has well-behaved characteristic curves similar to those of conventional (with gap) semiconductor materials. The introduced device has a clear off region and can be the prototype of devices suited for post-silicon nanoscale electron technology. The specific geometry overcomes the problems of triggering the minority charge current and gives a viable way for the design of electron devices based on large area monolayer graphene as substitute of standard semiconductors in the active area. The good field effect transistor behavior of the current versus the gate voltage makes the simulated device very promising and a challenging case for experimentalists.


title: Profit from Two Bitcoin Mining Tactics: Towing and Shutdown

abstract: Since Bitcoin's inception in 2008, it has became attractive investments for both trading and mining. To mine Bitcoins, a miner has to invest in computing power and pay for electricity to solve cryptographic puzzles for rewards, if it becomes the first to solve a puzzle, paid in Bitcoin. Given that mining is such a resource intensive effort, miners seek new strategies trying to make the mining process more profitable.


title: The Big Three: A Methodology to Increase Data Science ROI by Answering the Questions Companies Care About

abstract: Companies may be achieving only a third of the value they could be getting from data science in industry applications. In this paper, we propose a methodology for categorizing and answering 'The Big Three' questions (what is going on, what is causing it, and what actions can I take that will optimize what I care about) using data science. The applications of data science seem to be nearly endless in today's modern landscape, with each company jockeying for position in the new data and insights economy. Yet, data scientists seem to be solely focused on using classification, regression, and clustering methods to answer the question 'what is going on'. Answering questions about why things are happening or how to take optimal actions to improve metrics are relegated to niche fields of research and generally neglected in industry data science analysis. We survey technical methods to answer these other important questions, describe areas in which some of these methods are being applied, and provide a practical example of how to apply our methodology and selected methods to a real business use case.


title: 28nm Fully-Depleted SOI Technology: Cryogenic Control Electronics for Quantum Computing

abstract: This paper reports the first cryogenic characterization of 28nm Fully-Depleted-SOI CMOS technology. A comprehensive study of digital/analog performances and body-biasing from room to the liquid helium temperature is presented. Despite a cryogenic operation, effectiveness of body-biasing remains unchanged and provides an excellent $V_{TH}$ controllability. Low-temperature operation enables higher drive current and a largely reduced subthreshold swing (down to 7mV/dec). FDSOI can provide a valuable approach to cryogenic low-power electronics. Applications such as classical control hardware for quantum processors are envisioned.


title: Defect-Free Axially-Stacked GaAs/GaAsP Nanowire Quantum Dots with Strong Carrier Confinement

abstract: Axially-stacked quantum dots (QDs) in nanowires (NWs) have important applications in fabricating nanoscale quantum devices and lasers. Although their performances are very sensitive to crystal quality and structures, there is relatively little study on defect-free growth with Au-free mode and structure optimisation for achiving high performances. Here, we report a detailed study of the first self-catalyzed defect-free axially-stacked deep NWQDs. High structural quality is maintained when 50 GaAs QDs are placed in a single GaAsP NW. The QDs have very sharp interfaces (1.8~3.6 nm) and can be closely stacked with very similar structural properties. They exhibit the deepest carrier confinement (~90 meV) and largest exciton-biexciton splitting (~11 meV) among non-nitride III-V NWQDs, and can maintain good optical properties after being stored in ambient atmosphere for over 6 months due to excellent stability. Our study sets a solid foundation to build high-performance axially-stacked NWQD devices that are compatible with CMOS technologies.


title: A fungal mycelium templates the growth of aragonite needles

abstract: Fungi live within diverse environments and survive well under extreme conditions that are usually beyond the tolerance of most other organisms. In different environments fungi are known to induce precipitation of a wide range of minerals. Various species of fungi have been shown to facilitate calcium carbonate mineralization. Here, inspired by examples of needle-fiber calcite formed via fungus-induced biomineralization typically observed in soils and sediments, we utilized inactivated fungus to synthetically induce precipitation of CaCO3 needles. To our knowledge, the feasibility of growing aragonitic needles within fungal mycelium in vitro has not been previously demonstrated. The needles we obtained were curved, displayed hexagonal facets, and demonstrated high-aspect ratios close to 60. The size and shape of these synthetic needles matched those of the mycelium of the natural fungus. Utilizing high-resolution characterization techniques, we studied the morphology and the micro- and nanostructures of the aragonitic needles. Our findings showed that even inactivated fungal mycelium, if present in the crystallization environment, can serve as a template for the formation of high-aspect ratio fibers and can stabilize metastable polymorphs.


title: Efficient Electrical Spin-Splitter Based on Non-Relativistic Collinear Antiferromagnetism

abstract: Electrical spin-current generation is among the core phenomena driving the field of spintronics. Using {\em ab initio} calculations we show that a room-temperature metallic collinear antiferromagnet RuO$_2$ allows for highly efficient spin-current generation, arising from anisotropically-split bands with conserved up and down spins along the N\'eel vector axis. The zero net moment antiferromagnet acts as an electrical spin-splitter with a 34$^\circ$ propagation angle between spin-up and spin-down currents. Correspondingly, the spin-conductivity is a factor of three larger than the record value from a survey of 20,000 non-magnetic spin-Hall materials. We propose a versatile spin-splitter-torque concept utilizing antiferromagnetic RuO$_2$ films interfaced with a ferromagnet.


title: Multiplicity at any torus-fixed point in a Richardson variety in the Symplectic Grassmannian

abstract: We give a combinatorial description of the multiplicity at any torus fixed point on a Richardson variety in the Symplectic Grassmannian.


title: Role of interactions in a closed quenched system

abstract: We study the non-equilibrium steady states in a closed system consisting of interacting particles obeying exclusion principle with quenched hopping rate. Cluster mean field approach is utilized to theoretically analyze the system dynamics in terms of phase diagram, density profiles, current, etc. with respect to interaction energy $E$. It turns out that on increasing the interaction energy beyond a critical value, $E_c$, shock region shows non-monotonic behavior and contracts until another critical value $E_{c_1}$ is attained; a further increase leads to its expansion. Moreover, the phase diagram of an interacting system with specific set of parameters has a good agreement with its non-interacting analogue. For interaction energy below $E_c$, a new shock phase displaying features different from non-interacting version is observed leading to two distinct shock phases. We have also performed Monte Carlo simulations extensively to validate our theoretical findings.


title: Block-Approximated Exponential Random Graphs

abstract: An important challenge in the field of exponential random graphs (ERGs) is the fitting of non-trivial ERGs on large graphs. By utilizing fast matrix block-approximation techniques, we propose an approximative framework to such non-trivial ERGs that result in dyadic independence (i.e., edge independent) distributions, while being able to meaningfully model both local information of the graph (e.g., degrees) as well as global information (e.g., clustering coefficient, assortativity, etc.) if desired. This allows one to efficiently generate random networks with similar properties as an observed network, and the models can be used for several downstream tasks such as link prediction. Our methods are scalable to sparse graphs consisting of millions of nodes. Empirical evaluation demonstrates competitiveness in terms of both speed and accuracy with state-of-the-art methods -- which are typically based on embedding the graph into some low-dimensional space -- for link prediction, showcasing the potential of a more direct and interpretable probabalistic model for this task.


title: High energy scattering amplitude in the linearized gravitational theory

abstract: The asymptotic behavior of the elastic scattering amplitude by the exchange of graviton between two scalar particles at high energies and fixed momentum transfers is reconsidered in the Logunov-Tavkhelidze equation in the linearized gravitational theory. The corrections to the eikonal approximation in the quasi-potential approach of relative order $1/p$ is developed with the principal contributions at high energy. The eikonal expression of scattering amplitude and the formal first correction are derived. The Yukawa potential is applied to discuss the results.


title: A validated energy model of a solar dish-Stirling system considering the cleanliness of mirrors

abstract: Solar systems based on the coupling of parabolic concentrating collectors and thermal engines (i.e. dish-Stirling systems) are among the most efficient generators of solar power currently available. This study focuses on the modelling of functioning data from a 32 kWe dish-Stirling solar plant installed at a facility test site on the University of Palermo campus, in Southern Italy. The proposed model, based on real monitored data, the energy balance of the collector and the partial load efficiency of the Stirling engine, can be used easily to simulate the annual energy production of such systems, making use of the solar radiation database, with the aim of encouraging a greater commercialisation of this technology. Introducing further simplifying assumptions based on our experimental data, the model can be linearised providing a new analytical expression of the parameters that characterise the widely used Stine empirical model. The model was calibrated against data corresponding to the collector with clean mirrors and used to predict the net electric production of the dish-Stirling accurately. A numerical method for assessing the daily level of mirror soiling without the use of direct reflectivity measures was also defined. The proposed methodology was used to evaluate the history of mirror soiling for the observation period, which shows a strong correlation with the recorded sequence of rains and dust depositions. The results of this study emphasise how desert dust transport events, frequent occurrences in parts of the Mediterranean, can have a dramatic impact on the electric power generation of dish-Stirling plants.


title: The Cantor-Schr\"oder-Bernstein Theorem for $\infty$-groupoids

abstract: We show that the Cantor-Schr\"oder-Bernstein Theorem for homotopy types, or $\infty$-groupoids holds in the following form: For any two types, if each one is embedded into the other, then they are equivalent. The argument is developed in the language of homotopy type theory, or Voevodsky's univalent foundations (HoTT/UF), and requires classical logic. It follows that the theorem holds in any boolean $\infty$-topos.


title: The Probabilistic Model Checker Storm

abstract: We present the probabilistic model checker Storm. Storm supports the analysis of discrete- and continuous-time variants of both Markov chains and Markov decision processes. Storm has three major distinguishing features. It supports multiple input languages for Markov models, including the JANI and PRISM modeling languages, dynamic fault trees, generalized stochastic Petri nets, and the probabilistic guarded command language. It has a modular set-up in which solvers and symbolic engines can easily be exchanged. Its Python API allows for rapid prototyping by encapsulating Storm's fast and scalable algorithms. This paper reports on the main features of Storm and explains how to effectively use them. A description is provided of the main distinguishing functionalities of Storm. Finally, an empirical evaluation of different configurations of Storm on the QComp 2019 benchmark set is presented.


title: Higher Nash blowups of normal toric varieties in prime characteristic

abstract: We prove that the higher Nash blowup of a normal toric variety defined over a field of positive characteristic is an isomorphism if and only if it is non-singular. We also extend a result by R. Toh-Yama which shows that higher Nash blowups do not give a one-step resolution of the $A_3$-singularity. These results were previously known only in characteristic zero.


title: PCSGAN: Perceptual Cyclic-Synthesized Generative Adversarial Networks for Thermal and NIR to Visible Image Transformation

abstract: In many real world scenarios, it is difficult to capture the images in the visible light spectrum (VIS) due to bad lighting conditions. However, the images can be captured in such scenarios using Near-Infrared (NIR) and Thermal (THM) cameras. The NIR and THM images contain the limited details. Thus, there is a need to transform the images from THM/NIR to VIS for better understanding. However, it is non-trivial task due to the large domain discrepancies and lack of abundant datasets. Nowadays, Generative Adversarial Network (GAN) is able to transform the images from one domain to another domain. Most of the available GAN based methods use the combination of the adversarial and the pixel-wise losses (like $L_1$ or $L_2$) as the objective function for training. The quality of transformed images in case of THM/NIR to VIS transformation is still not up to the mark using such objective function. Thus, better objective functions are needed to improve the quality, fine details and realism of the transformed images. A new model for THM/NIR to VIS image transformation called Perceptual Cyclic-Synthesized Generative Adversarial Network (PCSGAN) is introduced to address these issues. The PCSGAN uses the combination of the perceptual (i.e., feature based) losses along with the pixel-wise and the adversarial losses. Both the quantitative and qualitative measures are used to judge the performance of the PCSGAN model over the WHU-IIP face and the RGB-NIR scene datasets. The proposed PCSGAN outperforms the state-of-the-art image transformation models, including Pix2pix, DualGAN, CycleGAN, PS2GAN, and PAN in terms of the SSIM, MSE, PSNR and LPIPS evaluation measures. The code is available at https://github.com/KishanKancharagunta/PCSGAN.


title: Covariant Quantum Mechanics and Quantum Spacetime

abstract: We present in the article the formulation of a version of Lorentz covariant quantum mechanics based on a group theoretical construction from a Heisenberg-Weyl symmetry with position and momentum operators transforming as Minkowski four-vectors under the Lorentz symmetry. The basic representation is identified as a coherent state representation, essentially an irreducible component of the regular representation, with the matching representation of an extension of the group $C^*$-algebra giving the algebra of observables. The key feature of the formulation is that it is not unitary but pseudo-unitary, exactly in the same sense as the Minkowski spacetime representation. Explicit wavefunction description is given without any restriction of the variable domains, yet with a finite integral inner product. The associated covariant harmonic oscillator Fock state basis has all the standard properties in exact analog to those of a harmonic oscillator with Euclidean position and momentum operators of any `dimension'. Galilean limit of the Lorentz symmetry and the classical limit of the Lorentz covariant framework are retrieved rigorously through appropriate symmetry contractions of the algebra and its representation, including the dynamics described through the symmetry of the phase space, given both in terms of real/complex number coordinates and noncommutative operator coordinates. The latter gives an explicit picture of the (projective) Hilbert space as a quantum/noncommutative spacetime.


title: Pompeiu's theorem and the moduli space of triangles

abstract: We introduce a kind of converse of Pompeiu's theorem. Fix an equilateral triangle $\triangle A_0B_0C_0$, then for any triangle $\triangle ABC$ there is a unique point $P$ inside the circumcircle $\Gamma_0$ of $\triangle A_0B_0C_0$ such that a triangle with edge lengths $PA_0, PB_0$, and $PC_0$ is similar to $\triangle ABC$. It follows that an open disc inside $\Gamma_0$ can be considered as a moduli space of similarity classes of triangles. We show that it is essentially equivalent to another moduli space based on a shape function of triangles which has been used in preceding studies.


title: Small-gain theorem for stability, cooperative control and distributed observation of infinite networks

abstract: Motivated by a paradigm shift towards a hyper-connected world, we develop a computationally tractable small-gain theorem for a network of infinitely many systems, termed as infinite networks. The proposed small-gain theorem addresses exponential input-to-state stability with respect to closed sets, which enables us to analyze diverse stability problems in a unified manner. The small-gain condition, expressed in terms of the spectral radius of a gain operator collecting all the information about the internal Lyapunov gains, can be numerically computed for a large class of systems in an efficient way. To demonstrate broad applicability of our small-gain theorem, we apply it to the stability analysis of infinite time-varying networks, to consensus in infinite-agent systems, as well as to the design of distributed observers for infinite networks.


title: Electrostatic Near-Limits Kinetic Energy Harvesting from Arbitrary Input Vibrations

abstract: The full architecture of an electrostatic kinetic energy harvester (KEH) based on the concept of near-limits KEH is reported. This concept refers to the conversion of kinetic energy to electric energy, from environmental vibrations of arbitrary forms, and at rates that target the physical limits set by the device's size and the input excitation characteristics. This is achieved thanks to the synthesis of particular KEH's mass dynamics, that maximize the harvested energy. Synthesizing these dynamics requires little hypotheses on the exact form of the input vibrations. In the proposed architecture, these dynamics are implemented by an adequate mechanical control which is synthesized by the electrostatic transducer. An interface circuit is proposed to carry out the necessary energy transfers between the transducer and the system's energy tank. A computation and finite-state automaton unit controls the interface circuit, based on the external input and on the system's mechanical state. The operation of the reported near-limits KEH is illustrated in simulations that demonstrate proof of concept of the proposed architecture. A figure of $68\%$ of the absolute limit of the KEH's input energy for the considered excitation is attained. This can be further improved by complete system optimization that takes into account the application constraints, the control law, the mechanical design of the transducer, the electrical interface design, and the sensing and computation blocks.


title: Graph Deconvolutional Generation

abstract: Graph generation is an extremely important task, as graphs are found throughout different areas of science and engineering. In this work, we focus on the modern equivalent of the Erdos-Renyi random graph model: the graph variational autoencoder (GVAE). This model assumes edges and nodes are independent in order to generate entire graphs at a time using a multi-layer perceptron decoder. As a result of these assumptions, GVAE has difficulty matching the training distribution and relies on an expensive graph matching procedure. We improve this class of models by building a message passing neural network into GVAE's encoder and decoder. We demonstrate our model on the specific task of generating small organic molecules


title: GRAPHITE: Generating Automatic Physical Examples for Machine-Learning Attacks on Computer Vision Systems

abstract: This paper investigates an adversary's ease of attack in generating adversarial examples for real-world scenarios. We address three key requirements for practical attacks for the real-world: 1) automatically constraining the size and shape of the attack so it can be applied with stickers, 2) transform-robustness, i.e., robustness of a attack to environmental physical variations such as viewpoint and lighting changes, and 3) supporting attacks in not only white-box, but also black-box hard-label scenarios, so that the adversary can attack proprietary models. In this work, we propose GRAPHITE, an efficient and general framework for generating attacks that satisfy the above three key requirements. GRAPHITE takes advantage of transform-robustness, a metric based on expectation over transforms (EoT), to automatically generate small masks and optimize with gradient-free optimization. GRAPHITE is also flexible as it can easily trade-off transform-robustness, perturbation size, and query count in black-box settings. On a GTSRB model in a hard-label black-box setting, we are able to find attacks on all possible 1,806 victim-target class pairs with averages of 77.8% transform-robustness, perturbation size of 16.63% of the victim images, and 126K queries per pair. For digital-only attacks where achieving transform-robustness is not a requirement, GRAPHITE is able to find successful small-patch attacks with an average of only 566 queries for 92.2% of victim-target pairs. GRAPHITE is also able to find successful attacks using perturbations that modify small areas of the input image against PatchGuard, a recently proposed defense against patch-based attacks.


title: 4D Semantic Cardiac Magnetic Resonance Image Synthesis on XCAT Anatomical Model

abstract: We propose a hybrid controllable image generation method to synthesize anatomically meaningful 3D+t labeled Cardiac Magnetic Resonance (CMR) images. Our hybrid method takes the mechanistic 4D eXtended CArdiac Torso (XCAT) heart model as the anatomical ground truth and synthesizes CMR images via a data-driven Generative Adversarial Network (GAN). We employ the state-of-the-art SPatially Adaptive De-normalization (SPADE) technique for conditional image synthesis to preserve the semantic spatial information of ground truth anatomy. Using the parameterized motion model of the XCAT heart, we generate labels for 25 time frames of the heart for one cardiac cycle at 18 locations for the short axis view. Subsequently, realistic images are generated from these labels, with modality-specific features that are learned from real CMR image data. We demonstrate that style transfer from another cardiac image can be accomplished by using a style encoder network. Due to the flexibility of XCAT in creating new heart models, this approach can result in a realistic virtual population to address different challenges the medical image analysis research community is facing such as expensive data collection. Our proposed method has a great potential to synthesize 4D controllable CMR images with annotations and adaptable styles to be used in various supervised multi-site, multi-vendor applications in medical image analysis.


title: Vortex cores in narrow thin-film strips

abstract: We study vortex current distributions in narrow thin-film superconducting strips. If one defines the vortex core ``boundary" as a curve where the current reaches the depairing value, intriguing features emerge. Our conclusions based on the London approach have only qualitative relevance since the approach breaks down near the core. Still, the main observation which might be useful is that the core size near the strip edges is smaller than in the rest of the strip. If so, the Bardeen-Stephen flux-flow resistivity should be reduced near the edges. Moreover, at elevated temperatures, when the depairing current is small, the vortex core may extend to the whole strip width, thus turning into an edge-to-edge phase-slip line.


title: Intelligent Reflecting Surface Aided Multiple Access: Capacity Region and Deployment Strategy

abstract: Intelligent reflecting surface (IRS) is a new promising technology that is able to manipulate the wireless propagation channel via smart and controllable signal reflection. In this paper, we investigate the capacity region of a multiple access channel (MAC) with two users sending independent messages to an access point (AP), aided by $M$ IRS reflecting elements. We consider two practical IRS deployment strategies that lead to different user-AP effective channels, namely, the distributed deployment where the $M$ reflecting elements form two IRSs, each deployed in the vicinity of one user, versus the centralized deployment where all the $M$ reflecting elements are deployed in the vicinity of the AP. For the distributed deployment, we derive the capacity region in closed-form; while for the centralized deployment, we derive a capacity region outer bound and propose an efficient rate-profile based method to characterize an achievable rate region (or capacity region inner bound). Furthermore, we compare the capacity regions of the two cases and draw useful insights into the optimal deployment of IRS in practical systems.


title: Wiener index, number of subtrees, and tree eccentric sequence

abstract: The eccentricity of a vertex $u$ in a connected graph $G$ is the distance between $u$ and a vertex farthest from it; the eccentric sequence of $G$ is the nondecreasing sequence of the eccentricities of $G$. In this paper, we determine the unique tree that minimises the Wiener index, i.e. the sum of distances between all unordered vertex pairs, among all trees with a given eccentric sequence. We show that the same tree maximises the number of subtrees among all trees with a given eccentric sequence, thus providing another example of negative correlation between the number of subtrees and the Wiener index of trees. Furthermore, we provide formulas for the corresponding extreme values of these two invariants in terms of the eccentric sequence. As a corollary to our results, we determine the unique tree that minimises the edge Wiener index, the vertex-edge Wiener index, the Schulz index (or degree distance), and the Gutman index among all trees with a given eccentric sequence.


title: Arrows of Time and Initial and Final Conditions in the Quantum Mechanics of Closed Systems Like the Universe

abstract: A model quantum cosmology is used to illustrate how arrows of time emerge in a universe governed by a time-neutral dynamical theory constrained by time asymmetric initial and final boundary conditions represented by initial and final density matrices. In a quantum universe universe arrows of time are described by the probabilities of appropriately coarse grained sets of histories of quantities like entropy that grow or decay. We show that the requirement of that these sets of histories decohere implies two things: (1) A time asymmetry between initial and final conditions that is a basis for arrows ot time. (2) How a final state of indifference that is represented by a final density matrix proportional to the unit density matrix is consistent with causality, and allows a finer-grained description of the model universe in terms of decoherent histories than any other final state.


title: A Divide and Conquer Algorithm of Bayesian Density Estimation

abstract: Data sets for statistical analysis become extremely large even with some difficulty of being stored on one single machine. Even when the data can be stored in one machine, the computational cost would still be intimidating. We propose a divide and conquer solution to density estimation using Bayesian mixture modeling including the infinite mixture case. The methodology can be generalized to other application problems where a Bayesian mixture model is adopted. The proposed prior on each machine or subsample modifies the original prior on both mixing probabilities as well as on the rest of parameters in the distributions being mixed. The ultimate estimator is obtained by taking the average of the posterior samples corresponding to the proposed prior on each subset. Despite the tremendous reduction in time thanks to data splitting, the posterior contraction rate of the proposed estimator stays the same (up to a log factor) as that of the original prior when the data is analyzed as a whole. Simulation studies also justify the competency of the proposed method compared to the established WASP estimator in the finite dimension case. In addition, one of our simulations is performed in a shape constrained deconvolution context and reveals promising results. The application to a GWAS data set reveals the advantage over a naive method that uses the original prior.


title: Product Subset Problem : Applications to number theory and cryptography

abstract: We consider applications of Subset Product Problem (SPP) in number theory and cryptography. We obtain a probabilistic algorithm that attack SPP and we analyze it with respect time/space complexity and success probability. In fact we provide an application to the problem of finding Carmichael numbers and an attack to Naccache-Stern knapsack cryptosystem, where we update previous results.


title: Visual Data Analysis and Simulation Prediction for COVID-19

abstract: The COVID-19 (formerly, 2019-nCoV) epidemic has become a global health emergency, as such, WHO declared PHEIC. China has taken the most hit since the outbreak of the virus, which could be dated as far back as late November by some experts. It was not until January 23rd that the Wuhan government finally recognized the severity of the epidemic and took a drastic measure to curtain the virus spread by closing down all transportation connecting the outside world. In this study, we seek to answer a few questions: How did the virus get spread from the epicenter Wuhan city to the rest of the country? To what extent did the measures, such as, city closure and community quarantine, help controlling the situation? More importantly, can we forecast any significant future development of the event had some of the conditions changed? By collecting and visualizing publicly available data, we first show patterns and characteristics of the epidemic development; we then employ a mathematical model of disease transmission dynamics to evaluate the effectiveness of some epidemic control measures, and more importantly, to offer a few tips on preventive measures.


title: Strong solutions of stochastic differential equations with coefficients in mixed-norm spaces

abstract: By studying parabolic equations in mixed-norm spaces, we prove the existence and uniqueness of strong solutions to stochastic differential equations driven by Brownian motion with coefficients in spaces with mixed-norm, which extends Krylov and R\"ockner's result in [11] and Zhang's result in [18].


title: Effective Capacity of NOMA with Finite Blocklength for Low-Latency Communications

abstract: In this paper, we investigate the link-layer rate of a two-users non-orthogonal multiple access (NOMA) network in finite blocklength (short packet communications) regime, where the two users are paired from a set of V users. The overall reliability consists of the transmission error probability and the queuing delay violation probability. The performance of the two-users NOMA network in finite blocklength is verified for achieving latency and reliability, using the effective capacity (EC) framework. Specifically, we derive closed-form expressions for the EC of the two-users NOMA network in finite blocklength regime, considering transmissions over Rayleigh fading channels. We also study a multiuser NOMA network and derive the total EC of two-users NOMA subsets, and show that the NOMA set with users having distinct channel conditions achieve maximum total EC. Focusing on a two-users NOMA network, we study the impact of the transmit signal-to-noise ratio, delay exponent, and transmission error probability, on the achievable EC of each user. The analysis shows that when the delay exponent is high, the delay violation probability does not improve below a certain value due to the dominant factor of the transmission error probability. The accuracy of the proposed closed-form expressions for the individual EC of a two-users NOMA is verified using the Monte-Carlo simulations.


title: The poset of morphism-extension classes of countable graphs

abstract: Let $\mathrm{XY_{L,T}}$ denote the class of countably infinite $L$-structures that satisfy the axioms $T$ and in which all homomorphisms of type X (these could be homomorphisms, monomorphisms, or isomorphisms) between finite substructures of $M$ are restrictions of an endomorphism of $M$ of type Y (for example, an automorphism or a surjective endomorphism). Lockett and Truss introduced 18 such morphism-extension classes for relational structures. For a given pair $L,T$, however, two or more morphism-extension properties may define the same class of structures. In this paper, we establish all equalities and inequalities between morphism-extension classes of countable (undirected, loopless) graphs.


title: Crisis contagion in the world trade network

abstract: We present a model of worldwide crisis contagion based on the Google matrix analysis of the world trade network obtained from the UN Comtrade database. The fraction of bankrupted countries exhibits an \textit{on-off} phase transition governed by a bankruptcy threshold $\kappa$ related to the trade balance of the countries. For $\kappa>\kappa_c$, the contagion is circumscribed to less than 10\% of the countries, whereas, for $\kappa<\kappa_c$, the crisis is global with about 90\% of the countries going to bankruptcy. We measure the total cost of the crisis during the contagion process. In addition to providing contagion scenarios, our model allows to probe the structural trading dependencies between countries. For different networks extracted from the world trade exchanges of the last two decades, the global crisis comes from the Western world. In particular, the source of the global crisis is systematically the Old Continent and The Americas (mainly US and Mexico). Besides the economy of Australia, those of Asian countries, such as China, India, Indonesia, Malaysia and Thailand, are the last to fall during the contagion. Also, the four BRIC are among the most robust countries to the world trade crisis.


title: Augmented Normalizing Flows: Bridging the Gap Between Generative Flows and Latent Variable Models

abstract: In this work, we propose a new family of generative flows on an augmented data space, with an aim to improve expressivity without drastically increasing the computational cost of sampling and evaluation of a lower bound on the likelihood. Theoretically, we prove the proposed flow can approximate a Hamiltonian ODE as a universal transport map. Empirically, we demonstrate state-of-the-art performance on standard benchmarks of flow-based generative modeling.


title: Stable manifolds of biholomorphisms in $\mathbb{C}^n$ asymptotic to formal curves

abstract: Given a germ of biholomorphism $F\in\mathrm{Diff}(\mathbb{C}^n,0)$ with a formal invariant curve $\Gamma$ such that the multiplier of the restricted formal diffeomorphism $F|_\Gamma$ is a root of unity or satisfies $|(F|_\Gamma)'(0)|<1$, we prove that either $\Gamma$ is contained in the set of periodic points of $F$ or there exists a finite family of stable manifolds of $F$ where all the orbits are asymptotic to $\Gamma$ and whose union eventually contains every orbit asymptotic to $\Gamma$. This result generalizes to the case where $\Gamma$ is a formal periodic curve.


title: UST branches, martingales, and multiple SLE(2)

abstract: We identify the local scaling limit of multiple boundary-to-boundary branches in a uniform spanning tree (UST) as a local multiple SLE(2), i.e., an SLE(2) process weighted by a suitable partition function. By recent results, this also characterizes the "global" scaling limit of the full collection of full curves. The identification is based on a martingale observable in the UST with $N$ branches, obtained by weighting the well-known martingale in the UST with one branch by the discrete partition functions of the models. The obtained weighting transforms of the discrete martingales and the limiting SLE processes, respectively, only rely on a discrete domain Markov property and (essentially) the convergence of partition functions. We illustrate their generalizability by sketching an analogous convergence proof for a boundary-visiting UST branch and a boundary-visiting SLE(2).


title: Structural, electronic, and magnetic properties of Vanadium-based Janus dichalcogenide monolayers: A first-principles study

abstract: The structural, electronic, and magnetic properties of VSSe, VSeTe, VSTe monolayers in both 2H and 1T phases are investigated via first-principles calculations. The 2H phase is energetically favorable in VSSe and VSeTe, whereas the 1T phase is lower in energy in VSTe. For V-based Janus monolayers in the 2H phase, calculations of the magnetic anisotropy show an easy-plane for the magnetic moment. As such, they should not exhibit a ferromagnetic phase transition, but instead, a Berezinskii-Kosterlitz-Thouless (BKT) transition. A classical XY model with nearest-neighbor coupling estimates critical temperatures (T$_{BKT}$) ranging from 106 K for VSSe to 46 K for VSTe.


title: Quantum Effects in Palatini Higgs Inflation

abstract: We study quantum effects in Higgs inflation in the Palatini formulation of gravity, in which the metric and connection are treated as independent variables. We exploit the fact that the cutoff, above which perturbation theory breaks down, is higher than the scale of inflation. Unless new physics above the cutoff leads to unnaturally large corrections, we can directly connect low-energy physics and inflation. On the one hand, the lower bound on the top Yukawa coupling due to collider experiments leads to an upper bound on the non-minimal coupling of the Higgs field to gravity: $\xi \lesssim 10^8$. On the other hand, the Higgs potential can only support successful inflation if $\xi \gtrsim 10^6$. This leads to a fairly strict upper bound on the top Yukawa coupling of $0.925$ (defined in the $\overline{\text{MS}}$-scheme at the energy scale $173.2\,\text{GeV}$) and constrains the inflationary prediction for the tensor-to-scalar ratio. Additionally, we compare our findings to metric Higgs inflation.


title: Controlling Computation versus Quality for Neural Sequence Models

abstract: Most neural networks utilize the same amount of compute for every example independent of the inherent complexity of the input. Further, methods that adapt the amount of computation to the example focus on finding a fixed inference-time computational graph per example, ignoring any external computational budgets or varying inference time limitations. In this work, we utilize conditional computation to make neural sequence models (Transformer) more efficient and computation-aware during inference. We first modify the Transformer architecture, making each set of operations conditionally executable depending on the output of a learned control network. We then train this model in a multi-task setting, where each task corresponds to a particular computation budget. This allows us to train a single model that can be controlled to operate on different points of the computation-quality trade-off curve, depending on the available computation budget at inference time. We evaluate our approach on two tasks: (i) WMT English-French Translation and (ii) Unsupervised representation learning (BERT). Our experiments demonstrate that the proposed Conditional Computation Transformer (CCT) is competitive with vanilla Transformers when allowed to utilize its full computational budget, while improving significantly over computationally equivalent baselines when operating on smaller computational budgets.


title: Demazure crystals for Kohnert polynomials

abstract: Kohnert polynomials are polynomials indexed by unit cell diagrams in the first quadrant defined earlier by the author and Searles that give a common generalization of Schubert polynomials and Demazure characters for the general linear group. Demazure crystals are certain truncations of normal crystals whose characters are Demazure characters. For each diagram satisfying a southwest condition, we construct a Demazure crystal whose character is the Kohnert polynomial for the given diagram, resolving an earlier conjecture of the author and Searles that these polynomials expand nonnegatively into Demazure characters. We give explicit formulas for the expansions with applications including a characterization of those diagrams for which the corresponding Kohnert polynomial is a single Demazure character.


title: The ultrafast onset of exciton formation in 2D semiconductors

abstract: The equilibrium and non-equilibrium optical properties of single-layer transition metal dichalcogenides (TMDs) are determined by strongly bound excitons. Exciton relaxation dynamics in TMDs have been extensively studied by time-domain optical spectroscopies. However, the formation dynamics of excitons following non-resonant photoexcitation of free electron-hole pairs have been challenging to directly probe because of their inherently fast timescales. Here we use extremely short optical pulses to non-resonantly excite an electron-hole plasma and show the formation of two-dimensional excitons in single-layer MoS2 on the timescale of 30 fs via the induced changes to photo-absorption. These formation dynamics are significantly faster than in conventional 2D quantum wells and are attributed to the intense Coulombic interactions present in 2D TMDs. A theoretical model of a coherent polarization that dephases and relaxes to an incoherent exciton population reproduces the experimental dynamics on the sub-100-fs timescale and sheds light into the underlying mechanism of how the lowest-energy excitons, which are the most important for optoelectronic applications, form from higher-energy excitations. Importantly, a phonon-mediated exciton cascade from higher energy states to the ground excitonic state is found to be the rate-limiting process. These results set an ultimate timescale of the exciton formation in TMDs and elucidate the exceptionally fast physical mechanism behind this process.


title: Velocity jump processes : an alternative to multi-timestep methods for faster and accurate molecular dynamics simulations

abstract: We propose a new route to accelerate molecular dynamics through the use of velocity jump processes allowing for an adaptive time-step specific to each atom-atom pair (2-body) interactions. We start by introducing the formalism of the new velocity jump molecular dynamics, ergodic with respect to the canonical measure. We then introduce the new BOUNCE integrator that allows for long-range forces to be evaluated at random and optimal time-steps, leading to strong savings in direct space. The accuracy and computational performances of a first BOUNCE implementation dedicated to classical (non-polarizable) force fields is tested in the cases of pure direct-space droplet-like simulations and of periodic boundary conditions (PBC) simulations using Smooth Particule Mesh Ewald. An analysis of the capability of BOUNCE to reproduce several condensed phase properties is provided. Since electrostatics and van der Waals 2-body contributions are evaluated much less often than with standard integrators using a 1fs timestep, up to a 400 % direct-space acceleration is observed. Applying the reversible reference system propagator algorithms (RESPA(1)) to reciprocal space (many-body) interactions allows BOUNCE-RESPA(1) to maintain large speedups in PBC while maintaining precision. Overall, we show that replacing the BAOAB integrator by the BOUNCE adaptive framework preserves a similar accuracy and leads to significant computational savings.


title: SMART U(1)$_X$ $-$ Standard Model with Axion, Right handed neutrinos, Two Higgs doublets and U(1)$_X$ gauge symmetry

abstract: To address five fundamental shortcomings of the Standard Model (SM) of particle physics and cosmology, we propose a SMART U(1)$_X$ model which is a $U(1)_X \times U(1)_{PQ}$ extension of the SM. The $U(1)_X$ gauge symmetry is a generalization of the well-known $U(1)_{B-L}$ symmetry and $U(1)_{PQ}$ is the global Peccie-Quinn (PQ) symmetry. Three right handed neutrinos are added to cancel $U(1)_X$ related anomalies, and they play a crucial role in understanding the observed neutrino oscillations and explaining the observed baryon asymmetry in the universe via leptogenesis. The PQ symmetry helps resolve the strong CP problem and also provides axion as a compelling dark matter (DM) candidate. The $U(1)_X$ gauge symmetry enables us to implement the inflection-point inflation scenario with $H_{inf} \lesssim 2 \times 10^{7}$ GeV, where $H_{inf}$ is the value of Hubble parameter during inflation. This allows us to overcome a potential axion domain wall problem as well as the axion isocurvature problem. The SMART U(1)$_X$ model can be merged with $SU(5)$ as we briefly show.


title: Targeted Forgetting and False Memory Formation in Continual Learners through Adversarial Backdoor Attacks

abstract: Artificial neural networks are well-known to be susceptible to catastrophic forgetting when continually learning from sequences of tasks. Various continual (or "incremental") learning approaches have been proposed to avoid catastrophic forgetting, but they are typically adversary agnostic, i.e., they do not consider the possibility of a malicious attack. In this effort, we explore the vulnerability of Elastic Weight Consolidation (EWC), a popular continual learning algorithm for avoiding catastrophic forgetting. We show that an intelligent adversary can bypass the EWC's defenses, and instead cause gradual and deliberate forgetting by introducing small amounts of misinformation to the model during training. We demonstrate such an adversary's ability to assume control of the model via injection of "backdoor" attack samples on both permuted and split benchmark variants of the MNIST dataset. Importantly, once the model has learned the adversarial misinformation, the adversary can then control the amount of forgetting of any task. Equivalently, the malicious actor can create a "false memory" about any task by inserting carefully-designed backdoor samples to any fraction of the test instances of that task. Perhaps most damaging, we show this vulnerability to be very acute; neural network memory can be easily compromised with the addition of backdoor samples into as little as 1% of the training data of even a single task.


title: Artificial Intelligence Forecasting of Covid-19 in China

abstract: BACKGROUND An alternative to epidemiological models for transmission dynamics of Covid-19 in China, we propose the artificial intelligence (AI)-inspired methods for real-time forecasting of Covid-19 to estimate the size, lengths and ending time of Covid-19 across China. METHODS We developed a modified stacked auto-encoder for modeling the transmission dynamics of the epidemics. We applied this model to real-time forecasting the confirmed cases of Covid-19 across China. The data were collected from January 11 to February 27, 2020 by WHO. We used the latent variables in the auto-encoder and clustering algorithms to group the provinces/cities for investigating the transmission structure. RESULTS We forecasted curves of cumulative confirmed cases of Covid-19 across China from Jan 20, 2020 to April 20, 2020. Using the multiple-step forecasting, the estimated average errors of 6-step, 7-step, 8-step, 9-step and 10-step forecasting were 1.64%, 2.27%, 2.14%, 2.08%, 0.73%, respectively. We predicted that the time points of the provinces/cities entering the plateau of the forecasted transmission dynamic curves varied, ranging from Jan 21 to April 19, 2020. The 34 provinces/cities were grouped into 9 clusters. CONCLUSIONS The accuracy of the AI-based methods for forecasting the trajectory of Covid-19 was high. We predicted that the epidemics of Covid-19 will be over by the middle of April. If the data are reliable and there are no second transmissions, we can accurately forecast the transmission dynamics of the Covid-19 across the provinces/cities in China. The AI-inspired methods are a powerful tool for helping public health planning and policymaking.


title: Handling Missing Annotations in Supervised Learning Data

abstract: Data annotation is an essential stage in supervised learning. However, the annotation process is exhaustive and time consuming, specially for large datasets. Activities of Daily Living (ADL) recognition is an example of systems that exploit very large raw sensor data readings. In such systems, sensor readings are collected from activity-monitoring sensors in a 24/7 manner. The size of the generated dataset is so huge that it is almost impossible for a human annotator to give a certain label to every single instance in the dataset. This results in annotation gaps in the input data to the adopting supervised learning system. The performance of the recognition system is negatively affected by these gaps. In this work, we propose and investigate three different paradigms to handle these gaps. In the first paradigm, the gaps are taken out by dropping all unlabeled readings. A single "Unknown" or "Do-Nothing" label is given to the unlabeled readings within the operation of the second paradigm. The last paradigm handles these gaps by giving every one of them a unique label identifying the encapsulating deterministic labels. Also, we propose a semantic preprocessing method of annotation gaps by constructing a hybrid combination of some of these paradigms for further performance improvement. The performance of the proposed three paradigms and their hybrid combination is evaluated using an ADL benchmark dataset containing more than $2.5\times 10^6$ sensor readings that had been collected over more than nine months. The evaluation results emphasize the performance contrast under the operation of each paradigm and support a specific gap handling approach for better performance.


title: Towards a More Complete and Accurate Experimental Nuclear Reaction Data Library (EXFOR): International Collaboration Between Nuclear Reaction Data Centres (NRDC)

abstract: The International Network of Nuclear Reaction Data Centres (NRDC) coordinated by the IAEA Nuclear Data Section (NDS) is successfully collaborating in the maintenance and development of the EXFOR library. As the scope of published data expands (e.g., to higher energy, to heavier projectile) to meet the needs from the frontier of sciences and applications, it becomes nowadays a hard and challenging task to maintain both completeness and accuracy of the whole EXFOR library. The paper describes evolution of the library with highlights on recent developments.


title: An SMT-Based Concolic Testing Tool for Logic Programs

abstract: Concolic testing mixes symbolic and concrete execution to generate test cases covering paths effectively. Its benefits have been demonstrated for more than 15 years to test imperative programs. Other programming paradigms, like logic programming, have received less attention. In this paper, we present a concolic-based test generation method for logic programs. Our approach exploits SMT-solving for constraint resolution. We then describe the implementation of a concolic testing tool for Prolog and validate it on some selected benchmarks.


title: A New Pricing Theory That Solves the St. Petersburg Paradox

abstract: The St. Petersburg Paradox, an important topic in probability theory, has not been solved in the last 280 years. Since Nicolaus Bernoulli proposed the St. Petersburg Paradox in 1738, many people had tried to solve it and had proposed various explanations, but all were not satisfactory. In this paper we propose a new pricing theory with several rules, which incidentally resolves this paradox. The new pricing theory states that so-called fair (reasonable) pricing should be judged by the seller and the buyer independently. Reasonable pricing for the seller may not be appropriate for the buyer. The seller cares about costs, while the buyer is concerned about the realistic prospect of returns.The pricing theory we proposed can be applied to financial markets to solve the confusion that financial asset return with fat tails distribution will cause the option pricing formula to fail, thus making up the theoretical defects of quantitative financial pricing theory.


title: Pricing Bitcoin Derivatives under Jump-Diffusion Models

abstract: In recent years cryptocurrency trading has captured the attention of practitioners and academics. The volume of the exchange with standard currencies has known a dramatic increasing of late. This paper addresses to the need of models describing a bitcoin-US dollar exchange dynamic and their use to evaluate European option having bitcoin as underlying asset.


title: Seeing Around Corners with Edge-Resolved Transient Imaging

abstract: Non-line-of-sight (NLOS) imaging is a rapidly growing field seeking to form images of objects outside the field of view, with potential applications in search and rescue, reconnaissance, and even medical imaging. The critical challenge of NLOS imaging is that diffuse reflections scatter light in all directions, resulting in weak signals and a loss of directional information. To address this problem, we propose a method for seeing around corners that derives angular resolution from vertical edges and longitudinal resolution from the temporal response to a pulsed light source. We introduce an acquisition strategy, scene response model, and reconstruction algorithm that enable the formation of 2.5-dimensional representations -- a plan view plus heights -- and a 180$^{\circ}$ field of view (FOV) for large-scale scenes. Our experiments demonstrate accurate reconstructions of hidden rooms up to 3 meters in each dimension.


title: Leakage detection for a transmon-based surface code

abstract: Leakage outside of the qubit computational subspace, present in many leading experimental platforms, constitutes a threatening error for quantum error correction (QEC) for qubits. We develop a leakage-detection scheme via Hidden Markov models (HMMs) for transmon-based implementations of the surface code. By performing realistic density-matrix simulations of the distance-3 surface code (Surface-17), we observe that leakage is sharply projected and leads to an increase in the surface-code defect probability of neighboring stabilizers. Together with the analog readout of the ancilla qubits, this increase enables the accurate detection of the time and location of leakage. We restore the logical error rate below the memory break-even point by post-selecting out leakage, discarding about 47% of the data. Leakage detection via HMMs opens the prospect for near-term QEC demonstrations, targeted leakage reduction and leakage-aware decoding and is applicable to other experimental platforms.


title: Fibration theorems \`a la Milnor for differentiable maps with non-isolated singularities

abstract: We prove fibration theorems \`a la Milnor for differentiable real maps with non isolated critical values. We study the situation for maps with linear discriminant, and prove that the concept of d-regularity is the key point for the existence of a Milnor fibration on the sphere. We also explain how one can modify the target space by homeomorphisms to linearize a general discriminant. Whenever the composed map is d-regular one has fibration on the sphere. Plenty of examples are discussed along the text.


title: Non-Archimedean Electrostatics

abstract: We introduce ensembles of repelling charged particles restricted to a ball in a non-archimedean field (such as the $p$-adic numbers) with interaction energy between pairs of particles proportional to the logarithm of the ($p$-adic) distance between them. In the {\em canonical ensemble}, a system of $N$ particles is put in contact with a heat bath at fixed inverse temperature $\beta$ and energy is allowed to flow between the system and the heat bath. Using standard axioms of statistical physics, the relative density of states is given by the $\beta$ power of the ($p$-adic) absolute value of the Vandermonde determinant in the locations of the particles. The partition function is the normalizing constant (as a function of $\beta$) of this ensemble, and we identify a recursion that allows this to be computed explicitly in finite time. Probabilities of interest, including the probabilities that specified subsets will have a prescribed occupation number of particles, and the conditional distribution of particles within a subset given a prescribed occupation number, are given explicitly in terms of the partition function. We then turn to the {\em grand canonical ensemble} where both the energy and number of particles are variable. We compute similar probabilities to those in the canonical ensemble and show how these probabilities can be given in terms the canonical and grand canonical partition functions. Finally, we briefly consider the multi-component ensemble where particles are allowed to take different integer charges, and we connect basic properties of this ensemble to the canonical and grand canonical ensembles.


title: Bayesian Structure Learning in Multi-layered Genomic Networks

abstract: Integrative network modeling of data arising from multiple genomic platforms provides insight into the holistic picture of the interactive system, as well as the flow of information across many disease domains including cancer. The basic data structure consists of a sequence of hierarchically ordered datasets for each individual subject, which facilitates integration of diverse inputs, such as genomic, transcriptomic, and proteomic data. A primary analytical task in such contexts is to model the layered architecture of networks where the vertices can be naturally partitioned into ordered layers, dictated by multiple platforms, and exhibit both undirected and directed relationships. We propose a multi-layered Gaussian graphical model (mlGGM) to investigate conditional independence structures in such multi-level genomic networks in human cancers. We implement a Bayesian node-wise selection (BANS) approach based on variable selection techniques that coherently accounts for the multiple types of dependencies in mlGGM; this flexible strategy exploits edge-specific prior knowledge and selects sparse and interpretable models. Through simulated data generated under various scenarios, we demonstrate that BANS outperforms other existing multivariate regression-based methodologies. Our integrative genomic network analysis for key signaling pathways across multiple cancer types highlights commonalities and differences of p53 integrative networks and epigenetic effects of BRCA2 on p53 and its interaction with T68 phosphorylated CHK2, that may have translational utilities of finding biomarkers and therapeutic targets.


title: Curves generating extremal rays in blowups of weighted projective planes

abstract: We consider blowups at a general point of weighted projective planes and, more generally, of toric surfaces with Picard number one. We give a unifying construction of negative curves on these blowups such that all previously known families appear as boundary cases of this. The classification consists of two classes of said curves, each depending on two parameters. Every curve in these two classes is algebraically related to other curves in both classes; this allows us to find their defining equations inductively. For each curve in our classification, we consider a family of blowups in which the curve defines an extremal class in the effective cone. We give a complete classification of these blowups into Mori Dream Spaces and non-Mori Dream Spaces. Our approach greatly simplifies previous proofs, avoiding positive characteristic methods and higher cohomology.


title: Jet quenching in a multi-stage Monte Carlo approach

abstract: We present a jet quenching model within a unified multi-stage framework and demonstrate for the first time a simultaneous description of leading hadrons, inclusive jets, and elliptic flow observables which spans multiple centralities and collision energies. This highlights one of the major successes of the JETSCAPE framework in providing a tool for setting up an effective parton evolution that includes a high-virtuality radiation dominated energy loss phase (MATTER), followed by a low-virtuality scattering dominated (LBT) energy loss phase. Measurements of jet and charged-hadron $R_{AA}$ set strong constraints on the jet quenching model. Jet-medium response is also included through a weakly-coupled transport description.


title: Agnostic Q-learning with Function Approximation in Deterministic Systems: Tight Bounds on Approximation Error and Sample Complexity

abstract: The current paper studies the problem of agnostic $Q$-learning with function approximation in deterministic systems where the optimal $Q$-function is approximable by a function in the class $\mathcal{F}$ with approximation error $\delta \ge 0$. We propose a novel recursion-based algorithm and show that if $\delta = O\left(\rho/\sqrt{\dim_E}\right)$, then one can find the optimal policy using $O\left(\dim_E\right)$ trajectories, where $\rho$ is the gap between the optimal $Q$-value of the best actions and that of the second-best actions and $\dim_E$ is the Eluder dimension of $\mathcal{F}$. Our result has two implications: 1) In conjunction with the lower bound in [Du et al., ICLR 2020], our upper bound suggests that the condition $\delta = \widetilde{\Theta}\left(\rho/\sqrt{\mathrm{dim}_E}\right)$ is necessary and sufficient for algorithms with polynomial sample complexity. 2) In conjunction with the lower bound in [Wen and Van Roy, NIPS 2013], our upper bound suggests that the sample complexity $\widetilde{\Theta}\left(\mathrm{dim}_E\right)$ is tight even in the agnostic setting. Therefore, we settle the open problem on agnostic $Q$-learning proposed in [Wen and Van Roy, NIPS 2013]. We further extend our algorithm to the stochastic reward setting and obtain similar results.


title: Pressure and temperature dependence of solubility and surface adsorption of nitrogen in the liquid hydrocarbon bodies on Titan

abstract: We have studied the pressure and temperature dependence of solubility of nitrogen in methane and ethane using vapor-liquid equilibrium simulations of binary mixtures of nitrogen in methane and ethane for a range of pressures between 1.5 atm and 3.5 atm and temperatures between 90 K and 110 K, thermodynamic conditions that may exist on the Saturn's moon, Titan. We find that the solubility of nitrogen in methane increases linearly with pressure while the solubility of nitrogen in ethane increases exponentially with pressure at temperature 90 K. Solubility of nitrogen in both methane and ethane exhibits an exponential decrease with temperature at a pressure of 3 atm. The solubility of nitrogen in methane is much larger compared to that in ethane in the range of pressure and temperature studied here. Our results are in quantitative agreement with the available experimental measurements of the solubility of nitrogen in methane and ethane. Furthermore, we find that the surface adsorption of nitrogen increases with increasing pressure at temperature 90 K, while the adsorption free energy increases with increasing pressure. Moreover, we find that the surface tension decreases linearly with pressure for both nitrogen-methane and nitrogen-ethane systems. The rate of decrease of surface tension with pressure for nitrogen-ethane system is much larger as compared to the nitrogen-methane system. Finally, we find that the absorption of a nitrogen molecule into the liquid-phase from the interface is diffusive and does not involve any appreciable energy barrier. Our results suggest that homogeneous nucleation of bubbles is unlikely on Titan and the bubble formation in the lakes on Titan must arise from heterogeneous nucleation of bubbles.


title: Compactifications of moduli of elliptic K3 surfaces: stable pair and toroidal

abstract: We describe two geometrically meaningful compactifications of the moduli space of elliptic K3 surfaces via stable slc pairs, for two different choices of a polarizing divisor, and show that their normalizations are two different toroidal compactifications of the moduli space, one for the ramification divisor and another for the rational curve divisor. In the course of the proof, we further develop the theory of integral affine spheres with 24 singularities. We also construct moduli of rational (generalized) elliptic stable slc surfaces of types ${\bf A_n}$ ($n\ge1$), ${\bf C_n}$ ($n\ge0$) and ${\bf E_n}$ ($n\ge0$).


title: Disease State Prediction From Single-Cell Data Using Graph Attention Networks

abstract: Single-cell RNA sequencing (scRNA-seq) has revolutionized biological discovery, providing an unbiased picture of cellular heterogeneity in tissues. While scRNA-seq has been used extensively to provide insight into both healthy systems and diseases, it has not been used for disease prediction or diagnostics. Graph Attention Networks (GAT) have proven to be versatile for a wide range of tasks by learning from both original features and graph structures. Here we present a graph attention model for predicting disease state from single-cell data on a large dataset of Multiple Sclerosis (MS) patients. MS is a disease of the central nervous system that can be difficult to diagnose. We train our model on single-cell data obtained from blood and cerebrospinal fluid (CSF) for a cohort of seven MS patients and six healthy adults (HA), resulting in 66,667 individual cells. We achieve 92 % accuracy in predicting MS, outperforming other state-of-the-art methods such as a graph convolutional network and a random forest classifier. Further, we use the learned graph attention model to get insight into the features (cell types and genes) that are important for this prediction. The graph attention model also allow us to infer a new feature space for the cells that emphasizes the differences between the two conditions. Finally we use the attention weights to learn a new low-dimensional embedding that can be visualized. To the best of our knowledge, this is the first effort to use graph attention, and deep learning in general, to predict disease state from single-cell data. We envision applying this method to single-cell data for other diseases.


title: The existence of minimizers for an isoperimetric problem with Wasserstein penalty term in unbounded domains

abstract: In this article, we consider the (double) minimization problem $$\min\left\{P(E;\Omega)+\lambda W_p(E,F):~E\subseteq\Omega,~F\subseteq \mathbb{R}^d,~\lvert E\cap F\rvert=0,~ \lvert E\rvert=\lvert F\rvert=1\right\},$$ where $p\geqslant 1$, $\Omega$ is a (possibly unbounded) domain in $\mathbb{R}^d$, $P(E;\Omega)$ denotes the relative perimeter of $E$ in $\Omega$ and $W_p$ denotes the $p$-Wasserstein distance. When $\Omega$ is unbounded and $d\geqslant 3$, it is an open problem proposed by Buttazzo, Carlier and Laborde in the paper ON THE WASSERSTEIN DISTANCE BETWEEN MUTUALLY SINGULAR MEASURES. We prove the existence of minimizers to this problem when $\frac{1}{p}+\frac{2}{d}>1$, $\Omega=\mathbb{R}^d$ and $\lambda$ is sufficiently small.


title: The rotating black hole interior: Insights from gravitational collapse in $AdS_3$ spacetime

abstract: We present results from a numerical study of rotating black hole formation in 3-dimensional asymptotically anti-de Sitter (AdS) spacetime, focusing on the structure of the black hole interior. While black holes in $AdS_3$ are of theoretical interest for a wide variety of reasons, we choose to study this system primarily as a toy model for astrophysical (4-dimensional) black holes formed from gravitational collapse. We investigate the effect of angular momentum on the geometry inside the event horizon, and see qualitative changes in the interior structure as a function of the spin parameter. For low spins, we find that a central spacelike curvature singularity forms, connecting to a singular, null Cauchy horizon. For spins above a threshold consistent with the linear analysis of Dias, Reall and Santos, curvature on the Cauchy horizon remains bounded, signaling a violation of the strong cosmic censorship conjecture. Further increasing the spin leads to a decrease in the relative size of the spacelike branch of the singularity, which vanishes completely above a second threshold. In these high-spin cases, the interior evolution is bounded by a regular Cauchy horizon, which extends all the way inward to a regular, timelike origin. We further explore the geodesic focusing ("gravitational shock-wave") effect predicted to occur along the outgoing branch of the inner horizon, first described by Marolf and Ori. Remarkably, we observe the effect at late times in all of the black holes we form, even those in which the inner apparent horizon collapses to zero radius early in their evolution.


title: A long neck principle for Riemannian spin manifolds with positive scalar curvature

abstract: We develop index theory on compact Riemannian spin manifolds with boundary in the case when the topological information is encoded by bundles which are supported away from the boundary. As a first application, we establish a "long neck principle" for a compact Riemannian spin $n$-manifold with boundary $X$, stating that if $\textrm{scal}(X)\geq n(n-1)$ and there is a nonzero degree map into the sphere $f\colon X\to S^n$ which is strictly area decreasing, then the distance between the support of $\textrm{d} f$ and the boundary of $X$ is at most $\pi/n$. This answers, in the spin setting and for strictly area decreasing maps, a question recently asked by Gromov. As a second application, we consider a Riemannian manifold $X$ obtained by removing $k$ pairwise disjoint embedded $n$-balls from a closed spin $n$-manifold $Y$. We show that if $\textrm{scal}(X)>\sigma>0$ and $Y$ satisfies a certain condition expressed in terms of higher index theory, then the radius of a geodesic collar neighborhood of $\partial X$ is at most $\pi \sqrt{(n-1)/(n\sigma)}$. Finally, we consider the case of a Riemannian $n$-manifold $V$ diffeomorphic to $N\times [-1,1]$, with $N$ a closed spin manifold with nonvanishing Rosenberg index. In this case, we show that if $\textrm{scal}(V)\geq\sigma>0$, then the distance between the boundary components of $V$ is at most $2\pi \sqrt{(n-1)/(n\sigma)}$. This last constant is sharp by an argument due to Gromov.


title: Blocking of martensitic transition at the nano-scale in the Ti2NiCu wedge

abstract: Shape memory effect associated with martensitic transformations is of the rapidly developing field in nanotechnologies, where industrial use of systems established on that effect provide greater flexibility on the nano-devices fabrication of various kind. And therefore it addresses questions to the phase transition phenomena at low-scale and its limitations and control. In this report, we studied the crystal structure of tapered plates of Ti$_2$NiCu alloy and the temperature $T_c$ at which the martensitic transition occurs. We demonstrated that $T_c$ has a strongly descending character as a function of the plate thickness $h$. The critical thickness value at which the transition completely suppressed is 20 nm. Moreover, the obtained results for $T_c(h)$ curves indicate the hysteretic nature of the transition. These findings open the pathway for size limits indication and regimes modulation where the alloy-based nano-mechanical devices can be tuned to operate more efficiently.


title: MOND-like Fractional Laplacian Theory

abstract: I provide a derivation of some characteristic effects of Milgrom's modified Newtonian dynamics (MOND) from a fractional version of Newton's theory based on the fractional Poisson equation. I employ the properties of the fractional Laplacian to investigate the features of the fundamental solution of the proposed model. The key difference between MOND and the fractional theory introduced here is that the latter is an inherently linear theory, featuring a characteristic length scale $\ell$, whilst the former is ultimately nonlinear in nature and it is characterized by an acceleration scale $a_0$. Taking advantage of the Tully-Fisher relation, as the fractional order $s$ approaches $3/2$, I then connect the length scale $\ell$, emerging from this modification of Newton's gravity, with the critical acceleration $a_0$ of MOND. Finally, implications for galaxy rotation curves of a variable-order version of the model are discussed.


title: Ramsey numbers of partial order graphs (comparability graphs) and implications in ring theory

abstract: For a partially ordered set $(A, \le)$, let $G_A$ be the simple, undirected graph with vertex set $A$ such that two vertices $a \neq b\in A$ are adjacent if either $a \le b$ or $b \le a$. We call $G_A$ the \emph{partial order graph} or \emph{comparability graph} of $A$. Further, we say that a graph $G$ is a partial order graph if there exists a partially ordered set $A$ such that $G = G_A$. For a class $\mathcal{C}$ of simple, undirected graphs and $n$, $m \ge 1$, we define the Ramsey number $\mathcal{R}_{\mathcal{C}}(m,n)$ with respect to $\mathcal{C}$ to be the minimal number of vertices $r$ such that every induced subgraph of an arbitrary partial order graph consisting of $r$ vertices contains either a complete $n$-clique $K_n$ or an independent set consisting of $m$ vertices. In this paper, we determine the Ramsey number with respect to some classes of partial order graphs. Furthermore, some implications of Ramsey numbers in ring theory are discussed.


title: On the Relationship between Software Complexity and Security

abstract: This work aims at discussing the complexity aspect of software while demonstrating its relationship with security. Complexity is an essential part of software; however, numerous studies indicate that they increase the vulnerability of the software systems and introduce bugs in the program. Many developers face difficulty when trying to understand the complex components of software. Complexity in software increases when objects in the software are used to design a more complex object while creating a hierarchical complexity in the system. However, it is necessary for the developers to strive for minimum complexity, as increased complexity introduces security risks in the software, which can cause severe monetary and reputational damage to a government or a private organization. It even causes bodily harm to human beings with various examples found in previous years where security breaches led to severe consequences. Hence it is vital to maintain low complexity and simple design of structure. Various developers tend to introduce deliberate complexities in the system so that they do not have to write the same program twice; however, it is getting problematic for the software organizations as the demands of security are continually increasing.


title: Precision Gating: Improving Neural Network Efficiency with Dynamic Dual-Precision Activations

abstract: We propose precision gating (PG), an end-to-end trainable dynamic dual-precision quantization technique for deep neural networks. PG computes most features in a low precision and only a small proportion of important features in a higher precision to preserve accuracy. The proposed approach is applicable to a variety of DNN architectures and significantly reduces the computational cost of DNN execution with almost no accuracy loss. Our experiments indicate that PG achieves excellent results on CNNs, including statically compressed mobile-friendly networks such as ShuffleNet. Compared to the state-of-the-art prediction-based quantization schemes, PG achieves the same or higher accuracy with 2.4$\times$ less compute on ImageNet. PG furthermore applies to RNNs. Compared to 8-bit uniform quantization, PG obtains a 1.2% improvement in perplexity per word with 2.7$\times$ computational cost reduction on LSTM on the Penn Tree Bank dataset. Code is available at: https://github.com/cornell-zhang/dnn-gating


title: Searching for heavy Higgs in supersymmetric final states at the LHC

abstract: In this work, we analyse and demonstrate possible strategies to explore extended Higgs sector of the Minimal Supersymmetric Standard Model (MSSM). In particular we concentrate on heavy Higgs decays to electroweakinos. We analyse the Higgs to electroweakino decays in the allowed MSSM parameter space after taking into account 13 TeV LHC searches for supersymmetric particles and phenomenological constraints such as flavour physics, Higgs measurements and dark matter constraints. We explore some novel aspects of these Higgs decays. The final states resulting from Higgs to electroweakino decays will have backgrounds arising from the Standard Model as well as direct electroweakino production at the LHC. We demonstrate explicit kinematical differences between Higgs to electroweakino decays and associated backgrounds. Furthermore, we demonstrate for a few specific example points, optimised analysis search strategies at the high luminosity LHC (HL-LHC) run. Finally, we comment on possible search strategies for heavy Higgs decays to exotic final states, where the lightest chargino is long lived and leads to a disappearing track at the LHC.


title: Pass-Efficient Randomized LU Algorithms for Computing Low-Rank Matrix Approximation

abstract: Low-rank matrix approximation is extremely useful in the analysis of data that arises in scientific computing, engineering applications, and data science. However, as data sizes grow, traditional low-rank matrix approximation methods, such as SVD and CPQR, are either prohibitively expensive or cannot provide sufficiently accurate results. A solution is to use randomized low-rank matrix approximation methods such as randomized SVD , and randomized LU decomposition on extremely large data sets. In this paper, we focus on the randomized LU decomposition method. First, we employ a reorthogonalization procedure to perform the power iteration of the existing randomized LU algorithm to compensate for the rounding errors caused by the power method. Then we propose a novel randomized LU algorithm, called PowerLU, for the fixed low-rank approximation problem. PowerLU allows for an arbitrary number of passes of the input matrix, $v \geq 2$. Recall that the existing randomized LU decomposition only allows an even number of passes. We prove the theoretical relationship between PowerLU and the existing randomized LU. Numerical experiments show that our proposed PowerLU is generally faster than the existing randomized LU decomposition, while remaining accurate. We also propose a version of PowerLU, called PowerLU_FP, for the fixed precision low-rank matrix approximation problem. PowerLU_FP is based on an efficient blocked adaptive rank determination Algorithm 4.1 proposed in this paper. We present numerical experiments that show that PowerLU_FP can achieve almost the same accuracy and is faster than the randomized blocked QB algorithm by Martinsson and Voronin. We finally propose a single-pass algorithm based on LU factorization. Tests show that the accuracy of our single-pass algorithm is comparable with the existing single-pass algorithms.


title: Propagation-assisted generation of intense few-femtosecond high-harmonic pulses

abstract: The ongoing development of intense high-harmonic generation (HHG) sources has recently enabled highly nonlinear ionization of atoms by the absorption of at least 10 extreme-ultraviolet (XUV) photons within a single atom [Senfftleben \textit{et al.}, arXiv1911.01375]. Here we investigate the role that reshaping of the fundamental, few-cycle, near-infrared (NIR) driving laser within the 30-cm-long HHG Xe medium plays in the generation of the intense HHG pulses. Using an incident NIR intensity that is higher than what is required for phase-matched HHG, signatures of reshaping are found by measuring the NIR blueshift and the fluorescence from the HHG medium along the propagation axis. These results are well reproduced by numerical calculations that show temporal compression of the NIR pulses in the HHG medium. The simulations predict that after refocusing an XUV beam waist radius of 320 nm and a clean attosecond pulse train can be obtained in the focal plane, with an estimated XUV peak intensity of 9x10^15 W/cm^2. Our results show that XUV intensities that were previously only available at large-scale facilities can now be obtained using moderately powerful table-top light sources.


title: On the eccentricity energy of complete mutipartite graph

abstract: The eccentricity (anti-adjacency) matrix $\varepsilon(G)$ of a graph $G$ is obtained from the distance matrix by retaining the eccentricities in each row and each column. The $\varepsilon$-eigenvalues of a graph $G$ are those of its eccentricity matrix $\varepsilon(G),$ and the eccentricity energy (or the $\varepsilon$-energy) of $G$ is the sum of the absolute values of $\varepsilon$-eigenvalues. In this paper, we establish some bounds for the $\varepsilon$-energy of the complete multipartite graph $K_{n_1, n_2, \ldots,n_p}$ of order $n= \sum_{i=1}^p n_i $ and characterize the extreme graphs. This partially answers the problem given in Wang {\em et al.} (2019). We finish the paper showing graphs that are not $\varepsilon$-cospectral with the same $\varepsilon$-energy.


title: Subset Sampling For Progressive Neural Network Learning

abstract: Progressive Neural Network Learning is a class of algorithms that incrementally construct the network's topology and optimize its parameters based on the training data. While this approach exempts the users from the manual task of designing and validating multiple network topologies, it often requires an enormous number of computations. In this paper, we propose to speed up this process by exploiting subsets of training data at each incremental training step. Three different sampling strategies for selecting the training samples according to different criteria are proposed and evaluated. We also propose to perform online hyperparameter selection during the network progression, which further reduces the overall training time. Experimental results in object, scene and face recognition problems demonstrate that the proposed approach speeds up the optimization procedure considerably while operating on par with the baseline approach exploiting the entire training set throughout the training process.


title: The continuum parabolic Anderson model with a half-Laplacian and periodic noise

abstract: We construct solutions of a renormalized continuum fractional parabolic Anderson model, formally given by $\partial_t u=-(-\Delta)^{1/2}u+\xi u$, where $\xi$ is a periodic spatial white noise. To be precise, we construct limits as $\varepsilon\to 0$ to solutions of $\partial_t u_\varepsilon=-(-\Delta)^{1/2}u_\varepsilon+(\xi_\varepsilon-C_\varepsilon)u_\varepsilon$, where $\xi_\varepsilon$ is a mollification of $\xi$ at scale $\varepsilon$ and $C_\varepsilon$ is a logarithmically diverging renormalization constant. We use a simple renormalization scheme based on that of Hairer and Labb\'e, "A simple construction of the continuum parabolic Anderson model on $\mathbf{R}^{2}$."