Neuro Symbolic Reinforcement Learning with First Order Logic. Deep reinforcement learning (RL) methods often require many trials before convergence, and no direct interpretability of trained policies is provided. In order to achieve fast convergence and interpretability for the policy in RL, we propose a novel RL method for text based games with a recent neuro symbolic framework called Logical Neural Network, which can learn symbolic and interpretable rules in their differentiable network. The method is first to extract first order logical facts from text observation and external word meaning network (ConceptNet), then train a policy in the network with directly interpretable logical operators. Our experimental results show RL training with the proposed method converges significantly faster than other state of the art neuro symbolic methods in a TextWorld benchmark.


Black Box Dissector Towards Erasing based Hard Label Model Stealing Attack. Previous studies have verified that the functionality of black box models can be stolen with full probability outputs. However, under the more practical hard label setting, we observe that existing methods suffer from catastrophic performance degradation. We argue this is due to the lack of rich information in the probability prediction and the overfitting caused by hard labels. To this end, we propose a novel hard label model stealing method termed emph black box dissector , which consists of two erasing based modules. One is a CAM driven erasing strategy that is designed to increase the information capacity hidden in hard labels from the victim model. The other is a random erasing based self knowledge distillation module that utilizes soft labels from the substitute model to mitigate overfitting. Extensive experiments on four widely used datasets consistently demonstrate that our method outperforms state of the art methods, with an improvement of at most 8.27 . We also validate the effectiveness and practical potential of our method on real world APIs and defense methods. Furthermore, our method promotes other downstream tasks, emph i.e. , transfer adversarial attacks.


A Combination of Multi Objective Genetic Algorithm and Deep Learning for Music Harmony Generation. Automatic Music Generation (AMG) has become an interesting research topic for many scientists in artificial intelligence, who are also interested in the music industry. One of the main challenges in AMG is that there is no clear objective evaluation criterion that can measure the music grammar, structural rules, and audience satisfaction. Also, original music contains different elements that should work together, such as melody, harmony, and rhythm but in the most of previous works, AMG works only for one element (e.g., melody). Therefore, in this paper, we propose a Multi Objective Genetic Algorithm (MO GA) to generate polyphonic music pieces, considering grammar and listener satisfaction. In this method, we use three objective functions. The first objective function is the accuracy of the generated music piece, based on music theory and the other two objective functions are modeled scores provided by music experts and ordinary listeners. The scoring of experts and listeners separately are modeled using Bi directional Long Short Term Memory (Bi LSTM) neural networks. The proposed music generation system tries to maximize mentioned objective functions to generate a new piece of music, including melody and harmony. The results show that the proposed method can generate pleasant pieces with desired styles and lengths, along with harmonic sounds that follow the grammar.


Forecasting blood sugar levels in Diabetes with univariate algorithms. AI procedures joined with wearable gadgets can convey exact transient blood glucose level forecast models. Also, such models can learn customized glucose insulin elements dependent on the sensor information gathered by observing a few parts of the physiological condition and every day movement of a person. Up to this point, the predominant methodology for creating information driven forecast models was to gather "however much information as could be expected" to help doctors and patients ideally change treatment. The goal of this work was to examine the base information assortment, volume, and speed needed to accomplish exact individual driven diminutive term expectation models. We built up a progression of these models utilizing distinctive AI time arrangement guaging strategies that are appropriate for execution inside a wearable processor. We completed a broad aloof patient checking concentrate in genuine conditions to fabricate a strong informational collection. The examination included a subset of type 1 diabetic subjects wearing a glimmer glucose checking framework. We directed a relative quantitative assessment of the presentation of the created information driven expectation models and comparing AI methods. Our outcomes show that precise momentary forecast can be accomplished by just checking interstitial glucose information over a brief timeframe and utilizing a low examining recurrence. The models created can anticipate glucose levels inside a 15 minute skyline with a normal mistake as low as 15.43 mg dL utilizing just 24 memorable qualities gathered inside a time of 6 hours, and by expanding the inspecting recurrence to incorporate 72 qualities, the normal blunder is limited to 10.15 mg dL. Our forecast models are reasonable for execution inside a wearable gadget, requiring the base equipment necessities while simultaneously accomplishing high expectation precision.


Linearizability A Typo. Linearizability is the de facto consistency condition for concurrent objects, widely used in theory and practice. Loosely speaking, linearizability classifies concurrent executions as correct if operations on shared objects appear to take effect instantaneously during the operation execution time. This paper calls attention to a somewhat neglected aspect of linearizability restrictions on how pending invocations are handled, an issue that has become increasingly important for software running on systems with non volatile main memory. Interestingly, the original published definition of linearizability includes a typo (a symbol is missing a prime) that concerns exactly this issue. In this paper we point out the typo and provide an amendment to make the definition complete. We believe that pointing this typo out rigorously and proposing a fix is important and timely.


Detecting Viruses in Contact Networks with Unreliable Detectors. This paper develops and analyzes optimization models for rapid detection of viruses in large contact networks. In the model, a virus spreads in a stochastic manner over an undirected connected graph, under various assumptions on the spread dynamics. A decision maker must place a limited number of detectors on a subset of the nodes in the graph in order to rapidly detect infection of the nodes by the virus. The objective is to determine the placement of these detectors so as to either maximize the probability of detection within a given time period or minimize the expected time to detection. Previous work in this area assumed that the detectors are perfectly reliable. In this work, it is assumed that the detectors may produce false negative results. In computational studies, the sample average approximation method is applied to solving the problem using a mixed integer program and a greedy heuristic. The heuristic is shown to be highly efficient and to produce high quality solutions. In addition, it is shown that the false negative effect can sometimes be ignored, without significant loss of solution quality, in the original optimization formulation.


Regulated Multiset Rewriting Systems. Multiset rewriting systems provide a formalism particularly suitable for the description of biological systems. We present an extension of this formalism with additional controls on the derivations as a tool for reducing possible non deterministic behaviour by providing additional knowledge about the system. We introduce several regulation mechanisms and compare their generative power.


On factorizations into coprime parts. Let f(n) and g(n) be the number of unordered and ordered factorizations of n into integers larger than one. Let F(n) and G(n) have the additional restriction that the factors are coprime. We establish the asymptotic bounds for the sums of F(n) beta and G(n) beta up to x for all real beta and the asymptotic bounds for f(n) beta and g(n) beta for all negative beta .


Ten Simple Rules for Success with HPC, i.e. Responsibly BASHing that Linux Cluster. High performance computing (HPC) clusters are widely used in house at scientific and academic research institutions. For some users, the transition from running their analyses on a single workstation to running them on a complex, multi tenanted cluster, usually employing some degree of parallelism, can be challenging, if not bewildering, especially for users whose role is not predominantly computational in nature. On the other hand, there are more experienced users, who can benefit from pointers on how to get the best from their use of HPC. This Ten Simple Rules guide is aimed at helping you identify ways to improve your utilisation of HPC, avoiding common pitfalls that can negatively impact other users and will also help ease the load (pun intended) on your HPC sysadmin. It is intended to provide technical advice common to the use of HPC platforms such as LSF, Slurm, PBS Torque, SGE, LoadLeveler and YARN, the scheduler used with Hadoop Spark platform.


Perspectives on Machine Learning from Psychology's Reproducibility Crisis. In the early 2010s, a crisis of reproducibility rocked the field of psychology. Following a period of reflection, the field has responded with radical reform of its scientific practices. More recently, similar questions about the reproducibility of machine learning research have also come to the fore. In this short paper, we present select ideas from psychology's reformation, translating them into relevance for a machine learning audience.


Spatiotemporal Data Mining A Survey on Challenges and Open Problems. Spatiotemporal data mining (STDM) discovers useful patterns from the dynamic interplay between space and time. Several available surveys capture STDM advances and report a wealth of important progress in this field. However, STDM challenges and problems are not thoroughly discussed and presented in articles of their own. We attempt to fill this gap by providing a comprehensive literature survey on state of the art advances in STDM. We describe the challenging issues and their causes and open gaps of multiple STDM directions and aspects. Specifically, we investigate the challenging issues in regards to spatiotemporal relationships, interdisciplinarity, discretisation, and data characteristics. Moreover, we discuss the limitations in the literature and open research problems related to spatiotemporal data representations, modelling and visualisation, and comprehensiveness of approaches. We explain issues related to STDM tasks of classification, clustering, hotspot detection, association and pattern mining, outlier detection, visualisation, visual analytics, and computer vision tasks. We also highlight STDM issues related to multiple applications including crime and public safety, traffic and transportation, earth and environment monitoring, epidemiology, social media, and Internet of Things.


Mechanistic determination of tear film thinning via fitting simplified models to tear breakup. Purpose To determine whether evaporation, tangential flow, or a combination of the two cause tear film breakup in a variety of instances to estimate related breakup parameters that cannot be measured in breakup during subject trials and to validate our procedure against previous work. Methods Five ordinary differential equation models for tear film thinning were designed that model evaporation, osmosis, and various types of flow. Eight tear film breakup instances of five healthy subjects that were identified in fluorescence images in previous work were fit with these five models. The fitting procedure used a nonlinear least squares optimization that minimized the difference of the computed theoretical fluorescent intensity from the models and the experimental fluorescent intensity from the images. The optimization was conducted over the evaporation rate and up to three flow rate parameters. The smallest norm of the difference was determined to correspond to the model that best explained the tear film dynamics. Results All of the breakup instances were best fit by models with time dependent flow. Our optimal parameter values and thinning rate and fluid flow profiles compare well with previous partial differential equation model results in most instances. Conclusion Our fitting procedure suggests that the combination of the Marangoni effect and evaporation cause most of the breakup instances. Comparison with results from previous work suggests that the simplified models can capture the essential tear film dynamics in most cases, thereby validating this procedure as one that could be used on many other instances.


