title: Duality between range and no-response tests and its application for inverse problems

abstract: In this paper we will show the duality between the range test (RT) and no-response test (NRT) for the inverse boundary value problem for the Laplace equation in $\Omega\setminus\overline D$ with an obstacle $D\Subset\Omega$ whose boundary $\partial D$ is visible from the boundary $\partial\Omega$ of $\Omega$ and a measurement is given as a set of Cauchy data on $\partial\Omega$. Here the Cauchy data is given by a unique solution $u$ of the boundary value problem for the Laplace equation in $\Omega\setminus\overline D$ with homogeneous and inhomogeneous Dirichlet boundary condition on $\partial D$ and $\partial\Omega$, respectively. These testing methods are domain sampling method to estimate the location of the obstacle using test domains and the associated indicator functions. Also both of these testing methods can test the analytic extension of $u$ to the exterior of a test domain. Since these methods are defined via some operators which are dual to each other, we could expect that there is a duality between the two methods. We will give this duality in terms of the equivalence of the pre-indicator functions associated to their indicator functions. As an application of the duality, the reconstruction of $D$ using the RT gives the reconstruction of $D$ using the NRT and vice versa. We will also give each of these reconstructions without using the duality if the Dirichlet data of the Cauchy data on $\partial\Omega$ is not identically zero and the solution to the associated forward problem does not have any analytic extension across $\partial D$. Moreover, we will show that these methods can still give the reconstruction of $D$ if $D$ is a convex polygon and it satisfies one of the following two properties: all of its corner angles are irrational and its diameter is less than its distance to $\partial\Omega$.


title: Non-Blocking Simultaneous Multithreading: Embracing the Resiliency of Deep Neural Networks

abstract: Deep neural networks (DNNs) are known for their inability to utilize underlying hardware resources due to hardware susceptibility to sparse activations and weights. Even in finer granularities, many of the non-zero values hold a portion of zero-valued bits that may cause inefficiencies when executed on hardware. Inspired by conventional CPU simultaneous multithreading (SMT) that increases computer resource utilization by sharing them across several threads, we propose non-blocking SMT (NB-SMT) designated for DNN accelerators. Like conventional SMT, NB-SMT shares hardware resources among several execution flows. Yet, unlike SMT, NB-SMT is non-blocking, as it handles structural hazards by exploiting the algorithmic resiliency of DNNs. Instead of opportunistically dispatching instructions while they wait in a reservation station for available hardware, NB-SMT temporarily reduces the computation precision to accommodate all threads at once, enabling a non-blocking operation. We demonstrate NB-SMT applicability using SySMT, an NB-SMT-enabled output-stationary systolic array (OS-SA). Compared with a conventional OS-SA, a 2-threaded SySMT consumes 1.4x the area and delivers 2x speedup with 33% energy savings and less than 1% accuracy degradation of state-of-the-art CNNs with ImageNet. A 4-threaded SySMT consumes 2.5x the area and delivers, for example, 3.4x speedup and 39% energy savings with 1% accuracy degradation of 40%-pruned ResNet-18.


title: Schiffer variations and the generic Torelli theorem for hypersurfaces

abstract: We show how to recover a general hypersurface in $\mathbb{P}^n$ of sufficiently large degree $d$ dividing $n+1$, from its finite order variation of Hodge structure. We also analyze the two other series of cases not covered by Donagi's generic Torelli theorem. Combined with Donagi's theorem, this shows that the generic Torelli theorem for hypersurfaces holds with finitely many exceptions.


title: Dynamic length scales in athermal, shear-driven, jamming of frictionless disks in two dimensions

abstract: We carry our numerical simulations of athermally sheared, bidisperse, frictionless disks in two dimensions. From an appropriately defined velocity correlation function, we determine that there are two diverging length scales, $\xi$ and $\ell$, as the jamming transition is approached. We analyze our results using a critical scaling ansatz for the correlation function, and argue that the more divergent length $\ell$ is a consequence of a dangerous irrelevant scaling variable, and that it is $\xi$ which is the correlation length that determines the divergence of the system viscosity as jamming is approached from below in the liquid phase. We find that $\xi\sim (\phi_J-\phi)^{-\nu}$ diverges with the critical exponent $\nu=1$. We provide evidence that $\xi$ measures the length scale of fluctuations in the rotation of the particle velocity field, while $\ell$ measures the length scale of fluctuations in the divergence of the velocity field.


title: A Novel Multi-Agent System for Complex Scheduling Problems

abstract: Complex scheduling problems require a large amount computation power and innovative solution methods. The objective of this paper is the conception and implementation of a multi-agent system that is applicable in various problem domains. Independent specialized agents handle small tasks, to reach a superordinate target. Effective coordination is therefore required to achieve productive cooperation. Role models and distributed artificial intelligence are employed to tackle the resulting challenges. We simulate a NP-hard scheduling problem to demonstrate the validity of our approach. In addition to the general agent based framework we propose new simulation-based optimization heuristics to given scheduling problems. Two of the described optimization algorithms are implemented using agents. This paper highlights the advantages of the agent-based approach, like the reduction in layout complexity, improved control of complicated systems, and extendability.


title: Efficient, arbitrarily high precision hardware logarithmic arithmetic for linear algebra

abstract: The logarithmic number system (LNS) is arguably not broadly used due to exponential circuit overheads for summation tables relative to arithmetic precision. Methods to reduce this overhead have been proposed, yet still yield designs with high chip area and power requirements. Use remains limited to lower precision or high multiply/add ratio cases, while much of linear algebra (near 1:1 multiply/add ratio) does not qualify. We present a dual-base approximate logarithmic arithmetic comparable to floating point in use, yet unlike LNS it is easily fully pipelined, extendable to arbitrary precision with $O(n^2)$ overhead, and energy efficient at a 1:1 multiply/add ratio. Compared to float32 or float64 vector inner product with FMA, our design is respectively 2.3x and 4.6x more energy efficient in 7 nm CMOS. It depends on exp and log evaluation 5.4x and 3.2x more energy efficient, at 0.23x and 0.37x the chip area for equivalent accuracy versus standard hyperbolic CORDIC using shift-and-add and approximated ODE integration in the style of Revol and Yakoubsohn. This technique is a novel design alternative for low power, high precision hardened linear algebra in computer vision, graphics and machine learning applications.


title: Confinement strategies in a simple SIR model

abstract: We propose a simple SIR model in order to investigate the impact of various confinement strategies on a most virulent epidemic. Our approach is motivated by the current COVID-19 pandemic. The main hypothesis is the existence of two populations of susceptible persons, one which obeys confinement and for which the infection rate does not exceed 1, and a population which, being non confined for various imperatives, can be substantially more infective. The model, initially formulated as a differential system, is discretised following a specific procedure, the discrete system serving as an integrator for the differential one. Our model is calibrated so as to correspond to what is observed in the COVID-19 epidemic. Several conclusions can be reached, despite the very simple structure of our model. First, it is not possible to pinpoint the genesis of the epidemic by just analysing data from when the epidemic is in full swing. It may well turn out that the epidemic has reached a sizeable part of the world months before it became noticeable. Concerning the confinement scenarios, a universal feature of all our simulations is that relaxing the lockdown constraints leads to a rekindling of the epidemic. Thus we sought the conditions for the second epidemic peak to be lower than the first one. This is possible in all the scenarios considered (abrupt, progressive or stepwise exit) but typically a progressive exit can start earlier than an abrupt one. However, by the time the progressive exit is complete, the overall confinement times are not too different. From our results, the most promising strategy is that of a stepwise exit. And in fact its implementation could be quite feasible, with the major part of the population (minus the fragile groups) exiting simultaneously but obeying rigorous distancing constraints.


title: Large deviations for a class of tempered subordinators and their inverse processes

abstract: We consider a class of tempered subordinators, namely a class of subordinators with one-dimensional marginal tempered distributions which belong to a family studied in [3]. The main contribution in this paper is a non-central moderate deviations result. More precisely we mean a class of large deviation principles that fill the gap between the (trivial) weak convergence of some non-Gaussian identically distributed random variables to their common law, and the convergence of some other related random variables to a constant. Some other minor results concern large deviations for the inverse of the tempered subordinators considered in this paper; actually, in some results, these inverse processes appear as random time-changes of other independent processes.


title: Degree-targeted cascades in modular, degree-heterogeneous networks

abstract: The dynamics of cascading activation, such as rapid changes in public opinion and the outbreak of disease epidemics, have a crucial dependence on the connectivity patterns among the agents. We study cascading dynamics in modular, degree-heterogeneous networks, and consider the impact of intra-module seeding strategy on inter-module spread. Specifically, we establish that although activating the highest-degree nodes is more effective than random selection at growing a cascade locally, there is a critical level of inter-module connectivity required for a cascade to cross from one module to another, irrespective of the seeding strategy. We present an analytical proof of this statement for the case that each module has the same degree distribution and all module pairs have the same inter-module connectivity, while our simulation results suggest its validity for more general situations, including a ring of modules. Interestingly, we find that on a network comprised of two modules, this critical level is primarily determined by the degree distribution of the \emph{alter} module, as opposed to the seed module. Our analytical approach extends a method developed by Gleeson, but is able to capture different seeding strategies using only one dynamical variable per module, namely the conditional exposure probability. Our work shows that the possibility of a global cascade depends sensitively on inter-module connectivity, and less on the intra-module seeding strategy. This suggests, for example, that slight changes to inter-module connectivity can be a feasible intervention strategy to promote or inhibit global cascades.


title: How Do Neural Networks Estimate Optical Flow? A Neuropsychology-Inspired Study

