title: Towards direct numerical simulation of turbulent co-current Taylor bubble flow

abstract: This paper present a simulation strategy for DNS of turbulent co-current Taylor bubble flow. This is a continuation of the work presented in [2] in which Large Eddy Simulation (LES) of co-current turbulent Taylor bubble flow was presented. It was observed that one of the main challenges is the physically accurate prediction of the behavior of the Taylor bubble skirt, and the related bubble shedding. An underestimation of the turbulent fluctuations in the wake of the Taylor bubble was observed in the LES results. It was suggested that this is related to over-prediction of the loss of void of the Taylor bubble due to bubble shedding induced by an LES mesh resolution which is not sufficient to capture the break-up and bubble formation accurately. To counter this, in the current work we present a DNS approach of co-current turbulent Taylor bubble flow called RK-Basilisk, based on the Basilisk code with local adaptive grid refinement. This strategy allows for very high mesh resolution near the bubble's interface while elsewhere the grid is allowed to be coarser. Each time step, the mesh is adapted based on a void fraction criterion. Basilisk's underlying data structure which is based on an 'octree' allows for much faster solution procedures and, therefore, a much greater number of grid points as compared to the LES simulations which were performed using the more general OpenFOAM code. We compare the results against experimental data of the same setting, as well as the co-current LES OpenFOAM results. The setting of turbulent Taylor bubble flow in co-current conditions allows for both lower order turbulence model development and the validation of more general two-phase modeling strategies. Taylor bubble flow in itself also bears relevance to specific two-phase flow situations. The current work contributes to an advancement in simulation capability for such situations.


title: Velocity averaging and H\"older regularity for kinetic Fokker-Planck equations with general transport operators and rough coefficients

abstract: This article addresses the local boundedness and H\"older continuity of weak solutions to kinetic Fokker-Planck equations with general transport operators and rough coefficients. These results are due to the mixing effect of diffusion and transport. Although the equation is parabolic only in the velocity variable, it has a hypoelliptic structure provided that the transport part $\partial_t+b(v)\cdot\nabla_x$ is nondegenerate in some sense. We achieve the results by revisiting the method, proposed by Golse, Imbert, Mouhot and Vasseur in the case $b(v)= v$, that combines the elliptic De Giorgi-Nash-Moser theory with velocity averaging lemmas.


title: CSTNet: A Dual-Branch Convolutional Network for Imaging of Reactive Flows using Chemical Species Tomography

abstract: Chemical Species Tomography (CST) has been widely used for in situ imaging of critical parameters, e.g. species concentration and temperature, in reactive flows. However, even with state-of-the-art computational algorithms the method is limited due to the inherently ill-posed and rank-deficient tomographic data inversion, and by high computational cost. These issues hinder its application for real-time flow diagnosis. To address them, we present here a novel CST-based convolutional neural Network (CSTNet) for high-fidelity, rapid, and simultaneous imaging of species concentration and temperature. CSTNet introduces a shared feature extractor that incorporates the CST measurement and sensor layout into the learning network. In addition, a dual-branch architecture is proposed for image reconstruction with crosstalk decoders that automatically learn the naturally correlated distributions of species concentration and temperature. The proposed CSTNet is validated both with simulated datasets, and with measured data from real flames in experiments using an industry-oriented sensor. Superior performance is found relative to previous approaches, in terms of robustness to measurement noise and millisecond-level computing time. This is the first time, to the best of our knowledge, that a deep learning-based algorithm for CST has been experimentally validated for simultaneous imaging of multiple critical parameters in reactive flows using a low-complexity optical sensor with severely limited number of laser beams.


title: A Combinatorial Characterization of Self-Stabilizing Population Protocols

abstract: We fully characterize self-stabilizing functions in population protocols for complete interaction graphs. In particular, we investigate self-stabilization in systems of $n$ finite state agents in which a malicious scheduler selects an arbitrary sequence of pairwise interactions under a global fairness condition. We show a necessary and sufficient condition for self-stabilization. Specifically we show that functions without certain set-theoretic conditions are impossible to compute in a self-stabilizing manner. Our main contribution is in the converse, where we construct a self-stabilizing protocol for all other functions that meet this characterization. Our positive construction uses Dickson's Lemma to develop the notion of the root set, a concept that turns out to fundamentally characterize self-stabilization in this model. We believe it may lend to characterizing self-stabilization in more general models as well.


title: Improved approximation ratios for two Euclidean maximum spanning tree problems

abstract: We study the following two maximization problems related to spanning trees in the Euclidean plane. It is not known whether or not these problems are NP-hard. We present approximation algorithms with better approximation ratios for both problems. The improved ratios are obtained mainly by employing the Steiner ratio, which has not been used in this context earlier. (i) Longest noncrossing spanning tree: Given a set of points in the plane, the goal is to find a maximum-length noncrossing spanning tree. Alon, Rajagopalan, and Suri (SoCG 1993) studied this problem for the first time and gave a $0.5$-approximation algorithm. Over the years, the approximation ratio has been successively improved to $0.502$, $0.503$, and to $0.512$ which is the current best ratio, due to Cabello et al.. We revisit this problem and improve the ratio further to $0.519$. The improvement is achieved by a collection of ideas, some from previous works and some new ideas (including the use of the Steiner ratio), along with a more refined analysis. (ii) Longest spanning tree with neighborhoods: Given a collection of regions (neighborhoods) in the plane, the goal is to select a point in each neighborhood so that the longest spanning tree on selected points has maximum length. We present an algorithm with approximation ratio $0.524$ for this problem. The previous best ratio, due to Chen and Dumitrescu, is $0.511$ which is in turn the first improvement beyond the trivial ratio $0.5$. Our algorithm is fairly simple, its analysis is relatively short, and it takes linear time after computing a diametral pair of points. The simplicity comes from the fact that our solution belongs to a set containing three stars and one double-star. The shortness and the improvement come from the use of the Steiner ratio.


title: Linear Breit-Wheeler pair production by high-energy bremsstrahlung photons colliding with an intense X-ray laser pulse

abstract: A possible setup for the experimental verification of linear Breit-Wheeler pair creation of electrons and positrons in photon-photon collisions is studied theoretically. It combines highly energetic bremsstrahlung photons, which are assumed to be generated by an incident beam of GeV electrons penetrating through a high-$Z$ target, with keV photons from an X-ray laser field, which is described as a focused Gaussian pulse. We discuss the dependencies of the pair yields on the incident electron energy, target thickness, laser parameters, and collision geometry. It is shown that, for suitable conditions which are nowadays in reach at X-ray laser facilities, the resulting number of created particles seems to be well accessible for enabling the first experimental observation of the linear Breit-Wheeler process $\gamma\gamma\to e^+e^-$.


title: Clinically Verified Hybrid Deep Learning System for Retinal Ganglion Cells Aware Grading of Glaucomatous Progression

abstract: Objective: Glaucoma is the second leading cause of blindness worldwide. Glaucomatous progression can be easily monitored by analyzing the degeneration of retinal ganglion cells (RGCs). Many researchers have screened glaucoma by measuring cup-to-disc ratios from fundus and optical coherence tomography scans. However, this paper presents a novel strategy that pays attention to the RGC atrophy for screening glaucomatous pathologies and grading their severity. Methods: The proposed framework encompasses a hybrid convolutional network that extracts the retinal nerve fiber layer, ganglion cell with the inner plexiform layer and ganglion cell complex regions, allowing thus a quantitative screening of glaucomatous subjects. Furthermore, the severity of glaucoma in screened cases is objectively graded by analyzing the thickness of these regions. Results: The proposed framework is rigorously tested on publicly available Armed Forces Institute of Ophthalmology (AFIO) dataset, where it achieved the F1 score of 0.9577 for diagnosing glaucoma, a mean dice coefficient score of 0.8697 for extracting the RGC regions and an accuracy of 0.9117 for grading glaucomatous progression. Furthermore, the performance of the proposed framework is clinically verified with the markings of four expert ophthalmologists, achieving a statistically significant Pearson correlation coefficient of 0.9236. Conclusion: An automated assessment of RGC degeneration yields better glaucomatous screening and grading as compared to the state-of-the-art solutions. Significance: An RGC-aware system not only screens glaucoma but can also grade its severity and here we present an end-to-end solution that is thoroughly evaluated on a standardized dataset and is clinically validated for analyzing glaucomatous pathologies.


title: A new wake-merging method for wind-farm power prediction in presence of heterogeneous background velocity fields

abstract: Many wind farms are placed near coastal regions or in proximity of orographic obstacles. The meso-scale gradients that develop in these zones make wind farms operating in velocity fields that are rarely uniform. However, all existing wake-merging methods in engineering wind-farm wake models assume a homogeneous background velocity field in and around the farm, relying on a single wind-speed value usually measured several hundreds of meters upstream of the first row of turbines. In this study, we derive a new momentum-conserving wake-merging method capable of superimposing the waked flow on a heterogeneous background velocity field. We couple the proposed wake-merging method with four different wake models, i.e. the Gaussian, super-Gaussian, double-Gaussian and Ishihara model, and we test its performance against LES data, dual-Doppler radar measurements and SCADA data from the Horns Rev, London Array, and Westermost Rough farm. Next to this, as an additional point of reference, the standard Jensen model with quadratic superposition is also included. Results show that the new method performs similarly to linear superposition of velocity deficits in homogeneous conditions but it shows better performance when a spatially varying background velocity is used. The most accurate estimates are obtained when the wake-merging method is coupled with the double-Gaussian and Gaussian single-wake model. The Ishihara model also shows good agreements with observations. In contrast to this, the Jensen and super-Gaussian wake model underestimate the farm power output for all wind speeds, wind directions and wind farms considered in our analysis.


title: Simulation of reactive flows using particle methods

abstract: We describe a new computational method for the numerically stable particle-based simulation of open-boundary flows, including volume conserving chemical reactions. The novel method is validated for the case of heterogeneous catalysis against a reliable reference simulation and is shown to deliver identical results while the computational efficiency is significantly increased.