Another proof that mathsf MM implies Woodin's axiom ( ) . Let mathsf MM ( kappa) state that the forcing axiom mathsf MM can be instantiated only for stationary set preserving posets of size at most kappa . We give a detailed account of Asper o and Schindler's proof that mathsf MM ( kappa) there are class many Woodin cardinals implies Woodin's axiom ( ) if Diamond kappa holds and kappa aleph 2 . Our presentation takes advantage of the notion of consistency property specifically we rephrase Asper o and Schindler's forcing as a specific instantiation of the notion of consistency property'' used by Makkai, Keisler, Mansfield and others in the study of infinitary logics. We also reorganize the order of presentation of the various parts of the proof. Taken aside these variations, our account is quite close to the original proof of Asper o and Schindler.


Efficient Combinatorial Optimization for Word level Adversarial Textual Attack. Over the past few years, various word level textual attack approaches have been proposed to reveal the vulnerability of deep neural networks used in natural language processing. Typically, these approaches involve an important optimization step to determine which substitute to be used for each word in the original input. However, current research on this step is still rather limited, from the perspectives of both problem understanding and problem solving. In this paper, we address these issues by uncovering the theoretical properties of the problem and proposing an efficient local search algorithm (LS) to solve it. We establish the first provable approximation guarantee on solving the problem in general cases.Extensive experiments involving 5 NLP tasks, 8 datasets and 26 NLP models show that LS can largely reduce the number of queries usually by an order of magnitude to achieve high attack success rates. Further experiments show that the adversarial examples crafted by LS usually have higher quality, exhibit better transferability, and can bring more robustness improvement to victim models by adversarial training.


Club Stationary Reflection and the Special Aronszajn Tree Property. We prove that it is consistent that Club Stationary Reflection and the Special Aronszajn Tree Property simultaneously hold on omega 2 , thereby contributing to the study of the tension between compactness and incompactness in set theory. The poset which produces the final model follows the collapse of a weakly compact cardinal first with an iteration of club adding (with anticipation) and second with an iteration specializing Aronszajn trees. In the first part of the paper, we prove a general theorem about specializing Aronszajn trees after forcing with what we call mathcal F WC Strongly Proper posets. This type of poset, of which the Levy collapse is a degenerate example, uses systems of exact residue functions to create many strongly generic conditions. We prove a new result about stationary set preservation by quotients of this kind of poset as a corollary, we show that the original Laver Shelah model satisfies a strong stationary reflection principle, though it fails to satisfy the full Club Stationary Reflection. In the second part, we show that the composition of collapsing and club adding (with anticipation) is an mathcal F WC Strongly Proper poset. After proving a new result about Aronszajn tree preservation, we show how to obtain the final model.


Continuous Time Video Generation via Learning Motion Dynamics with Neural ODE. In order to perform unconditional video generation, we must learn the distribution of the real world videos. In an effort to synthesize high quality videos, various studies attempted to learn a mapping function between noise and videos, including recent efforts to separate motion distribution and appearance distribution. Previous methods, however, learn motion dynamics in discretized, fixed interval timesteps, which is contrary to the continuous nature of motion of a physical body. In this paper, we propose a novel video generation approach that learns separate distributions for motion and appearance, the former modeled by neural ODE to learn natural motion dynamics. Specifically, we employ a two stage approach where the first stage converts a noise vector to a sequence of keypoints in arbitrary frame rates, and the second stage synthesizes videos based on the given keypoints sequence and the appearance noise vector. Our model not only quantitatively outperforms recent baselines for video generation, but also demonstrates versatile functionality such as dynamic frame rate manipulation and motion transfer between two datasets, thus opening new doors to diverse video generation applications.


CDA a Cost Efficient Content based Multilingual Web Document Aligner. We introduce a Content based Document Alignment approach (CDA), an efficient method to align multilingual web documents based on content in creating parallel training data for machine translation (MT) systems operating at the industrial level. CDA works in two steps (i) projecting documents of a web domain to a shared multilingual space then (ii) aligning them based on the similarity of their representations in such space. We leverage lexical translation models to build vector representations using TF IDF. CDA achieves performance comparable with state of the art systems in the WMT 16 Bilingual Document Alignment Shared Task benchmark while operating in multilingual space. Besides, we created two web scale datasets to examine the robustness of CDA in an industrial setting involving up to 28 languages and millions of documents. The experiments show that CDA is robust, cost effective, and is significantly superior in (i) processing large and noisy web data and (ii) scaling to new and low resourced languages.


Towards Indirect Top Down Road Transport Emissions Estimation. Road transportation is one of the largest sectors of greenhouse gas (GHG) emissions affecting climate change. Tackling climate change as a global community will require new capabilities to measure and inventory road transport emissions. However, the large scale and distributed nature of vehicle emissions make this sector especially challenging for existing inventory methods. In this work, we develop machine learning models that use satellite imagery to perform indirect top down estimation of road transport emissions. Our initial experiments focus on the United States, where a bottom up inventory was available for training our models. We achieved a mean absolute error (MAE) of 39.5 kg CO 2 of annual road transport emissions, calculated on a pixel by pixel (100 m 2 ) basis in Sentinel 2 imagery. We also discuss key model assumptions and challenges that need to be addressed to develop models capable of generalizing to global geography. We believe this work is the first published approach for automated indirect top down estimation of road transport sector emissions using visual imagery and represents a critical step towards scalable, global, near real time road transportation emissions inventories that are measured both independently and objectively.


Absolute retracts for finite distributive lattices and slim semimodular lattices. We describe the absolute retracts for the following classes of finite lattices (1) slim semimodular lattices, (2) finite distributive lattices, and for each positive integer n , (3) at most n dimensional finite distributive lattices. Although the singleton lattice is the only absolute retract for the first class, this result has paved the way to some other classes. For the second class, we prove that the absolute retracts are exactly the finite boolean lattices this generalizes a 1979 result of J. Schmid. For the third class, the absolute retracts are the finite boolean lattices of dimension at most n and the direct products of n nontrivial finite chains. Also, we point out that in each of these classes, the algebraically closed lattices and the strongly algebraically closed lattices are the same as the absolute retracts. Slim (and necessarily planar) semimodular lattices were introduced by G. Gr "atzer and E. Knapp in 2007, and they have been intensively studied since then. Algebraically closed and strongly algebraically closed lattices have been investigated by J. Schmid and, in several papers, by A. Molkhasi.


Type Based Termination for Futures. In sequential functional languages, sized types enable termination checking of programs with complex patterns of recursion in the presence of mixed inductive coinductive types. In this paper, we adapt sized types and their metatheory to the concurrent setting. We extend the semi axiomatic sequent calculus, a subsuming paradigm for futures based functional concurrency, and its underlying operational semantics with recursion and arithmetic refinements. The latter enables a new and highly general sized type scheme we call sized type refinements. As a widely applicable technical device, we type recursive programs with infinitely deep typing derivations that unfold all recursive calls. Then, we observe that certain such derivations can be made infinitely wide but finitely deep. The resulting trees serve as the induction target of our strong normalization result, which we develop via a novel logical relations argument.


Counting Tripods on the Torus. Motivated by the problem of counting finite BPS webs, we count certain immersed metric graphs, tripods, on the flat torus. Classical Euclidean geometry turns this into a lattice point counting problem in mathbb C 2 , and we give an asymptotic counting result using lattice point counting techniques.


On Gorenstein Fano threefolds. We classify the non toric, Q factorial, Gorenstein, log terminal Fano threefolds of Picard number one that admit an effective action of a two dimensional algebraic torus.


TL SDD A Transfer Learning Based Method for Surface Defect Detection with Few Samples. Surface defect detection plays an increasingly important role in manufacturing industry to guarantee the product quality. Many deep learning methods have been widely used in surface defect detection tasks, and have been proven to perform well in defects classification and location. However, deep learning based detection methods often require plenty of data for training, which fail to apply to the real industrial scenarios since the distribution of defect categories is often imbalanced. In other words, common defect classes have many samples but rare defect classes have extremely few samples, and it is difficult for these methods to well detect rare defect classes. To solve the imbalanced distribution problem, in this paper we propose TL SDD a novel Transfer Learning based method for Surface Defect Detection. First, we adopt a two phase training scheme to transfer the knowledge from common defect classes to rare defect classes. Second, we propose a novel Metric based Surface Defect Detection (M SDD) model. We design three modules for this model (1) feature extraction module containing feature fusion which combines high level semantic information with low level structural information. (2) feature reweighting module transforming examples to a reweighting vector that indicates the importance of features. (3) distance metric module learning a metric space in which defects are classified by computing distances to representations of each category. Finally, we validate the performance of our proposed method on a real dataset including surface defects of aluminum profiles. Compared to the baseline methods, the performance of our proposed method has improved by up to 11.98 for rare defect classes.