abstract: End-to-end trained convolutional neural networks have led to a breakthrough in optical flow estimation. The most recent advances focus on improving the optical flow estimation by improving the architecture and setting a new benchmark on the publicly available MPI-Sintel dataset. Instead, in this article, we investigate how deep neural networks estimate optical flow. A better understanding of how these networks function is important for (i) assessing their generalization capabilities to unseen inputs, and (ii) suggesting changes to improve their performance. For our investigation, we focus on FlowNetS, as it is the prototype of an encoder-decoder neural network for optical flow estimation. Furthermore, we use a filter identification method that has played a major role in uncovering the motion filters present in animal brains in neuropsychological research. The method shows that the filters in the deepest layer of FlowNetS are sensitive to a variety of motion patterns. Not only do we find translation filters, as demonstrated in animal brains, but thanks to the easier measurements in artificial neural networks, we even unveil dilation, rotation, and occlusion filters. Furthermore, we find similarities in the refinement part of the network and the perceptual filling-in process which occurs in the mammal primary visual cortex.


title: Non-linear interlinkages and key objectives amongst the Paris Agreement and the Sustainable Development Goals

abstract: The United Nations' ambitions to combat climate change and prosper human development are manifested in the Paris Agreement and the Sustainable Development Goals (SDGs), respectively. These are inherently inter-linked as progress towards some of these objectives may accelerate or hinder progress towards others. We investigate how these two agendas influence each other by defining networks of 18 nodes, consisting of the 17 SDGs and climate change, for various groupings of countries. We compute a non-linear measure of conditional dependence, the partial distance correlation, given any subset of the remaining 16 variables. These correlations are treated as weights on edges, and weighted eigenvector centralities are calculated to determine the most important nodes. We find that SDG 6, clean water and sanitation, and SDG 4, quality education, are most central across nearly all groupings of countries. In developing regions, SDG 17, partnerships for the goals, is strongly connected to the progress of other objectives in the two agendas whilst, somewhat surprisingly, SDG 8, decent work and economic growth, is not as important in terms of eigenvector centrality.


title: Residual-type a posteriori error analysis of HDG methods for Neumann boundary control problems

abstract: We study a posteriori error analysis of linear-quadratic boundary control problems under bilateral box constraints on the control which acts through a Neumann type boundary condition. We adopt the hybridizable discontinuous Galerkin method as discretization technique, and the flux variables, the scalar variables and the boundary trace variables are all approximated by polynomials of degree k. As for the control variable, it is discretized by the variational discretization concept. Then an efficient and reliable a posteriori error estimator is introduced, and we prove that the error estimator provides an upper bound and a lower bound for the error. Finally, numerical results are presented to illustrate the performance of the obtained a posteriori error estimator.


title: Quantization Guided JPEG Artifact Correction

abstract: The JPEG image compression algorithm is the most popular method of image compression because of its ability for large compression ratios. However, to achieve such high compression, information is lost. For aggressive quantization settings, this leads to a noticeable reduction in image quality. Artifact correction has been studied in the context of deep neural networks for some time, but the current state-of-the-art methods require a different model to be trained for each quality setting, greatly limiting their practical application. We solve this problem by creating a novel architecture which is parameterized by the JPEG files quantization matrix. This allows our single model to achieve state-of-the-art performance over models trained for specific quality settings.


title: Combining multimodal information for Metal Artefact Reduction: An unsupervised deep learning framework

abstract: Metal artefact reduction (MAR) techniques aim at removing metal-induced noise from clinical images. In Computed Tomography (CT), supervised deep learning approaches have been shown effective but limited in generalisability, as they mostly rely on synthetic data. In Magnetic Resonance Imaging (MRI) instead, no method has yet been introduced to correct the susceptibility artefact, still present even in MAR-specific acquisitions. In this work, we hypothesise that a multimodal approach to MAR would improve both CT and MRI. Given their different artefact appearance, their complementary information can compensate for the corrupted signal in either modality. We thus propose an unsupervised deep learning method for multimodal MAR. We introduce the use of Locally Normalised Cross Correlation as a loss term to encourage the fusion of multimodal information. Experiments show that our approach favours a smoother correction in the CT, while promoting signal recovery in the MRI.


title: Protecting a Bosonic Qubit with Autonomous Quantum Error Correction

abstract: To build a universal quantum computer from fragile physical qubits, effective implementation of quantum error correction (QEC) is an essential requirement and a central challenge. Existing demonstrations of QEC are based on a schedule of discrete error syndrome measurements and adaptive recovery operations. These active routines are hardware intensive, prone to introducing and propagating errors, and expected to consume a vast majority of the processing power in a large-scale quantum computer. In principle, QEC can be realized autonomously and continuously by tailoring dissipation within the quantum system, but this strategy has remained challenging so far. Here we encode a logical qubit in Schr\"odinger cat-like multiphoton states of a superconducting cavity, and demonstrate a corrective dissipation process that directly stabilizes an error syndrome operator: the photon number parity. Implemented with continuous-wave control fields only, this passive protocol realizes autonomous correction against single-photon loss and boosts the coherence time of the multiphoton qubit by over a factor of two. Notably, QEC is realized in a modest hardware setup with neither high-fidelity readout nor fast digital feedback, in contrast to the technological sophistication required for prior QEC demonstrations. Compatible with other error suppression and phase stabilization techniques, our experiment suggests reservoir engineering as a resource-efficient alternative or supplement to active QEC in future quantum computing architectures.


title: Locality of Interatomic Interactions in Self-Consistent Tight Binding Models

abstract: A key starting assumption in many classical interatomic potential models for materials is a site energy decomposition of the potential energy surface into contributions that only depend on a small neighbourhood. Under a natural stability condition, we construct such a spatial decomposition for self-consistent tight binding models, extending recent results for linear tight binding models to the non-linear setting.


title: Constraint-based Wavevector- and Frequency-dependent Exchange-Correlation Kernel of the Uniform Electron Gas

abstract: According to time-dependent density functional theory, the exact exchange-correlation kernel f$_{xc}$(n, q, $\omega$) determines not only the ground-state energy but also the excited-state energies/lifetimes and time-dependent linear density response of an electron gas of uniform density n $=$ 3/(4$\pi$r$^3_s$). Here we propose a parametrization of this function based upon the satisfaction of exact constraints. For the static ($\omega$ = 0) limit, we modify the model of Constantin and Pitarke at small wavevector q to recover the known second-order gradient expansion, plus other changes. For all frequencies $\omega$ at q $=$ 0, we use the model of Gross, Kohn, and Iwamoto. A Cauchy integral extends this model to complex $\omega$ and implies the standard Kramers-Kronig relations. A scaling relation permits closed forms for not only the imaginary but also the real part of f$_{xc}$ for real $\omega$. We then combine these ingredients by damping out the $\omega$ dependence at large q in the same way that the q dependence is damped. Away from q $=$ 0 and $\omega$ $=$ 0, the correlation contribution to the kernel becomes dominant over exchange, even at r$_s$ $=$ 4, the valence electron density of metallic sodium. The resulting correlation energy from integration over imaginary $\omega$ is essentially exact. The plasmon pole of the density response function is found by analytic continuation of f$_{xc}$ to $\omega$ just below the real axis, and the resulting plasmon lifetime first decreases from infinity and then increases as q grows from 0 toward the electron-hole continuum. A static charge-density wave is found for r$_s$ $>$ 69, and shown to be associated with softening of the plasmon mode. The exchange-only version of our static kernel confirms Overhauser's 1968 prediction that correlation enhances the charge-density wave.


title: Boundary amenability and measure equivalence rigidity among two-dimensional Artin groups of hyperbolic type

abstract: We study $2$-dimensional Artin groups of hyperbolic type from the viewpoint of measure equivalence, and establish rigidity theorems. We first prove that they are boundary amenable. So is every group acting discretely by simplicial isometries on a connected piecewise hyperbolic $\mathrm{CAT}(-1)$ simplicial complex with countably many simplices in finitely many isometry types, assuming that vertex stabilizers are boundary amenable. Consequently, they satisfy the Novikov conjecture. We then show that measure equivalent $2$-dimensional Artin groups of hyperbolic type have isomorphic fixed set graphs -- an analogue of the curve graph, introduced by Crisp. This yields classification results. We obtain strong rigidity theorems. Let $G=G_\Gamma$ be a $2$-dimensional Artin group of hyperbolic type, with $\mathrm{Out}(G)$ finite. When the automorphism groups of the fixed set graph and of the Cayley complex $\mathfrak{C}$ coincide, every countable group $H$ which is measure equivalent to $G$, is commensurable to a lattice in $\mathrm{Aut}(\mathfrak{C})$. This happens whenever $\Gamma$ is triangle-free with all labels at least $3$ -- unless $G$ is commensurable to the direct sum of $\mathbb{Z}$ and a free group. When $\Gamma$ satisfies an additional star-rigidity condition, then $\mathrm{Aut}(\mathfrak{C})$ is countable, and $H$ is almost isomorphic to $G$. This has applications to orbit equivalence rigidity, and rigidity results for von Neumann algebras associated to ergodic actions of Artin groups. We also derive a rigidity statement regarding possible lattice envelopes of certain Artin groups, and a cocycle superrigidity theorem from higher-rank lattices to $2$-dimensional Artin groups of hyperbolic type.


title: Nielsen equivalence in closed 2-orbifold groups

abstract: We prove that any generating tuple of the fundamental group of a sufficiently large 2-dimensional orbifold is represented by an almost orbifold covering. As a corollary we obtain a generalization of Louder's Theorem which asserts that any two generating tuples of the fundamental group of a closed surface are Nielsen equivalent.


title: Tracemax: A Novel Single Packet IP Traceback Strategy for Data-Flow Analysis

abstract: The identification of the exact path that packets are routed on in the network is quite a challenge. This paper presents a novel, efficient traceback strategy named Tracemax in context of a defense system against distributed denial of service (DDoS) attacks. A single packet can be directly traced over many more hops than the current existing techniques allow. In combination with a defense system it differentiates between multiple connections. It aims to letting non-malicious connections pass while bad ones get thwarted. The novel concept allows detailed analyses of the traffic and the transmission path through the network. The strategy can effectively reduce the effect of common bandwidth and resource consumption attacks, foster early warning and prevention as well as higher the availability of the network services for the wanted customers.