title: Supervised learning of few dirty bosons with variable particle number

abstract: We investigate the supervised machine learning of few interacting bosons in optical speckle disorder via artificial neural networks. The learning curve shows an approximately universal power-law scaling for different particle numbers and for different interaction strengths. We introduce a network architecture that can be trained and tested on heterogeneous datasets including different particle numbers. This network provides accurate predictions for all system sizes included in the training set and, by design, is suitable to attempt extrapolations to (computationally challenging) larger sizes. Notably, a novel transfer-learning strategy is implemented, whereby the learning of the larger systems is substantially accelerated and made consistently accurate by including in the training set many small-size instances.


title: Phase behaviour of Lennard-Jones particles in two dimensions

abstract: The phase diagram of the prototypical two-dimensional Lennard-Jones system, while extensively investigated, is still debated. In particular, there are controversial results in the literature as concern the existence of the hexatic phase and the melting scenario. Here, we study the phase behaviour of 2D LJ particles via large-scale numerical simulations. We demonstrate that at high temperature, when the attraction in the potential plays a minor role, melting occurs via a continuous solid-hexatic transition followed by a first-order hexatic-fluid transition. As the temperature decreases, the density range where the hexatic phase occurs shrinks so that at low-temperature melting occurs via a first-order liquid-solid transition. The temperature where the hexatic phase disappears is well above the liquid-gas critical temperature. The evolution of the density of topological defects confirms this scenario.


title: The X-ray properties of Eta Carinae during its 2020 X-ray minimum

abstract: The massive binary system Eta Carinae is characterized by intense colliding winds that form shocks and emit X-rays. The system is highly eccentric ($e\simeq0.9$), resulting in modulated X-ray emission during its 5.54 year orbit. The X-ray flux increases in the months prior to periastron passage, exhibiting strong flares, then rapidly declines to a flat minimum lasting a few weeks, followed by a gradual recovery. We present Neutron Star Interior Composition Explorer (NICER) telescope spectra obtained before, during, and after the 2020 X-ray minimum, and perform spectral analysis to establish the temporal behavior of X-ray flux and X-ray-absorbing column density ($N_{\rm H}(t)$) for the 2-10 keV and 5-10 keV energy ranges. The latter range is dominated by the stellar wind collision region and, therefore, these spectral parameters - in particular, $N_{\rm H}(t)$ - serves as a potentially stringent constraint on the binary orientation. We compare the observed $N_{\rm H}(t)$ results to the behavior predicted by a simple geometrical model in an attempt to ascertain which star is closer to us at periastron: the more massive primary ($\omega \simeq 240$-$270^\circ$), or the secondary ($\omega \simeq 90^\circ$). We find that the variations in column density, both far from periastron and around periastron passage, support the latter configuration ($\omega \simeq 90^\circ$). The 2020 X-ray minimum showed the fastest recovery among the last five minima, providing additional evidence for a recent weakening of the primary star's wind.


title: Silicon Flexes Muscles: Giant Electrochemical Actuation in a Nanoporous Silicon-Polypyrrole Hybrid Material

abstract: The absence of piezoelectricity in silicon makes direct electro-mechanical applications of this mainstream semiconductor impossible. Integrated electrical control of the silicon mechanics, however, would open up new perspectives for on-chip actuorics. Here, we combine wafer-scale nanoporosity in single-crystalline silicon with polymerization of an artificial muscle material inside pore space to synthesize a composite that shows macroscopic electrostrain in aqueous electrolyte. The voltage-strain coupling is 3 orders of magnitude larger than the best-performing ceramics in terms of piezoelectric actuation. We trace this huge electroactuation to the concerted action of 100 billions of nanopores per square centimetre cross-section and to potential-dependent pressures of up to 150 atmospheres at the single-pore scale. The exceptionally small operation voltages (0.4-0.9 V) along with the sustainable and biocompatible base materials make this hybrid promising for bio-actuator applications.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.soft physics.app-ph physics.chem-ph


title: Floquet higher order topological insulator in a periodically driven bipartite lattice

abstract: Floquet higher order topological insulators (FHOTIs) are a novel topological phase that can occur in periodically driven lattices. An appropriate experimental platform to realize FHOTIs has not yet been identified. We introduce a periodically-driven bipartite (two-band) system that hosts FHOTI phases, and predict that this lattice can be realized in experimentally-realistic optical waveguide arrays, similar to those previously used to study anomalous Floquet insulators. The model exhibits interesting phase transitions from first-order to second-order topological matter by tuning a coupling strength parameter, without breaking lattice symmetry. In the FHOTI phase, the lattice hosts corner modes at eigenphase $0$ or $\pi$, which are robust against disorder in the individual couplings.


title: A Co-Interactive Transformer for Joint Slot Filling and Intent Detection

abstract: Intent detection and slot filling are two main tasks for building a spoken language understanding (SLU) system. The two tasks are closely related and the information of one task can be utilized in the other task. Previous studies either model the two tasks separately or only consider the single information flow from intent to slot. None of the prior approaches model the bidirectional connection between the two tasks simultaneously. In this paper, we propose a Co-Interactive Transformer to consider the cross-impact between the two tasks. Instead of adopting the self-attention mechanism in vanilla Transformer, we propose a co-interactive module to consider the cross-impact by building a bidirectional connection between the two related tasks. In addition, the proposed co-interactive module can be stacked to incrementally enhance each other with mutual features. The experimental results on two public datasets (SNIPS and ATIS) show that our model achieves the state-of-the-art performance with considerable improvements (+3.4% and +0.9% on overall acc). Extensive experiments empirically verify that our model successfully captures the mutual interaction knowledge.


title: Large Product Key Memory for Pretrained Language Models

abstract: Product key memory (PKM) proposed by Lample et al. (2019) enables to improve prediction accuracy by increasing model capacity efficiently with insignificant computational overhead. However, their empirical application is only limited to causal language modeling. Motivated by the recent success of pretrained language models (PLMs), we investigate how to incorporate large PKM into PLMs that can be finetuned for a wide variety of downstream NLP tasks. We define a new memory usage metric, and careful observation using this metric reveals that most memory slots remain outdated during the training of PKM-augmented models. To train better PLMs by tackling this issue, we propose simple but effective solutions: (1) initialization from the model weights pretrained without memory and (2) augmenting PKM by addition rather than replacing a feed-forward network. We verify that both of them are crucial for the pretraining of PKM-augmented PLMs, enhancing memory utilization and downstream performance. Code and pretrained weights are available at https://github.com/clovaai/pkm-transformers.


title: Characteristic Time Variability of Gravitational-Wave and Neutrino Signals from Three-dimensional Simulations of Non-Rotating and Rapidly Rotating Stellar Core-Collapse

abstract: We present results from full general relativistic three-dimensional hydrodynamics simulations of stellar core collapse of a 70 M$_\odot$ star with spectral neutrino transport. To investigate the impact of rotation on non-axisymmetric instabilities, we compute three models by parametrically changing the initial strength of rotation. The most rapidly rotating model exhibits a transient development of the low-$T/|W|$ instability with one-armed spiral flow at the early postbounce phase. Subsequently, the two-armed spiral flow appears, which persists during the simulation time. The moderately rotating model also shows the growth of the low-$T/|W|$ instability, but only with the two-armed spiral flow. In the nonrotating model, a vigorous activity of the standing accretion-shock instability (SASI) is only observed. The SASI is first dominated by the sloshing mode, which is followed by the spiral SASI until the black hole formation. We present a spectrogram analysis of the gravitational waves (GWs) and neutrinos, focusing on the time correlation. Our results show that characteristic time modulations in the GW and neutrino signals can be linked to the growth of the non-axisymmetric instabilities. We find that the degree of the protoneutron star (PNS) deformation, depending upon which modes of the non-axisymmetric instabilities develop, predominantly affects the characteristic frequencies of the correlated GW and neutrino signals. We point out that these signals would be simultaneously detectable by the current-generation detectors up to $\sim10$ kpc. Our findings suggest that the joint observation of GWs and neutrinos is indispensable for extracting information on the PNS evolution preceding the black hole formation.


title: MHD stability of JT-60SA operation scenarios driven by passing energetic particles for a hot Maxwellian model

abstract: We analyze the effects of the passing energetic particles on the resistive ballooning modes (RBM) and the energetic particle driven modes in JT-60SA plasma, which leads to the prediction of the stability in N-NBI heated plasma. The analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species assuming an averaged Maxwellian EP distribution fitted to the slowing down distribution, including the effect of the acoustic modes. The simulations show the possible destabilization of a 3/2-4/2 TAE with a frequency (f) of 115 kHz, a 6/4-7/4 TAE with f=98 kHz and a 6/4 or 7/4 BAE with f=57 kHz in the ITER-like inductive scenario. If the energetic particle beta increases, Beta induced AEs (BAE), Toroidal AEs (TAE) and Elliptical AEs (EAE) are destabilized between the inner-middle plasma region, leading to the overlapping of AE of different toroidal families. If these instabilities coexist in the non-linear saturation phase the EP transport could be enhanced leading to a lower heating efficiency. For a hypothetical configuration based on the ITER-like inductive scenario but an center peaked EP profile, the EP beta threshold increases and several BAEs are destabilized in the inner plasma region, indicating an improved AE stability with respect to the off-axis peaked EP profile. In addition, the analysis of a hypothetical JT-60SA scenario with a resonant q=1 in the inner plasma region shows the destabilization of fishbones-like instabilities by the off-axis peaked EP profile. Also, the EPs have an stabilizing effect on the RBM, stronger as the population of EP with low energies (below 250 keV) increases at the plasma pedestal.


title: Lattice Bounded Distance Equivalence for 1D Delone Sets with Finite Local Complexity