CInC Flow Characterizable Invertible 3x3 Convolution. Normalizing flows are an essential alternative to GANs for generative modelling, which can be optimized directly on the maximum likelihood of the dataset. They also allow computation of the exact latent vector corresponding to an image since they are composed of invertible transformations. However, the requirement of invertibility of the transformation prevents standard and expressive neural network models such as CNNs from being directly used. Emergent convolutions were proposed to construct an invertible 3 times 3 CNN layer using a pair of masked CNN layers, making them inefficient. We study conditions such that 3 times 3 CNNs are invertible, allowing them to construct expressive normalizing flows. We derive necessary and sufficient conditions on a padded CNN for it to be invertible. Our conditions for invertibility are simple, can easily be maintained during the training process. Since we require only a single CNN layer for every effective invertible CNN layer, our approach is more efficient than emerging convolutions. We also proposed a coupling method, Quad coupling. We benchmark our approach and show similar performance results to emergent convolutions while improving the model's efficiency.


Bounded generation for congruence subgroups of rm Sp 4(R) . This paper describes a bounded generation result concerning the minimal natural number K such that for Q(C 2,2R) A varepsilon phi (2x)A 1 x in R,A in rm Sp 4(R), phi in C 2 , one has N C 2,2R X 1 cdots X K forall 1 leq i leq K X i in Q(C 2,2R) for rings of algebraic integers R and the principal congruence subgroup N C 2,2R in rm Sp 4(R). This gives an explicit version of an abstract bounded generation result of a similar type as presented by Morris. Furthermore, the result presented does not depend on several number theoretic quantities unlike Morris' result. Using this bounded generation result, we further give explicit bounds for the strong boundedness of rm Sp 4(R) for certain examples of rings R, thereby giving explicit versions of results in an earlier paper. We further give a classification of normally generating subsets of rm Sp 4(R) for R a ring of algebraic integers.


Deep Spherical Manifold Gaussian Kernel for Unsupervised Domain Adaptation. Unsupervised Domain adaptation is an effective method in addressing the domain shift issue when transferring knowledge from an existing richly labeled domain to a new domain. Existing manifold based methods either are based on traditional models or largely rely on Grassmannian manifold via minimizing differences of single covariance matrices of two domains. In addition, existing pseudo labeling algorithms inadequately consider the quality of pseudo labels in aligning the conditional distribution between two domains. In this work, a deep spherical manifold Gaussian kernel (DSGK) framework is proposed to map the source and target subspaces into a spherical manifold and reduce the discrepancy between them by embedding both extracted features and a Gaussian kernel. To align the conditional distributions, we further develop an easy to hard pseudo label refinement process to improve the quality of the pseudo labels and then reduce categorical spherical manifold Gaussian kernel geodesic loss. Extensive experimental results show that DSGK outperforms state of the art methods, especially on challenging cross domain learning tasks.


When standard network measures fail to rank journals A theoretical and empirical analysis. Journal rankings are widely used and are often based on citation data in combination with a network perspective. We argue that some of these network based rankings can produce misleading results. From a theoretical point of view, we show that the standard network modelling approach of citation data at the journal level (i.e., the projection of paper citations onto journals) introduces fictitious relations among journals. To overcome this problem, we propose a citation path perspective, and empirically show that rankings based on the network and the citation path perspective are very different. Based on our theoretical and empirical analysis, we highlight the limitations of standard network metrics, and propose a method to overcome these limitations and compute journal rankings.


Localization of extriangulated categories. In this article, we show that the localization of an extriangulated category by a multiplicative system satisfying mild assumptions can be equipped with a natural, universal structure of an extriangulated category. This construction unifies the Serre quotient of abelian categories and the Verdier quotient of triangulated categories. Indeed we give such a construction for a bit wider class of morphisms, so that it covers several other localizations appeared in the literature, such as Rump's localization of exact categories by biresolving subcategories, localizations of extriangulated categories by means of Hovey twin cotorsion pairs, and the localization of exact categories by two sided admissibly percolating subcategories.


Fast Payoff Matrix Sparsification Techniques for Structured Extensive Form Games. The practical scalability of many optimization algorithms for large extensive form games is often limited by the games' huge payoff matrices. To ameliorate the issue, Zhang and Sandholm (2020) recently proposed a sparsification technique that factorizes the payoff matrix mathbf A into a sparser object mathbf A hat mathbf A mathbf U mathbf V top , where the total combined number of nonzeros of hat mathbf A , mathbf U , and mathbf V is significantly smaller. Such a factorization can be used in place of the original payoff matrix in many optimization algorithm, such as interior point and second order methods, thus increasing the size of games that can be handled. Their technique significantly sparsifies poker (end)games, standard benchmarks used in computational game theory, AI, and more broadly. We show that the existence of extremely sparse factorizations in poker games can be tied to their particular emph Kronecker product structure. We clarify how such structure arises and introduce the connection between that structure and sparsification. By leveraging such structure, we give two ways of computing strong sparsifications of poker games (as well as any other game with a similar structure) that are i) orders of magnitude faster to compute, ii) more numerically stable, and iii) produce a dramatically smaller number of nonzeros than the prior technique. Our techniques enable for the first time effective computation of high precision Nash equilibria and strategies subject to constraints on the amount of allowed randomization. Furthermore, they significantly speed up parallel first order game solving algorithms we show state of the art speed on a GPU.


Monoid extensions and the Grothendieck construction. In category theory circles it is well known that the Schreier theory of group extensions can be understood in terms of the Grothendieck construction on indexed categories. However, it is seldom discussed how this relates to extensions of monoids. We provide an introduction to the generalised Grothendieck construction and apply it to recover classifications of certain classes of monoid extensions (including Schreier and weakly Schreier extensions in particular).


Real Time Predictive Maintenance using Autoencoder Reconstruction and Anomaly Detection. Rotary machine breakdown detection systems are outdated and dependent upon routine testing to discover faults. This is costly and often reactive in nature. Real time monitoring offers a solution for detecting faults without the need for manual observation. However, manual interpretation for threshold anomaly detection is often subjective and varies between industrial experts. This approach is ridged and prone to a large number of false positives. To address this issue, we propose a Machine Learning (ML) approach to model normal working operation and detect anomalies. The approach extracts key features from signals representing known normal operation to model machine behaviour and automatically identify anomalies. The ML learns generalisations and generates thresholds based on fault severity. This provides engineers with a traffic light system were green is normal behaviour, amber is worrying and red signifies a machine fault. This scale allows engineers to undertake early intervention measures at the appropriate time. The approach is evaluated on windowed real machine sensor data to observe normal and abnormal behaviour. The results demonstrate that it is possible to detect anomalies within the amber range and raise alarms before machine failure.


Comparison of Poisson structures on moduli spaces. Let X be a complex irreducible smooth projective curve, and let mathbb L be an algebraic line bundle on X with a nonzero section sigma 0 . Let mathcal M denote the moduli space of stable Hitchin pairs (E, , theta) , where E is an algebraic vector bundle on X of fixed rank r and degree delta , and theta , in , H 0(X, , End(E) otimes K X otimes mathbb L ) . Associating to every stable Hitchin pair its spectral data, an isomorphism of mathcal M with a moduli space mathcal P of stable sheaves of pure dimension one on the total space of K X otimes mathbb L is obtained. Both the moduli spaces mathcal P and mathcal M are equipped with algebraic Poisson structures, which are constructed using sigma 0 . Here we prove that the above isomorphism between mathcal P and mathcal M preserves the Poisson structures.


Iterative Training Finding Binary Weight Deep Neural Networks with Layer Binarization. In low latency or mobile applications, lower computation complexity, lower memory footprint and better energy efficiency are desired. Many prior works address this need by removing redundant parameters. Parameter quantization replaces floating point arithmetic with lower precision fixed point arithmetic, further reducing complexity. Typical training of quantized weight neural networks starts from fully quantized weights. Quantization creates random noise. As a way to compensate for this noise, during training, we propose to quantize some weights while keeping others in floating point precision. A deep neural network has many layers. To arrive at a fully quantized weight network, we start from one quantized layer and then quantize more and more layers. We show that the order of layer quantization affects accuracies. Order count is large for deep neural networks. A sensitivity pre training is proposed to guide the layer quantization order. Recent work in weight binarization replaces weight input matrix multiplication with additions. We apply the proposed iterative training to weight binarization. Our experiments cover fully connected and convolutional networks on MNIST, CIFAR 10 and ImageNet datasets. We show empirically that, starting from partial binary weights instead of from fully binary ones, training reaches fully binary weight networks with better accuracies for larger and deeper networks. Layer binarization in the forward order results in better accuracies. Guided layer binarization can further improve that. The improvements come at a cost of longer training time.


On canonical radial Kaehler metrics. We prove that a radial Kaehler metric g is Kaehler Einstein if and only if one of the following conditions is satisfied 1. g is extremal and it is associated to a Kaehler Ricci soliton 2. two different generalized scalar curvatures of g are constant 3. g is extremal (not cscK) and one of its generalized scalar curvature is constant.


Large time behavior of finite state mean field systems with multi classes. We study in this paper the large time asymptotics of the empirical vector associated with a family of finite state mean field systems with multi classes. The empirical vector is composed of local empirical measures characterizing the different classes within the system. As the number of particles in the system goes to infinity, the empirical vector process converges towards the solution to a McKean Vlasov system. First, we investigate the large deviations principles of the invariant distribution from the limiting McKean Vlasov system. Then, we examine the metastable phenomena arising at a large scale and large time. Finally, we estimate the rate of convergence of the empirical vector process to its invariant measure. Given the local homogeneity in the system, our results are established in a product space.