title: Gaussian Decomposition of {\lambda}21-cm Interstellar HI profiles

abstract: Following an established protocol of science, that results must be reproducible, we examine the Gaussian fits to Galactic 21-cm emission profiles obtained by two seemingly complementary methods: the semi-automated approach based on the method used by Verschuur (2004) and the automated technique of Nidever et al. (2008). Both methods use data from the Leiden/Argentine/Bonn all-sky survey. The appeal of an automated routine is great, if for no other reason than the time saved over semi-automated fits. The pitfalls, however, are often unanticipated, and the most important aspect of any algorithm is the reproducibility of the results. The comparisons led to the identification of four problems with the Nidever et al. (2008) analysis: (1) different methods of calculating the reduced chi-squared measuring the goodness of fit; (2) an ultra-broad component found bridging the gap between low and intermediate velocity gas; (3) the lack of an imposed spatial coherence allowing different components to appear and disappear in profiles separated by a fraction of a beam width; and (4) multiple, fundamentally different solutions for the profiles at both the North and South Galactic Poles. A two-step method would improve the algorithm, where an automated fit is followed by a quality-assurance, visual inspection. Confirming evidence emerges from this study of a pervasive component with a line width of order 34 km/s, which may be explained by the Critical Ionization Velocity (CIV) of helium. Since the Nidever et al. (2008) paper contains the only result in the refereed literature that contradicts the CIV model, it is important to understand the flaws in the analysis that let to this contradiction.


title: Robust Partial Matching for Person Search in the Wild

abstract: Various factors like occlusions, backgrounds, etc., would lead to misaligned detected bounding boxes , e.g., ones covering only portions of human body. This issue is common but overlooked by previous person search works. To alleviate this issue, this paper proposes an Align-to-Part Network (APNet) for person detection and re-Identification (reID). APNet refines detected bounding boxes to cover the estimated holistic body regions, from which discriminative part features can be extracted and aligned. Aligned part features naturally formulate reID as a partial feature matching procedure, where valid part features are selected for similarity computation, while part features on occluded or noisy regions are discarded. This design enhances the robustness of person search to real-world challenges with marginal computation overhead. This paper also contributes a Large-Scale dataset for Person Search in the wild (LSPS), which is by far the largest and the most challenging dataset for person search. Experiments show that APNet brings considerable performance improvement on LSPS. Meanwhile, it achieves competitive performance on existing person search benchmarks like CUHK-SYSU and PRW.


title: Convex Analysis and Duality

abstract: Convexity is an important notion in non linear optimization theory as well as in infinite dimensional functional analysis. As will be seen below, very simple and powerful tools will be derived from elementary duality arguments (which are byproducts of the Moreau-Fenchel transform and Hahn Banach Theorem). We will emphasize on applications to a large range of variational problems. Some arguments of measure theory will be skipped.


title: Mott transition in two-band fermion model with on-site Coulomb repulsion

abstract: We provide analytical and numerical solution of the two band fermion model with on-site Coulomb at half filling. In limiting cases for generate bands and one flat band, the model reduces to the Hubbard and Falicov-Kimball models, respectively. We have shown that the insulator state emerges at half filling due to hybridization of fermions of different bands with momenta k and k$+\pi$. Such hybridization breaks the conservation of the number of particles in each band, the Mott transition is a consequence of spontaneous symmetry breaking. A gap in the spectrum is calculated depending on the magnitude of on-site Coulomb repulsion and the width of the band for the chain, as well as for square and cubic lattices. The proposed approach allows us to describe the formation of the gap in the fermion spectra in the Hubbard and Falicov-Kimball models within the framework of the same mechanism for an arbitrary dimension of the system.


title: Topological Quantum Walk with Discrete Time-Glide Symmetry

abstract: Discrete quantum walks are periodically driven systems with discrete time evolution. In contrast to ordinary Floquet systems, no microscopic Hamiltonian exists, and the one-period time evolution is given directly by a series of unitary operators. Regarding each constituent unitary operator as a discrete time step, we formulate discrete space-time symmetry in quantum walks and evaluate the corresponding symmetry protected topological phases. In particular, we study chiral and/or time-glide symmetric topological quantum walks in this formalism. Due to discrete nature of time evolution,the topological classification is found to be different from that in conventional Floquet systems. As a concrete example, we study a two-dimensional quantum walk having both chiral and time-glide symmetries, and identify the anomalous edge states protected by these symmetries.


title: On Eagleson's theorem in the non-stationary setup

abstract: The classical Eagleson's theorem states that if appropriately normalized Birkhoff sums generated by a measurable function and a probability preserving transformation converge in distribution, then they also converge in distribution with respect to any probability measure which is absolutely continuous with respect to the invariant one. In this short note we prove several versions of Eagleson's theorem for some classes of non-stationary stochastic processes which satisfy certain type of decay of correlations.


title: Double- and simple-layer potentials for generalized singular elliptic equations and their applications to the solving the Dirichlet problem

abstract: Potentials play an important role in solving boundary value problems for elliptic equations. In the middle of the last century, a potential theory was constructed for a two-dimensional elliptic equation with one singular coefficient. In the study of potentials, the properties of the fundamental solutions of the given equation are essentially and fruitfully used. At the present time, fundamental solutions of a multidimensional elliptic equation with several singular coefficients are already known. In this paper, we investigate the double- and simple-layer potentials for this kind of elliptic equations. Results from potential theory allow us to represent the solution of the boundary value problems in integral equation form. By using a decomposition formula and other identities for the Lauricella's hypergeometric function in many variables, we prove limiting theorems and derive integral equations concerning a densities of the double- and simple-layer potentials. The obtained results are applied to find an explicit solution of the Dirichlet problem for the generalized singular elliptic equation in the some part of the multidimensional ball.


title: Continuous-Discrete Filtering and Smoothing on Submanifolds of Euclidean Space

abstract: In this paper the issue of filtering and smoothing in continuous discrete time is studied when the state variable evolves in some submanifold of Euclidean space, which may not have the usual Lebesgue measure. Formal expressions for prediction and smoothing problems are derived, which agree with the classical results except that the formal adjoint of the generator is different in general. For approximate filtering and smoothing the projection approach is taken, where it turns out that the prediction and smoothing equations are the same as in the case when the state variable evolves in Euclidean space. The approach is used to develop projection filters and smoothers based on the von Mises-Fisher distribution.


title: A brief and personal history of stochastic partial differential equations

abstract: We trace the evolution of the theory of stochastic partial differential equations from the foundation to its development, until the recent solution of long-standing problems on well-posedness of the KPZ equation and the stochastic quantization in dimension three.


title: Design-controlled Synthesis of IrO$_2$ sub-monolayers on Au Nanodendrites: Marrying Plasmonic and Electrocatalytic Properties

abstract: We develop herein plasmonic-catalytic Au-IrO$_2$ nanostructures with a morphology optimized for efficient light harvesting and catalytic surface area; the nanoparticles have a dendritic morphology, with closely spaced Au branches all partially covered by an ultrathin (1 nm) IrO$_2$ shell. This nanoparticle architecture optimizes optical features due to the interactions of closely spaced plasmonic branches forming electromagnetic hot spots, and the ultra-thin IrO$_2$ layer maximizes efficient use of this expensive catalyst. This concept was evaluated towards the enhancement of the electrocatalytic performances towards the oxygen evolution reaction (OER) as a model transformation. The OER can play a central role in meeting future energy demands but the performance of conventional electrocatalysts in this reaction is limited by the sluggish OER kinetics. We demonstrate an improvement of the OER performance for one of the most active OER catalysts, IrO$_2$, by harvesting plasmonic effects from visible light illumination in multimetallic nanoparticles. We find that the OER activity for the Au-IrO$_2$ nanodendrites can be improved under LSPR excitation, matching best properties reported in the literature. Our simulations and electrocatalytic data demonstrate that the enhancement in OER activities can be attributed to an electronic interaction between Au and IrO$_2$ and to the activation of Ir-O bonds by LSPR excited hot holes, leading to a change in the reaction mechanism (rate-determinant step) under visible light illumination.


title: Augmented Curation of Unstructured Clinical Notes from a Massive EHR System Reveals Specific Phenotypic Signature of Impending COVID-19 Diagnosis

abstract: Understanding the temporal dynamics of COVID-19 patient phenotypes is necessary to derive fine-grained resolution of pathophysiology. Here we use state-of-the-art deep neural networks over an institution-wide machine intelligence platform for the augmented curation of 15.8 million clinical notes from 30,494 patients subjected to COVID-19 PCR diagnostic testing. By contrasting the Electronic Health Record (EHR)-derived clinical phenotypes of COVID-19-positive (COVIDpos, n=635) versus COVID-19-negative (COVIDneg, n=29,859) patients over each day of the week preceding the PCR testing date, we identify anosmia/dysgeusia (37.4-fold), myalgia/arthralgia (2.6-fold), diarrhea (2.2-fold), fever/chills (2.1-fold), respiratory difficulty (1.9-fold), and cough (1.8-fold) as significantly amplified in COVIDpos over COVIDneg patients. The specific combination of cough and diarrhea has a 3.2-fold amplification in COVIDpos patients during the week prior to PCR testing, and along with anosmia/dysgeusia, constitutes the earliest EHR-derived signature of COVID-19 (4-7 days prior to typical PCR testing date). This study introduces an Augmented Intelligence platform for the real-time synthesis of institutional knowledge captured in EHRs. The platform holds tremendous potential for scaling up curation throughput, with minimal need for retraining underlying neural networks, thus promising EHR-powered early diagnosis for a broad spectrum of diseases.