abstract: Spectra of suitably chosen Pisot-Vijayaraghavan numbers represent non-trivial examples of self-similar Delone point sets of finite local complexity, indispensable in quasicrystal modeling. For the case of quadratic Pisot units we characterize, dependingly on digits in the corresponding numeration systems, the spectra which are bounded distance to an average lattice. Our method stems in interpretation of the spectra in the frame of the cut-and-project method. Such structures are coded by an infinite word over a finite alphabet which enables us to exploit combinatorial notions such as balancedness, substitutions and the spectrum of associated incidence matrices.


title: Partial wave correlations with uncertainties in the Two Nucleon System

abstract: The NN scattering problem is usually analyzed in terms of partial waves and the corresponding coupled channel phase-shifts and mixing angles, but the available experiments induce correlations amongst the corresponding channels with different quantum numbers. Based on the Granada-2013 database we analyze the meaning and impact of those correlations taking into account both the purely statistical ones reflecting the primary experimental data uncertainties as well as the systematic ones exhibiting the ambiguities in the form of the potential representing the unknown nuclear force for distances below 3fm. We also ponder on the implications of fixing nuclear potentials by direct fits to the phase-shifts stemming from a full fledged partial wave analysis (PWA) inferred from experimental data.


title: Rapid and sensitive detection of SARS-CoV-2 with functionalized magnetic nanoparticles

abstract: The outbreak of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) threatens global medical systems and economies, and rules our daily living life. Controlling the outbreak of SARS-CoV-2 has become one of the most important and urgent strategies throughout the whole world. As of October, 2020, there have not yet been any medicines or therapies to be effective against SARS-CoV-2. Thus, rapid and sensitive diagnostics is the most important measures to control the outbreak of SARS-CoV-2. Homogeneous biosensing based on magnetic nanoparticles (MNPs) is one of the most promising approaches for rapid and highly sensitive detection of biomolecules. This paper proposes an approach for rapid and sensitive detection of SARS-CoV-2 with functionalized MNPs via the measurement of their magnetic response in an ac magnetic field. Experimental results demonstrate that the proposed approach allows the rapid detection of mimic SARS-CoV-2 with a limit of detection of 0.084 nM (5.9 fmole). The proposed approach has great potential for designing a low-cost and point-of-care device for rapid and sensitive diagnostics of SARS-CoV-2.


title: Quantum error mitigation as a universal error-minimization technique: applications from NISQ to FTQC eras

abstract: In the early years of fault-tolerant quantum computing (FTQC), it is expected that the available code distance and the number of magic states will be restricted due to the limited scalability of quantum devices and the insufficient computational power of classical decoding units. Here, we integrate quantum error correction and quantum error mitigation into an efficient FTQC architecture that effectively increases the code distance and $T$-gate count at the cost of constant sampling overheads in a wide range of quantum computing regimes. For example, while we need $10^4$ to $10^{10}$ logical operations for demonstrating quantum advantages from optimistic and pessimistic points of view, we show that we can reduce the required number of physical qubits by $80\%$ and $45\%$ in each regime. From another perspective, when the achievable code distance is up to about 11, our scheme allows executing $10^3$ times more logical operations. This scheme will dramatically alleviate the required computational overheads and hasten the arrival of the FTQC era.


title: How Does the Symmetry of S1 Influence Exciton Transport in Conjugated Polymers?

abstract: Many optoelectronic devices based on organic materials require rapid and long-range singlet exciton transport. Key factors that control the transport of singlet excitons includes the electronic structure of the material, disorder and exciton-phonon coupling. An important parameter whose influence on exciton transport has not been explored is the symmetry of the singlet electronic state (S1). Here, we employ femtosecond transient absorption spectroscopy and microscopy to reveal the relationship between the symmetry of S1 and exciton transport in highly aligned, near-disorder free, one-dimensional conjugated polymers based on polydiacetylene.

categories: cond-mat.mes-hall cond-mat.mtrl-sci physics.app-ph physics.chem-ph quant-ph


title: Cornering the axion with CP-violating interactions

abstract: Besides CP-preserving interactions, axions and axion-like particles may also have small CP-violating scalar Yukawa interactions with nucleons and electrons. Any such interaction will generate macroscopic monopole-dipole forces which can be searched for experimentally. When the best experimental limits on scalar interactions are combined with stellar energy-loss arguments constraining pseudoscalar interactions, strong bounds can be set on CP-violating axion couplings which almost intersect the expectation for QCD models. Over the years, both astrophysical and laboratory tests have improved. We provide a much-needed up-to-date compilation of these constraints, showing improvements in some regions of parameter space by factors between 40 and 130. We advocate experimental opportunities, without astrophysical or dark-matter assumptions, to track down the axion in the lesser-explored corners of its parameter space.


title: On the boundedness of infinite matrix products with alternating factors from two sets of matrices

abstract: We consider the question of the boundedness of matrix products $A_{n}B_{n}\cdots A_{1}B_{1}$ with factors from two sets of matrices, $A_{i}\in\mathscr{A}$ and $B_{i}\in\mathscr{B}$, due to an appropriate choice of matrices $\{B_{i}\}$. It is assumed that for any sequence of matrices $\{A_{i}\}$ there is a sequence of matrices $\{B_{i}\}$ for which the sequence of matrix products $\{A_{n}B_{n}\cdots A_{1}B_{1}\}_{n=1}^{\infty}$ is norm bounded. Some situations are described in which in this case the norms of matrix products $A_{n}B_{n}\cdots A_{1}B_{1}$ are uniformly bounded, that is, $\|A_{n}B_{n}\cdots A_{1}B_{1}\|\le C$ for all natural numbers $n$, where $C>0$ is some constant independent of the sequence $\{A_{i}\}$ and the corresponding sequence $\{B_{i}\}$. In the general case, the question of the validity of the corresponding statement remains open.


title: Conditional Goodness-of-Fit Tests for Discrete Distributions

abstract: In this paper, we address the problem of testing goodness-of-fit for discrete distributions, where we focus on the geometric distribution. We define new likelihood-based goodness-of-fit tests using the beta-geometric distribution and the type I discrete Weibull distribution as alternative distributions. The tests are compared in a simulation study, where also the classical goodness-of-fit tests are considered for comparison. Throughout the paper we consider conditional testing given a minimal sufficient statistic under the null hypothesis, which enables the calculation of exact p-values. For this purpose, a new method is developed for drawing conditional samples from the geometric distribution and the negative binomial distribution. We also explain briefly how the conditional approach can be modified for the binomial, negative binomial and Poisson distributions. It is finally noted that the simulation method may be extended to other discrete distributions having the same sufficient statistic, by using the Metropolis-Hastings algorithm.


title: Modeling of the ECCD injection effect on the Heliotron J and LHD plasma stability

abstract: The aim of the study is to analyze the stability of the Energetic Particle Modes (EPM) and Alfven Eigenmodes (AE) in Helitron J and LHD plasma if the electron cyclotron current drive (ECCD) is applied. The analysis is performed using the code FAR3d that solves the reduced MHD equations describing the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moments for the energetic particle (EP) species, including the effect of the acoustic modes. The Landau damping and resonant destabilization effects are added via the closure relation. The simulation results show that the n=1 EPM and n=2 Global AE (GAE) in Heliotron J plasma can be stabilized if the magnetic shear is enhanced at the plasma periphery by an increase (co-ECCD injection) or decrease (ctr-ECCD injection) of the rotational transform at the magnetic axis iota0. In the ctr-ECCD simulations, the EPM/AE growth rate decreases only below a given iota0, similar to the ECCD intensity threshold observed in the experiments. In addition, ctr-ECCD simulations show an enhancement of the continuum damping. The simulations of the LHD discharges with ctr-ECCD injection indicate the stabilization of the n=1 EPM, n=2 Toroidal AE (TAE) and n=3 TAE, caused by an enhancement of the continuum damping in the inner plasma leading to a higher EP beta threshold with respect to the co- and no-ECCD simulations.


title: High-intensity scaling in UV-modified QED

abstract: QED perturbation theory in a background field has been conjectured to break down for sufficiently high field intensity. The high-intensity behavior of a field theory is however intertwined with its high-energy (UV) behavior. Here we show that a UV modification of QED changes the high-intensity behaviour of observables. Specifically we study non-linear Compton scattering in a constant crossed field in QED with an additional Pauli term. In the UV modified theory the cross section exhibits a faster power-law scaling with intensity than in ordinary QED.


title: Intrinsic Hierarchical Clustering Behavior Recovers Higher Dimensional Shape Information

abstract: We show that specific higher dimensional shape information of point cloud data can be recovered by observing lower dimensional hierarchical clustering dynamics. We generate multiple point samples from point clouds and perform hierarchical clustering within each sample to produce dendrograms. From these dendrograms, we take cluster evolution and merging data that capture clustering behavior to construct simplified diagrams that record the lifetime of clusters akin to what zero dimensional persistence diagrams do in topological data analysis. We compare differences between these diagrams using the bottleneck metric, and examine the resulting distribution. Finally, we show that statistical features drawn from these bottleneck distance distributions detect artefacts of, and can be tapped to recover higher dimensional shape characteristics.


title: Zhang-Zhang Polynomials of Ribbons

abstract: We report a closed-form formula for the Zhang-Zhang polynomial (aka ZZ polynomial or Clar covering polynomial) of an important class of elementary pericondensed benzenoids $Rb\left(n_{1},n_{2},m_{1},m_{2}\right)$ usually referred to as ribbons. A straightforward derivation is based on the recently developed interface theory of benzenoids [Langner and Witek, MATCH Commun. Math. Comput. Chem. 84, 143--176 (2020)]. The discovered formula provides compact expressions for various topological invariants of $Rb\left(n_{1},n_{2},m_{1},m_{2}\right)$: the number of Kekul\'e structures, the number of Clar covers, its Clar number, and the number of Clar structures. The last two classes of elementary benzenoids, for which closed-form ZZ polynomial formulas remain to be found, are hexagonal flakes $O\left(k,m,n\right)$ and oblate rectangles $Or\left(m,n\right)$.