Understanding and Creating Art with AI Review and Outlook. Technologies related to artificial intelligence (AI) have a strong impact on the changes of research and creative practices in visual arts. The growing number of research initiatives and creative applications that emerge in the intersection of AI and art, motivates us to examine and discuss the creative and explorative potentials of AI technologies in the context of art. This paper provides an integrated review of two facets of AI and art 1) AI is used for art analysis and employed on digitized artwork collections 2) AI is used for creative purposes and generating novel artworks. In the context of AI related research for art understanding, we present a comprehensive overview of artwork datasets and recent works that address a variety of tasks such as classification, object detection, similarity retrieval, multimodal representations, computational aesthetics, etc. In relation to the role of AI in creating art, we address various practical and theoretical aspects of AI Art and consolidate related works that deal with those topics in detail. Finally, we provide a concise outlook on the future progression and potential impact of AI technologies on our understanding and creation of art.


Distributed Resilient Submodular Action Selection in Adversarial Environments. In this letter, we consider a distributed submodular maximization problem for multi robot systems when attacked by adversaries. One of the major challenges for multi robot systems is to increase resilience against failures or attacks. This is particularly important for distributed systems under attack as there is no central point of command that can detect, mitigate, and recover from attacks. Instead, a distributed multi robot system must coordinate effectively to overcome adversarial attacks. In this work, our distributed submodular action selection problem models a broad set of scenarios where each robot in a multi robot system has multiple action selections that may fulfill a global objective, such as exploration or target tracking. To increase resilience in this context, we propose a fully distributed algorithm to guide each robot's action selection when the system is attacked. The proposed algorithm guarantees performance in a worst case scenario where up to a portion of the robots malfunction due to attacks. Importantly, the proposed algorithm is also consistent, as it is shown to converge to the same solution as a centralized method. Finally, a distributed resilient multi robot exploration problem is presented to confirm the performance of the proposed algorithm.


Classification of simple Harish Chandra modules for basic classical map superalgebras. We obtain a classification of simple modules with finite weight multiplicities over basic classical map superalgebras. Any such module is parabolic induced from a simple cuspidal bounded module over a cuspidal map superalgebra. Further on, any simple cuspidal bounded module is isomorphic to an evaluation module. As an application, we obtain a classification of all simple Harish Chandra modules for basic classical loop superalgebras. Extending these results to affine Kac Moody Lie superalgebras obtained by adding the degree derivation, we construct a family of bounded simple modules of level zero, and conjecture that all bounded simple cuspidal modules belong to this family. Finally, we show that for affine Kac Moody Lie superalgebras of type I the Kac induction functor reduces the classification of all simple bounded modules to the classification of the same class of modules over the even part, whose classification is claimed by Dimitrov and Grantcharov.


Identities for a parametric Weyl algebra over a ring. In 2013 Benkart, Lopes and Ondrus introduced and studied in a series of papers the infinite dimensional unital associative algebra A h generated by elements x,y, which satisfy the relation yx xy h for some 0 neq h in FF x . We generalize this construction to A h( B) by working over the fixed FF algebra B instead of FF . We describe the polynomial identities for A h( B) over the infinite field FF in case h in B x satisfies certain restrictions.


Agile Methods in Higher Education Adapting and Using eduScrum with Real World Projects. This Innovative Practice Full Paper presents our learnings of the process to perform a Master of Science class with eduScrum integrating real world problems as projects. We prepared, performed, and evaluated an agile educational concept for the new Master of Science program Digital Transformation organized and provided by the department of business computing at the University of Applied Sciences and Arts Hochschule Hannover in Germany. The course deals with innovative methodologies of agile project management and is attended by 25 students. We performed the class due the summer term in 2019 and 2020 as a teaching pair. The eduScrum method has been used in different educational contexts, including higher education. During the approach preparation, we decided to use challenges, problems, or questions from the industry. Thus, we acquired four companies and prepared in coordination with them dedicated project descriptions. Each project description was refined in the form of a backlog (list of requirements). We divided the class into four eduScrum teams, one team for each project. The subdivision of the class was done randomly. Since we wanted to integrate realistic projects into industry partners' implementation, we decided to adapt the eduScrum approach. The eduScrum teams were challenged with different projects, e.g., analyzing a dedicated phenomenon in a real project or creating a theoretical model for a company's new project management approach. We present our experiences of the whole process to prepare, perform and evaluate an agile educational approach combined with projects from practice. We found, that the students value the agile method using real world problems. The paper contributes to the distribution of methods for higher education teaching in the classroom and distance learning.


Revisiting the representation theorem of finite distributive lattices with principal congruences. A classical result of R. ,P. Dilworth states that every finite distributive lattice D can be represented as the congruence lattice of a finite lattice L . A sharper form was published in G. Gr "atzer and E. ,T. Schmidt in 1962, adding the requirement that all congruences in L be principal. Another variant, published in 1998 by the authors and E. ,T. Schmidt, constructs a planar semimodular lattice L . In this paper, we merge these two results we construct L as a planar semimodular lattice in which all congruences are principal. This paper relies on the techniques developed by the authors and E. ,T. Schmidt in the 1998 paper.


Multiscale mechanical model based on patient specific geometry application to early keratoconus development. Keratoconus is a pathology of the cornea associated with a tissue thinning and a weakening of its mechanical properties. However, it remains elusive which aspect is the leading cause of the disease. To investigate this question, we combined a multiscale model with a patient geometry in order to simulate the mechanical response of healthy and pathological corneas under intraocular pressure. The constitutive behavior of the cornea is described through an energy function which takes into account the isotropic matrix of the cornea, the geometric structure of collagen lamellae and the quasi incompressibility of the tissue. A micro sphere description is implemented to take into account the typical features of the collagen lamellae as obtained experimentally, namely their orientation, their stiffness and their dispersion, as well as the their unfolding stretch, at which they start to provide a significant force. A set of reference parameters is obtained to fit experimental inflation data of the literature. We show that the most sensitive parameter is the unfolding stretch, as a small variation of this parameter induces a major change in the corneal apex displacement. The keratoconus case is then studied by separating the impact of the geometry and the one of the mechanics. We computed the evolution of the SimK (a clinical indicator of cornea curvature) and elevation maps we were able to reproduce the reported changes of SimK with pressure only by a mechanical weakening, and not by a change in geomtry. More specifically, the weakening has to target the lamellae and not the matrix. The mechanical weakening leads to elevations close to early stage keratoconus, but our model lacks the remodeling component to couple the change in mechanics with changes in geometry.


Swivel Hardening WebAssembly against Spectre. We describe Swivel, a new compiler framework for hardening WebAssembly (Wasm) against Spectre attacks. Outside the browser, Wasm has become a popular lightweight, in process sandbox and is, for example, used in production to isolate different clients on edge clouds and function as a service platforms. Unfortunately, Spectre attacks can bypass Wasm's isolation guarantees. Swivel hardens Wasm against this class of attacks by ensuring that potentially malicious code can neither use Spectre attacks to break out of the Wasm sandbox nor coerce victim code another Wasm client or the embedding process to leak secret data. We describe two Swivel designs, a software only approach that can be used on existing CPUs, and a hardware assisted approach that uses extension available in Intel 11th generation CPUs. For both, we evaluate a randomized approach that mitigates Spectre and a deterministic approach that eliminates Spectre altogether. Our randomized implementations impose under 10.3 overhead on the Wasm compatible subset of SPEC 2006, while our deterministic implementations impose overheads between 3.3 and 240.2 . Though high on some benchmarks, Swivel's overhead is still between 9x and 36.3x smaller than existing defenses that rely on pipeline fences.


Rigidity of saddle loops. A saddle loop is a germ of a holomorphic foliation near a homoclinic saddle connection. We prove that they are classied by their Poincar 'e rst return map. We also prove that they are formally rigid when the Poincar 'e map is multivalued. Finally, we provide a list of all analytic classes of Liouville integrable saddle loops.


Deep Learning Models in Detection of Dietary Supplement Adverse Event Signals from Twitter. Objective The objective of this study is to develop a deep learning pipeline to detect signals on dietary supplement related adverse events (DS AEs) from Twitter. Material and Methods We obtained 247,807 tweets ranging from 2012 to 2018 that mentioned both DS and AE. We annotated biomedical entities and relations on 2,000 randomly selected tweets. For the concept extraction task, we compared the performance of traditional word embeddings with SVM, CRF and LSTM CRF classifiers to BERT models. For the relation extraction task, we compared GloVe vectors with CNN classifiers to BERT models. We chose the best performing models in each task to assemble an end to end deep learning pipeline to detect DS AE signals and compared the results to the known DS AEs from a DS knowledge base (i.e., iDISK). Results In both tasks, the BERT based models outperformed traditional word embeddings. The best performing concept extraction model is the BioBERT model that can identify supplement, symptom, and body organ entities with F1 scores of 0.8646, 0.8497, and 0.7104, respectively. The best performing relation extraction model is the BERT model that can identify purpose and AE relations with F1 scores of 0.8335 and 0.7538, respectively. The end to end pipeline was able to extract DS indication and DS AEs with an F1 score of 0.7459 and 0,7414, respectively. Comparing to the iDISK, we could find both known and novel DS AEs. Conclusion We have demonstrated the feasibility of detecting DS AE signals from Twitter with a BioBERT based deep learning pipeline.