title: Vacua in novel 4D Einstein-Gauss-Bonnet Gravity: pathology and instability?

abstract: We show an inconsistence of the novel 4D Einstein-Gauss-Bonnet gravity by considering a quantum tunneling process of vacua. Using standard semi-classical techniques, we analytically study the vacuum decay rate for all allowed cases in the parameter space. It turns out, without exception, that the theory either encounters a disastrous divergence of vacuum decay rate, or exhibits a confusing complex value of vacuum decay rate, or involves an infinitely instability (a large vacuum mixing). These suggest a strong possibility that the theory, at least the vacuum of the theory, is either unphysical or unstable, or has no well-defined limit as $D\rightarrow 4$.


title: On the Evaluation of Military Simulations: Towards A Taxonomy of Assessment Criteria

abstract: In the area of military simulations, a multitude of different approaches is available. Close Combat Tactical Trainer, Joint Tactical Combat Training System, Battle Force Tactical Training or Warfighter's Simulation 2000 are just some examples within the history of the large DoD Development Program in Modelling and Simulation, representing just a small piece of the variety of diverse solutions. Very often, individual simulators are very unique and so it is often difficult to classify military simulations even for experienced users. This circumstance is further boosted due to the fact that in the field of military simulations - unlike in other areas - no general classification for military simulations exists. To address this shortcoming, this publication is dedicated to the idea of providing a first contribution to the development of a commonly accepted taxonomy in the area of military simulations. To this end, the problem field is structured into three main categories (general functional requirements for simulators, special military requirements for simulators and non-functional requirements for simulators). Based upon that, individual categories are provided with appropriate classes. For a better understanding, the taxonomy is also applied to a concrete example (NetLogo Rebellion).


title: Uniform H\"older-norm bounds for finite element approximations of second-order elliptic equations

abstract: We develop a discrete counterpart of the De Giorgi-Nash-Moser theory, which provides uniform H\"older-norm bounds on continuous piecewise affine finite element approximations of second-order linear elliptic problems of the form $-\nabla \cdot(A\nabla u)=f-\nabla\cdot F$ with $A\in L^\infty(\Omega;\mathbb{R}^{n\times n})$ a uniformly elliptic matrix-valued function, $f\in L^{q}(\Omega)$, $F\in L^p(\Omega;\mathbb{R}^n)$, with $p > n$ and $q > n/2$, on $A$-nonobtuse shape-regular triangulations, which are not required to be quasi-uniform, of a bounded polyhedral Lipschitz domain $\Omega \subset \mathbb{R}^n$.


title: Hot-Starting the Ac Power Flow with Convolutional Neural Networks

abstract: Obtaining good initial conditions to solve the Newton-Raphson (NR) based ac power flow (ACPF) problem can be a very difficult task. In this paper, we propose a framework to obtain the initial bus voltage magnitude and phase values that decrease the solution iterations and time for the NR based ACPF model, using the dc power flow (DCPF) results and one dimensional convolutional neural networks (1D CNNs). We generate the dataset used to train the 1D CNNs by sampling from a distribution of load demands, and by computing the DCPF and ACPF results for each sample. Experiments on the IEEE 118-bus and \textsc{Pegase} 2869-bus study systems show that we can achieve 33.56\% and 30.06\% reduction in solution time, and 66.47% and 49.52% reduction in solution iterations per case, respectively. We include the 1D CNN architectures and the hyperparameters used, which can be expanded on by the future studies on this topic.


title: Compatibility of weak approximation for zero-cycles on products of varieties

abstract: Zero-cycles are conjectured to satisfy weak approximation with Brauer-Manin obstruction for proper smooth varieties defined over number fields. Roughly speaking, we prove that the conjecture is compatible for products of rationally connected varieties, K3 surfaces, Kummer varieties, and one curve.


title: New veto hodoscope ANTI-0 for the NA62 experiment at CERN

abstract: The NA62 experiment is a fixed-target experiment at CERN SPS. The main goal of the experiment is to measure the branching ratio of the ultra-rare kaon decay $K^{+}\to\pi^{+}\nu\bar\nu$. The NA62 detector allows also to study other rare kaon decays and to search for very weakly coupled particles of MeV-GeV mass-scale. The new ANTI-0 hodoscope is proposed and designed to veto events with charged halo particles entering the decay volume. It is now being assembled at CERN. The detector design, the performance simulation and the results of measurements with cosmic rays and test beams for the individual elements are presented. The commissioning and the first run of data-taking is scheduled for after the end of LS2 long shutdown (April 2021).


title: Deep Unfolded Multicast Beamforming

abstract: Multicast beamforming is a promising technique for multicast communication. Providing an efficient and powerful beamforming design algorithm is a crucial issue because multicast beamforming problems such as a max-min-fair problem are NP-hard in general. Recently, deep learning-based approaches have been proposed for beamforming design. Although these approaches using deep neural networks exhibit reasonable performance gain compared with conventional optimization-based algorithms, their scalability is an emerging problem for large systems in which beamforming design becomes a more demanding task. In this paper, we propose a novel deep unfolded trainable beamforming design with high scalability and efficiency. The algorithm is designed by expanding the recursive structure of an existing algorithm based on projections onto convex sets and embedding a constant number of trainable parameters to the expanded network, which leads to a scalable and stable training process. Numerical results show that the proposed algorithm can accelerate its convergence speed by using unsupervised learning, which is a challenging training process for deep unfolding.


title: Model of a spin-$1/2$ electric charge in $F\left( B^2\right)$ modified Weyl gravity

abstract: Within $F\left( B^2\right)$ modified Weyl gravity, we consider a model of a spin-$1/2$ electric charge consisting of interior and exterior regions. The interior region is determined by quantum gravitational effects whose approximate description is carried out using Weyl gravity nonminimally coupled to a massless Dirac spinor field. The interior region is embedded in exterior Minkowski spacetime, and the joining surface is a two-dimensional torus. It is shown that mass, electric charge, and spin of the object suggested may be the same as those for a real electron.


title: End-to-End Whisper to Natural Speech Conversion using Modified Transformer Network

abstract: Machine recognition of an atypical speech like whispered speech, is a challenging task. We introduce whisper-to-natural-speech conversion using sequence-to-sequence approach by proposing enhanced transformer architecture, which uses both parallel and non-parallel data. We investigate different features like Mel frequency cepstral coefficients and smoothed spectral features. The proposed networks are trained end-to-end using supervised approach for feature-to-feature transformation. Further, we also investigate the effectiveness of embedded auxillary decoder used after N encoder sub-layers, trained with the frame-level objective function for identifying source phoneme labels. We show results on opensource wTIMIT and CHAINS datasets by measuring word error rate using end-to-end ASR and also BLEU scores for the generated speech. Alternatively, we also propose a novel method to measure spectral shape of it by measuring formant distributions w.r.t. reference speech, as formant divergence metric. We have found whisper-to-natural converted speech formants probability distribution is similar to the groundtruth distribution. To the authors' best knowledge, this is the first time enhanced transformer has been proposed, both with and without auxiliary decoder for whisper-to-natural-speech conversion and vice versa.


title: Propagative oscillations in co-directional polariton waveguide couplers

abstract: We report on novel exciton-polariton routing devices created to study and purposely guide light-matter particles in their condensate phase. In a co-directional coupling device, two waveguides are connected by a partially etched section which facilitates tunable coupling of the adjacent channels. This evanescent coupling of the two macroscopic wavefunctions in each waveguide reveals itself in real space oscillations of the condensate. This Josephson-like oscillation has only been observed in coupled polariton traps so far. Here, we report on a similar coupling behavior in a controllable, propagative waveguide-based design. By controlling the gap width, channel length or the propagation energy, the exit port of the polariton flow can be chosen. This co-directional polariton device is a passive and scalable coupler element that can serve in compact, next generation logic architectures.


title: Microrheological approach for the viscoelastic response of gels

abstract: In this paper I present a simple and self-consistent framework based on microrheology that allows one to obtain the mechanical response of viscoelastic fluids and gels from the motion of probe particles immersed on it. By considering a non-markovian Langevin equation, I obtain general expressions for the mean-squared displacement and the time-dependent diffusion coefficient that are directly related to the memory kernels and the response function, and which allow one to obtain estimates for the complex shear modulus and the complex viscosity of the material. The usefulness of such approach is demonstrated by applying it to describe experimental data on chemically cross-linked polyacrylamide through its sol-gel transition.


title: MorphStore: Analytical Query Engine with a Holistic Compression-Enabled Processing Model

abstract: In this paper, we present MorphStore, an open-source in-memory columnar analytical query engine with a novel holistic compression-enabled processing model. Basically, compression using lightweight integer compression algorithms already plays an important role in existing in-memory column-store database systems, but mainly for base data. In particular, during query processing, these systems only keep the data compressed until an operator cannot process the compressed data directly, whereupon the data is decompressed, but not recompressed. Thus, the full potential of compression during query processing is not exploited. To overcome that, we developed a novel compression-enabled processing model as presented in this paper. As we are going to show, the continuous usage of compression for all base data and all intermediates is very beneficial to reduce the overall memory footprint as well as to improve the query performance.


title: Geometry optimisation of a transparent axisymmetric ion trap for the MORA project