title: A practical guide towards agile test-driven development for scientific software projects

abstract: Software testing has received much attention over the last years and has reached such critical importance that agile software development practices put software testing at its core. Agile software development is successfully applied in large-scale industrial software developments but due to its granular responsibilities with roles assigned to various members of the development team, these practices may not be applicable to scientific code development, especially in an academic environment, where it is not uncommon that the codebase is developed, maintained and used by a single person. Even for collaborative scientific software development, financed through external grants, the end-users are typically still part of the development team. This is in contrast to how software is developed in many industries, where the development team and end-users are two separate entities. There are, however, many good code development practices that can be adopted for scientific software projects. Specifically, the intention of this article is to take the centrepiece of agile software development and tailor it to scientific and academic, single-user code development. In this study, a c++ starter project is developed and made available, based on the meson build system, which provides native support for software testing. It is used to show how a simple linear algebra application, found in many scientific and academic applications, can be developed and how simple unit, integration and system tests can be created that are managed through the meson build system. In this way, we are able to minimise software defects and reduce the risk to interpret incorrect data generated by erroneous software that may result in the wrong conclusions to be drawn. Each layer of testing presents one additional layer of protection and we will explore how these may be incorporated with minimum overhead to produce bug-free software.


title: BGM: Building a Dynamic Guidance Map without Visual Images for Trajectory Prediction

abstract: Visual images usually contain the informative context of the environment, thereby helping to predict agents' behaviors. However, they hardly impose the dynamic effects on agents' actual behaviors due to the respectively fixed semantics. To solve this problem, we propose a deterministic model named BGM to construct a guidance map to represent the dynamic semantics, which circumvents to use visual images for each agent to reflect the difference of activities in different periods. We first record all agents' activities in the scene within a period close to the current to construct a guidance map and then feed it to a Context CNN to obtain their context features. We adopt a Historical Trajectory Encoder to extract the trajectory features and then combine them with the context feature as the input of the social energy based trajectory decoder, thus obtaining the prediction that meets the social rules. Experiments demonstrate that BGM achieves state-of-the-art prediction accuracy on the two widely used ETH and UCY datasets and handles more complex scenarios.


title: Consistent Specification Test of the Quantile Autoregression

abstract: This paper proposes a test for the joint hypothesis of correct dynamic specification and no omitted latent factors for the Quantile Autoregression. If the composite null is rejected we proceed to disentangle the cause of rejection, i.e., dynamic misspecification or an omitted variable. We establish the asymptotic distribution of the test statistics under fairly weak conditions and show that factor estimation error is negligible. A Monte Carlo study shows that the suggested tests have good finite sample properties. Finally, we undertake an empirical illustration of modelling GDP growth and CPI inflation in the United Kingdom, where we find evidence that factor augmented models are correctly specified in contrast with their non-augmented counterparts when it comes to GDP growth, while also exploring the asymmetric behaviour of the growth and inflation distributions.


title: Population Based Training for Data Augmentation and Regularization in Speech Recognition

abstract: Varying data augmentation policies and regularization over the course of optimization has led to performance improvements over using fixed values. We show that population based training is a useful tool to continuously search those hyperparameters, within a fixed budget. This greatly simplifies the experimental burden and computational cost of finding such optimal schedules. We experiment in speech recognition by optimizing SpecAugment this way, as well as dropout. It compares favorably to a baseline that does not change those hyperparameters over the course of training, with an 8% relative WER improvement. We obtain 5.18% word error rate on LibriSpeech's test-other.


title: Surface-tension-driven coarsening in mass-conserved reaction-diffusion systems

abstract: Mass conservation in chemical species appears in a broad class of reaction-diffusion systems (RDs) and is known to bring about coarsening of the pattern in chemical concentration. Recent theoretical studies on RDs with mass conservation (MCRDs) reported that the interfacial curvature between two states contributes to the coarsening process, reminiscent of phase separation phenomena. However, since MCRDs do not presuppose a variational principle, it is largely unknown whether description of surface tension is operative or not. In this study, we numerically and theoretically explore the coarsening process of patterns in MCRDs in two and three dimensions. We identify the parameter regions where the homogeneous steady state becomes stable, unstable, and metastable. In the unstable region, pattern formation is triggered by usual Turing instability, whereas in the metastable region, nucleation-growth-type pattern formation is observed. In the later stage, spherical droplet patterns are observed in both regions, where they obey a relation similar to the Young-Laplace law and coarsen following the evaporation-condensation mechanism. These results demonstrate that in the presence of a conserved variable, a physical quantity similar to surface tension is relevant to MCRDs, which provides new insight into molecular self-assembly driven by chemical reactions.


title: Two-dimensional Fourier Continuation and applications

abstract: This paper presents a "two-dimensional Fourier Continuation" method (2D-FC) for construction of bi-periodic extensions of smooth non-periodic functions defined over general two-dimensional smooth domains. The approach can be directly generalized to domains of any given dimensionality, and even to non-smooth domains, but such generalizations are not considered here. The 2D-FC extensions are produced in a two-step procedure. In the first step the one-dimensional Fourier Continuation method is applied along a discrete set of outward boundary-normal directions to produce, along such directions, continuations that vanish outside a narrow interval beyond the boundary. Thus, the first step of the algorithm produces "blending-to-zero along normals" for the given function values. In the second step, the extended function values are evaluated on an underlying Cartesian grid by means of an efficient, high-order boundary-normal interpolation scheme. A Fourier Continuation expansion of the given function can then be obtained by a direct application of the two-dimensional FFT algorithm. Algorithms of arbitrarily high order of accuracy can be obtained by this method. The usefulness and performance of the proposed two-dimensional Fourier Continuation method are illustrated with applications to the Poisson equation and the time-domain wave equation within a bounded domain. As part of these examples the novel "Fourier Forwarding" solver is introduced which, propagating plane waves as they would in free space and relying on certain boundary corrections, can solve the time-domain wave equation and other hyperbolic partial differential equations within general domains at computing costs that grow sublinearly with the size of the spatial discretization.


title: IRX-1D: A Simple Deep Learning Architecture for Remote Sensing Classifications

abstract: We proposes a simple deep learning architecture combining elements of Inception, ResNet and Xception networks. Four new datasets were used for classification with both small and large training samples. Results in terms of classification accuracy suggests improved performance by proposed architecture in comparison to Bayesian optimised 2D-CNN with small training samples. Comparison of results using small training sample with Indiana Pines hyperspectral dataset suggests comparable or better performance by proposed architecture than nine reported works using different deep learning architectures. In spite of achieving high classification accuracy with limited training samples, comparison of classified image suggests different land cover classes are assigned to same area when compared with the classified image provided by the model trained using large training samples with all datasets.


title: Injecting Word Information with Multi-Level Word Adapter for Chinese Spoken Language Understanding

abstract: In this paper, we improve Chinese spoken language understanding (SLU) by injecting word information. Previous studies on Chinese SLU do not consider the word information, failing to detect word boundaries that are beneficial for intent detection and slot filling. To address this issue, we propose a multi-level word adapter to inject word information for Chinese SLU, which consists of (1) sentence-level word adapter, which directly fuses the sentence representations of the word information and character information to perform intent detection and (2) character-level word adapter, which is applied at each character for selectively controlling weights on word information as well as character information. Experimental results on two Chinese SLU datasets show that our model can capture useful word information and achieve state-of-the-art performance.


title: Simulation of catalytic reactions in open-cell foam structures

abstract: We describe a technique for particle-based simulations of heterogeneous catalysis in open-cell foam structures, which is based on isotropic Stochastic Rotation Dynamics (iSRD) together with Constructive Solid Geometry (CSG). The approach is validated by means of experimental results for the low temperature water-gas shift reaction in an open-cell foam structure modeled as inverse sphere packing. Considering the relation between Sherwood and Reynolds number, we find two distinct regimes meeting approximately at the strut size Reynolds number 10. For typical parameters from the literature, we find that the catalyst density in the washcoat can be reduced considerably without a notable loss of conversion efficiency. We vary the porosity to determine optimum open-cell foam structures, which combine low flow resistance with high conversion efficiency and find large porosity values to be favorable not only in the mass transfer limited regime but also in the intermediate regime.


title: HLT-NUS Submission for NIST 2019 Multimedia Speaker Recognition Evaluation

abstract: This work describes the speaker verification system developed by Human Language Technology Laboratory, National University of Singapore (HLT-NUS) for 2019 NIST Multimedia Speaker Recognition Evaluation (SRE). The multimedia research has gained attention to a wide range of applications and speaker recognition is no exception to it. In contrast to the previous NIST SREs, the latest edition focuses on a multimedia track to recognize speakers with both audio and visual information. We developed separate systems for audio and visual inputs followed by a score level fusion of the systems from the two modalities to collectively use their information. The audio systems are based on x-vector based speaker embedding, whereas the face recognition systems are based on ResNet and InsightFace based face embeddings. With post evaluation studies and refinements, we obtain an equal error rate (EER) of 0.88% and an actual detection cost function (actDCF) of 0.026 on the evaluation set of 2019 NIST multimedia SRE corpus.