Revisit of macroscopic dynamics for non equilibrium chemical reactions from a Hamiltonian viewpoint. Most biochemical reactions in living cells are open systems interacting with environment through chemostats to exchange both energy and materials. At a mesoscopic scale, the number of each species in those biochemical reactions can be modeled by a random time changed Poisson processes. To characterize macroscopic behaviors in the large volume limit, the law of large numbers in the path space determines a mean field limit nonlinear reaction rate equation describing the dynamics of the concentration of species, while the WKB expansion for the chemical master equation yields a Hamilton Jacobi equation (HJE) and the Lagrangian gives the good rate function in the large deviation principle. We decompose a general macroscopic reaction rate equation into a conservative part and a dissipative part in terms of the stationary solution to HJE. This stationary solution is used to determine the energy landscape and thermodynamics for general chemical reactions, which particularly maintains a positive entropy production rate at a non equilibrium steady state. The associated energy dissipation law is proved together with the passage from the mesoscopic to macroscopic one. Furthermore, we use a reversible Hamiltonian to study a class of non equilibrium enzyme reactions, which identifies a new concept of balance within the same reaction vector due to flux grouping degeneracy. This macroscopic reversibility, brought by the reversibility of the chemical reaction jumping process, gives an Onsager type strong gradient flow. The reversible Hamiltonian also yields a time reversal symmetry for the corresponding Lagrangian. Thus a modified time reversed least action path serves as the transition paths with associated path affinities and energy barriers.


A preliminary study on evaluating Consultation Notes with Post Editing. Automatic summarisation has the potential to aid physicians in streamlining clerical tasks such as note taking. But it is notoriously difficult to evaluate these systems and demonstrate that they are safe to be used in a clinical setting. To circumvent this issue, we propose a semi automatic approach whereby physicians post edit generated notes before submitting them. We conduct a preliminary study on the time saving of automatically generated consultation notes with post editing. Our evaluators are asked to listen to mock consultations and to post edit three generated notes. We time this and find that it is faster than writing the note from scratch. We present insights and lessons learnt from this experiment.


Internal tide generation from isolated seamounts and continental shelves. We model linear, inviscid non hydrostatic internal tides generated by the interaction of a barotropic tide with variable topography in two dimensions. We first derive an asymptotic solution for the nonuniform barotropic flow over the topography that serves as forcing for the baroclinic equations. The resulting internal tide generation problem is reformulated as a Coupled Mode System (CMS) by means of a series decomposition of the baroclinic stream function in terms of vertical basis functions. We solve this CMS numerically and also provide a method for estimating the sea surface signature of internal tides. We consider several seamounts and shelf profiles and perform calculations for a wide range of (topographic) heights and slopes. For subcritical topographies, the energy flux as a function of height exhibits local maxima, separated by cases of weakly or even non radiating topographies. For supercritical topographies, the energy flux generally increases with height and criticality. Our calculations agree with the Weak Topography Approximation only for very small heights. Perhaps more surprisingly, they agree with the Knife Edge model only for moderately supercritical topographies. We also compare the effect of the adjusted barotropic tide on the energy flux and the local properties of the baroclinic field with other semi analytical methods based on a uniform barotropic tide. We observe significant differences in the flow field near the topographies only.


E Commerce in Turkey and SAP Integrated E Commerce System. E commerce is a kind of shopping by use of the internet. E commerce, very different from the usual shopping concept, is compatible with today's economic dynamics. E commerce is becoming an indispensable method with the increase of internet usage. With the use of E commerce, there are also a number of advantages for companies. On the other hand, SAP is a pioneer and leader in the company resource planning software sector. SAP is very important for large scale companies. They manage all their processes on SAP and its integration is very important with other related software. In this article, we give brief information on some important aspects of e commerce and propose a solution for ERP integration of an e commerce system.


Extending LIME for Business Process Automation. AI business process applications automate high stakes business decisions where there is an increasing demand to justify or explain the rationale behind algorithmic decisions. Business process applications have ordering or constraints on tasks and feature values that cause lightweight, model agnostic, existing explanation methods like LIME to fail. In response, we propose a local explanation framework extending LIME for explaining AI business process applications. Empirical evaluation of our extension underscores the advantage of our approach in the business process setting.


Sobolev embeddings for Fractional Haj l asz Sobolev spaces in the setting of rearrangement invariant spaces. We obtain symmetrization inequalities in the context of Fractional Hajlasz Sobolev spaces in the setting of rearrangement invariant spaces and prove that for a large class of measures our symmetrization inequalities are equivalent to the lower bound of the measure.


Retrieve Fast, Rerank Smart Cooperative and Joint Approaches for Improved Cross Modal Retrieval. Current state of the art approaches to cross modal retrieval process text and visual input jointly, relying on Transformer based architectures with cross attention mechanisms that attend over all words and objects in an image. While offering unmatched retrieval performance, such models 1) are typically pretrained from scratch and thus less scalable, 2) suffer from huge retrieval latency and inefficiency issues, which makes them impractical in realistic applications. To address these crucial gaps towards both improved and efficient cross modal retrieval, we propose a novel fine tuning framework that turns any pretrained text image multi modal model into an efficient retrieval model. The framework is based on a cooperative retrieve and rerank approach which combines 1) twin networks (i.e., a bi encoder) to separately encode all items of a corpus, enabling efficient initial retrieval, and 2) a cross encoder component for a more nuanced (i.e., smarter) ranking of the retrieved small set of items. We also propose to jointly fine tune the two components with shared weights, yielding a more parameter efficient model. Our experiments on a series of standard cross modal retrieval benchmarks in monolingual, multilingual, and zero shot setups, demonstrate improved accuracy and huge efficiency benefits over the state of the art cross encoders.


Full color nanorouter for high resolution imaging. Pixel scaling effects have been a major issue for the development of high resolution color image sensors due to the reduced photoelectric signal and the color crosstalk. Various structural color techniques have been proposed and demonstrated the large freedom in color manipulation by the structure design. However, the optical efficiency and the color distortion limit the practical applications due to its intrinsic filtering mechanism. Instead, on chip full color routing is quite charming for improving the signal to noise ratio. In this paper, a single layer quick response code like nanorouter is proposed for the full color light routing in a pixel level of image sensors. It shows much higher routing efficiency than various planar lens schemes for signal wavelength focusing. Moreover, over 60 signal enhancement with robust polarization insensitivity is obtained in all three primary color bands with a same nanorouter by a multi objective optimization method. Negligible color distortion is observed from the reconstructed color image. Such a simple nanorouter structure is promising for the development of image sensor, photovoltaics and display.