abstract: In the frame of the project MORA (Matter's Origin from the Radio Activity of trapped and oriented ions), a transparent axially symmetric radio-frequency ion trap (MORATrap) was designed in order to measure the triple correlation parameter $D$ in nuclear $\beta-$decay of laser-polarised ions. The trap design was inspired from the LPCTrap geometry, operated at GANIL from 2005 to 2013. In a real (non-ideal) Paul trap, the quadrupole electric potential is not perfect leading to instabilities in ion motion and therefore affecting the overall trapping efficiency. This paper presents a numerical method aiming to optimise the geometry of a trap. It is applied to MORATrap in order to improve the trapping efficiency and to enlarge the axial transparent solid angle compared to LPCTrap. In the whole optimisation process, numerical computation of electric potential and field was carried out using an electrostatic solver based on boundary element method (BEM). The optimisation consisted in minimising an objective function (fitness function) depending on higher order multipoles of the potential. Finally, systematic changes of trap dimensions and electrode displacements were applied to investigate geometrical effects on the potential quality.


title: Smart Radio Environments Empowered by Reconfigurable Intelligent Surfaces: How it Works, State of Research, and Road Ahead

abstract: What is a reconfigurable intelligent surface? What is a smart radio environment? What is a metasurface? How do metasurfaces work and how to model them? How to reconcile the mathematical theories of communication and electromagnetism? What are the most suitable uses and applications of reconfigurable intelligent surfaces in wireless networks? What are the most promising smart radio environments for wireless applications? What is the current state of research? What are the most important and challenging research issues to tackle? These are a few of the many questions that we investigate in this short opus, which has the threefold objective of introducing the emerging research field of smart radio environments empowered by reconfigurable intelligent surfaces, putting forth the need of reconciling and reuniting C. E. Shannon's mathematical theory of communication with G. Green's and J. C. Maxwell's mathematical theories of electromagnetism, and reporting pragmatic guidelines and recipes for employing appropriate physics-based models of metasurfaces in wireless communications.


title: Irrationality of growth constants associated with polynomial recursions

abstract: We consider integer sequences that satisfy a recursion of the form $x_{n+1} = P(x_n)$ for some polynomial $P$ of degree $d > 1$. If such a sequence tends to infinity, then it satisfies an asymptotic formula of the form $x_n \sim A \alpha^{d^n}$, but little can be said about the constant $\alpha$. In this paper, we show that $\alpha$ is always irrational or an integer. In fact, we prove a stronger statement: if a sequence $G_n$ satisfies an asymptotic formula of the form $G_n = A \alpha^n + B + O(\alpha^{-\epsilon n})$, where $A,B$ are algebraic and $\alpha > 1$, and the sequence contains infinitely many integers, then $\alpha$ is irrational or an integer.


title: The Flow Physics of COVID-19

abstract: Flow physics plays a key role in nearly every facet of the COVID-19 pandemic. This includes the generation and aerosolization of virus-laden respiratory droplets from a host, its airborne dispersion and deposition on surfaces, as well as the subsequent inhalation of these bioaerosols by unsuspecting recipients. Fluid dynamics is also key to preventative measures such as the use of face masks, hand-washing, ventilation of indoor environments, and even social distancing. This article summarizes what we know, and more importantly, what we need to learn about the science underlying these issues so that we are better prepared to tackle the next outbreak of COVID-19 or a similar disease.


title: Stable pair invariants of local Calabi-Yau 4-folds

abstract: In 2008, Klemm-Pandharipande defined Gopakumar-Vafa type invariants of a Calabi-Yau 4-fold $X$ using Gromov-Witten theory. Recently, Cao-Maulik-Toda proposed a conjectural description of these invariants in terms of stable pair theory. When $X$ is the total space of the sum of two line bundles over a surface $S$, and all stable pairs are scheme theoretically supported on the zero section, we express stable pair invariants in terms of intersection numbers on Hilbert schemes of points on $S$. As an application, we obtain new verifications of the Cao-Maulik-Toda conjectures for low degree curve classes and find connections to Carlsson-Okounkov numbers. Some of our verifications involve genus zero Gopakumar-Vafa type invariants recently determined in the context of the log-local principle by Bousseau-Brini-van Garrel. Finally, using the vertex formalism, we provide a few more verifications of the Cao-Maulik-Toda conjectures when thickened curves contribute and also for the case of local $\mathbb{P}^3$.


title: Searching for obscured AGN in z $\sim$ 2 submillimetre galaxies

abstract: Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $\mu$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales.


title: The Manin-Peyre conjecture for smooth spherical Fano varieties of semisimple rank one

abstract: The Manin-Peyre conjecture is established for a class of smooth spherical Fano varieties of semisimple rank one. This includes all smooth spherical Fano threefolds of type T as well as some higher-dimensional smooth spherical Fano varieties.


title: Fourier decay of self-similar measures and self-similar sets of uniqueness

abstract: In this paper, we investigate the Fourier transform of self-similar measures on R. We provide quantitative decay rates of Fourier transform of some self-similar measures. Our method is based on random walks on lattices and Diophantine approximation in number fields. We also completely identify all self-similar sets which are sets of uniqueness. This generalizes a classical result of Salem and Zygmund.


title: A 4-6 GHz Radio Recombination Line Survey in the Milky Way

abstract: We performed a radio recombination line (RRL) survey to construct a high-mass star-forming region (HMSFR) sample in the Milky Way based on the all-sky Wide-Field Infrared Survey Explorer ($\textit{All-WISE}$) point source catalog. The survey was observed with the Shanghai 65m Tianma radio telescope (TMRT) covering 10 hydrogen RRL transitions ranging from H98$\alpha$ to H113$\alpha$ (corresponding to the rest frequencies of 4.5$-$6.9 GHz) simultaneously. Out of 3348 selected targets, we identified an HMSFR sample consisting of 517 sources traced by RRLs, a large fraction of this sample (486) locate near the Galactic plane ($|$$\textit{b}$$|$ $<$ 2 deg). In addition to the hydrogen RRLs, we also detected helium and carbon RRLs towards 49 and 23 sources respectively. We cross-match the RRL detections with the 6.7 methanol maser sources built up in previous works for the same target sample, as a result, 103 HMSFR sources were found to harbor both emissions. In this paper, we present the HMSFR catalog accompanied by the measured RRL line properties and a correlation with our methanol maser sample, which is believed to tracer massive stars at earlier stages. The construction of an HMSFR sample consisting of sources in various evolutionary stages indicated by different tracers is fundamental for future studies of high-mass star formation in such regions.


title: Use of a Generalized Energy Mover's Distance in the Search for Rare Phenomena at Colliders

abstract: In this paper, we expand on the previously proposed concept of Energy Mover's Distance. The resulting observables are shown to provide a way of identifying rare processes in proton-proton collider experiments. It is shown that different processes are grouped together differently and that this can contribute to the improvement of experimental analyses. The $t\bar{t}Z$ production at the Large Hadron Collider is used as a benchmark to illustrate the applicability of the method. Furthermore, we study the use of these observables as new features which can be used in the training of Deep Neural Networks.


title: Stress mirror polishing for future large lightweight mirrors: design using shape optimization

abstract: This study proposes a new way to manufacture large lightweight aspherics for space telescopes using Stress Mirror Polishing (SMP). This technique is well known to allow reaching high quality optical surfaces in a minimum time period, thanks to a spherical full-size polishing tool. To obtain the correct surface' aspheric shape, it is necessary to define precisely the thickness distribution of the mirror to be deformed, according to the manufacturing parameters. We first introduce active optics and Stress Mirror Polishing techniques, and then, we describe the process to obtain the appropriate thickness mirror distribution allowing to generate the required aspheric shape during polishing phase. Shape optimization procedure using PYTHON programing and NASTRAN optimization solver using Finite Element Model (FEM) is developed and discussed in order to assist this process. The main result of this paper is the ability of the shape optimization process to support SMP technique to generate a peculiar aspherical shape from a spherical optical surface thanks to a thickness distribution reshaping. This paper is primarily focused on a theoretical framework with numerical simulations as the first step before the manufacturing of a demonstrator. This two-steps approach was successfully used for previous projects.


title: A Generalization of the Allreduce Operation

abstract: Allreduce is one of the most frequently used MPI collective operations, and thus its performance attracts much attention in the past decades. Many algorithms were developed with different properties and purposes. We present a novel approach to communication description based on the permutations inspired by the mathematics of a Rubik's cube where the moves form a mathematical structure called group. Similarly, cyclic communication patterns between a set of $P$ processes may be described by a permutation group. This new approach allows constructing a generalization of the widely used Allreduce algorithms such as Ring, Recursive Doubling and Recursive Halving. Using the developed approach we build an algorithm that successfully solves the well-known problem of the non-power-of-two number of processes which breaks down the performance of many existing algorithms. The proposed algorithm provides a general solution for any number of processes with the dynamically changing amount of communication steps between $\lceil \log{P} \rceil$ for the latency-optimal version and $2 \cdot \lceil \log{P} \rceil$ for the bandwidth-optimal case.


title: Deep-COVID: Predicting COVID-19 From Chest X-Ray Images Using Deep Transfer Learning

abstract: The COVID-19 pandemic is causing a major outbreak in more than 150 countries around the world, having a severe impact on the health and life of many people globally. One of the crucial step in fighting COVID-19 is the ability to detect the infected patients early enough, and put them under special care. Detecting this disease from radiography and radiology images is perhaps one of the fastest ways to diagnose the patients. Some of the early studies showed specific abnormalities in the chest radiograms of patients infected with COVID-19. Inspired by earlier works, we study the application of deep learning models to detect COVID-19 patients from their chest radiography images. We first prepare a dataset of 5,000 Chest X-rays from the publicly available datasets. Images exhibiting COVID-19 disease presence were identified by board-certified radiologist. Transfer learning on a subset of 2,000 radiograms was used to train four popular convolutional neural networks, including ResNet18, ResNet50, SqueezeNet, and DenseNet-121, to identify COVID-19 disease in the analyzed chest X-ray images. We evaluated these models on the remaining 3,000 images, and most of these networks achieved a sensitivity rate of 98% ($\pm$ 3%), while having a specificity rate of around 90%. Besides sensitivity and specificity rates, we also present the receiver operating characteristic (ROC) curve, precision-recall curve, average prediction, and confusion matrix of each model. We also used a technique to generate heatmaps of lung regions potentially infected by COVID-19 and show that the generated heatmaps contain most of the infected areas annotated by our board certified radiologist. While the achieved performance is very encouraging, further analysis is required on a larger set of COVID-19 images, to have a more reliable estimation of accuracy rates. The dataset, model implementations (in PyTorch), and evaluations, are all made publicly available for research community at https://github.com/shervinmin/DeepCovid.git


title: Design and construction of MuTe: a hybrid Muon Telescope to study Colombian Volcanoes

abstract: We present a hybrid Muon Telescope, MuTe, designed and built for imaging active Colombian volcanoes. The MuTe has a resolution of tens of meters, low power consumption, robustness and transportability making it suitable for use in difficult access zones where active volcanoes usually are. The main feature of MuTe is the implementation of a hybrid detection technique combining two scintillation panels for particle tracking and a Water Cherenkov Detector for filtering background sources due to the electromagnetic component of extended air showers and multiple particle events. MuTe incorporates particle-identification techniques for reducing the background noise sources and discrimination of fake events by a picosecond Time-of-Flight system. We also describe the mechanical behavior of the MuTe during typical tremors and wind conditions at the observation place, as well as the frontend electronics design and power consumption.


title: A simple proof of regularity for $C^{1,\alpha}$ interface transmission problems

abstract: We give a simple proof of a recent result in [1] by Caffarelli, Soria-Carro, and Stinga about the $C^{1,\alpha}$ regularity of weak solutions to transmission problems with $C^{1,\alpha}$ interfaces. Our proof does not use the mean value property or the maximum principle, and also works for more general elliptic systems. Some extensions to $C^{1,\text{Dini}}$ interfaces and to domains with multiple sub-domains are also discussed.


title: R package SamplingStrata: new developments and extension to Spatial Sampling

abstract: The R package SamplingStrata was developed in 2011 as an instrument to optimize the design of stratified samples. The optimization is performed by considering the stratification variables available in the sampling frame, and the precision constraints on target estimates of the survey (Ballin & Barcaroli, 2014). The genetic algorithm at the basis of the optimization step explores the universe of the possible alternative stratifications determining for each of them the best allocation, that is the one of minumum total size that allows to satisfy the precision constraints: the final optimal solution is the one that ensures the global minimum sample size. One fundamental requirement to make this approach feasible is the possibility to estimate the variability of target variables in generated strata; in general, as target variable values are not available in the frame, but only proxy ones, anticipated variance is calculated by modelling the relations between target and proxy variables. In case of spatial sampling, it is important to consider not only the total model variance, but also the co-variance derived by the spatial auto-correlation. The last release of SamplingStrata enables to consider both components of variance, thus allowing to harness spatial auto-correlation in order to obtain more efficient samples.


title: ClovaCall: Korean Goal-Oriented Dialog Speech Corpus for Automatic Speech Recognition of Contact Centers

abstract: Automatic speech recognition (ASR) via call is essential for various applications, including AI for contact center (AICC) services. Despite the advancement of ASR, however, most publicly available call-based speech corpora such as Switchboard are old-fashioned. Also, most existing call corpora are in English and mainly focus on open domain dialog or general scenarios such as audiobooks. Here we introduce a new large-scale Korean call-based speech corpus under a goal-oriented dialog scenario from more than 11,000 people, i.e., ClovaCall corpus. ClovaCall includes approximately 60,000 pairs of a short sentence and its corresponding spoken utterance in a restaurant reservation domain. We validate the effectiveness of our dataset with intensive experiments using two standard ASR models. Furthermore, we release our ClovaCall dataset and baseline source codes to be available via https://github.com/ClovaAI/ClovaCall.


title: Awareness of crash risk improves Kelly strategies in simulated financial time series

abstract: We simulate a simplified version of the price process including bubbles and crashes proposed in Kreuser and Sornette (2018). The price process is defined as a geometric random walk combined with jumps modelled by separate, discrete distributions associated with positive (and negative) bubbles. The key ingredient of the model is to assume that the sizes of the jumps are proportional to the bubble size. Thus, the jumps tend to efficiently bring back excess bubble prices close to a normal or fundamental value (efficient crashes). This is different from existing processes studied that assume jumps that are independent of the mispricing. The present model is simplified compared to Kreuser and Sornette (2018) in that we ignore the possibility of a change of the probability of a crash as the price accelerates above the normal price. We study the behaviour of investment strategies that maximize the expected log of wealth (Kelly criterion) for the risky asset and a risk-free asset. We show that the method behaves similarly to Kelly on Geometric Brownian Motion in that it outperforms other methods in the long-term and it beats classical Kelly. As a primary source of outperformance, we determine knowledge about the presence of crashes, but interestingly find that knowledge of only the size, and not the time of occurrence, already provides a significant and robust edge. We then perform an error analysis to show that the method is robust with respect to variations in the parameters. The method is most sensitive to errors in the expected return.


title: A Massive Molecular Torus inside a Gas-Poor Cirnumnuclear Disk in the Radio Galaxy NGC 1052 Discovered with ALMA

abstract: We report ALMA observations of NGC 1052 to quest mass accretion in a gas-poor active galactic nucleus (AGN). We detected CO emission representing a rotating ring-like circumnuclear disk (CND) seen edge-on with the gas mass of $5.3 \times 10^{5}$ M$_{\odot}$. The CND has smaller gas mass than that in typical Seyfert galaxies with circumnuclear star formation and is too gas-poor to drive mass accretion onto the central engine. The continuum emission casts molecular absorption features of CO, HCN, HCO$^+$, SO, SO$_2$, CS, CN, and H$_2$O, with H$^{13}$CN and HC$^{15}$N and vibrationally-excited (v$_2 = 1$) HCN and HCO$^+$. Broader absorption line widths than CND emission line widths imply presence of a geometrically thick molecular torus with a radius of $2.4 \pm 1.3$ pc and a thickness ratio of $0.7 \pm 0.3$. We estimate the H$_2$ column density of $(3.3 \pm 0.7) \times 10^{25}$ cm$^{-2}$ using H$^{12}$CN, H$^{13}$CN, and HCO$^{+}$ absorption features and adopting abundance ratio of $^{12}$C-to-$^{13}$C and a HCO$^{+}$-to-H$_2$, and derived the torus gas mass of $(1.3 \pm 0.3) \times 10^7$ M$_{\odot}$, which is $\sim 9$\% of the central black-hole mass. The molecular gas in the torus is clumpy with the estimated covering factor of $0.17^{+0.06}_{-0.03}$. The gas density of clumps inside the torus is inferred to be $(6.4 \pm 1.3) \times 10^7$ cm$^{-3}$, which meets the excitation conditions of H$_2$O maser. Specific angular momentum in the torus exceeds a flat-rotation curve extrapolated from that of the CND, indicating a Keplerian rotation inside a 14.4-pc sphere of influence.


title: Learning Ising models from one or multiple samples

abstract: There have been two separate lines of work on estimating Ising models: (1) estimating them from multiple independent samples under minimal assumptions about the model's interaction matrix; and (2) estimating them from one sample in restrictive settings. We propose a unified framework that smoothly interpolates between these two settings, enabling significantly richer estimation guarantees from one, a few, or many samples. Our main theorem provides guarantees for one-sample estimation, quantifying the estimation error in terms of the metric entropy of a family of interaction matrices. As corollaries of our main theorem, we derive bounds when the model's interaction matrix is a (sparse) linear combination of known matrices, or it belongs to a finite set, or to a high-dimensional manifold. In fact, our main result handles multiple independent samples by viewing them as one sample from a larger model, and can be used to derive estimation bounds that are qualitatively similar to those obtained in the afore-described multiple-sample literature. Our technical approach benefits from sparsifying a model's interaction network, conditioning on subsets of variables that make the dependencies in the resulting conditional distribution sufficiently weak. We use this sparsification technique to prove strong concentration and anti-concentration results for the Ising model, which we believe have applications beyond the scope of this paper.


title: Complex phase masks for fabricating OH-suppression filters for astronomy

abstract: The design of a complex phase mask (CPM) for inscribing aperiodic filters in fibers is presented. The complex structure of the mask surface relief consists of discrete aperiodic phase-steps at periodic intervals. We show that the diffraction of the inscribing laser beam from the phase-step locations of the CPM produces periodically located half phase-steps along the fiber. The accumulated phase, along with controlled index modulation, generates the desired aperiodic reflection spectrum. Compared to a complex 'running-light' interferometry based inscription method, CPM offers the well known convenience and reproducibility of the standard phase mask inscription technique. The complexity of an aperiodic grating is encoded into the structure of the CPM. Complex filters fabricated with CPM can be used for suppressing the near infrared (NIR) OH- emission lines generated in the upper atmosphere, improving the performance of ground based telescopes.


title: COVID-19 Should be Suppressed by Mixed Constraints -- from Simulations on Constrained Scale-Free Networks

abstract: The spreading of virus infection is here simulated over artificial human networks. Here, the real-space urban life of people is modeled as a scale-free network with constraints. A scale-free network has been adopted for modeling on-line communities so far but is employed here for the aim to represent peoples' social behaviors where the generated communities are restricted reflecting the spatiotemporal constraints in the real life. As a result, three findings and a policy proposal have been obtained. First, the height of the peaks in the time sequence of the number of infection cases tends to get reduced corresponding to the upper bound of the size of groups where all members meet. Second, if we adopt the constraint on m0, the number of all other people one meets separately each at a time, to the range between 2 and 8, its effect on the suppression of infections may be weak as far as we allow group meetings of size W of 8 or larger. Third, such a moderate constraint may temporarily seem to work for the reduction of infections in the early stage but it may turn out to be just a delay of peaks. Based on these results, a policy is proposed here: for quickly suppressing the number of infections, restrict W to less than 4 if the constraint to make m0 at most 1 is too strict. If W is set to less than 4, setting m0 to 4 or less works for quick reduction of infections according to the result.


title: Network-inspired versus Kozeny-Carman based permeability-porosity relations applied to Biot's poroelasticity model

abstract: Water injection in the aquifer induces deformations in the soil. These mechanical deformations give rise to a change in porosity and permeability, which results in non-linearity of the mathematical problem. Assuming that the deformations are very small, the model provided by Biot's theory of linear poroelasticity is used to determine the local displacement of the skeleton of a porous medium, as well as the fluid flow through the pores. In this continuum scale model, the Kozeny-Carman equation is commonly used to determine the permeability of the porous medium from the porosity. The Kozeny-Carman relation states that flow through the pores is possible at a certain location as long as the porosity is larger than zero at this location in the aquifer. However, from network models it is known that percolation thresholds exist, indicating that the permeability will be equal to zero if the porosity becomes smaller than these thresholds. In this paper, the relationship between permeability and porosity is investigated. A new permeability-porosity relation, based on the percolation theory, is derived and compared with the Kozeny-Carman relation. The strongest feature of the new approach is related to its capability to give a good description of the permeability in case of low porosities. However, with this network-inspired approach small values of the permeability are more likely to occur. Since we show that the solution of Biot's model converges to the solution of a saddle point problem for small time steps and low permeability, we need stabilisation in the finite element approximation.


title: Complex-Object Visual Inspection via Multiple Lighting Configurations

abstract: The design of an automatic visual inspection system is usually performed in two stages. While the first stage consists in selecting the most suitable hardware setup for highlighting most effectively the defects on the surface to be inspected, the second stage concerns the development of algorithmic solutions to exploit the potentials offered by the collected data. In this paper, first, we present a novel illumination setup embedding four illumination configurations to resemble diffused, dark-field, and front lighting techniques. Second, we analyze the contributions brought by deploying the proposed setup in training phase only - mimicking the scenario in which an already developed visual inspection system cannot be modified on the customer site - and in evaluation phase. Along with an exhaustive set of experiments, in this paper, we demonstrate the suitability of the proposed setup for effective illumination of complex-objects, defined as manufactured items with variable surface characteristics that cannot be determined a priori. Moreover, we discuss the importance of multiple light configurations availability during training and their natural boosting effect which, without the need to modify the system design in evaluation phase, lead to improvements in the overall system performance.


title: Stochastic Peierls-Nabarro Model for Dislocations in High Entropy Alloys

abstract: High entropy alloys (HEAs) are single phase crystals that consist of random solid solutions of multiple elements in approximately equal proportions. This class of novel materials have exhibited superb mechanical properties, such as high strength combined with other desired features. The strength of crystalline materials is associated with the motion of dislocations. In this paper, we derive a stochastic continuum model based on the Peierls-Nabarro framework for inter-layer dislocations in a bilayer HEA from an atomistic model that incorporates the atomic level randomness. We use asymptotic analysis and limit theorem in the convergence from the atomistic model to the continuum model. The total energy in the continuum model consists of a stochastic elastic energy in the two layers, and a stochastic misfit energy that accounts for the inter-layer nonlinear interaction. The obtained continuum model can be considered as a stochastic generalization of the classical, deterministic Peierls-Nabarro model for the dislocation core and related properties. This derivation also validates the stochastic model adopted by Zhang et al. (Acta Mater. 166, 424-434, 2019).


title: Conditional-UNet: A Condition-aware Deep Model for Coherent Human Activity Recognition From Wearables

abstract: Recognizing human activities from multi-channel time series data collected from wearable sensors is ever more practical. However, in real-world conditions, coherent activities and body movements could happen at the same time, like moving head during walking or sitting. A new problem, so-called "Coherent Human Activity Recognition (Co-HAR)", is more complicated than normal multi-class classification tasks since signals of different movements are mixed and interfered with each other. On the other side, we consider such Co-HAR as a dense labelling problem that classify each sample on a time step with a label to provide high-fidelity and duration-varied support to applications. In this paper, a novel condition-aware deep architecture "Conditional-UNet" is developed to allow dense labeling for Co-HAR problem. We also contribute a first-of-its-kind Co-HAR dataset for head movement recognition under walk or sit condition for future research. Experiments on head gesture recognition show that our model achieve overall 2%-3% performance gain of F1 score over existing state-of-the-art deep methods, and more importantly, systematic and comprehensive improvements on real head gesture classes.


title: An Euler-Maclaurin formula for polygonal sums

abstract: We prove an Euler-Maclaurin formula for double polygonal sums and, as a corollary, we obtain approximate quadrature formulas for integrals of smooth functions over polygons with integer vertices. Our Euler-Maclaurin formula is in the spirit of Pick's theorem on the number of integer points in an integer polygon and involves weighted Riemann sums, using tools from Harmonic analysis. Finally, we also exhibit a classical trick, dating back to Huygens and Newton, to accelerate convergence of these Riemann sums.


title: Exploring the superwind mechanism for generating ultrahigh-energy cosmic rays using large-scale modeling of starbursts

abstract: The Pierre Auger Collaboration has provided a compelling indication for a possible correlation between the arrival directions of ultrahigh-energy cosmic rays and nearby starburst galaxies. Herein we show how the latest large-scale modeling of starburst galaxies is compatible with the cosmic rays producing the anisotropy signal being accelerated at the terminal shock of superwinds.


title: Criteria Based Evaluation Framework for Service Oriented Methodologies

abstract: Service Oriented Software Engineering is based on concepts and principles for constructing complex enterprise systems in which services as building block of the system, are distributed in large networks. The main goal of the service oriented methodologies is to define a process for development and maintenance of service based systems. Most of the Service Oriented methodologies are not mature enough compared with traditional software development methodologies such as Object Oriented or Component-Based. Hence, defining an evaluation framework will be useful for comparing methodologies for identifying their strengths and weaknesses, defining new methodologies or extending existing Service Oriented methodologies. At the time being, there is no complete evaluation framework for evaluating Service Oriented methodologies. The principal objective of this paper is to introduce a comprehensive evaluation framework for evaluating Service-Oriented methodologies. This evaluation tool is appropriate for methodology engineers to develop new methodologies, as well as project managers to select an appropriate methodology at a specific project.


title: A Procedure for Extracting Software Development Process Patterns

abstract: Process patterns represent well-structured and successful recurring activities of Software Development Methodologies. They are able to form a library of reusable building blocks that can be utilized in Situational Method Engineering for constructing a custom SDM or enhancing an existing one to fit specific project situation. Recently, some researchers have subjectively extracted process patterns from existing SDMs based on cumulative experience in various domains; however, how to objectively extract process patterns from SDMs by adopting a systematic procedure has remained as question. In this regard, this paper is concerned with a procedure aiming to take process patterns out of existing SDMs. An example illustrates applicability of the proposed procedure for extracting process patterns in a specific context.


title: Process Patterns for Service Oriented Development

abstract: Software systems development nowadays has moved towards dynamic composition of services that run on distributed infrastructures aligned with continuous changes in the system requirements. Consequently, software developers need to tailor project specific methodologies to fit their methodology requirements. Process patterns present a suitable solution by providing reusable method chunks of software development methodologies for constructing methodologies to fit specific requirements. In this paper, we propose a set of high-level service-oriented process patterns that can be used for constructing and enhancing situational service-oriented methodologies. We show how these patterns are used to construct a specific service-oriented methodology for the development of a sample system. Keywords. Service-Oriented Software Development Methodologies, Process Patterns, Process Meta-Model, Situational Method Engineering


title: Open Clusters as Tracers on Radial Migration of the Galactic Disk

abstract: Radial migration is an important process in the Galactic disk. A few open clusters show some evidence on this mechanism but there is no systematic study. In this work, we investigate the role of radial migration on the Galactic disk based on a large sample of 146 open clusters with homogeneous metallicity and age from Netopil et al. and kinematics calculated from Gaia DR2. The birth site Rb, guiding radius Rg and other orbital parameters are calculated, and the migration distance |Rg-Rb| is obtained, which is a combination of metallicity, kinematics and age information. It is found that 44% open clusters have |Rg-Rb|< 1 kpc, for which radial migration (churning) is not significant. Among the remaining 56% open clusters with |Rg-Rb|> 1 kpc, young ones with t<1.0 Gyr tend to migrate inward, while older clusters usually migrate outward. Different mechanisms of radial migration between young and old clusters are suggested based on their different migration rates, Galactic locations and orbital parameters. For the old group, we propose a plausible way to estimate migration rate and obtain a reasonable value of 1.5(+-0.5) kpc/Gyr based on ten intermediate-age clusters at the outer disk, where the existence of several special clusters implies its complicate formation history.


title: Escaping points of commuting meromorphic functions with finitely many poles

abstract: Let $f$ and $g$ be commuting meromorphic functions with finitely many poles. By studying the behaviour of Fatou components under this commuting relation, we prove that $f$ and $g$ have the same Julia set whenever $f$ and $g$ have no simply connected fast-escaping wandering domains. By combining this with a recent result of Tsantaris', we obtain the strongest statement (to date) regarding the Julia sets of commuting meromorphic functions. In order to highlight the difference to the entire case, we show that transcendental meromorphic functions with finitely many poles have orbits that alternate between approaching a pole and escaping to infinity at strikingly fast rates.


title: Suppressed cosmic growth in coupled vector-tensor theories

abstract: We study a coupled dark energy scenario in which a massive vector field $A_{\mu}$ with broken $U(1)$ gauge symmetry interacts with the four-velocity $u_c^{\mu}$ of cold dark matter (CDM) through the scalar product $Z=-u_c^{\mu} A_{\mu}$. This new coupling corresponds to the momentum transfer, so that the background vector and CDM continuity equations do not have explicit interacting terms analogous to the energy exchange. Hence the observational preference of uncoupled generalized Proca theories over the $\Lambda$CDM model can be still maintained at the background level. Meanwhile, the same coupling strongly affects the evolution of cosmological perturbations. While the effective sound speed of CDM vanishes, the propagation speed and no-ghost condition of a longitudinal scalar of $A_{\mu}$ and the CDM no-ghost condition are subject to nontrivial modifications by the $Z$ dependence in the Lagrangian. We propose a concrete dark energy model and show that the gravitational interaction on scales relevant to the linear growth of large-scale structures can be smaller than the Newton constant at low redshifts. This leads to the suppression of growth rates of both CDM and total matter density perturbations, so our model allows an interesting possibility for reducing the tension of matter density contrast $\sigma_8$ between high- and low-redshift measurements.


title: Linear Expected Complexity for Directional and Multiplicative Voronoi Diagrams

abstract: While the standard unweighted Voronoi diagram in the plane has linear worst-case complexity, many of its natural generalizations do not. This paper considers two such previously studied generalizations, namely multiplicative and semi Voronoi diagrams. These diagrams both have quadratic worst-case complexity, though here we show that their expected complexity is linear for certain natural randomized inputs. Specifically, we argue that the expected complexity is linear for: (1) semi Voronoi diagrams when the visible direction is randomly sampled, and (2) for multiplicative diagrams when either weights are sampled from a constant-sized set, or the more challenging case when weights are arbitrary but locations are sampled from a square.


title: Anomalous decay rate of quasinormal modes in Schwarzschild-dS and Schwarzschild-AdS black holes

abstract: Recently an anomalous decay rate of the quasinormal modes of a massive scalar field in Schwarzschild black holes backgrounds was reported in which the longest-lived modes are the ones with higher angular number, for a scalar field mass smaller than a critical value, while that beyond this value the behaviour is inverted. In this work, we extend the study to other asymptotic geometries, such as, Schwarzschild-de Sitter and Schwarzschild-AdS black holes. Mainly, we found that such behaviour and the critical mass are present in the Schwarzschild-de Sitter background. Also, we found that the value of the critical mass increases when the cosmological constant increases and also when the overtone number is increasing. On the other hand, despite the critical mass is not present in Schwarzschild-AdS black holes backgrounds, the decay rate of the quasinormal modes always exhibits an anomalous behaviour.


title: Spectral theory of fluctuations in time-average statistical mechanics of reversible and driven systems

abstract: We present a spectral-theoretic approach to time-average statistical mechanics for general, non-equilibrium initial conditions. We consider the statistics of bounded, local additive functionals of reversible as well as irreversible ergodic stochastic dynamics with continuous or discrete state-space. We derive exact results for the mean, fluctuations and correlations of time average observables from the eigenspectrum of the underlying generator of Fokker-Planck or master equation dynamics, and discuss the results from a physical perspective. Feynman-Kac formulas are re-derived using It\^o calculus and combined with non-Hermitian perturbation theory. The emergence of the universal central limit law in a spectral representation is shown explicitly on large deviation time-scales. For reversible dynamics with equilibrated initial conditions we derive a general upper bound to fluctuations of occupation measures in terms of an integral of the return probability. Simple, exactly solvable examples are analyzed to demonstrate how to apply the theory. As a biophysical example we revisit the Berg-Purcell problem on the precision of concentration measurements by a single receptor. Our results are directly applicable to a diverse range of phenomena underpinned by time-average observables and additive functionals in physical, chemical, biological, and economical systems.

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: MixPUL: Consistency-based Augmentation for Positive and Unlabeled Learning

abstract: Learning from positive and unlabeled data (PU learning) is prevalent in practical applications where only a couple of examples are positively labeled. Previous PU learning studies typically rely on existing samples such that the data distribution is not extensively explored. In this work, we propose a simple yet effective data augmentation method, coined~\algo, based on \emph{consistency regularization} which provides a new perspective of using PU data. In particular, the proposed~\algo~incorporates supervised and unsupervised consistency training to generate augmented data. To facilitate supervised consistency, reliable negative examples are mined from unlabeled data due to the absence of negative samples. Unsupervised consistency is further encouraged between unlabeled datapoints. In addition,~\algo~reduces margin loss between positive and unlabeled pairs, which explicitly optimizes AUC and yields faster convergence. Finally, we conduct a series of studies to demonstrate the effectiveness of consistency regularization. We examined three kinds of reliable negative mining methods. We show that~\algo~achieves an averaged improvement of classification error from 16.49 to 13.09 on the CIFAR-10 dataset across different positive data amount.


title: Multigrid preconditioners for the hybridized Discontinuous Galerkin discretisation of the shallow water equations

abstract: Numerical climate- and weather-prediction requires the fast solution of the equations of fluid dynamics. Discontinuous Galerkin (DG) discretisations have several advantageous properties. They can be used for arbitrary domains and support a structured data layout, which is important on modern chip architectures. For smooth solutions, higher order approximations can be particularly efficient since errors decrease exponentially in the polynomial degree. Due to the wide separation of timescales in atmospheric dynamics, semi-implicit time integrators are highly efficient, since the implicit treatment of fast waves avoids tight constraints on the time step size, and can therefore improve overall efficiency. However, if implicit-explicit (IMEX) integrators are used, a large linear system of equations has to be solved in every time step. A particular problem for DG discretisations of velocity-pressure systems is that the normal Schur-complement reduction to an elliptic system for the pressure is not possible since the numerical fluxes introduce artificial diffusion terms. For the shallow water equations, which form an important model system, hybridised DG methods have been shown to overcome this issue. However, no attention has been paid to the efficient solution of the resulting linear system of equations. In this paper we address this issue and show that the elliptic system for the flux unknowns can be solved efficiently with a non-nested multigrid algorithm. The method is implemented in the Firedrake library and we demonstrate the excellent performance of the algorithm both for an idealised stationary flow problem in a flat domain and for non-stationary setups in spherical geometry from the Williamson et al. testsuite. In the latter case the performance of our bespoke multigrid preconditioner (although itself not highly optimised) is comparable to that of a highly optimised direct solver.

categories: physics.comp-ph cs.NA math.NA physics.flu-dyn


title: The three body first post-Newtonian effects on the secular dynamics of a compact binary near a spinning supermassive black hole

abstract: The binary black holes (BBHs) formed near the supermassive black holes (SMBHs) in the galactic nuclei would undergo eccentricity excitation due to the gravitational perturbations from the SMBH and therefore merger more efficiently. In this paper, we study the coupling of the three body 1st post-Newtonian (PN) effects with the spin effects from the SMBH in the hierarchical triple system. We extend previous work by including the coupling between the de Sitter precession and the Lense-Thirring precession from the SMBH spin. This coupling includes both the precessions of the inner orbit angular momentum and the Runge-Lenz vector around the outer orbit angular momentum in a general reference frame. We find the change of the (maximal) eccentricity in the neighboring Kozai-Lidov cycles due to spin effects is detectable by LISA in the future. Our general argument on the coupling of the three body 1PN effects in three body systems could be extended to any other situation as long as the outer orbital plane evolves.


title: Mapping the working of environmental effects in A963

abstract: We qualitatively assess and map the relative contribution of pre-processing and cluster related processes to the build-up of A963, a massive cluster at z=0.2 showing an unusually high fraction of star forming galaxies in its interior. We use Voronoi binning of positions of cluster members on the plane of the sky in order to map the 2D variations of galaxy properties in the centre and infall region of A963. We map four galaxy parameters (fraction of star forming galaxies, specific star formation rate, HI deficiency and age of the stellar population) based on full SED fitting, 21cm imaging and optical spectroscopy. We find an extended region dominated by passive galaxies along a north-south axis crossing the cluster centre, possibly associated with known filaments of the large-scale structure. There are signs that the passive galaxies in this region were quenched long before their arrival in the vicinity of the cluster. Contrary to that, to the east and west of the cluster centre lie regions of recent accretion dominated by gas rich, actively star forming galaxies not associated with any substructure or filament. The few passive galaxies in this region appear to be recently quenched, and some gas rich galaxies show signs of ongoing ram-pressure stripping. We report the first tentative observations at 21cm of ongoing ram-pressure stripping at z=0.2, as well as observed inflow of low-entropy gas into the cluster along filaments of the large-scale structure. The observed galaxy content of A963 is a result of strongly anisotropic accretion of galaxies with different properties. Gas rich, star forming galaxies are being accreted from the east and west of the cluster and these galaxies are being quenched at r<R200, while the bulk of the accretion, containing multiple groups, happens along the north-south axis and brings mostly passive galaxies.


title: A non-cooperative meta-modeling game for automated third-party calibrating, validating, and falsifying constitutive laws with parallelized adversarial attacks

abstract: The evaluation of constitutive models, especially for high-risk and high-regret engineering applications, requires efficient and rigorous third-party calibration, validation and falsification. While there are numerous efforts to develop paradigms and standard procedures to validate models, difficulties may arise due to the sequential, manual and often biased nature of the commonly adopted calibration and validation processes, thus slowing down data collections, hampering the progress towards discovering new physics, increasing expenses and possibly leading to misinterpretations of the credibility and application ranges of proposed models. This work attempts to introduce concepts from game theory and machine learning techniques to overcome many of these existing difficulties. We introduce an automated meta-modeling game where two competing AI agents systematically generate experimental data to calibrate a given constitutive model and to explore its weakness, in order to improve experiment design and model robustness through competition. The two agents automatically search for the Nash equilibrium of the meta-modeling game in an adversarial reinforcement learning framework without human intervention. By capturing all possible design options of the laboratory experiments into a single decision tree, we recast the design of experiments as a game of combinatorial moves that can be resolved through deep reinforcement learning by the two competing players. Our adversarial framework emulates idealized scientific collaborations and competitions among researchers to achieve a better understanding of the application range of the learned material laws and prevent misinterpretations caused by conventional AI-based third-party validation.