title: M{\"o}bius-invariant self-avoidance energies for non-smooth sets in arbitrary dimensions

abstract: In the present paper we investigate generalizations of O'Hara's M\"obius energy on curves \cite{ohara_1991a}, to M\"obius-invariant energies on non-smooth subsets of $\R^n$ of arbitrary dimension and co-dimension. In particular, we show under mild assumptions on the local flatness of an admissible possibly unbounded set $\Sigma\subset \R^n$ that locally finite energy implies that $\Sigma$ is, in fact, an embedded Lipschitz submanifold of $\R^n$ -- sometimes even smoother (depending on the a priorily given additional regularity of the admissible set). We also prove, on the other hand, that a local graph structure of low fractional Sobolev regularity on a set $\Sigma$ is already sufficient to guarantee finite energy of $\Sigma$. This type of Sobolev regularity is exactly what one would expect in view of Blatt's characterization \cite{blatt_2012a} of the correct energy space for the M\"obius energy on closed curves. Our results hold in particular for Kusner and Sullivan's cosine energy $E_\textnormal{KS}$ \cite{kusner-sullivan_1997} since one of the energies considered here is equivalent to $E_\textnormal{KS}$.


title: Classification of Speech with and without Face Mask using Acoustic Features

abstract: The understanding and interpretation of speech can be affected by various external factors. The use of face masks is one such factors that can create obstruction to speech while communicating. This may lead to degradation of speech processing and affect humans perceptually. Knowing whether a speaker wears a mask may be useful for modeling speech for different applications. With this motivation, finding whether a speaker wears face mask from a given speech is included as a task in Computational Paralinguistics Evaluation (ComParE) 2020. We study novel acoustic features based on linear filterbanks, instantaneous phase and long-term information that can capture the artifacts for classification of speech with and without face mask. These acoustic features are used along with the state-of-the-art baselines of ComParE functionals, bag-of-audio-words, DeepSpectrum and auDeep features for ComParE 2020. The studies reveal the effectiveness of acoustic features, and their score level fusion with the ComParE 2020 baselines leads to an unweighted average recall of 73.50% on the test set.


title: On generators of transition semigroups associated to semilinear stochastic partial differential equations

abstract: Let $\mathcal{X}$ be a real separable Hilbert space. Let $Q$ be a linear, self-adjoint, positive, trace class operator on $\mathcal{X}$, let $F:\mathcal{X}\rightarrow\mathcal{X}$ be a (smooth enough) function and let $\{W(t)\}_{t\geq 0}$ be a $\mathcal{X}$-valued cylindrical Wiener process. For $\alpha\in [0,1/2]$ we consider the operator $A:=-(1/2)Q^{2\alpha-1}:Q^{1-2\alpha}(\mathcal{X})\subseteq\mathcal{X}\rightarrow\mathcal{X}$. We are interested in the mild solution $X(t,x)$ of the semilinear stochastic partial differential equation \begin{gather} \left\{\begin{array}{ll} dX(t,x)=\big(AX(t,x)+F(X(t,x))\big)dt+ Q^{\alpha}dW(t), & t>0;\\ X(0,x)=x\in \mathcal{X}, \end{array} \right. \end{gather} and in its associated transition semigroup \begin{align} P(t)\varphi(x):=E[\varphi(X(t,x))], \qquad \varphi\in B_b(\mathcal{X}),\ t\geq 0,\ x\in \mathcal{X}; \end{align} where $B_b(\mathcal{X})$ is the space of the real-valued, bounded and Borel measurable functions on $\mathcal{X}$. In this paper we study the behavior of the semigroup $P(t)$ in the space $L^2(\mathcal{X},\nu)$, where $\nu$ is the unique invariant probability measure of \eqref{Tropical}, when $F$ is dissipative and has polynomial growth. Then we prove the logarithmic Sobolev and the Poincar\'e inequalities and we study the maximal Sobolev regularity for the stationary equation \[\lambda u-N_2 u=f,\qquad \lambda>0,\ f\in L^2(\mathcal{X},\nu);\] where $N_2$ is the infinitesimal generator of $P(t)$ in $L^2(\mathcal{X},\nu)$.


title: Emotion Invariant Speaker Embeddings for Speaker Identification with Emotional Speech

abstract: Emotional state of a speaker is found to have significant effect in speech production, which can deviate speech from that arising from neutral state. This makes identifying speakers with different emotions a challenging task as generally the speaker models are trained using neutral speech. In this work, we propose to overcome this problem by creation of emotion invariant speaker embedding. We learn an extractor network that maps the test embeddings with different emotions obtained using i-vector based system to an emotion invariant space. The resultant test embeddings thus become emotion invariant and thereby compensate the mismatch between various emotional states. The studies are conducted using four different emotion classes from IEMOCAP database. We obtain an absolute improvement of 2.6% in accuracy for speaker identification studies using emotion invariant speaker embedding against average speaker model based framework with different emotions.


title: An Experimental Analysis of Indoor Spatial Queries: Modeling, Indexing, and Processing

abstract: Indoor location-based services (LBS), such as POI search and routing, are often built on top of typical indoor spatial queries. To support such queries and indoor LBS, multiple techniques including model/indexes and search algorithms have been proposed. In this work, we conduct an extensive experimental study on existing proposals for indoor spatial queries. We survey five model/indexes, compare their algorithmic characteristics, and analyze their space and time complexities. We also design an in-depth benchmark with real and synthetic datasets, evaluation tasks and performance metrics. Enabled by the benchmark, we obtain and report the performance results of all model/indexes under investigation. By analyzing the results, we summarize the pros and cons of all techniques and suggest the best choice for typical scenarios.


title: Accretion and inter-cycle variations in the PMS interacting binary AK Sco

abstract: There are only a handful of known short-period pre-main sequence spectroscopic binaries with significant accretion rates (Class II sources). AK Sco stands out in this list because the system is composed of two equal mass F5 stars in a highly eccentric orbit thus both stars get as close as 11 stellar radii at periastron passage. This configuration is optimal for accretion studies because enhance accretion events can be precisely timed at periastron passage. In this work, we present the results from the monitoring of the AK Sco system with Hubble during three consecutive periastron passages. These data provide a unique dataset to spectroscopically characterize accretion and evaluate the inter-cycle variability of the system. Clear evidence of accretion rate enhancement was observed in cycle1 and 3: the blueing of the near UV continuum, the sudden flux increase of important accretion tracers, such as the N V, Si IV and C IV lines, and also of neutral/singly ionized species such as O I and C II. Also, variations in the Si III]/C III] ratio reveals an enhancement of the electron density by an order of magnitude during the periastron passage. Moreover, in cycle 3, the spectral resolution of the observations obtained with the Cosmic Origins Spectrograph enabled discerning that the flow was channelled preferentially into one of the two components. The most remarkable feature in the cycle-to-cycle variations was the detection of a notable increase of the UV flux from cycle 1 to 2 that was not accompanied by enhanced accretion signatures.


title: Top pair production at ultra-high energies

abstract: The top quark, the heaviest quark and, indeed, the heaviest elementary particle known today, constitutes a novel probe of the long-lived medium in quark-gluon phase which, as expected, can be produced even in light nuclei collisions at ultra-high energies. Some distinctive features are considered for particle production in the top sector in ultra-high energy domain. The antitop-top pair production is studied within the quantum chromodynamics and effective field theory approach used for calculations of total partonic cross sections. Predictions for all observables are computed at NNLO in quantum chromodynamics and at LO in effective field theory. These quantitative results can be important for both the future collider experiments at center-of-mass energy frontier and the improvement of the phenomenological models for development of the cosmic ray cascades in ultra-high energy domain. Thus the study allows the better understanding of heavy particle production and emphasizes the exciting interrelation between the high-energy physics on accelerators and ultra-high energy cosmic ray measurements.


title: Frames of group-sets and their application in bundle theory

abstract: We study fiber bundles where the fibers are not a group $G$, but a free $G$-space with disjoint orbits. These bundles closely resemble principal bundles, hence we call them semi-principal bundles. The study of such bundles is facilitated by defining the notion of a basis of a $G$-set, in analogy with a basis of a vector space. The symmetry group of these bases is a wreath product. Similar to vector bundles, using the notion of a basis induces a frame bundle construction, which in this case results in a principal bundle with the wreath product as structure group. This construction can be formalized in the language of a functor, which retracts the semi-principal bundles to the principal bundles. In addition, semi-principal bundles support parallel transport just like principal bundles, and this carries over to the frame bundle.


title: Graphical Calculi and their Conjecture Synthesis

abstract: Categorical Quantum Mechanics, and graphical calculi in particular, has proven to be an intuitive and powerful way to reason about quantum computing. This work continues the exploration of graphical calculi, inside and outside of the quantum computing setting, by investigating the algebraic structures with which we label diagrams. The initial aim for this was Conjecture Synthesis; the algorithmic process of creating theorems. To this process we introduce a generalisation step, which itself requires the ability to infer and then verify parameterised families of theorems. This thesis introduces such inference and verification frameworks, in doing so forging novel links between graphical calculi and fields such as Algebraic Geometry and Galois Theory. These frameworks inspired further research into the design of graphical calculi, and we introduce two important new calculi here. First is the calculus RING, which is initial among ring-based qubit graphical calculi, and in turn inspired the introduction and classification of phase homomorphism pairs also presented here. The second is the calculus ZQ, an edge-decorated calculus which naturally expresses arbitrary qubit rotations, eliminating the need for non-linear rules such as (EU) of ZX. It is expected that these results will be of use to those creating optimisation schemes and intermediate representations for quantum computing, to those creating new graphical calculi, and for those performing conjecture synthesis.


title: The effect of ion irradiation on dephasing of coherent optical phonons in GaP

abstract: The dephasing of coherent longitudinal optical (LO) phonons in ion-irradiated GaP has been investigated with a femtosecond pump-probe technique based on electro-optic sampling. The dephasing time of the coherent LO phonon is found to be dramatically prolonged by the introduction of a small amount of defects by means of Ga-ion irradiation. The maximum dephasing time observed at room temperature is 9.1 ps at a Ga$^{+}$ ion dose of 10$^{13}$/cm$^{2}$, which is significantly longer than the value of 8.3 ps for GaP before ion irradiation. The longer dephasing time is explained in terms of the suppression of electron-LO-phonon scattering by the presence of defect-induced deep levels.


title: Nuclear Pairing Gaps and Neutron Star Cooling