On some geometric properties of sets related to the d Hausorff content and new Whitney type coverings. Let S subset mathbb R n be a d thick set for some d in (0,n) , i.e., there exists a number lambda in (0,1 such that the d Hausdorff content mathcal H d infty (Q cap S) geq lambda (l(Q)) d for all cubes Q subset mathbb R n centered at x in S with side lengths l(Q) in (0,1 . Given a cube Q Q(x,l) with x in mathbb R n and l in (0,1 , we show that if the d Hausdorff content mathcal H d infty (Q cap S) is sufficiently small then the cube Q is S porous, i.e., there is a cube Q' subset Q setminus S with l(Q') approx l(Q) . In the case of an arbitrary set E subset mathbb R n , given a cube Q Q(x,l) with x in mathbb R n and l in (0,1 we show that if the d Hausdorff content mathcal H d infty (Q cap E) is sufficiently small then the cube Q is S hollow i.e., there exists a Borel set U subset Q setminus E with mathcal H n (U) approx mathcal H n (Q) . Furthermore, we introduce the special (d, lambda) thick distance rho E,d, lambda and show that U can be chosen to lie in the complement of the delta l neighborhood of E for some delta delta(n,d, lambda) in (0,1) . As an application we obtain new Whitney type coverings. The sharpness of the results is illustrated by several examples.


Summarizing and Analyzing the Privacy Preserving Techniques in Bitcoin and other Cryptocurrencies. Bitcoin and many other similar Cryptocurrencies have been in existence for over a decade, prominently focusing on decentralized, pseudo anonymous ledger based transactions. Many protocol improvements and changes have resulted in new variants of Cryptocurrencies that are known for their peculiar characteristics. For instance, Storjcoin is a Proof of Storage based Cryptocurrency that incentivizes its peers based on the amount of storage owned by them. Cryptocurrencies like Monero strive for user privacy by using privacy centric cryptographic algorithms. While Cryptocurrencies strive to maintain peer transparency by making the transactions and the entire ledger public, user privacy is compromised at times. Monero and many other privacy centric Cryptocurrencies have significantly improved from the original Bitcoin protocol after several problems were found in the protocol. Most of these deficiencies were related to the privacy of users. Even though Bitcoin claims to have pseudo anonymous user identities, many attacks have managed to successfully de anonymize users. In this paper, we present some well known attacks and analysis techniques that have compromised the privacy of Bitcoin and many other similar Cryptocurrencies. We also analyze and study different privacy preserving algorithms and the problems these algorithms manage to solve. Lastly, we touch upon the ethics, impact, legality, and acceptance of imposing these privacy algorithms.


Multicomponent Heat Bath Configuration Interaction with the Perturbative Correction for the Calculation of Protonic Excited States. In this study, we extend the multicomponent heat bath configuration interaction (HCI) method to excited states. Previous multicomponent HCI studies have been performed using only the variational stage of the HCI algorithm as they have largely focused on the calculation of protonic densities. Because this study focuses on energetic quantities, a second order perturbative correction after the variational stage is essential. Therefore, this study implements the second order Epstein Nesbet correction to the variational stage of multicomponent HCI for the first time. Additionally, this study introduces a new procedure for calculating reference excitation energies for multicomponent methods using the Fourier grid Hamiltonian (FGH) method, which should allow the one particle electronic basis set errors to be better isolated from errors arising from an incomplete description of electron proton correlation. The excited state multicomponent HCI method is benchmarked by computing protonic excitations of the HCN and FHF molecules and is shown to be of similar accuracy to previous excited state multicomponent methods such as multicomponent time dependent density functional theory and equation of motion coupled cluster theory relative to the new FGH reference values.


Functional Collection Programming with Semi Ring Dictionaries. This paper introduces semi ring dictionaries, a powerful class of compositional and purely functional collections that subsume other collection types such as sets, multisets, arrays, vectors, and matrices. We developed SDQL, a statically typed language that can express relational algebra with aggregations, linear algebra, and functional collections over data such as relations and matrices using semi ring dictionaries. Furthermore, thanks to the algebraic structure behind these dictionaries, SDQL unifies a wide range of optimizations commonly used in databases (DB) and linear algebra (LA). As a result, SDQL enables efficient processing of hybrid DB and LA workloads, by putting together optimizations that are otherwise confined to either DB systems or LA frameworks. We show experimentally that a handful of DB and LA workloads can take advantage of the SDQL language and optimizations. SDQL can be competitive with or outperforms a host of systems that are state of the art in their own domain in memory DB systems Typer and Tectorwise for (flat, not nested) relational data SciPy for LA workloads sparse tensor compiler taco the Trance nested relational engine and the in database machine learning engines LMFAO and Morpheus for hybrid DB LA workloads over relational data.


SCIMAT Science and Mathematics Dataset. In this work, we announce a comprehensive well curated and opensource dataset with millions of samples for pre college and college level problems in mathematicsand science. A preliminary set of results using transformer architecture with character to character encoding is shown. The dataset identifies some challenging problem and invites research on better architecture search


A Geometric Classification of the Holomorphic Vertex Operator Algebras of Central Charge 24. We associate with a generalised deep hole of the Leech lattice vertex operator algebra a generalised hole diagram. We show that this Dynkin diagram determines the generalised deep hole up to conjugacy and that there are exactly 70 such diagrams. In an earlier work we proved a bijection between the generalised deep holes and the strongly rational, holomorphic vertex operator algebras of central charge 24 with non trivial weight 1 space. Hence, we obtain a new, geometric classification of these vertex operator algebras.


Generalized torsion for knots with arbitrarily high genus. In a group, a non trivial element is called a generalized torsion element if some non empty finite product of its conjugates equals to the identity. We say that a knot has generalized torsion if its knot group admits such an element. For a (2, 2q 1) torus knot K, we demonstrate that there are infinitely many unknots c such that p twisting K about c yields a twist family, which consists of hyperbolic knots with generalized torsion whenever p 3. This gives a new infinite class of hyperbolic knots having generalized torsion. In particular, each class contains knots with arbitrarily high genus. We also show that some twisted torus knots, including the ( 2, 3, 7) pretzel knot, have generalized torsion. Since generalized torsion is an obstruction for having bi order, these knots have non bi orderable knot groups.


FakeNewsLab Experimental Study on Biases and Pitfalls Preventing us from Distinguishing True from False News. Misinformation posting and spreading in Social Media is ignited by personal decisions on the truthfulness of news that may cause wide and deep cascades at a large scale in a fraction of minutes. When individuals are exposed to information, they usually take a few seconds to decide if the content (or the source) is reliable, and eventually to share it. Although the opportunity to verify the rumour is often just one click away, many users fail to make a correct evaluation. We studied this phenomenon by implementing a web based questionnaire that was compiled by 7,298 different volunteers. Participants were asked to mark 20 news as true or false. Interestingly, false news is correctly identified more frequently than true news, but showing the full article instead of just the title, surprisingly, does not increase general accuracy. Also, displaying the original source of the news may contribute to mislead the user in some cases, while the wisdom of the crowd can positively assist individuals' ability to classify correctly. Furthermore, participants that autonomously opened an additional browser tab while compiling the survey show higher accuracy than users that did not. This suggests a parallel fact checking activity that support users' decisions. Finally, users that declare themselves as young adults are also those who open a new tab more often, suggesting more familiarity with the Web.


Word Level Coreference Resolution. Recent coreference resolution models rely heavily on span representations to find coreference links between word spans. As the number of spans is O(n 2) in the length of text and the number of potential links is O(n 4) , various pruning techniques are necessary to make this approach computationally feasible. We propose instead to consider coreference links between individual words rather than word spans and then reconstruct the word spans. This reduces the complexity of the coreference model to O(n 2) and allows it to consider all potential mentions without pruning any of them out. We also demonstrate that, with these changes, SpanBERT for coreference resolution will be significantly outperformed by RoBERTa. While being highly efficient, our model performs competitively with recent coreference resolution systems on the OntoNotes benchmark.


What's wrong with this video? Comparing Explainers for Deepfake Detection. Deepfakes are computer manipulated videos where the face of an individual has been replaced with that of another. Software for creating such forgeries is easy to use and ever more popular, causing serious threats to personal reputation and public security. The quality of classifiers for detecting deepfakes has improved with the releasing of ever larger datasets, but the understanding of why a particular video has been labelled as fake has not kept pace. In this work we develop, extend and compare white box, black box and model specific techniques for explaining the labelling of real and fake videos. In particular, we adapt SHAP, GradCAM and self attention models to the task of explaining the predictions of state of the art detectors based on EfficientNet, trained on the Deepfake Detection Challenge (DFDC) dataset. We compare the obtained explanations, proposing metrics to quantify their visual features and desirable characteristics, and also perform a user survey collecting users' opinions regarding the usefulness of the explainers.


On the optimal controllability time for linear hyperbolic systems with time dependent coefficients. The optimal time for the controllability of linear hyperbolic systems in one dimensional space with one side controls has been obtained recently for time independent coefficients in our previous works. In this paper, we consider linear hyperbolic systems with time varying zero order terms. We show the possibility that the optimal time for the null controllability becomes significantly larger than the one of the time invariant setting even when the zero order term is indefinitely differentiable. When the analyticity with respect to time is imposed for the zero order term, we also establish that the optimal time is the same as in the time independent setting.


In what sense space dimensionality can be used to cast light into cultural anthropology?. Humans have always constructed spaces, through Mythos and Logos, as part of an aspiration to capture the essence of the changing world. This has been a permanent endeavour since the invention of language. By doing this, in fact, Humankind started constructing itself we are beings in constant evolutionary process in real and imaginary spaces. Our concepts of Space and our anthropological ideas, specially the fundamental concepts of subject and subjectivity, are intertwined and intimately connected. We believe that the great narratives about Humanity, which ultimately define our view of ourselves, are entangled with those concepts that Cassirer identified as the cornerstones of culture space, time, and number. To explore these ideas, the authors wrote an essay, in 2017, in a book format, in which the fundamental role of real and imaginary spaces (and especially of their dimensionalities) in the History of Culture was discussed. This book, titled "O Livro, o Espa c c o e a Natureza Ensaio Sobre a Leitura do Mundo, as Muta c c oes da Cultura e do Sujeito", has a preface written by Francisco Antonio Doria. As many of the issues treated there are among his multiple interests, it was decided to revisit here the problems of subjectivity and subject's relationship with the dimensionality of space including the question of the architecture of books and other writing supports.


The Bou 'e Dupuis formula and the exponential hypercontractivity in the Gaussian space. This paper concerns a variational representation formula for Wiener functionals. Let B B t t ge 0 be a standard d dimensional Brownian motion. Bou 'e and Dupuis (1998) showed that, for any bounded measurable functional F(B) of B up to time 1 , the expectation mathbb E ! left e F(B) right admits a variational representation in terms of drifted Brownian motions. In this paper, with a slight modification of insightful reasoning by Lehec (2013) allowing also F(B) to be a functional of B over the whole time interval, we prove that the Bou 'e Dupuis formula holds true provided that both e F(B) and F(B) are integrable, relaxing conditions in earlier works. We also show that the formula implies the exponential hypercontractivity of the Ornstein Uhlenbeck semigroup in mathbb R d , and hence, due to their equivalence, implies the logarithmic Sobolev inequality in the d dimensional Gaussian space.


Geometry of fitness landscapes Peaks, shapes and universal positive epistasis. Darwinian evolution is driven by random mutations, genetic recombination (gene shuffling) and selection that favors genotypes with high fitness. For systems where each genotype can be represented as a bitstring of length L , an overview of possible evolutionary trajectories is provided by the L cube graph with nodes labeled by genotypes and edges directed toward the genotype with higher fitness. Peaks (sinks in the graphs) are important since a population can get stranded at a suboptimal peak. The fitness landscape is defined by the fitness values of all genotypes in the system. Some notion of curvature is necessary for a more complete analysis of the landscapes, including the effect of recombination. The shape approach uses triangulations (shapes) induced by fitness landscapes. The main topic for this work is the interplay between peak patterns and shapes. Because of constraints on the shapes for L 3 imposed by peaks, there are in total 25 possible combinations of peak patterns and shapes. Similar constraints exist for higher L . Specifically, we show that the constraints induced by the staircase triangulation can be formulated as a condition of emph universal positive epistasis , an order relation on the fitness effects of arbitrary sets of mutations that respects the inclusion relation between the corresponding genetic backgrounds. We apply the concept to a large protein fitness landscape for an immunoglobulin binding protein expressed in Streptococcal bacteria.


Safe reopening of university campuses is possible with COVID 19 vaccination. We construct an agent based SEIR model to simulate COVID 19 spread at a 16000 student mostly non residential urban university during the Fall 2021 Semester. We find that mRNA vaccine coverage above 80 makes it possible to safely reopen to in person instruction. If vaccine coverage is 100 , then our model indicates that facemask use is not necessary. Our simulations with vaccine coverage below 70 exhibit a right skew for total infections over the semester, which suggests that high levels of infection are not exceedingly rare with campus social connections the main transmission route. Less effective vaccines or incidence of new variants may require additional intervention such as screening testing to reopen safely.


Information Driven Modeling of Biomolecular Complexes. Proteins play crucial roles in every cellular process by interacting with each other, with nucleic acids, metabolites, and other molecules. The resulting assemblies can be very large and intricate and pose challenges to experimental methods. In the current era of integrative modeling, it is often only by a combination of various experimental techniques and computations that 3D models of those molecular machines can be obtained. Among the various computational approaches available, molecular docking is often the method of choice when it comes to predicting 3D structures of complexes. Docking can generate particularly accurate models when taking into account the available information on the complex of interest. We review here the use of experimental and bioinformatics data in protein protein docking, describing recent software developments and highlighting applications for the modeling of antibody antigen complexes and membrane protein complexes, and the use of evolutionary and shape information.


Fresnel reflection boundary for radiative transport lattice Boltzmann methods in highly scattering volume. With its roots in kinetic theory, the lattice Boltzmann method (LBM) cannot only be used to solve complex fluid flows but also radiative transport in volume. The present work derives a novel Fresnel boundary scheme for radiative transport LBM, based on Fresnel's equation, which depicts the partly reflected radiation on surfaces. Driven from a boundary modeling and discussion on the microscopic level, incorporating Fresnel's equation, it is developed a boundary model for the mesoscopic radiative transport LBM. At an intermediate step, the Fresnel's equation is related to well known partial differential (Robin) equations, based on a bottom up approach where the P1 Approximation is deployed. To connect the novel boundary scheme to the so derived target equation, a Chapman Enskog expansion is examined in addition. Both techniques together, point out how to interpret microscopic modeling by the means of macroscopic expressions and as a consequence how, to chose simulation parameters according to the specific boundary. The numerical tests suggest that the proposed boundary is first order convergent. The paper closes with a showcase, where the novel boundary method for radiative transport LBM is applied to a setup with multiple LED spots.


A Method to Analyze Multiple Social Identities in Twitter Bios. Twitter users signal social identity in their profile descriptions, or bios, in a number of important but complex ways that are not well captured by existing characterizations of how identity is expressed in language. Better ways of defining and measuring these expressions may therefore be useful both in understanding how social identity is expressed in text, and how the self is presented on Twitter. To this end, the present work makes three contributions. First, using qualitative methods, we identify and define the concept of a personal identifier, which is more representative of the ways in which identity is signaled in Twitter bios. Second, we propose a method to extract all personal identifiers expressed in a given bio. Finally, we present a series of validation analyses that explore the strengths and limitations of our proposed method. Our work opens up exciting new opportunities at the intersection between the social psychological study of social identity and the study of how we compose the self through markers of identity on Twitter and in social media more generally.


Enjoy the Salience Towards Better Transformer based Faithful Explanations with Word Salience. Pretrained transformer based models such as BERT have demonstrated state of the art predictive performance when adapted into a range of natural language processing tasks. An open problem is how to improve the faithfulness of explanations (rationales) for the predictions of these models. In this paper, we hypothesize that salient information extracted a priori from the training data can complement the task specific information learned by the model during fine tuning on a downstream task. In this way, we aim to help BERT not to forget assigning importance to informative input tokens when making predictions by proposing SaLoss an auxiliary loss function for guiding the multi head attention mechanism during training to be close to salient information extracted a priori using TextRank. Experiments for explanation faithfulness across five datasets, show that models trained with SaLoss consistently provide more faithful explanations across four different feature attribution methods compared to vanilla BERT. Using the rationales extracted from vanilla BERT and SaLoss models to train inherently faithful classifiers, we further show that the latter result in higher predictive performance in downstream tasks.


Optimizing Bayesian Recurrent Neural Networks on an FPGA based Accelerator. Neural networks have demonstrated their outstanding performance in a wide range of tasks. Specifically recurrent architectures based on long short term memory (LSTM) cells have manifested excellent capability to model time dependencies in real world data. However, standard recurrent architectures cannot estimate their uncertainty which is essential for safety critical applications such as in medicine. In contrast, Bayesian recurrent neural networks (RNNs) are able to provide uncertainty estimation with improved accuracy. Nonetheless, Bayesian RNNs are computationally and memory demanding, which limits their practicality despite their advantages. To address this issue, we propose an FPGA based hardware design to accelerate Bayesian LSTM based RNNs. To further improve the overall algorithmic hardware performance, a co design framework is proposed to explore the most fitting algorithmic hardware configurations for Bayesian RNNs. We conduct extensive experiments on healthcare applications to demonstrate the improvement of our design and the effectiveness of our framework. Compared with GPU implementation, our FPGA based design can achieve up to 10 times speedup with nearly 106 times higher energy efficiency. To the best of our knowledge, this is the first work targeting acceleration of Bayesian RNNs on FPGAs.


deepstruct linking deep learning and graph theory. deepstruct connects deep learning models and graph theory such that different graph structures can be imposed on neural networks or graph structures can be extracted from trained neural network models. For this, deepstruct provides deep neural network models with different restrictions which can be created based on an initial graph. Further, tools to extract graph structures from trained models are available. This step of extracting graphs can be computationally expensive even for models of just a few dozen thousand parameters and poses a challenging problem. deepstruct supports research in pruning, neural architecture search, automated network design and structure analysis of neural networks.


Dense e e plasma formation in magnetic dipole wave vacuum breakdown by 10 PW class lasers. When studying the interaction of matter with extreme fields using multipetawatt lasers, there are two limiting cases maximizing either the electric field or the magnetic field. Here, the main attention is paid to the optimal configuration of laser beams in the form of an m dipole wave, which maximizes the magnetic field, and the corresponding production of pair plasma via a QED cascade using 10 PW class lasers. We show that the threshold of vacuum breakdown with respect to avalanche like pair generation is about 10 PW. Using 3D PIC modeling in the specified fields, we go deeper into the physics of vacuum breakdown, i.e. we examined in detail the individual trajectories of particles produced in inhomogeneous electric and magnetic fields, the space time distributions of pair densities on the avalanche stage, and the energy distributions of charged particles and gamma photons. The forming plasma structures represent concentric rings around the central magnetic axis, which can result in significant change of laser plasma interaction in comparison with the case of an e dipole wave.


Market Making with Stochastic Liquidity Demand Simultaneous Order Arrival and Price Change Forecasts. We provide an explicit characterization of the optimal market making strategy in a discrete time Limit Order Book (LOB). In our model, the number of filled orders during each period depends linearly on the distance between the fundamental price and the market maker's limit order quotes, with random slope and intercept coefficients. The high frequency market maker (HFM) incurs an end of the day liquidation cost resulting from linear price impact. The optimal placement strategy incorporates in a novel and parsimonious way forecasts about future changes in the asset's fundamental price. We show that the randomness in the demand slope reduces the inventory management motive, and that a positive correlation between demand slope and investors' reservation prices leads to wider spreads. Our analysis reveals that the simultaneous arrival of buy and sell market orders (i) reduces the shadow cost of inventory, (ii) leads the HFM to reduce price pressures to execute larger flows, and (iii) introduces patterns of nonlinearity in the intraday dynamics of bid and ask spreads. Our empirical study shows that the market making strategy outperforms those which ignores randomness in demand, simultaneous arrival of buy and sell market orders, and local drift in the fundamental price.


Fast n fold Boolean Convolution via Additive Combinatorics. We consider the problem of computing the Boolean convolution (with wraparound) of n vectors of dimension m , or, equivalently, the problem of computing the sumset A 1 A 2 ldots A n for A 1, ldots,A n subseteq mathbb Z m . Boolean convolution formalizes the frequent task of combining two subproblems, where the whole problem has a solution of size k if for some i the first subproblem has a solution of size i and the second subproblem has a solution of size k i . Our problem formalizes a natural generalization, namely combining solutions of n subproblems subject to a modular constraint. This simultaneously generalises Modular Subset Sum and Boolean Convolution (Sumset Computation). Although nearly optimal algorithms are known for special cases of this problem, not even tiny improvements are known for the general case. We almost resolve the computational complexity of this problem, shaving essentially a factor of n from the running time of previous algorithms. Specifically, we present a emph deterministic algorithm running in emph almost linear time with respect to the input plus output size k . We also present a emph Las Vegas algorithm running in emph nearly linear expected time with respect to the input plus output size k . Previously, no deterministic or randomized o(nk) algorithm was known. At the heart of our approach lies a careful usage of Kneser's theorem from Additive Combinatorics, and a new deterministic almost linear output sensitive algorithm for non negative sparse convolution. In total, our work builds a solid toolbox that could be of independent interest.


The pigenhole principle and multicolor Ramsey numbers. For integers k,r geq 2 , the diagonal Ramsey number R r(k) is the minimum N in mathbb N such that every r coloring of the edges of a complete graph on N vertices yields on a monochromatic subgraph on k vertices. Here we make a careful effort of extracting explicit upper bounds for R r(k) from the pigeonhole principle alone. Our main term improves on previously documented explicit bounds for r geq 3 , and we also consider an often ignored secondary term, which allows us to subtract a uniformly bounded below positive proportion of the main term. Asymptotically, we give a self contained proof that R r(k) leq left( frac 3 e 2 right) frac (r(k 2))! ((k 2)!) r (1 o r to infty (1)), and we conclude by noting that our methods combine with previous estimates on R r(3) to improve the constant frac 3 e 2 to frac 3 e 2 frac d 48 , where d 66 R 4(3) geq 4 . We also compare our formulas, and previously documented formulas, to some collected numerical data.


Provably efficient, succinct, and precise explanations. We consider the problem of explaining the predictions of an arbitrary blackbox model f given query access to f and an instance x , output a small set of x 's features that in conjunction essentially determines f(x) . We design an efficient algorithm with provable guarantees on the succinctness and precision of the explanations that it returns. Prior algorithms were either efficient but lacked such guarantees, or achieved such guarantees but were inefficient. We obtain our algorithm via a connection to the problem of sl implicitly learning decision trees. The implicit nature of this learning task allows for efficient algorithms even when the complexity of f necessitates an intractably large surrogate decision tree. We solve the implicit learning problem by bringing together techniques from learning theory, local computation algorithms, and complexity theory. Our approach of "explaining by implicit learning" shares elements of two previously disparate methods for post hoc explanations, global and local explanations, and we make the case that it enjoys advantages of both.


Global Research Trends in the Modern Language Journal from 1999 to 2018 A Data Driven Analysis. The present study conducts a scientometric study of the Modern Language Journal literature from 1999 to 2018 based on the database of Web of Science, 2018. A total of 2564 items resulted from the publication name using "Modern Language Journal" as the search term. Based on the number of publications during the study period no consistent growth is observed in the research activities pertaining to the journal. The annual distribution of publications, number of authors, institution productivity, country wise publications and Citations are analyzed. Highly productive authors, institutions, and countries are identified. The results reveal that the maximum number of papers 179 is published in the year 1999. It was also observed that Byrnes H is the most productive, contributed 51 publications and Kramsch C is most cited author in the field having 543 global citations. The highest number (38.26 ) of publications, contributed from USA and the foremost productive establishment was University of Iowa.


Typical structure of hereditary properties of binary matroids. We prove an arithmetic analogue of the typical structure theorem for graph hereditary properties due to Alon, Balogh, Bollob 'as and Morris.


Adversarial Motion Modelling helps Semi supervised Hand Pose Estimation. Hand pose estimation is difficult due to different environmental conditions, object and self occlusion as well as diversity in hand shape and appearance. Exhaustively covering this wide range of factors in fully annotated datasets has remained impractical, posing significant challenges for generalization of supervised methods. Embracing this challenge, we propose to combine ideas from adversarial training and motion modelling to tap into unlabeled videos. To this end we propose what to the best of our knowledge is the first motion model for hands and show that an adversarial formulation leads to better generalization properties of the hand pose estimator via semi supervised training on unlabeled video sequences. In this setting, the pose predictor must produce a valid sequence of hand poses, as determined by a discriminative adversary. This adversary reasons both on the structural as well as temporal domain, effectively exploiting the spatio temporal structure in the task. The main advantage of our approach is that we can make use of unpaired videos and joint sequence data both of which are much easier to attain than paired training data. We perform extensive evaluation, investigating essential components needed for the proposed framework and empirically demonstrate in two challenging settings that the proposed approach leads to significant improvements in pose estimation accuracy. In the lowest label setting, we attain an improvement of 40 in absolute mean joint error.


Mereological Emptiness for the Signed Number Problem. It is well known (at least in the education research literature) that primary school students face considerable difficulties in the understanding of negative integers (and numbers), related operations and their visualizations. In the present research, a proposal based on mereological principles is advanced towards solving aspects of the problem by the present author. Her proposal builds on the actual semantic content of what is taught expressed in algebraic terms, and natural types that can potentially be an explicit habit.


On the inhomogeneous NLS with inverse square potential. We consider the inhomogeneous nonlinear Schr "odinger equation with inverse square potential in mathbb R N i u t mathcal L a u lambda x b u alpha u 0, mathcal L a Delta frac a x 2 , where lambda pm1 , alpha,b 0 and a frac (N 2) 2 4 . We first establish sufficient conditions for global existence and blow up in H 1 a( mathbb R N) for lambda 1 , using a Gagliardo Nirenberg type estimate. In the sequel, we study local and global well posedness in H 1 a( mathbb R N) in the H 1 subcritical case, applying the standard Strichartz estimates combined with the fixed point argument. The key to do that is to establish good estimates on the nonlinearity. Making use of these estimates, we also show a scattering criterion and construct a wave operator in H 1 a( mathbb R N) , for the mass supercritical and energy subcritical case.


Three Operator Splitting with Subgradients, Stochastic Gradients, and Adaptive Learning Rates. Three Operator Splitting (TOS) (Davis Yin, 2017) can minimize the sum of multiple convex functions effectively when an efficient gradient oracle or proximal operator is available for each term. This requirement often fails in machine learning applications (i) instead of full gradients only stochastic gradients may be available and (ii) instead of proximal operators, using subgradients to handle complex penalty functions may be more efficient and realistic. Motivated by these concerns, we analyze three potentially valuable extensions of TOS. The first two permit using subgradients and stochastic gradients, and are shown to ensure a mathcal O (1 sqrt t ) convergence rate. The third extension AdapTOS endows TOS with adaptive step sizes. For the important setting of optimizing a convex loss over the intersection of convex sets AdapTOS attains universal convergence rates, i.e., the rate adapts to the unknown smoothness degree of the objective. We compare our proposed methods with competing methods on various applications.


Computing and Authentication Practices in Global Oil and Gas Fields. Oil and gas fields are a critical part of our infrastructure, and vulnerable to attack by powerful adversaries. In addition, these are often difficult work environments, with constraints on space, clothing, and more. Yet there is little research on the technology practices and constraints of workers in these environments. We present what we believe is the first survey of oil and gas field workers located around the world. We establish the presence and status of a variety of computing devices and of the security practices that govern their use. We also determine the working conditions (such as personal protective equipment) under which these devices are used, which impacts usable security aspects like feasible forms of authentication. We extend these basic insights with additional information from a small number of in depth interviews. Our preliminary work suggests many directions for improving security in this critical sector.


Conservative Policy Construction Using Variational Autoencoders for Logged Data with Missing Values. In high stakes applications of data driven decision making like healthcare, it is of paramount importance to learn a policy that maximizes the reward while avoiding potentially dangerous actions when there is uncertainty. There are two main challenges usually associated with this problem. Firstly, learning through online exploration is not possible due to the critical nature of such applications. Therefore, we need to resort to observational datasets with no counterfactuals. Secondly, such datasets are usually imperfect, additionally cursed with missing values in the attributes of features. In this paper, we consider the problem of constructing personalized policies using logged data when there are missing values in the attributes of features in both training and test data. The goal is to recommend an action (treatment) when Xt , a degraded version of Xb with missing values, is observed. We consider three strategies for dealing with missingness. In particular, we introduce the textit conservative strategy where the policy is designed to safely handle the uncertainty due to missingness. In order to implement this strategy we need to estimate posterior distribution p( Xb Xt) , we use variational autoencoder to achieve this. In particular, our method is based on partial variational autoencoders (PVAE) which are designed to capture the underlying structure of features with missing values.


Euclidean domains in complex manifolds. In this paper we find big Euclidean domains in complex manifolds. We consider open neighbourhoods of sets of the form K cup M in a complex manifold X , where K is a compact mathscr O(U) convex set in an open Stein neighbourhood U of K , M is an embedded Stein submanifold of X , and K cap M is compact and mathscr O(M) convex. We prove a Docquier Grauert type theorem concerning biholomorphic equivalence of neighbourhoods of such sets, and we give sufficient conditions for the existence of Stein neighbourhoods of K cup M , biholomorphic to domains in mathbb C n with n dim X , such that M is mapped onto a closed complex submanifold of mathbb C n .


Infrasound of a Wind Turbine Reanalyzed as Power Spectrum and Power Spectral Density. The infrasound levels due to the blade tower interaction generated by a wind turbine in the publication by Pilger and Ceranna (JSV, Vol. 388, pp. 188 200, 2017) have to be corrected to be interpreted as sound pressure level. Also, the electrical power of the wind turbine should be corrected for the high wind case to 660 kW. We provide a reanalysis of the measured data with a power spectrum showing levels for the low frequency signal of the wind turbine about 34 dB below the original work. All measured levels at a distance of 200 m from the wind turbine's infrasound signal are well below the hearing threshold.


HS integral and Eisenstein integral mixed circulant graphs. A mixed graph is called emph second kind hermitian integral (or emph HS integral ) if the eigenvalues of its Hermitian adjacency matrix of second kind are integers. A mixed graph is called emph Eisenstein integral if the eigenvalues of its (0, 1) adjacency matrix are Eisenstein integers. We characterize the set S for which a mixed circulant graph text Circ ( mathbb Z n, S) is HS integral. We also show that a mixed circulant graph is Eisenstein integral if and only if it is HS integral. Further, the eigenvalues and the HS eigenvalues of some oriented circulant graphs are expressed in terms of generalized M ddot text o bius function.


A positivity and monotonicity preserving nonlinear operator splitting approach for approximating solutions to quenching combustion semilinear partial differential equations. In recent years, there has been a large increase in interest in numerical algorithms which preserve various qualitative features of the original continuous problem. Herein, we propose and investigate a numerical algorithm which preserves qualitative features of so called quenching combustion partial differential equations (PDEs). Such PDEs are often used to model solid fuel ignition processes or enzymatic chemical reactions and are characterized by their singular nonlinear reaction terms and the exhibited positivity and monotonicity of their solutions on their time intervals of existence. In this article, we propose an implicit nonlinear operator splitting algorithm which allows for the natural preservation of these features. The positivity and monotonicity of the algorithm is rigorously proven. Furthermore, the convergence analysis of the algorithm is carried out and the explicit dependence on the singularity is quantified in a nonlinear setting.