abstract: We study the cooling of isolated neutron stars with particular regard to the importance of nuclear pairing gaps. A microscopic nuclear equation of state derived in the Brueckner-Hartree-Fock approach is used together with compatible neutron and proton pairing gaps. We then study the effect of modifying the gaps on the final deduced neutron star mass distributions. We find that a consistent description of all current cooling data can be achieved and a reasonable neutron star mass distribution can be predicted employing the (slightly reduced by about 40\%) proton 1S0 Bardeen-Cooper-Schrieffer (BCS) gaps and no neutron 3P2 pairing.


title: AT-MFCGA: An Adaptive Transfer-guided Multifactorial Cellular Genetic Algorithm for Evolutionary Multitasking

abstract: Transfer Optimization is an incipient research area dedicated to solving multiple optimization tasks simultaneously. Among the different approaches that can address this problem effectively, Evolutionary Multitasking resorts to concepts from Evolutionary Computation to solve multiple problems within a single search process. In this paper we introduce a novel adaptive metaheuristic algorithm to deal with Evolutionary Multitasking environments coined as Adaptive Transfer-guided Multifactorial Cellular Genetic Algorithm (AT-MFCGA). AT-MFCGA relies on cellular automata to implement mechanisms in order to exchange knowledge among the optimization problems under consideration. Furthermore, our approach is able to explain by itself the synergies among tasks that were encountered and exploited during the search, which helps us to understand interactions between related optimization tasks. A comprehensive experimental setup is designed to assess and compare the performance of AT-MFCGA to that of other renowned evolutionary multitasking alternatives (MFEA and MFEA-II). Experiments comprise 11 multitasking scenarios composed of 20 instances of 4 combinatorial optimization problems, yielding the largest discrete multitasking environment solved to date. Results are conclusive in regard to the superior quality of solutions provided by AT-MFCGA with respect to the rest of the methods, which are complemented by a quantitative examination of the genetic transferability among tasks throughout the search process.


title: Near-Optimal Comparison Based Clustering

abstract: The goal of clustering is to group similar objects into meaningful partitions. This process is well understood when an explicit similarity measure between the objects is given. However, far less is known when this information is not readily available and, instead, one only observes ordinal comparisons such as "object i is more similar to j than to k." In this paper, we tackle this problem using a two-step procedure: we estimate a pairwise similarity matrix from the comparisons before using a clustering method based on semi-definite programming (SDP). We theoretically show that our approach can exactly recover a planted clustering using a near-optimal number of passive comparisons. We empirically validate our theoretical findings and demonstrate the good behaviour of our method on real data.


title: Elementary Integral Series for Heun Functions. With an Application to Black-Hole Perturbation Theory

abstract: Heun differential equations are the most general second order Fuchsian equations with four regular singularities. An explicit integral series representation of Heun functions involving only elementary integrands has hitherto been unknown and noted as an important open problem in a recent review. We provide explicit integral representations of the solutions of all equations of the Heun class: general, confluent, bi-confluent, doubly-confluent and triconfluent, with integrals involving only rational functions and exponential integrands. All the series are illustrated with concrete examples of use. These results stem from the technique of path-sums, which we use to evaluate the path-ordered exponential of a variable matrix chosen specifically to yield Heun functions. We demonstrate the utility of the integral series by providing the first representation of the solution to the Teukolsky radial equation governing the metric perturbations of rotating black holes that is convergent everywhere from the black hole horizon up to spatial infinity.


title: Predicting Typological Features in WALS using Language Embeddings and Conditional Probabilities: \'UFAL Submission to the SIGTYP 2020 Shared Task

abstract: We present our submission to the SIGTYP 2020 Shared Task on the prediction of typological features. We submit a constrained system, predicting typological features only based on the WALS database. We investigate two approaches. The simpler of the two is a system based on estimating correlation of feature values within languages by computing conditional probabilities and mutual information. The second approach is to train a neural predictor operating on precomputed language embeddings based on WALS features. Our submitted system combines the two approaches based on their self-estimated confidence scores. We reach the accuracy of 70.7% on the test data and rank first in the shared task.


title: Is the atomic metal vapor a dielectric state?

abstract: We propose a simple method for calculating the metal vapor conductivity at the critical point and near-critical isotherms. This method's base is the hypothesis of an electron jellium's existence as an origin of the conduction band in metal vapor's gaseous phase. Satisfactory agreement with the experimental data for alkali metals (Cs, Rb) concludes that the "cold ionization" mechanism is possible in the critical point's vicinity. The liquid metal-dielectric transition, as previously thought, does not occur; but the liquid metal-gaseous metal transition occurs probably. The gaseous metal is not a dielectric state of the metal.


title: Minimization of mean-CVaR evacuation time of a crowd using rescue guides: a scenario-based approach

abstract: In case of a threat in a public space, the crowd in it should be moved to a shelter or evacuated without delays. Risk management and evacuation planning in public spaces should also take into account uncertainties in the traffic patterns of crowd flow. One way to account for the uncertainties is to make use of safety staff, or guides, that lead the crowd out of the building according to an evacuation plan. Nevertheless, solving the minimum time evacuation plan is a computationally demanding problem. In this paper, we model the evacuating crowd and guides as a multi-agent system with the social force model. To represent uncertainty, we construct probabilistic scenarios. The evacuation plan should work well both on average and also for the worst-performing scenarios. Thus, we formulate the problem as a bi-objective scenario optimization problem, where the mean and conditional value-at-risk (CVaR) of the evacuation time are objectives. A solution procedure combining numerical simulation and genetic algorithm is presented. We apply it to the evacuation of a fictional passenger terminal. In the mean-optimal solution, guides are assigned to lead the crowd to the nearest exits, whereas in the CVaR-optimal solution the focus is on solving the physical congestion occurring in the worst-case scenario. With one guide positioned behind each agent group near each exit, a plan that minimizes both objectives is obtained.


title: VECMAtk: A Scalable Verification, Validation and Uncertainty Quantification Toolkit for Scientific Simulations

abstract: We present the VECMA toolkit (VECMAtk), a flexible software environment for single and multiscale simulations that introduces directly applicable and reusable procedures for verification, validation (V&V), sensitivity analysis (SA) and uncertainty quantification (UQ). It enables users to verify key aspects of their applications, systematically compare and validate the simulation outputs against observational or benchmark data, and run simulations conveniently on any platform from the desktop to current multi-petascale computers. In this sequel to our paper on VECMAtk which we presented last year, we focus on a range of functional and performance improvements that we have introduced, cover newly introduced components, and applications examples from seven different domains such as conflict modelling and environmental sciences. We also present several implemented patterns for UQ/SA and V&V, and guide the reader through one example concerning COVID-19 modelling in detail.


title: Recycled pulsars with multipolar magnetospheres from accretion-induced magnetic burial

abstract: Many millisecond pulsars are thought to be old neutron stars spun up (`recycled') during an earlier accretion phase. They typically have relatively weak ($\lesssim 10^{9} \text{ G}$) dipole field strengths, consistent with accretion-induced magnetic burial. Recent data from the Neutron Star Interior Composition Explorer indicate that hot spots atop the recycled pulsar PSR J0030$+$0451 are not antipodal, so that the magnetic field cannot be that of a centered dipole. In this paper it is shown that multipolarity is naturally expected in the burial scenario because of equatorial field line compression. Grad-Shafranov equilibria are constructed to show how magnetic multipole moments can be calculated in terms of various properties, such as the amount of accreted mass and the crustal equation of state.


title: Active Dipolar Colloids in Three Dimensions: Strings, Sheets, Labyrinthine Textures and Crystals

abstract: Active matter exhibits striking behaviour reminiscent of living matter and molecular fluids, and has promising applications in drug delivery or mixing at the micron scale. Active colloidal systems provide important models with simple and controllable interactions amenable to theory and computer simulation. Experimental work is dominated by (quasi) two--dimensional (2d) systems at relatively low concentration, and rather less is known of the 3d case at concentrations pertinent to motility induced phase separation or to mimic morphogenesis. Here we investigate a 3d experimental system of active colloids in a dense suspension up to volume fractions of 0.5. The particles in our system are self-propelled in the lateral plane under an AC electric field. The field in addition induces an electric dipole, and the competition between activity and both steric and dipolar interactions gives rise to phase behaviour ranging from an active gas to a dynamic labyrinthine phase as well as dense tetragonal and hexagonal crystals. Intermediate volume fractions are characterised by two--dimensional sheets with large fluctuations reminiscent of active membranes. These active sheets break symmetry in a direction perpendicular to the applied field. Moreover, the relationship between electric field and the particle dynamics depends in a complex and unexpected manner upon the position in the state diagram.

categories: cond-mat.soft cond-mat.mtrl-sci cond-mat.stat-mech physics.chem-ph


title: Charge and hydration structure of dendritic polyelectrolytes: molecular simulations of polyglycerol sulphate

abstract: Macromolecules based on dendritic or hyperbranched polyelectrolytes have been emerging as high potential candidates for biomedical applications. Here we study the charge and solvation structure of dendritic polyglycerol sulphate (dPGS) of generations 0 to 3 in aqueous sodium chloride solution by explicit-solvent molecular dynamics computer simulations. We characterize dPGS by calculating several important properties such as relevant dPGS radii, molecular distributions, the solvent accessible surface area, and the partial molecular volume. In particular, as the dPGS exhibits high charge renormalization effects, we address the challenges of how to obtain a well-defined effective charge and surface potential of dPGS for practical applications. We compare implicit- and explicit-solvent approaches in our all-atom simulations with the coarse-grained simulations from our previous work. We find consistent values for the effective electrostatic size (i.e., the location of the effective charge of a Debye--H\"{u}ckel sphere) within all the approaches, deviating at most by the size of a water molecule. Finally, the excess chemical potential of water insertion into dPGS and its thermodynamic signature are presented and rationalized.


title: Rough-Surface Shadowing of Self-Affine Random Rough Surfaces

abstract: Light scattering from self-affine homogeneous isotropic random rough surfaces is studied using the ray-optics approximation. Numerical methods are developed to accelerate the first-order scattering simulations from surfaces represented as single-connected single-valued random fields, and to store the results of the simulations into a numerical reflectance model. Horizon mapping and marching methods are developed to accelerate the simulation. Emphasis is given to the geometric shadowing and masking effects as a function of surface roughness, especially, to the azimuthal rough-surface shadowing effect.


title: Phenomenological constraints on the transport properties of QCD matter with data-driven model averaging

abstract: Using combined data from the Relativistic Heavy Ion and Large Hadron Colliders, we constrain the shear and bulk viscosities of quark-gluon plasma (QGP) at temperatures of ${\sim\,}150{-}350$ MeV. We use Bayesian inference to translate experimental and theoretical uncertainties into probabilistic constraints for the viscosities. With Bayesian Model Averaging we account for the irreducible model ambiguities in the transition from a fluid description of the QGP to hadronic transport in the final evolution stage, providing the most reliable phenomenological constraints to date on the QGP viscosities.


title: Response theory for nonequilibrium steady-states of open quantum systems

abstract: We introduce a response theory for open quantum systems within nonequilibrium steady-states subject to a Hamiltonian perturbation. Working in the weak system-bath coupling regime, our results are derived within the Lindblad-Gorini-Kossakowski-Sudarshan formalism. We find that the response of the system to a small perturbation is not simply related to a correlation function within the system, unlike traditional linear response theory in closed systems or expectations from the fluctuation-dissipation theorem. In limiting cases, when the perturbation is small relative to the coupling to the surroundings or when it does not lead to a change of the eigenstructure of the system, a perturbative expansion exists where the response function is related to a sum of a system correlation functions and additional forces induced by the surroundings. Away from these limiting regimes however, the secular approximation results in a singular response that cannot be captured within the traditional approach, but can be described by reverting to a microscopic Hamiltonian description. These findings are illustrated by explicit calculations in coupled qubits and anharmonic oscillators in contact with bosonic baths at different temperatures.


title: Solving non-linear integral equations for laser pulse retrieval with Newton's method

abstract: We present an algorithm based on numerical techniques that have become standard for solving nonlinear integral equations: Newton's method, homotopy continuation, the multilevel method and random projection to solve the inversion problem that appears when retrieving the electric field of an ultrashort laser pulse from a 2-dimensional intensity map measured with Frequency-resolved optical gating (FROG), dispersion-scan or amplitude-swing experiments. Here we apply the solver to FROG and specify the necessary modifications for similar integrals. Unlike other approaches we transform the integral and work in time-domain where the integral can be discretised as an over-determined polynomial system and evaluated through list auto-correlations. The solution curve is partially continues and partially stochastic, consisting of small linked path segments and enables the computation of optimal solutions in the presents of noise. Interestingly, this is a novel method to find real solutions of polynomial systems which are notoriously difficult to find. We show how to implement adaptive Tikhonov-type regularization to smooth the solution when dealing with noisy data, we compare the results for noisy test data with a least-squares solver and propose the L-curve method to fine-tune the regularization parameter.


title: A symmetry-preserving second-order time-accurate PISO-based method

abstract: A new conservative symmetry-preserving second-order time-accurate PISO-based pressure-velocity coupling for solving the incompressible Navier-Stokes equations on unstructured collocated grids is presented in this paper. This new method for implicit time stepping is an extension of the conservative symmetry-preserving incremental-pressure projection method for explicit time stepping and unstructured collocated meshes of Trias et al. (2014). In order to assess and compare both methods, we have implemented them within one unified solver in the open source code OpenFOAM. We combine both methods with a Butcher tableau for a family of explicit and diagonally implicit Runge-Kutta temporal schemes. We assess the energy conservation properties of the implemented discretisation methods and the temporal consistency of the selected Runge-Kutta schemes using Taylor-Green vortex and lid-driven cavity flow test cases. Although both implemented methods are based on a symmetry-preserving discretisation, we show that both methods still produce a small amount of numerical dissipation when the total pressure is directly solved from a Poisson equation. This numerical dissipation is mainly caused by the corresponding pressure error which is of $O(\Delta t \Delta h^2)$. When an incremental-pressure approach is used, where a pressure correction is solved from a Poisson equation, the pressure error reduces to $O(\Delta t^2 \Delta h^2)$, yielding better conservation properties: both methods are then effectively fully-conservative. Furthermore, we conclude that all selected explicit and implicit higher order temporal schemes suffer from a reduction of the temporal order to approximately one when the pressure Poisson equation is based on the total pressure due to the presence of a pressure error of $O(\Delta t \Delta h^2)$.


title: Fraunhofer oscillations of the critical current at a varying Zeeman field in a spin-orbit coupled Josephson junction

abstract: The Zeeman interaction results in spontaneous current through a Josephson contact with a spin-orbit coupled normal metal, even in the absence of any voltage, or phase bias. In the case of the Rashba spin orbit coupling of electrons in a two-dimensional (2D) electron gas this effect takes place for the Zeeman field which is parallel to the 2D system and to superconducting contacts. At the same time, the spontaneous current is absent when this field is perpendicular to the contacts. It is shown that in the latter case it may manifest itself in oscillations of the critical Josephson current at the varying Zeeman energy. These oscillations have a form of the Fraunhofer diffraction pattern. The Josephson current under the phase bias was calculated based on the semiclassical Green functions for a disordered 2D electron gas with the strong spin orbit coupling, as well as for surface electrons of a three dimensional topological insulator. In the latter case the diffraction pattern was found to be most pronounced, while in the Rashba gas the oscillations of the critical current are weaker.


title: Assessing Classifier Fairness with Collider Bias

abstract: The increasing application of machine learning techniques in everyday decision-making processes has brought concerns about the fairness of algorithmic decision-making. This paper concerns the problem of collider bias which produces spurious associations in fairness assessment and develops theorems to guide fairness assessment avoiding the collider bias. We consider a real-world application of auditing a trained classifier by an audit agency. We propose an unbiased assessment algorithm by utilising the developed theorems to reduce collider biases in the assessment. Experiments and simulations show the proposed algorithm reduces collider biases significantly in the assessment and is promising in auditing trained classifiers.


title: Prioritized Level Replay

abstract: Environments with procedurally generated content serve as important benchmarks for testing systematic generalization in deep reinforcement learning. In this setting, each level is an algorithmically created environment instance with a unique configuration of its factors of variation. Training on a prespecified subset of levels allows for testing generalization to unseen levels. What can be learned from a level depends on the current policy, yet prior work defaults to uniform sampling of training levels independently of the policy. We introduce Prioritized Level Replay (PLR), a general framework for selectively sampling the next training level by prioritizing those with higher estimated learning potential when revisited in the future. We show TD-errors effectively estimate a level's future learning potential and, when used to guide the sampling procedure, induce an emergent curriculum of increasingly difficult levels. By adapting the sampling of training levels, PLR significantly improves sample efficiency and generalization on Procgen Benchmark--matching the previous state-of-the-art in test return--and readily combines with other methods. Combined with the previous leading method, PLR raises the state-of-the-art to over 76% improvement in test return relative to standard RL baselines.


title: Extending C++ for Heterogeneous Quantum-Classical Computing

abstract: We present qcor - a language extension to C++ and compiler implementation that enables heterogeneous quantum-classical programming, compilation, and execution in a single-source context. Our work provides a first-of-its-kind C++ compiler enabling high-level quantum kernel (function) expression in a quantum-language agnostic manner, as well as a hardware-agnostic, retargetable compiler workflow targeting a number of physical and virtual quantum computing backends. qcor leverages novel Clang plugin interfaces and builds upon the XACC system-level quantum programming framework to provide a state-of-the-art integration mechanism for quantum-classical compilation that leverages the best from the community at-large. qcor translates quantum kernels ultimately to the XACC intermediate representation, and provides user-extensible hooks for quantum compilation routines like circuit optimization, analysis, and placement. This work details the overall architecture and compiler workflow for qcor, and provides a number of illuminating programming examples demonstrating its utility for near-term variational tasks, quantum algorithm expression, and feed-forward error correction schemes.


title: Cinema Darkroom: A Deferred Rendering Framework for Large-Scale Datasets

abstract: This paper presents a framework that fully leverages the advantages of a deferred rendering approach for the interactive visualization of large-scale datasets. Geometry buffers (G-Buffers) are generated and stored in situ, and shading is performed post hoc in an interactive image-based rendering front end. This decoupled framework has two major advantages. First, the G-Buffers only need to be computed and stored once---which corresponds to the most expensive part of the rendering pipeline. Second, the stored G-Buffers can later be consumed in an image-based rendering front end that enables users to interactively adjust various visualization parameters---such as the applied color map or the strength of ambient occlusion---where suitable choices are often not known a priori. This paper demonstrates the use of Cinema Darkroom on several real-world datasets, highlighting CD's ability to effectively decouple the complexity and size of the dataset from its visualization.


title: Direct reaction contribution in the very low energy $^{19}$F(p,$\alpha_0$) reaction

abstract: The direct reaction component of the $^{19}$F(p,$\alpha_{0}$) reaction at E$_{cm}$=180-600 keV is studied for the data that has become very recently available. This component has been found to be significant in the present work using the direct pickup model in framework of the DWBA formalism and indicate the strong cluster structure of $^{19}$F.


title: Self-consistent ladder D$\Gamma$A approach

abstract: We present and implement a self-consistent D$\Gamma$A approach for multi-orbital models and ab initio materials calculations. It is applied to the one-band Hubbard model at various interaction strengths with and without doping, to the two-band Hubbard model with two largely different bandwidths, and to SrVO$_3$. The self-energy feedback reduces critical temperatures compared to dynamical mean-field theory, even to zero temperature in two-dimensions. Compared to a one-shot, non-self-consistent calculation the non-local correlations are significantly reduced when they are strong. In case non-local correlations are weak to moderate as for SrVO$_3$, one-shot calculations are sufficient.


title: Mixing for the primitive equations under bounded non-degenerate noise

abstract: We study the stochastic 3D primitive equations of the atmospheric mechanics. We consider them under a bounded and non-degenerate noise, which is statistically periodic in time with period $1$. In such a case we prove that the associated integer-time Markov chain is mixing, which means that there exists a unique stationary measure to which converge the laws all trajectories of this Markov chain.


title: Phase coexistence and negative thermal expansion in the triple perovskite iridate Ba$_{3}$CoIr$_{2}$O$_{9}$

abstract: The anomalous thermal expansion in a layered 3$d$-5$d$ based triple perovskite iridate Ba$_{3}$CoIr$_{2}$O$_{9}$ is investigated using high resolution synchrotron diffraction. Below the magneto-structural transition at 107\,K, the onset of antiferromagnetic order is associated with a monoclinic distortion of the hexagonal structure. Deeper within the magnetically ordered state, a part of the monoclinic phase distorts even further, and both these structural phases co-exist down to the lowest measured temperatures. We observe negative thermal expansion in this phase co-existence regime, which appears to be intimately connected to the temperature driven relative fractions of these monoclinic phases. The significant NTE observed in this system could be driven by magnetic exchange striction, and is of relevance to a number of systems with pronounced spin orbit interactions.


title: Scattering from rough thin films: DDA-simulations

abstract: We investigate the wave-optical light scattering properties of deformed thin circular films of constant thickness using the discrete-dipole approximation. Effects on the intensity distribution of the scattered light due to different statistical roughness models, model dependant roughness parameters, and uncorrelated random small-scale porosity of the inhomogeneous medium are studied. The usability of discrete-dipole approximation to rough-surface scattering problems is evaluated by considering thin films as computationally feasible rough-surface analogs. The effects due to small-scale inhomogeneity of the scattering medium are compared with the analytic approximation by Maxwell Garnett and the results are found to agree with the approximation.


title: Excited-state dynamics of structurally characterized crystal of SnxSb1-x

abstract: The topological behavior of heavy metal alloys opens a vast area for incredible research and future technology. Here, we extend our previous report about the superconducting properties of Sn0.4Sb0.6 along with the compositional variation of Sn and Sb in SnxSb1-x (with (X=0.5 and 0.6)) to study the detailed optical properties. Structural and morphological details of grown crystal are carried from the previous study. Further, the samples are excited by a pump of 2.61 eV with a broad probe of 0.77-1.54 eV in the NIR regime for transient reflectance ultrafast studies (TRUS) measurements. The differential reflectance profile shows an unprecedented negative magnitude, and the average power-dependent analysis of this negative trend has been analyzed. This article not only provides evidence of band filling phenomenon in the samples but also shows that with the variation of average power, there is a definite increase in the excited charge carriers, and thereby enhancing the band filling response. The estimated value of the bandgap between the band filled states and valence state is also determined from these studies. The nonlinear properties and bandgap analysis of the studied topological alloys and similar materials help in the advancement of various nonlinear optical applications.


title: More on $\mathcal{N} =2$ S-folds

abstract: We carry out a systematic study of 4d $\mathcal{N} = 2$ preserving S-folds of F-theory 7-branes and the worldvolume theories on D3-branes probing them. They consist of two infinite series of theories, which we denote following the original papers by $\mathcal{S}^{(r)}_{G,\ell}$ for $\ell = 2,3,4$ and $\mathcal{T}^{(r)}_{G,\ell}$ for $\ell = 2,3,4,5,6$. Their distinction lies in the discrete torsion carried by the S-fold and in the difference in the asymptotic holonomy of the gauge bundle on the 7-brane. We study various properties of these theories, using diverse field theoretical and string theoretical methods.


title: The soluble radical and orbits of certain maps on finite groups

abstract: For each element $u$ in a finite group $G$ define a map $\theta_u\colon G\to G$ by $\theta_u(g)=[g^{-u},g]$ and set $\Theta_G(u)=\{g\in G\mid \theta_u^n(g)=g \hbox{ for some } n>0\}$. Then $\theta_u$ induces a permutation of $\Theta_G(u)$; let $\beta_G(u)$ be the number of orbits apart from $\{1\}$. Building on work of J.N. Bray, R.A. Wilson and the second author, we show that the index of the soluble radical of a finite group $G$ is bounded in terms of the values of $\beta_G(u)$ for $2$-elements $u$.


title: Interplay between coherent and incoherent decay processes in chaotic systems: the role of quantum interference

abstract: The population decay due to a small opening in an otherwise closed cavity supporting chaotic classical dynamics displays a quantum correction on top of the classical exponential form, a pure manifestation of quantum coherence that acquires a universal form and can be explained by path interference. Being coherent, such enhancement is prone to decoherence effects due to the coupling of the system to an external environment. We study this interplay between incoherent and coherent quantum corrections to decay by evaluating, within a Caldeira-Legget scenario, off-diagonal contributions to the decoherence functional coming from pairs of correlated classical paths in the time regime where dissipative effects are neglected and decoherence does not affect the classical dynamics, but quantum interference must be accounted for. We find that the competing effects of interference and decoherence lead to a universal non-monotonous form for the survival probability depending only on the universality class, coupling strength, and macroscopic parameters of the cavity.


title: Exponential corrected thermodynamics of black holes

abstract: Recently, it is reported that thermal fluctuations, which are interpreted as quantum effects, modify the black holes entropy by an exponential term. We now find the effect of such modification on the black hole mass and other thermodynamics quantities. We find that the Schwarzschild black hole mass decreased by thermal fluctuations. Hence, we study exponential corrected thermodynamics and statistics of black holes by computing the partition function. We obtain the special condition on the event horizon radius to satisfy Smarr-Gibbs-Duhem relation in the presence of quantum correction. As we know, the Schwarzschild black hole is unstable, while the effect of exponential correction is the stability of $4D$ Schwarzschild black hole as well as the Schwarzschild-AdS black hole at a small area. On the other hand, a $5D$ Schwarzschild black hole is completely unstable. The effect of the quantum correction on the Reissner-Nordstr\"{o}m black hole is instability at quantum scales. Finally, we consider the most general case of charged AdS black hole and study the corrected thermodynamics.


title: How to form a wormhole

abstract: We provide a simple but very useful description of the process of wormhole formation. We place two massive objects in two parallel universes (modeled by two branes). Gravitational attraction between the objects competes with the resistance coming from the brane tension. For sufficiently strong attraction, the branes are deformed, objects touch and a wormhole is formed. Our calculations show that more massive and compact objects are more likely to fulfill the conditions for wormhole formation. This implies that we should be looking for wormholes either in the background of black holes and compact stars, or massive microscopic relics. Our formation mechanism applies equally well for a wormhole connecting two objects in the same universe.


title: Artificial intelligence supported anemia control system (AISACS) to prevent anemia in maintenance hemodialysis patients

abstract: Anemia, for which erythropoiesis-stimulating agents (ESAs) and iron supplements (ISs) are used as preventive measures, presents important difficulties for hemodialysis patients. Nevertheless, the number of physicians able to manage such medications appropriately is not keeping pace with the rapid increase of hemodialysis patients. Moreover, the high cost of ESAs imposes heavy burdens on medical insurance systems. An artificial-intelligence-supported anemia control system (AISACS) trained using administration direction data from experienced physicians has been developed by the authors. For the system, appropriate data selection and rectification techniques play important roles. Decision making related to ESAs poses a multi-class classification problem for which a two-step classification technique is introduced. Several validations have demonstrated that AISACS exhibits high performance with correct classification rates of 72-87% and clinically appropriate classification rates of 92-98%.


title: An Optimum Algorithm for Quantum Search

abstract: This paper discusses an improvement to Grover's algorithm for searches where target states are Hamming weight eigenstates and search space is not ordered. It is shown that under these conditions search efficiency depends on the smaller number of 0's and 1's, not the total length, of binary string of target state, and that Grover's algorithm can be improved whenever number of 0's and number 1's are not equal. In particular, improvement can be exponential when number of 0's or number of 1's is very small relative to binary string length. One interesting application is that Dicke state preparation, which in Grover's algorithm is P on average, can be made poly-efficient in all cases. For decision making process, this improvement won't improve computation efficiency, but can make implementation much simpler.


title: Reward Machines: Exploiting Reward Function Structure in Reinforcement Learning

abstract: Reinforcement learning (RL) methods usually treat reward functions as black boxes. As such, these methods must extensively interact with the environment in order to discover rewards and optimal policies. In most RL applications, however, users have to program the reward function and, hence, there is the opportunity to make the reward function visible -- to show the reward function's code to the RL agent so it can exploit the function's internal structure to learn optimal policies in a more sample efficient manner. In this paper, we show how to accomplish this idea in two steps. First, we propose reward machines, a type of finite state machine that supports the specification of reward functions while exposing reward function structure. We then describe different methodologies to exploit this structure to support learning, including automated reward shaping, task decomposition, and counterfactual reasoning with off-policy learning. Experiments on tabular and continuous domains, across different tasks and RL agents, show the benefits of exploiting reward structure with respect to sample efficiency and the quality of resultant policies. Finally, by virtue of being a form of finite state machine, reward machines have the expressive power of a regular language and as such support loops, sequences and conditionals, as well as the expression of temporally extended properties typical of linear temporal logic and non-Markovian reward specification.


title: MolDesigner: Interactive Design of Efficacious Drugs with Deep Learning

abstract: The efficacy of a drug depends on its binding affinity to the therapeutic target and pharmacokinetics. Deep learning (DL) has demonstrated remarkable progress in predicting drug efficacy. We develop MolDesigner, a human-in-the-loop web user-interface (UI), to assist drug developers leverage DL predictions to design more effective drugs. A developer can draw a drug molecule in the interface. In the backend, more than 17 state-of-the-art DL models generate predictions on important indices that are crucial for a drug's efficacy. Based on these predictions, drug developers can edit the drug molecule and reiterate until satisfaction. MolDesigner can make predictions in real-time with a latency of less than a second.