Neglecting discretization corrections in regularized singular or nearly singular integrals. A method for computing singular or nearly singular integrals on closed surfaces was presented by J. T. Beale, W. Ying, and J. R. Wilson Comm. Comput. Phys. 20 (2016), 733 753, arXiv 1508.00265 and applied to single and double layer potentials for harmonic functions. It uses regularized kernels, a straightforward quadrature rule, and corrections added for smoothing and discretization errors. In this note we give estimates for the discretization corrections which show that they can reasonably be neglected with proper choice of numerical parameters.


Fine Grained Visual Classification with Efficient End to end Localization. The term fine grained visual classification (FGVC) refers to classification tasks where the classes are very similar and the classification model needs to be able to find subtle differences to make the correct prediction. State of the art approaches often include a localization step designed to help a classification network by localizing the relevant parts of the input images. However, this usually requires multiple iterations or passes through a full classification network or complex training schedules. In this work we present an efficient localization module that can be fused with a classification network in an end to end setup. On the one hand the module is trained by the gradient flowing back from the classification network. On the other hand, two self supervised loss functions are introduced to increase the localization accuracy. We evaluate the new model on the three benchmark datasets CUB200 2011, Stanford Cars and FGVC Aircraft and are able to achieve competitive recognition performance.


Giant dipole oscillations and ionization of heavy atoms by intense electromagnetic fields. Binding energy of heavy atoms is estimated by means of the Thomas Fermi model, giant dipole oscillation are highlighted and ionization is discussed.


On cascade of kinetic energy in compressible hydrodynamic turbulence. Properties of the turbulent cascade of kinetic energy are studied using direct numerical simulations of three dimensional hydrodynamic decaying turbulence with a moderate Reynolds number and the initial Mach number M 1 . Compressible and incompressible versions of the Karman Howarth Monin (KHM) and low pass filtering coarse graining approaches are compared. In the simulation the total energy is well conserved the scale dependent KHM and coarse grained energy equations are also well conserved the two approaches show similar results, the system does not have an inertial range for the cascade of kinetic energy, the region where this cascade dominates also have a non negligible contribution of the kinetic energy decay, dissipation, and pressure dilatation effects. While the two approaches give semi quantitatively similar results for the kinetic energy cascade, dissipation and pressure dilatation rates, they differ in the increment separation and filtering scales these scales are not simply related. The two approaches may be used to find the inertial range and to determine the cascade dissipation rate of the kinetic energy.


Optimizing Memory Access Patterns for Deep Learning Accelerators. Deep learning (DL) workloads are moving towards accelerators for faster processing and lower cost. Modern DL accelerators are good at handling the large scale multiply accumulate operations that dominate DL workloads however, it is challenging to make full use of the compute power of an accelerator since the data must be properly staged in a software managed scratchpad memory. Failing to do so can result in significant performance loss. This paper proposes a systematic approach which leverages the polyhedral model to analyze all operators of a DL model together to minimize the number of memory accesses. Experiments show that our approach can substantially reduce the impact of memory accesses required by common neural network models on a homegrown AWS machine learning inference chip named Inferentia, which is available through Amazon EC2 Inf1 instances.


The Multi Agent Programming Contest A r 'esum 'e. The Multi Agent Programming Contest, MAPC, is an annual event organized since 2005 out of Clausthal University of Technology. Its aim is to investigate the potential of using decentralized, autonomously acting intelligent agents, by providing a complex scenario to be solved in a competitive environment. For this we need suitable benchmarks where agent based systems can shine. We present previous editions of the contest and also its current scenario and results from its use in the 2019 MAPC with a special focus on its suitability. We conclude with lessons learned over the years.


Demonstration of Two Color X ray Free Electron Laser Pulses with a Sextupole Magnet. We present measurements of two color X ray free electron laser (FEL) pulses generated with a novel scheme utilizing a sextupole magnet. The sextupole, in combination with a standard orbit control tool, is used to suppress the radiation from the bunch core, while keeping the head and the tail of the beam lasing, each at a different photon energy. The method is simple, cost effective and applicable at any repetition rate. We demonstrate the tunability of the scheme and discuss its advantages and practical limitations.


A two stage fourth order gas kinetic CPR method for the Navier Stokes equations on triangular meshes. A highly efficient gas kinetic scheme with fourth order accuracy in both space and time is developed for the Navier Stokes equations on triangular meshes. The scheme combines an efficient correction procedure via reconstruction (CPR) framework with a robust gas kinetic flux formula, which computes both the flux and its time derivative. The availability of the flux time derivative makes it straightforward to adopt an efficient two stage temporal discretization to achieve fourth order time accuracy. In addition, through the gas kinetic evolution model, the inviscid and viscous fluxes are coupled and computed uniformly without any separate treatment for the viscous fluxes. As a result, the current scheme is more efficient than traditional explicit CPR methods with a separate treatment for viscous fluxes, and a fourth order Runge Kutta approach. Furthermore, a robust and accurate subcell finite volume (SCFV) limiting procedure is extended to the CPR framework for troubled cells, resulting in subcell resolution of flow discontinuities. Numerical tests demonstrate the high accuracy, efficiency and robustness of the current scheme in a wide range of inviscid and viscous flow problems from subsonic to supersonic speeds.


A Domain Specific Language for Verifying Software Requirement Constraints. Software requirement analysis can certainly benefit from prevention and early detection of failures, in particular by some kind of automatic analysis. Formal methods offer means to represent and analyze requirements with rigorous tools, avoiding ambiguities and allowing automatic verification of requirement consistency. However, formalisms often clash in the culture or lack of skills of software analysts, making them challenging to apply. In this article, we propose a Domain Specific Language (DSL) based on Set Theory for requirement analysts. The Graphical InvaRiant Language (GIRL) can be used to specify software requirement structural invariants, with entities and their relationships. Those invariants can then have their consistency evaluated by the Alloy Analyzer, based on a mapping semantics we provide for transforming GIRL models into Alloy specifications with no user intervention. With a prototypical language editor and transformations implemented into an Eclipse plugin, we carried out a qualitative study with requirement analysts working for a government software company in Brazil, to evaluate usability and effectiveness of the GIRL based analysis of real software requirements. The participants were able to effectively use the underlying formal analysis, since 79 out of 80 assigned invariants were correctly modeled. While participants perceived as low the complexity of learning and using GIRL's simplest, set based structures and relationships, the most complex logical structures, such as quantification and implication, were challenging. Furthermore, almost all post study evaluations from the participants were positive, especially as a tool for discovering requirement inconsistencies.


Structured strong boldsymbol ell ifications for structured matrix polynomials in the monomial basis. In the framework of Polynomial Eigenvalue Problems, most of the matrix polynomials arising in applications are structured polynomials (namely (skew )symmetric, (skew )Hermitian, (anti )palindromic, or alternating). The standard way to solve Polynomial Eigenvalue Problems is by means of linearizations. The most frequently used linearizations belong to general constructions, valid for all matrix polynomials of a fixed degree, known as em companion linearizations . It is well known, however, that is not possible to construct companion linearizations that preserve any of the previous structures for matrix polynomials of even degree. This motivates the search for more general companion forms, in particular em companion ell ifications . In this paper, we present, for the first time, a family of (generalized) companion ell ifications that preserve any of these structures, for matrix polynomials of degree k (2d 1) ell . We also show how to construct sparse ell ifications within this family. Finally, we prove that there are no structured companion quadratifications for quartic matrix polynomials.


Open Domain Event Extraction Using Neural Latent Variable Models. We consider open domain event extraction, the task of extracting unconstraint types of events from news clusters. A novel latent variable neural model is constructed, which is scalable to very large corpus. A dataset is collected and manually annotated, with task specific evaluation metrics being designed. Results show that the proposed unsupervised model gives better performance compared to the state of the art method for event schema induction.


Randomly connected networks generate emergent selectivity and predict decoding properties of large populations of neurons. Advances in neural recording methods enable sampling from populations of thousands of neurons during the performance of behavioral tasks, raising the question of how recorded activity relates to the theoretical models of computations underlying performance. In the context of decision making in rodents, patterns of functional connectivity between choice selective cortical neurons, as well as broadly distributed choice information in both excitatory and inhibitory populations, were recently reported 1 . The straightforward interpretation of these data suggests a mechanism relying on specific patterns of anatomical connectivity to achieve selective pools of inhibitory as well as excitatory neurons. We investigate an alternative mechanism for the emergence of these experimental observations using a computational approach. We find that a randomly connected network of excitatory and inhibitory neurons generates single cell selectivity, patterns of pairwise correlations, and indistinguishable excitatory and inhibitory readout weight distributions, as observed in recorded neural populations. Further, we make the readily verifiable experimental predictions that, for this type of evidence accumulation task, there are no anatomically defined sub populations of neurons representing choice, and that choice preference of a particular neuron changes with the details of the task. This work suggests that distributed stimulus selectivity and patterns of functional organization in population codes could be emergent properties of randomly connected networks.


A Modified Fourier Mellin Approach for Source Device Identification on Stabilized Videos. To decide whether a digital video has been captured by a given device, multimedia forensic tools usually exploit characteristic noise traces left by the camera sensor on the acquired frames. This analysis requires that the noise pattern characterizing the camera and the noise pattern extracted from video frames under analysis are geometrically aligned. However, in many practical scenarios this does not occur, thus a re alignment or synchronization has to be performed. Current solutions often require time consuming search of the realignment transformation parameters. In this paper, we propose to overcome this limitation by searching scaling and rotation parameters in the frequency domain. The proposed algorithm tested on real videos from a well known state of the art dataset shows promising results.


Continuity and Discontinuity of Seminorms on Infinite Dimensional Vector Spaces. Let S be a seminorm on an infinite dimensional real or complex vector space X . Our purpose in this note is to study the continuity and discontinuity properties of S with respect to certain norm topologies on X .


Ulam stability of Jensen functional inequality in groups. In this paper we establish the Ulam stability of Jensen functional inequality in some classes of groups.


Self Supervised Learning for Stereo Reconstruction on Aerial Images. Recent developments established deep learning as an inevitable tool to boost the performance of dense matching and stereo estimation. On the downside, learning these networks requires a substantial amount of training data to be successful. Consequently, the application of these models outside of the laboratory is far from straight forward. In this work we propose a self supervised training procedure that allows us to adapt our network to the specific (imaging) characteristics of the dataset at hand, without the requirement of external ground truth data. We instead generate interim training data by running our intermediate network on the whole dataset, followed by conservative outlier filtering. Bootstrapped from a pre trained version of our hybrid CNN CRF model, we alternate the generation of training data and network training. With this simple concept we are able to lift the completeness and accuracy of the pre trained version significantly. We also show that our final model compares favorably to other popular stereo estimation algorithms on an aerial dataset.


A spectral method to incidence balance of oriented hypergraphs and induced signed hypergraphs. An oriented hypergraph is a hypergraph together with an incidence orientation such that each edge vertex incidence is given a label of 1 or 1 . An oriented hypergraph is called incidence balanced if there exists a bipartition of the vertex set such that every edge intersects one part of the bipartition in positively incident vertices with the edge and other part in negatively incident vertices with the edge. In this paper, we investigate the incidence balance of oriented hypergraphs and induced signed hypergraphs by means of the eigenvalues of associated matrices and tensors, and provide a spectral method to characterize the incidence balanced oriented hypergraphs and their induced signed hypergraphs.


Identification and Remediation of Self Admitted Technical Debt in Issue Trackers. Technical debt refers to taking shortcuts to achieve short term goals, which might negatively influence software maintenance in the long term. There is increasing attention on technical debt that is admitted by developers in source code comments (termed as self admitted technical debt or SATD). But SATD in issue trackers is relatively unexplored. We performed a case study, where we manually examined 500 issues from two open source projects (i.e. Hadoop and Camel), which contained 152 SATD items. We found that 1) eight types of technical debt are identified in issues, namely architecture, build, code, defect, design, documentation, requirement, and test debt 2) developers identify technical debt in issues in three different points in time, and a small part is identified by its creators 3) the majority of technical debt is paid off, 4) mostly by those who identified it or created it 5) the median time and average time to repay technical debt are 872.3 and 25.0 hours respectively.


Differential Projective Modules over Differential Rings, II. Differential modules over a commutative differential ring R which are finitely generated projective as ring modules, with differential homomorphisms, form an additive category, so their isomorphism classes form a monoid. We study the quotient monoid of this monoid by the submonoid of isomorphism classes of free modules with component wise derivation. This quotient monoid has the reduced K group of R (ignoring the derivation) as an image and contains the reduced K group of the constants of R as its subgroup of units. This monoid provides a description of the isomorphism classes of differential projective R modules up to an equivalence.


On 2 local nonlinear surjective isometries on normed spaces and C algebras. We prove that, if the closed unit ball of a normed space X has sufficiently many extreme points, then every mapping Phi from X into itself with the following property is affine For any pair of points in X , there exists a (not necessarily linear) surjective isometry on X that coincides with Phi at the two points. We also consider surjectivity of such a mapping in some special cases including C algebras.


The star avoidance game. Let n, k be positive integers. The (k 1) star avoidance game on K n is played as follows. Two players take it in turn to claim a (previously unclaimed) edge of the complete graph on n vertices. The first player to claim all edges of a subgraph isomorphic to a (k 1) star loses. Equivalently, each player must keep all degrees in the subgraph formed by his edges at most k . If all edges have been chosen and neither player has lost, the game is declared a draw. We prove that, for each fixed k , the game is a win for the second player for all n sufficiently large.


Early life imprints the hierarchy of T cell clone sizes. The adaptive immune system responds to pathogens by selecting clones of cells with specific receptors. While clonal selection in response to particular antigens has been studied in detail, it is unknown how a lifetime of exposures to many antigens collectively shape the immune repertoire. Here, through mathematical modeling and statistical analyses of T cell receptor sequencing data we demonstrate that clonal expansions during a perinatal time window leave a long lasting imprint on the human T cell repertoire. We demonstrate how the empirical scaling law relating the rank of the largest clones to their size can emerge from clonal growth during repertoire formation. We statistically identify early founded clones and find that they are indeed highly enriched among the largest clones. This enrichment persists even after decades of human aging, in a way that is quantitatively predicted by a model of fluctuating clonal selection. Our work presents a quantitative theory of human T cell dynamics compatible with the statistical laws of repertoire organization and provides a mechanism for how early clonal dynamics imprint the hierarchy of T cell clone sizes with implications for pathogen defense and autoimmunity.


Full Vectorial Modeling of Second Harmonic Generation in III V on insulator Nanowires. We model second harmonic generation in subwavelength III V on insulator waveguides. The large index contrast induces strong longitudinal electric field components that play an important role in the nonlinear conversion. We show that many different waveguide dimensions are suitable for efficient conversion of a fundamental quasi TE pump mode around the 1550 nm telecommunication wavelength to a higher order second harmonic mode.


Impedance boundary conditions for acoustic time harmonic wave propagation in viscous gases. We present Helmholtz or Helmholtz like equations for the approximation of the time harmonic wave propagation in gases with small viscosity, which are completed with local boundary conditions on rigid walls. We derived approximative models based on the method of multiple scales for the pressure and the velocity separately, both up to order 2. Approximations to the pressure are described by the Helmholtz equations with impedance boundary conditions, which relate its normal derivative to the pressure itself. The boundary conditions from first order on are of Wentzell type and include a second tangential derivative of the pressure proportional to the square root of the viscosity, and take thereby absorption inside the viscosity boundary layer of the underlying velocity into account. The velocity approximations are described by Helmholtz like equations for the velocity, where the Laplace operator is replaced by nabla operatorname div , and the local boundary conditions relate the normal velocity component to its divergence. The velocity approximations are for the so called far field and do not exhibit a boundary layer. Including a boundary corrector, the so called near field, the velocity approximation is accurate even up to the domain boundary. The boundary conditions are stable and asymptotically exact, which is justified by a complete mathematical analysis. The results of some numerical experiments are presented to illustrate the theoretical foundation.


The talented monoid of a Leavitt path algebra. There is a tight relation between the geometry of a directed graph and the algebraic structure of a Leavitt path algebra associated to it. In this note, we show a similar connection between the geometry of the graph and the structure of a certain monoid associated to it. This monoid is isomorphic to the positive cone of the graded K0 group of the Leavitt path algebra which is naturally equipped with a Z action. As an example, we show that a graph has a cycle without an exit if and only if the monoid has a periodic element. Consequently a graph has Condition (L) if and only if the group Z acts freely on the monoid. We go on to show that the algebraic structure of Leavitt path algebras (such as simplicity, purely infinite simplicity, or the lattice of ideals) can be described completely via this monoid. Therefore an isomorphism between the monoids (or graded K0's) of two Leavitt path algebras implies that the algebras have similar algebraic structures. These all confirm that the graded Grothendieck group could be a sought after complete invariant for the classification of Leavitt path algebras.


Retirement decision and optimal consumption investment under addictive habit persistence. This paper studies the retirement decision, optimal investment and consumption strategies under habit persistence for an agent with the opportunity to design the retirement time. The optimization problem is formulated as an interconnected optimal stopping and stochastic control problem (Stopping Control Problem) in a finite time horizon. The problem contains three state variables wealth x , habit level h and wage rate w . We aim to derive the retirement boundary of this "wealth habit wage" triplet (x,h,w) . The complicated dual relation is proposed and proved to convert the original problem to the dual one. We obtain the retirement boundary of the dual variables based on an obstacle type free boundary problem. Using dual relation we find the retirement boundary of primal variables and feed back forms of optimal strategies. We show that if the so called "de facto wealth" exceeds a critical proportion of wage, it will be optimal for the agent to choose to retire immediately. In numerical applications, we show how "de facto wealth" determines the retirement decisions and optimal strategies. Moreover, we observe discontinuity at retirement boundary investment proportion always jumps down upon retirement, while consumption may jump up or jump down, depending on the change of marginal utility. We also find that the agent with higher standard of life tends to work longer.


Photophoretic Levitation of Macroscopic Nanocardboard Plates. Scaling down miniature rotorcraft and flapping wing flyers to sub centimeter dimensions is challenging due to complex electronics requirements, manufacturing limitations, and the increase in viscous damping at low Reynolds numbers. Photophoresis, or light driven fluid flow, was previously used to levitate solid particles without any moving parts but only at the scale of 1 20 micrometers. Here, we leverage stiff architectured plates with 50 nm thickness to realize photophoretic levitation at the millimeter to centimeter scales. Instead of creating lift through conventional rotors or wings, our levitation occurs due to light induced thermal transpiration through micro channels within our plates. At atmospheric pressure, the plates hover above a solid substrate at heights of 0.5 mm by creating an air cushion beneath the plate. Moreover, at reduced pressures (10 1000 Pa), the increased speed of thermal transpiration through the plate's channels creates an air jet that enables mid air levitation and allows the plates to carry small payloads heavier than the plates themselves. This new propulsion mechanism for macroscopic structures creates the prospect of future nanotechnology enabled flying vehicles without any moving parts in the Earth's upper atmosphere and at the surface of other planets.


SimUSR A Simple but Strong Baseline for Unsupervised Image Super resolution. In this paper, we tackle a fully unsupervised super resolution problem, i.e., neither paired images nor ground truth HR images. We assume that low resolution (LR) images are relatively easy to collect compared to high resolution (HR) images. By allowing multiple LR images, we build a set of pseudo pairs by denoising and downsampling LR images and cast the original unsupervised problem into a supervised learning problem but in one level lower. Though this line of study is easy to think of and thus should have been investigated prior to any complicated unsupervised methods, surprisingly, there are currently none. Even more, we show that this simple method outperforms the state of the art unsupervised method with a dramatically shorter latency at runtime, and significantly reduces the gap to the HR supervised models. We submitted our method in NTIRE 2020 super resolution challenge and won 1st in PSNR, 2nd in SSIM, and 13th in LPIPS. This simple method should be used as the baseline to beat in the future, especially when multiple LR images are allowed during the training phase. However, even in the zero shot condition, we argue that this method can serve as a useful baseline to see the gap between supervised and unsupervised frameworks.


Dynamic characterization of beating sinusoidal wave oscillations in a thin slice solid state laser with coupled orthogonally polarized transverse modes. Dynamic characterization of self induced beating sinusoidal wave (BSW) oscillations and related phenomena, which were previously explored in a thin slice dual polarization solid state laser operating in the regime of incomplete mode locking of orthogonally polarized transverse eigenmodes, namely, a quasi locked state Laser Physics Letters, 15 (2018) 075001 , is performed in terms of the amplitude correlation coefficient and intensity circulation among various orthogonally polarized mode pairs. Polarization dependent symmetry breaking featuring a nonreciprocal intensity flow for different orthogonally polarized mode pairs operating in the BSW state are found. The singular state is hidden in a nonlinear system where a perfect mode locked BSW state with a large amplitude correlation coefficient featuring reciprocal intensity flow is established for a transverse mode pair possessing particular polarization directions that are related to the intensity ratio of orthogonally polarized transverse eigenmodes. Self organized chaos synchronization is shown to take place for such a singular transverse mode pair subjected to self mixing modulations.


Fully quantum embedding with density functional theory for full configuration interaction quantum Monte Carlo. In common with many high accuracy electronic structure methods, the initiator adaptation of full configuration interaction quantum Monte Carlo (i FCIQMC) has difficulty treating realistic systems with large numbers of electrons. This barrier has prevented the application of i FCIQMC to questions of catalysis that, even for the simplest of models, require high accuracy modeling of several features of the electronic structure, such as strong and dynamic correlation, and localized vs. delocalized bonding. We here present a fully quantum embedded version of i FCIQMC , which we apply to calculate the bond dissociation energy of an ionic bond (LiH) and a covalent bond (HF) physisorbed to a benzene molecule. The embedding is performed using a recently developed Huzinaga projection operator approach, which affords good synergy with i FCIQMC by minimizing the number of orbitals in the calculation. We find that, without embedding, i FCIQMC struggles to converge these calculations due to their substantial system sizes and a lack of error cancellation between reactants and products. With embedding, the i FCIQMC calculation converges straightforwardly to CCSD(T) benchmarks. Our results suggest that embedded i FCIQMC will be able treat system sizes well beyond our current reach (even though embedding introduces an error). We discuss how embedding might be improved (and thus the introduced error reduced) using i FCIQMC energies as benchmarks.


On the validity of memristor modeling in the neural network literature. An analysis of the literature shows that there are two types of non memristive models that have been widely used in the modeling of so called "memristive" neural networks. Here, we demonstrate that such models have nothing in common with the concept of memristive elements they describe either non linear resistors or certain bi state systems, which all are devices without memory. Therefore, the results presented in a significant number of publications are at least questionable, if not completely irrelevant to the actual field of memristive neural networks.


Near optimal sampling strategies for multivariate function approximation on general domains. In this paper, we address the problem of approximating a multivariate function defined on a general domain in d dimensions from sample points. We consider weighted least squares approximation in an arbitrary finite dimensional space P from independent random samples taken according to a suitable measure. In general, least squares approximations can be inaccurate and ill conditioned when the number of sample points M is close to N dim(P) . To counteract this, we introduce a novel method for sampling in general domains which leads to provably accurate and well conditioned approximations. The resulting sampling measure is discrete, and therefore straightforward to sample from. Our main result shows near optimal sample complexity for this procedure specifically, M mathcal O (N log(N)) samples suffice for a well conditioned and accurate approximation. Numerical experiments on polynomial approximation in general domains confirm the benefits of this method over standard sampling.


Unsupervised Surgical Instrument Segmentation via Anchor Generation and Semantic Diffusion. Surgical instrument segmentation is a key component in developing context aware operating rooms. Existing works on this task heavily rely on the supervision of a large amount of labeled data, which involve laborious and expensive human efforts. In contrast, a more affordable unsupervised approach is developed in this paper. To train our model, we first generate anchors as pseudo labels for instruments and background tissues respectively by fusing coarse handcrafted cues. Then a semantic diffusion loss is proposed to resolve the ambiguity in the generated anchors via the feature correlation between adjacent video frames. In the experiments on the binary instrument segmentation task of the 2017 MICCAI EndoVis Robotic Instrument Segmentation Challenge dataset, the proposed method achieves 0.71 IoU and 0.81 Dice score without using a single manual annotation, which is promising to show the potential of unsupervised learning for surgical tool segmentation.


The Origins and Prevalence of Texture Bias in Convolutional Neural Networks. Recent work has indicated that, unlike humans, ImageNet trained CNNs tend to classify images by texture rather than by shape. How pervasive is this bias, and where does it come from? We find that, when trained on datasets of images with conflicting shape and texture, CNNs learn to classify by shape at least as easily as by texture. What factors, then, produce the texture bias in CNNs trained on ImageNet? Different unsupervised training objectives and different architectures have small but significant and largely independent effects on the level of texture bias. However, all objectives and architectures still lead to models that make texture based classification decisions a majority of the time, even if shape information is decodable from their hidden representations. The effect of data augmentation is much larger. By taking less aggressive random crops at training time and applying simple, naturalistic augmentation (color distortion, noise, and blur), we train models that classify ambiguous images by shape a majority of the time, and outperform baselines on out of distribution test sets. Our results indicate that apparent differences in the way humans and ImageNet trained CNNs process images may arise not primarily from differences in their internal workings, but from differences in the data that they see.


ProvMark A Provenance Expressiveness Benchmarking System. System level provenance is of widespread interest for applications such as security enforcement and information protection. However, testing the correctness or completeness of provenance capture tools is challenging and currently done manually. In some cases there is not even a clear consensus about what behavior is correct. We present an automated tool, ProvMark, that uses an existing provenance system as a black box and reliably identifies the provenance graph structure recorded for a given activity, by a reduction to subgraph isomorphism problems handled by an external solver. ProvMark is a beginning step in the much needed area of testing and comparing the expressiveness of provenance systems. We demonstrate ProvMark's usefuless in comparing three capture systems with different architectures and distinct design philosophies.


Surface magnetic structure investigation of a nanolaminated Mn 2 GaC thin film using a magnetic field microscope based on Nitrogen Vacancy centers. This work presents a magnetic field imaging method based on color centers in diamond crystal applied to a thin film of a nanolaminated Mn 2 GaC MAX phase. Magnetic properties of the surface related structures have been described around the first order transition at 214 K by performing measurements in the temperature range between 200 K and 235 K with the surface features fading out by increasing temperature above the transition temperature. The results presented here demonstrate how Nitrogen Vacancy center based magnetic microscopy can supplement the traditionally used set of experimental techniques, giving additional information of microscopic scale magnetic field features, and allowing to investigate the temperature dependent magnetic behavior. The additional information acquired in this way is relevant to applications where surface magnetic properties are of essence.


Network memory in the movement of hospital patients carrying drug resistant bacteria. Hospitals constitute highly interconnected systems that bring into contact an abundance of infectious pathogens and susceptible individuals, thus making infection outbreaks both common and challenging. In recent years, there has been a sharp incidence of antimicrobial resistance amongst healthcare associated infections, a situation now considered endemic in many countries. Here we present network based analyses of a data set capturing the movement of patients harbouring drug resistant bacteria across three large London hospitals. We show that there are substantial memory effects in the movement of hospital patients colonised with drug resistant bacteria. Such memory effects break first order Markovian transitive assumptions and substantially alter the conclusions from the analysis, specifically on node rankings and the evolution of diffusive processes. We capture variable length memory effects by constructing a lumped state memory network, which we then use to identify overlapping communities of wards. We find that these communities of wards display a quasi hierarchical structure at different levels of granularity which is consistent with different aspects of patient flows related to hospital locations and medical specialties.


Approximation and Convergence of Large Atomic Congestion Games. We consider the question of whether, and in what sense, Wardrop equilibria provide a good approximation for Nash equilibria in atomic unsplittable congestion games with a large number of small players. We examine two different definitions of small players. In the first setting, we consider games where each player's weight is small. We prove that when the number of players goes to infinity and their weights to zero, the random flows in all (mixed) Nash equilibria for the finite games converge in distribution to the set of Wardrop equilibria of the corresponding nonatomic limit game. In the second setting, we consider an increasing number of players with a unit weight that participate in the game with a decreasingly small probability. In this case, the Nash equilibrium flows converge in total variation towards Poisson random variables whose expected values are Wardrop equilibria of a different nonatomic game with suitably defined costs. The latter can be viewed as symmetric equilibria in a Poisson game in the sense of Myerson, establishing a plausible connection between the Wardrop model for routing games and the stochastic fluctuations observed in real traffic. In both settings we provide explicit approximation bounds, and we study the convergence of the price of anarchy. Beyond the case of congestion games, we prove a general result on the convergence of large games with random players towards Poisson games.


Self Training Autonomous Driving Agent. Intrinsically, driving is a Markov Decision Process which suits well the reinforcement learning paradigm. In this paper, we propose a novel agent which learns to drive a vehicle without any human assistance. We use the concept of reinforcement learning and evolutionary strategies to train our agent in a 2D simulation environment. Our model's architecture goes beyond the World Model's by introducing difference images in the auto encoder. This novel involvement of difference images in the auto encoder gives better representation of the latent space with respect to the motion of vehicle and helps an autonomous agent to learn more efficiently how to drive a vehicle. Results show that our method requires fewer (96 less) total agents, (87.5 less) agents per generations, (70 less) generations and (90 less) rollouts than the original architecture while achieving the same accuracy of the original.


On Rayner structures. In this note, we study substructures of generalised power series fields induced by families of well ordered subsets of the group of exponents. We characterise the set theoretic and algebraic properties of the induced substructures in terms of conditions on the families. We extend the work of Rayner by giving both emph necessary and sufficient conditions to obtain truncation closed subgroups, subrings and subfields.


Hardware In the Loop for Connected Automated Vehicles Testing in Real Traffic. We present a hardware in the loop (HIL) simulation setup for repeatable testing of Connected Automated Vehicles (CAVs) in dynamic, real world scenarios. Our goal is to test control and planning algorithms and their distributed implementation on the vehicle hardware and, possibly, in the cloud. The HIL setup combines PreScan for perception sensors, road topography, and signalized intersections Vissim for traffic micro simulation ETAS DESK LABCAR a dynamometer for vehicle and powertrain dynamics and on board electronic control units for CAV real time control. Models of traffic and signalized intersections are driven by real world measurements. To demonstrate this HIL simulation setup, we test a Model Predictive Control approach for maximizing energy efficiency of CAVs in urban environments.


Generalized L p norm joint inversion of gravity and magnetic data using cross gradient constraint. A generalized unifying approach for L p norm joint inversion of gravity and magnetic data using the cross gradient constraint is presented. The presented framework incorporates stabilizers that use L 0 , L 1 , and L 2 norms of the model parameters, and or the gradient of the model parameters. Furthermore, the formulation is developed from standard approaches for independent inversion of single data sets, and, thus, also facilitates the inclusion of necessary model and data weighting matrices that provide, for example, depth weighting and imposition of hard constraint data. The developed efficient algorithm can, therefore, be employed to provide physically relevant smooth, sparse, or blocky target(s) which are relevant to the geophysical community. Here, the nonlinear objective function, that describes the inclusion of all stabilizing terms and the fit to data measurements, is minimized iteratively by imposing stationarity on the linear equation that results from applying linearization of the objective function about a starting model. To numerically solve the resulting linear system, at each iteration, the conjugate gradient algorithm is used. The general framework is then validated for three dimensional synthetic models for both sparse and smooth reconstructions, and the results are compared with those of individual gravity and magnetic inversions. It is demonstrated that the presented joint inversion algorithm is practical and significantly improves reconstructed models obtained by independent inversion.


Flat Sampled Data Systems, an Algorithmic Approach. Flatness of sampled data systems can be characterized by a simple property. They must admit the transformation to special representations, which are the series or partial series connection of a Brunovsky normal form and a complement. It is shown that this property follows from an integrability condition, which must be met by the flat parametrization. The series connections admit a simplification by reduction of the original problem to a simpler one, which allows us to develop two algorithms, where the first version delivers the flat outputs directly, but one has to solve linear PDEs or nonlinear ODEs. But one gets the flat outputs directly. The second version overcomes this problem, but one can only test the existence of flat outputs.


Overfullness of critical class 2 graphs with a small core degree. Let G be a simple graph, and let n , Delta(G) and chi' (G) be the order, the maximum degree and the chromatic index of G , respectively. We call G overfull if E(G) lfloor n 2 rfloor Delta(G) , and critical if chi'(H) chi'(G) for every proper subgraph H of G . Clearly, if G is overfull then chi'(G) Delta(G) 1 . The core of G , denoted by G Delta , is the subgraph of G induced by all its maximum degree vertices. Hilton and Zhao conjectured that for any critical class 2 graph G with Delta(G) ge 4 , if the maximum degree of G Delta is at most two, then G is overfull, which in turn gives Delta(G) n 2 1 . We show that for any critical class 2 graph G , if the minimum degree of G Delta is at most two and Delta(G) n 2 1 , then G is overfull.


A Stock Options Metaphor for Content Delivery Networks. The concept of Stock Options is used to address the scarcity of resources, not adequately addressed by the previous tools of our Prediction Mechanism. Using a Predictive Reservation Scheme, network and disk resources are being monitored through well established techniques (Kernel Regression Estimators) in a given time frame. Next, an Secondary Market mechanism significantly improves the efficiency and robustness of our Predictive Reservation Scheme by allowing the fast exchange of unused (remaining) resources between the Origin Servers (CDN Clients). This exchange can happen, either by implementing socially optimal practices or by allowing automatic electronic auctions at the end of the day or at shorter time intervals. Finally, we further enhance our Prediction Mechanism Stock Options are obtained and exercised, depending on the lack of resources at the end of day. As a result, Origin Servers may acquire resources (if required) at a normal price. The effectiveness of our mechanism further improves.


A simple proof of the L 'evy Khintchine formula for subordinators. We present a relatively simple and mostly elementary proof of the L 'evy Khintchine formula for subordinators. The main idea is to study the Poisson process time changed by the subordinator. The technical tools used are conditional expectations, probability generating function and convergence of discrete random variables.


Routing and Wavelength Assignment with Protection A Quadratic Unconstrained Binary Optimization Approach. The routing and wavelength assignment with protection is an important problem in telecommunications. Given an optical network and incoming connection requests, a commonly studied variant of the problem aims to grant maximum number of requests by assigning lightpaths at minimum network resource usage level, while ensuring the provided services remain functional in case of a single link failure through dedicated path protection. We consider a practically relevant version where alternative lightpaths for requests are assumed to be given as a precomputed set, and show that it is NP hard. We formulate the problem as an integer programming (IP) model, and also use it as a foundation to develop a novel quadratic unconstrained binary optimization (QUBO) model, which can be both directly solved by a state of the art solver like GUROBI. We present necessary and sufficient conditions on objective function parameters to prioritize request granting objective over wavelength link usage for both models, and a sufficient condition to ensure the exactness of the QUBO model. Moreover, we implement a problem specific branch and cut algorithm for the IP model, and employ a new quantum inspired technology, Digital Annealer (DA), for the QUBO model. We conduct computational experiments on a large suite of instances that are hard to optimally solve in order to assess the efficiency and efficacy of all of these approaches as well as a problem specific heuristic. The results show that the emerging technology DA outperforms the considered established techniques coupled with GUROBI, in finding mostly significantly better or as good solutions in only two minutes compared to two hours of run time, whereas the problem specific heuristic fails to be competitive.


On Learning Semantic Representations for Million Scale Free Hand Sketches. In this paper, we study learning semantic representations for million scale free hand sketches. This is highly challenging due to the domain unique traits of sketches, e.g., diverse, sparse, abstract, noisy. We propose a dual branch CNNRNN network architecture to represent sketches, which simultaneously encodes both the static and temporal patterns of sketch strokes. Based on this architecture, we further explore learning the sketch oriented semantic representations in two challenging yet practical settings, i.e., hashing retrieval and zero shot recognition on million scale sketches. Specifically, we use our dual branch architecture as a universal representation framework to design two sketch specific deep models (i) We propose a deep hashing model for sketch retrieval, where a novel hashing loss is specifically designed to accommodate both the abstract and messy traits of sketches. (ii) We propose a deep embedding model for sketch zero shot recognition, via collecting a large scale edge map dataset and proposing to extract a set of semantic vectors from edge maps as the semantic knowledge for sketch zero shot domain alignment. Both deep models are evaluated by comprehensive experiments on million scale sketches and outperform the state of the art competitors.


From the 1 2 3 Conjecture to the Riemann Hypothesis. This survey presents some combinatorial problems with number theoretic flavor. Our journey starts from a simple graph coloring question, but at some point gets close to a dangerous territory of the Riemann Hypothesis. We will mostly focus on open problems, but we will also provide some simple proofs, just for adorning.


Bounds for rm GL 2 times GL 2 L functions in depth aspect. Let f and g be holomorphic or Maass cusp forms for rm SL 2( mathbb Z ) and let chi be a primitive Dirichlet character of prime power conductor mathfrak q p kappa with p prime and kappa 12 . A subconvex bound for the central values of the Rankin Selberg L functions L(s,f otimes g otimes chi) is proved in the depth aspect L left( frac 1 2 ,f otimes g otimes chi right) ll f,g, varepsilon p 3 4 mathfrak q 15 16 varepsilon .


On the permutative equivalence of squares of unconditional bases. We prove that if the squares of two unconditional bases are equivalent up to a permutation, then the bases themselves are permutatively equivalent. This settles a twenty year old question raised by Casazza and Kalton in Uniqueness of unconditional bases in Banach spaces, Israel J. Math. 103 (1998), 141 175 . Solving this problem provides a new paradigm to study the uniqueness of unconditional basis in the general framework of quasi Banach spaces. Multiple examples are given to illustrate how to put in practice this theoretical scheme. Among the main applications of this principle we obtain the uniqueness of unconditional basis up to permutation of finite sums of quasi Banach spaces with this property.


Vocabulary Learning via Optimal Transport for Neural Machine Translation. The choice of token vocabulary affects the performance of machine translation. This paper aims to figure out what is a good vocabulary and whether one can find the optimal vocabulary without trial training. To answer these questions, we first provide an alternative understanding of the role of vocabulary from the perspective of information theory. Motivated by this, we formulate the quest of vocabularization finding the best token dictionary with a proper size as an optimal transport (OT) problem. We propose VOLT, a simple and efficient solution without trial training. Empirical results show that VOLT outperforms widely used vocabularies in diverse scenarios, including WMT 14 English German and TED's 52 translation directions. For example, VOLT achieves almost 70 vocabulary size reduction and 0.5 BLEU gain on English German translation. Also, compared to BPE search, VOLT reduces the search time from 384 GPU hours to 30 GPU hours on English German translation. Codes are available at https github.com Jingjing NLP VOLT .


Asymptotic upper bound for tangential speed of parabolic semigroups of holomorphic self maps in the unit disc. We show that the tangential speed of a parabolic semigroup of holomorphic self maps in the unit disc is asymptotically bounded from above by (1 2)logt, proving a conjecture by Bracci. In order to show the proof we need a result of "asymptotical monotonicity" of the tangential speed for proper pairs of parabolic semigroups with positive hyperbolic step.


Rigidity for Bach flat metrics on manifolds with boundary and applications. In the article we consider Bach flat metrics on four manifolds with boundary, with conformally invariant boundary conditions. We show that such metrics arise naturally as critical points of the Weyl energy under a constraint. We then prove a rigidity result if a Yamabe metric associated to a critical metric when restricted to the boundary is isometric to the round three sphere, then the critical metric must be isometric to the standard upper hemisphere.


Modal Group Velocity Mismatch Induced Intermodal Modulation Instability in Step index Fiber. We present detailed experimental study on noise seeded intermodal modulation instability (IM MI) in normal dispersion region of a conventional step index fiber. The sharp refractive index contrast between core and cladding leads to large group velocity mismatch between the spatial modes, coaxing to efficient IM MI and generation of multiple spectral peaks along with Raman peaks. Evolution of the spectrum with pump powers and fiber lengths are observed. Experimental findings are well supported with the theoretical framework based on bimodal MI model considering the distinct dispersion parameters of the participating modes.


A family of multi parameterized proximal point algorithms. In this paper, a multi parameterized proximal point algorithm combining with a relaxation step is developed for solving convex minimization problem subject to linear constraints. We show its global convergence and sublinear convergence rate from the prospective of variational inequality. Preliminary numerical experiments on testing a sparse minimization problem from signal processing indicate that the proposed algorithm performs better than some well established methods


Mathematical aspects relative to the fluid statics of a self gravitating perfect gas isothermal sphere. In the present paper we analyze and discuss some mathematical aspects of the fluid static configurations of a self gravitating perfect gas enclosed in a spherical solid shell. The mathematical model we consider is based on the well known Lane Emden equation, albeit under boundary conditions that differ from those usually assumed in the astrophysical literature. The existence of multiple solutions requires particular attention in devising appropriate numerical schemes apt to deal with and catch the solution multiplicity as efficiently and accurately as possible. In sequence, we describe some analytical properties of the model, the two algorithms used to obtain numerical solutions, and the numerical results for two selected cases.


GeoRocket A scalable and cloud based data store for big geospatial files. We present GeoRocket, a software for the management of very large geospatial datasets in the cloud. GeoRocket employs a novel way to handle arbitrarily large datasets by splitting them into chunks that are processed individually. The software has a modern reactive architecture and makes use of existing services including Elasticsearch and storage back ends such as MongoDB or Amazon S3. GeoRocket is schema agnostic and supports a wide range of heterogeneous geospatial file formats. It is also format preserving and does not alter imported data in any way. The main benefits of GeoRocket are its performance, scalability, and usability, which make it suitable for a number of scientific and commercial use cases dealing with very high data volumes, complex datasets, and high velocity (Big Data). GeoRocket also provides many opportunities for further research in the area of geospatial data management.


Knowledge Distillation for Multi task Learning. Multi task learning (MTL) is to learn one single model that performs multiple tasks for achieving good performance on all tasks and lower cost on computation. Learning such a model requires to jointly optimize losses of a set of tasks with different difficulty levels, magnitudes, and characteristics (e.g. cross entropy, Euclidean loss), leading to the imbalance problem in multi task learning. To address the imbalance problem, we propose a knowledge distillation based method in this work. We first learn a task specific model for each task. We then learn the multi task model for minimizing task specific loss and for producing the same feature with task specific models. As the task specific network encodes different features, we introduce small task specific adaptors to project multi task features to the task specific features. In this way, the adaptors align the task specific feature and the multi task feature, which enables a balanced parameter sharing across tasks. Extensive experimental results demonstrate that our method can optimize a multi task learning model in a more balanced way and achieve better overall performance.


Generating Question Answer Hierarchies. The process of knowledge acquisition can be viewed as a question answer game between a student and a teacher in which the student typically starts by asking broad, open ended questions before drilling down into specifics (Hintikka, 1981 Hakkarainen and Sintonen, 2002). This pedagogical perspective motivates a new way of representing documents. In this paper, we present SQUASH (Specificity controlled Question Answer Hierarchies), a novel and challenging text generation task that converts an input document into a hierarchy of question answer pairs. Users can click on high level questions (e.g., "Why did Frodo leave the Fellowship?") to reveal related but more specific questions (e.g., "Who did Frodo leave with?"). Using a question taxonomy loosely based on Lehnert (1978), we classify questions in existing reading comprehension datasets as either "general" or "specific". We then use these labels as input to a pipelined system centered around a conditional neural language model. We extensively evaluate the quality of the generated QA hierarchies through crowdsourced experiments and report strong empirical results.


Using an Automated Heterogeneous Robotic System for Radiation Surveys. During missions involving radiation exposure, unmanned robotic platforms may embody a valuable tool, especially thanks to their capability of replacing human operators in certain tasks to eliminate the health risks associated with such an environment. Moreover, rapid development of the technology allows us to increase the automation rate, making the human operator generally less important within the entire process. This article presents a multi robotic system designed for highly automated radiation mapping and source localization. Our approach includes a three phase procedure comprising sequential deployment of two diverse platforms, namely, an unmanned aircraft system (UAS) and an unmanned ground vehicle (UGV), to perform aerial photogrammetry, aerial radiation mapping, and terrestrial radiation mapping. The central idea is to produce a sparse dose rate map of the entire study site via the UAS and, subsequently, to perform detailed UGV based mapping in limited radiation contaminated regions. To accomplish these tasks, we designed numerous methods and data processing algorithms to facilitate, for example, digital elevation model (DEM) based terrain following for the UAS, automatic selection of the regions of interest, obstacle map based UGV trajectory planning, and source localization. The overall usability of the multi robotic system was demonstrated by means of a one day, authentic experiment, namely, a fictitious car accident including the loss of several radiation sources. The ability of the system to localize radiation hotspots and individual sources has been verified.


On the minimality of Keplerian arcs with fixed negative energy. We revisit a classical result by Jacobi on the local minimality, as critical points of the corresponding energy functional, of fixed energy solutions of the Kepler equation joining two distinct points with the same distance from the origin. Our proof relies on the Morse index theorem, together with a characterization of the conjugate points as points of geodesic bifurcation.


Context Aware Personality Inference in Dyadic Scenarios Introducing the UDIVA Dataset. This paper introduces UDIVA, a new non acted dataset of face to face dyadic interactions, where interlocutors perform competitive and collaborative tasks with different behavior elicitation and cognitive workload. The dataset consists of 90.5 hours of dyadic interactions among 147 participants distributed in 188 sessions, recorded using multiple audiovisual and physiological sensors. Currently, it includes sociodemographic, self and peer reported personality, internal state, and relationship profiling from participants. As an initial analysis on UDIVA, we propose a transformer based method for self reported personality inference in dyadic scenarios, which uses audiovisual data and different sources of context from both interlocutors to regress a target person's personality traits. Preliminary results from an incremental study show consistent improvements when using all available context information.


A pressure robust staggered discontinuous Galerkin method for the Stokes equations. In this paper we propose a pressure robust staggered discontinuous Galerkin method for the Stokes equations on general polygonal meshes by using piecewise constant approximations. We modify the right hand side of the body force in the discrete formulation by exploiting divergence preserving velocity reconstruction operator, which is the crux for pressure independent velocity error estimates. The optimal convergence for velocity gradient, velocity and pressure are proved. In addition, we are able to prove the superconvergence of velocity approximation by the incorporation of divergence preserving velocity reconstruction operator in the dual problem, which is also an important contribution of this paper. Finally, several numerical experiments are carried out to confirm the theoretical findings.


Timed Automata Benchmark Description. This report contains the descriptions of the timed automata (models) and the properties (specifications) that are used as the "benchmark examples in Data structure choices for on the fly model checking of real time systems" and "The power of proofs New algorithms for timed automata model checking." The four models from those sources are CSMA, FISCHER, LEADER, and GRC. Additionally we include in this report two additional models FDDI and PATHOS. These six models are often used to benchmark timed automata model checker speed throughout timed automata model checking papers.


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


Phase lag predicts nonlinear response maxima in liquid sloshing experiments. Mass spring models are essential for the description of sloshing resonances in engineering. By experimentally measuring the liquid's centre of mass in a horizontally oscillated rectangular tank, we show that low amplitude sloshing obeys the Duffing equation. A bending of the response curve in analogy to a softening spring is observed, with growing hysteresis as the driving amplitude increases. At large amplitudes, complex wave patterns emerge (including wave breaking and run up at the tank walls), competition between flow states is observed and the dynamics departs progressively from Duffing. We also provide a quantitative comparison of wave shapes and response curves to the predictions of a multimodal model based on potential flow theory (Faltinsen Timokha 2009) and show that it systematically overestimates the sloshing amplitudes and the hysteresis. We find that the phase lag between the liquid's centre of mass and the forcing is the key predictor of the nonlinear response maxima. The phase lag reflects precisely the onset of deviations from Duffing dynamics and most importantly at resonance the sloshing motion always lags the driving by 90 deg (independently of the wave pattern). This confirms the theoretical 90 deg phase lag criterion (Cenedese Haller 2020).


Exploring Content Selection in Summarization of Novel Chapters. We present a new summarization task, generating summaries of novel chapters using summary chapter pairs from online study guides. This is a harder task than the news summarization task, given the chapter length as well as the extreme paraphrasing and generalization found in the summaries. We focus on extractive summarization, which requires the creation of a gold standard set of extractive summaries. We present a new metric for aligning reference summary sentences with chapter sentences to create gold extracts and also experiment with different alignment methods. Our experiments demonstrate significant improvement over prior alignment approaches for our task as shown through automatic metrics and a crowd sourced pyramid analysis. We make our data collection scripts available at https github.com manestay novel chapter dataset .


Exploring Longitudinal Effects of Session based Recommendations. Session based recommendation is a problem setting where the task of a recommender system is to make suitable item suggestions based only on a few observed user interactions in an ongoing session. The lack of long term preference information about individual users in such settings usually results in a limited level of personalization, where a small set of popular items may be recommended to many users. This repeated exposure of such a subset of the items through the recommendations may in turn lead to a reinforcement effect over time, and to a system which is not able to help users discover new content anymore to the desirable extent. In this work, we investigate such potential longitudinal effects of session based recommendations in a simulation based approach. Specifically, we analyze to what extent algorithms of different types may lead to concentration effects over time. Our experiments in the music domain reveal that all investigated algorithms both neural and heuristic ones may lead to lower item coverage and to a higher concentration on a subset of the items. Additional simulation experiments however also indicate that relatively simple re ranking strategies, e.g., by avoiding too many repeated recommendations in the music domain, may help to deal with this problem.


A Robust Biped Locomotion Based on Linear Quadratic Gaussian Controller and Divergent Component of Motion. Generating robust locomotion for a humanoid robot in the presence of disturbances is difficult because of its high number of degrees of freedom and its unstable nature. In this paper, we used the concept of Divergent Component of Motion (DCM) and propose an optimal closed loop controller based on Linear Quadratic Gaussian to generate a robust and stable walking for humanoid robots. The biped robot dynamics has been approximated using the Linear Inverted Pendulum Model (LIPM). Moreover, we propose a controller to adjust the landing location of the swing leg to increase the withstanding level of the robot against a severe external push. The performance and also the robustness of the proposed controller is analyzed and verified by performing a set of simulations using mbox MATLAB . The simulation results showed that the proposed controller is capable of providing a robust walking even in the presence of disturbances and in challenging situations.


A Kingman like Theorem. We provide a Kingman like Theorem for arbitrary finite measures and a version of Birkhoff's Theorem for bounded observable. As an application, we show that Birkhoff's limit exists for some continuous observable, in an example of Bowen.


Measuring Impact of Climate Change on Tree Species analysis of JSDM on FIA data. One of the first beings affected by changes in the climate are trees, one of our most vital resources. In this study tree species interaction and the response to climate in different ecological environments is observed by applying a joint species distribution model to different ecological domains in the United States. Joint species distribution models are useful to learn inter species relationships and species response to the environment. The climates' impact on the tree species is measured through species abundance in an area. We compare the model's performance across all ecological domains and study the sensitivity of the climate variables. With the prediction of abundances, tree species populations can be predicted in the future and measure the impact of climate change on tree populations.


Fields of dimension one algebraic over a global or local field need not be of type C 1 . Let (K, v) be a Henselian discrete valued field with a quasifinite residue field. This paper proves the existence of an algebraic extension E K satisfying the following (i) E has dimension dim (E) le 1 , i.e. the Brauer group Br (E prime ) is trivial, for every algebraic extension E prime E (ii) finite extensions of E are not C 1 fields. This, applied to the maximal algebraic extension K of the field mathbb Q of rational numbers in the field mathbb Q p of p adic numbers, for a given prime p , proves the existence of an algebraic extension E p mathbb Q , such that dim (E p ) le 1 , E p is not a C 1 field, and E p has a Henselian valuation of residual characteristic p .


On the importance of system view centric validation for the design and operation of a crypto based digital economy. Ubiquitous connectivity, networked computation, open technologies and advances in intelligent web approaches enable the third generation web, Web3.0. Not least due to advancements in crypto technologies adoption in real world applications, it is now possible to convert the internet of things (IoT) into an economy of things (EoT), which basically refers to a heterogeneous digital economy of everything (dEoE). Based on the realization of the deficiencies of Web2.0 and the loss of confidence in a central intermediaries based economy, a societal desire to overcome these deficiencies has emerged. In the beginning, the crypto movement has been enthusiastic about the ability to encode all necessities for a digital economy in algorithms and code. However, in the mean time, insight has gained ground that on the one hand, designing the EoT is a quite complex endeavor and the resulting system might fail to achieve the targeted ideals, even when no intermediaries are present. We emphasize this position by arguing that the EoT in fact corresponds to a complex system operated in open contexts, for which a sophisticated design approach and a set of measures are inevitable, in order to be able to achieve a valid design and operation. In addition, we argue that the high level guiding principles necessary for a valid and effective EoT align well with the ideals of the democratized Web3.0 movement, and are no naive rapture, but rather adequate adoption provided can lead to socioeconomic efficiency. Over and above, we postulate that in this kind of EoT, commercial companies not only play a role, but can benefit alongside the society. In more game theoretical jargon, one could state that it is rational for the society and enterprises to undertake the effort to establish a valid EoT based on these ideals.


Online information on medical cannabis may rise unrealistic expectations and downplay potential side effects. There is a growing literature on the potential medical uses of Cannabis sativa and cannabinoid compounds. Although these have only been approved by regulatory agencies for few indications, there is a hype about their possible benefits in a variety of conditions and a large market in the wellness industry. As in many cases patients search for information on cannabis products online, we have analyzed the information on medical cannabis available on the Internet. Analyzing 176 webpages returned by a search engine, we found that more than half of them were news websites. Pain, epilepsy and multiple sclerosis were the most frequently therapeutic areas mentioned by the webpages, which did not always match those for which there is regulatory approval. Information was also incomplete, with only 22 of the webpages mentioning potential side effects. Health portal websites provided the most complete information. On average, 80 of webpages had a neutral stance on the potential benefits of medical cannabis, with commercial websites having more frequently a positive stance (67 ). We conclude that the information that can be found online could raise unrealistic expectations regarding therapeutic areas for which science based evidence is often still weak.


Convection under internal waves in a mountain lake. Turbulent mixing processes in deep alpine Lake Garda (I) have not extensively been observed. Knowledge about drivers of turbulent fluxes are important for insights in the transport of matter, nutrients and pollutants, in the lake and in natural water bodies in general. In this paper, the occurrence of internal wave induced turbulent convection is addressed as opposed to the more common shear induced turbulence in a density stratified environment. Observations are analyzed from a dedicated yearlong mooring holding 100 high resolution temperature sensors at 1.5 m intervals under a single current meter in the deeper half of the 344 m deep lake center. Episodically, the weakly density stratified waters in the lower 50 m above the lake floor show spectral slope and coherence evidence of short term (15 to 30 minutes) convective motions under internal waves that are supported by the stronger stratified waters above. The near homogeneous conditions are not attributable to frictional Ekman dynamics, but to large scale internal wave crests.


Fluorescence detected Fourier transform electronic spectroscopy by phase tagged photon counting. Fluorescence detected Fourier transform (FT) spectroscopy is a technique in which the relative paths of an optical interferometer are controlled to excite a material sample, and the ensuing fluorescence is detected as a function of the interferometer path delay and relative phase. A common approach to enhance the signal to noise ratio in these experiments is to apply a continuous phase sweep to the relative optical path, and to detect the resulting modulated fluorescence using a phase sensitive lock in amplifier. In many important situations, the fluorescence signal is too weak to be measured using a lock in amplifier, so that photon counting techniques are preferred. Here we introduce an approach to low signal fluorescence detected FT spectroscopy, in which individual photon counts are assigned to a modulated interferometer phase ('phase tagged photon counting,' or PTPC), and the resulting data are processed to construct optical spectra. We studied the fluorescence signals of a molecular sample excited resonantly by a pulsed coherent laser over a range of photon flux and visibility levels. We compare the performance of PTPC to standard lock in detection methods and establish the range of signal parameters over which meaningful measurements can be carried out. We find that PTPC generally outperforms the lock in detection method, with the dominant source of measurement uncertainty being associated with the statistics of the finite number of samples of the photon detection rate.


Strong laws of large numbers for a growth fragmentation process with bounded cell sizes. Growth fragmentation processes model systems of cells that grow continuously over time and then fragment into smaller pieces. Typically, on average, the number of cells in the system exhibits asynchronous exponential growth and, upon compensating for this, the distribution of cell sizes converges to an asymptotic profile. However, the long term stochastic behaviour of the system is more delicate, and its almost sure asymptotics have been so far largely unexplored. In this article, we study a growth fragmentation process whose cell sizes are bounded above, and prove the existence of regimes with differing almost sure long term behaviour.


Monogenic trinomials with non squarefree discriminant. For each integer n ge 2 , we identify new infinite families of monogenic trinomials f(x) x n Ax m B with non squarefree discriminant, many of which have small Galois group. These families are thus different from many previous examinations of specific trinomial forms in the literature. Moreover, in certain situations when A B ge 2 with fixed n and m , we produce asymptotics on the number of such trinomials with A le X .


A Fast Dictionary Learning Method for Coupled Feature Space Learning. In this letter, we propose a novel computationally efficient coupled dictionary learning method that enforces pairwise correlation between the atoms of dictionaries learned to represent the underlying feature spaces of two different representations of the same signals, e.g., representations in different modalities or representations of the same signals measured with different qualities. The jointly learned correlated feature spaces represented by coupled dictionaries are used in sparse representation based classification, recognition and reconstruction tasks. The presented experimental results show that the proposed coupled dictionary learning method has a significantly lower computational cost. Moreover, the visual presentation of jointly learned dictionaries shows that the pairwise correlations between the corresponding atoms are ensured.


Sparse Solutions for Inverse Problems in Reproducing Kernel Hilbert Spaces. A linear operator in a Hilbert space defined through the form of Riesz representation naturally introduces a reproducing kernel Hilbert space structure over the range space. Such formulation, called H HK formulation in this paper, possesses a built in mechanism to solve some basic type problems in the formulation by using the basis method, that include identification of the range space, the inversion problem, and the Moore Penrose pseudo (generalized) inversion problem. After a quick survey of the existing theory, the aim of the article is to establish connection between this formulation with sparse series representation, and in particular with one called pre orthogonal adaptive Fourier decomposition (POAFD), the latter being one, most recent and well developed, with great efficiency and wide and deep connections with traditional analysis. Within the matching pursuit methodology the optimality of POAFD is theoretically guaranteed. In practice POAFD offers fast converging numerical solutions.


The mobile phone as a free rotation laboratory. Modern mobile phones contain a three axis microelectromechanical system (MEMS) gyroscope, capable of taking accurate measurements of the angular velocity along the three principal axes of the phone with a sampling rate of 100 Hz or better. If the phone is tossed in the air, then, neglecting air resistance, it is in free rotation (rotation in the absence of a torque) with respect to its centre of mass, and the phone's gyroscope can be used to record the rotational dynamics. This enables experimental investigation of free rotation. In this paper, we use a mobile phone to demonstrate the steady states for rotation of the phone about two of its principal axes, and the instability in rotation about the third corresponding to the intermediate moment of inertia. We also show the approximate conservation of angular momentum and rotational kinetic energy during motion in the air, and compare the data with numerical solution of Euler's equations for free rotation. Our results demonstrate the capability of smartphones for investigating free rotation, and should be of interest to college and university teachers developing "at home" physics labs for remote learning.


Renewal reward perspective on linear switching diffusion systems. In many biological systems, the movement of individual agents is commonly characterized as having multiple qualitatively distinct behaviors that arise from various biophysical states. This is true for vesicles in intracellular transport, micro organisms like bacteria, or animals moving within and responding to their environment. For example, in cells the movement of vesicles, organelles and other cargo are affected by their binding to and unbinding from cytoskeletal filaments such as microtubules through molecular motor proteins. A typical goal of theoretical or numerical analysis of models of such systems is to investigate the effective transport properties and their dependence on model parameters. While the effective velocity of particles undergoing switching diffusion is often easily characterized in terms of the long time fraction of time that particles spend in each state, the calculation of the effective diffusivity is more complicated because it cannot be expressed simply in terms of a statistical average of the particle transport state at one moment of time. However, it is common that these systems are regenerative, in the sense that they can be decomposed into independent cycles marked by returns to a base state. Using decompositions of this kind, we calculate effective transport properties by computing the moments of the dynamics within each cycle and then applying renewal reward theory. This method provides a useful alternative large time analysis to direct homogenization for linear advection reaction diffusion partial differential equation models. Moreover, it applies to a general class of semi Markov processes and certain stochastic differential equations that arise in models of intracellular transport. Applications of the proposed framework are illustrated for case studies such as mRNA transport in developing oocytes and processive cargo movement by teams of motor proteins.


Shuffle algebra realization of quantum affine superalgebra U v ( hat mathfrak D (2,1 theta)) . We give shuffle algebra realization of positive part of quantum affine superalgebra U v ( widehat mathfrak D (2,1 theta)) associated to any simple root systems. We also determine the shuffle algebra associated to widehat mathfrak sl (2 1) with odd root system when v is a primitive root of unity of even order, generalizing results in cite FJMMT03 .


Interfacing Graphene Based Materials With Neural Cells. The scientific community has witnessed an exponential increase in the applications of graphene and graphene based materials in a wide range of fields. For what concerns neuroscience, the interest raised by these materials is two fold. On one side, nanosheets made of graphene or graphene derivatives (graphene oxide, or its reduced form) can be used as carriers for drug delivery. Here, an important aspect is to evaluate their toxicity, which strongly depends on flake composition, chemical functionalization and dimensions. On the other side, graphene can be exploited as a substrate for tissue engineering. In this case, conductivity is probably the most relevant amongst the various properties of the different graphene materials, as it may allow to instruct and interrogate neural networks, as well as to drive neural growth and differentiation. In this review, we try to give a comprehensive view of the accomplishments and new challenges of the field, as well as which in our view are the most exciting directions to take in the immediate future. These include the need to engineer multifunctional nanoparticles able to cross the blood brain barrier to reach neural cells, and to achieve on demand delivery of specific drugs. We describe the state of the art in the use of graphene materials to engineer three dimensional scaffolds to drive neuronal growth and regeneration in vivo, and the possibility of using graphene as a component of hybrid composites multi layer organic electronics devices. Last but not least, we address the need of an accurate theoretical modeling of the interface between graphene and biological material, by modeling the interaction of graphene with proteins and cell membranes at the nanoscale, and describing the physical mechanism(s) of charge transfer by which the various graphene materials can influence the excitability and physiology of neural cells.


Nonlinear inviscid damping for a class of monotone shear flows in finite channel. We prove the nonlinear inviscid damping for a class of monotone shear flows in T times 0,1 for initial perturbation in Gevrey 1 s ( s 2 ) class with compact support. The main idea of the proof is to use the wave operator of a slightly modified Rayleigh operator in a well chosen coordinate system.


KeepAugment A Simple Information Preserving Data Augmentation Approach. Data augmentation (DA) is an essential technique for training state of the art deep learning systems. In this paper, we empirically show data augmentation might introduce noisy augmented examples and consequently hurt the performance on unaugmented data during inference. To alleviate this issue, we propose a simple yet highly effective approach, dubbed emph KeepAugment , to increase augmented images fidelity. The idea is first to use the saliency map to detect important regions on the original images and then preserve these informative regions during augmentation. This information preserving strategy allows us to generate more faithful training examples. Empirically, we demonstrate our method significantly improves on a number of prior art data augmentation schemes, e.g. AutoAugment, Cutout, random erasing, achieving promising results on image classification, semi supervised image classification, multi view multi camera tracking and object detection.


Towards better substitution based word sense induction. Word sense induction (WSI) is the task of unsupervised clustering of word usages within a sentence to distinguish senses. Recent work obtain strong results by clustering lexical substitutes derived from pre trained RNN language models (ELMo). Adapting the method to BERT improves the scores even further. We extend the previous method to support a dynamic rather than a fixed number of clusters as supported by other prominent methods, and propose a method for interpreting the resulting clusters by associating them with their most informative substitutes. We then perform extensive error analysis revealing the remaining sources of errors in the WSI task. Our code is available at https github.com asafamr bertwsi.


Depth from a polarisation RGB stereo pair. In this paper, we propose a hybrid depth imaging system in which a polarisation camera is augmented by a second image from a standard digital camera. For this modest increase in equipment complexity over conventional shape from polarisation, we obtain a number of benefits that enable us to overcome longstanding problems with the polarisation shape cue. The stereo cue provides a depth map which, although coarse, is metrically accurate. This is used as a guide surface for disambiguation of the polarisation surface normal estimates using a higher order graphical model. In turn, these are used to estimate diffuse albedo. By extending a previous shape from polarisation method to the perspective case, we show how to compute dense, detailed maps of absolute depth, while retaining a linear formulation. We show that our hybrid method is able to recover dense 3D geometry that is superior to state of the art shape from polarisation or two view stereo alone.


A KAM Theorem for finitely differentiable Hamiltonian systems. Given l 2 nu 2d geq 4 , we prove the persistence of a Cantor family of KAM tori of measure O( varepsilon 1 2 nu l ) for any non degenerate nearly integrable Hamiltonian system of class C l( mathscr D times mathbb T d) , where mathscr D subset mathbb R d is a bounded domain, provided that the size varepsilon of the perturbation is sufficiently small. This extends a result by D. Salamon in cite salamon2004kolmogorov according to which we do have the persistence of a single KAM torus in the same framework. Moreover, it is well known that, for the persistence of a single torus, the regularity assumption can not be improved.


Evidence for a multi level trophic organization of the human gut microbiome. The human gut microbiome is a complex ecosystem, in which hundreds of microbial species and metabolites coexist, in part due to an extensive network of cross feeding interactions. However, both the large scale trophic organization of this ecosystem, and its effects on the underlying metabolic flow, remain unexplored. Here, using a simplified model, we provide quantitative support for a multi level trophic organization of the human gut microbiome, where microbes consume and secrete metabolites in multiple iterative steps. Using a manually curated set of metabolic interactions between microbes, our model suggests about four trophic levels, each characterized by a high level to level metabolic transfer of byproducts. It also quantitatively predicts the typical metabolic environment of the gut (fecal metabolome) in approximate agreement with the real data. To understand the consequences of this trophic organization, we quantify the metabolic flow and biomass distribution, and explore patterns of microbial and metabolic diversity in different levels. The hierarchical trophic organization suggested by our model can help mechanistically establish causal links between the abundances of microbes and metabolites in the human gut.


Orbital angular momentum tuning using a phase only parallel aligned LCoS display. We propose two schemes for achieving continuous and complete orbital angular momentum tuning using a simple device based on a phase only parallel aligned liquid crystal on silicon display, where the tuning is managed by changing the input beam state of polarization. We use an optical tweezers toolbox to simulate the optical force field of the tuned beams applied to dielectric spherical particles in order to prove the applicability of the method.


Rinascimento Optimising Statistical Forward Planning Agents for Playing Splendor. Game based benchmarks have been playing an essential role in the development of Artificial Intelligence (AI) techniques. Providing diverse challenges is crucial to push research toward innovation and understanding in modern techniques. Rinascimento provides a parameterised partially observable multiplayer card based board game, these parameters can easily modify the rules, objectives and items in the game. We describe the framework in all its features and the game playing challenge providing baseline game playing AIs and analysis of their skills. We reserve to agents' hyper parameter tuning a central role in the experiments highlighting how it can heavily influence the performance. The base line agents contain several additional contribution to Statistical Forward Planning algorithms.


AnonymousNet Natural Face De Identification with Measurable Privacy. With billions of personal images being generated from social media and cameras of all sorts on a daily basis, security and privacy are unprecedentedly challenged. Although extensive attempts have been made, existing face image de identification techniques are either insufficient in photo reality or incapable of balancing privacy and usability qualitatively and quantitatively, i.e., they fail to answer counterfactual questions such as "is it private now?", "how private is it?", and "can it be more private?" In this paper, we propose a novel framework called AnonymousNet, with an effort to address these issues systematically, balance usability, and enhance privacy in a natural and measurable manner. The framework encompasses four stages facial attribute estimation, privacy metric oriented face obfuscation, directed natural image synthesis, and adversarial perturbation. Not only do we achieve the state of the arts in terms of image quality and attribute prediction accuracy, we are also the first to show that facial privacy is measurable, can be factorized, and accordingly be manipulated in a photo realistic fashion to fulfill different requirements and application scenarios. Experiments further demonstrate the effectiveness of the proposed framework.


Global Nash Kuiper theorem for compact manifolds. We obtain global extensions of the celebrated Nash Kuiper theorem for C 1, theta isometric immersions of compact manifolds with optimal H "older exponent. In particular for the Weyl problem of isometrically embedding a convex compact surface in 3 space, we show that the Nash Kuiper non rigidity prevails upto exponent theta 1 5 . This extends previous results on embedding 2 discs as well as higher dimensional analogues.


Notes on k rainbow independent domination in graphs. The k rainbow independent domination number of a graph G , denoted gamma rm rik (G) , is the cardinality of a smallest set consisting of two vertex disjoint independent sets V 1 and V 2 for which every vertex in V(G) setminus (V 1 cup V 2) has neighbors in both V 1 and V 2 . This domination invariant was proposed by v S umenjak, Rall and Tepeh in (Applied Mathematics and Computation 333(15), 2018 353 361), which allows to reduce the problem of computing the independent domination number of the generalized prism G Box K k to an integer labeling problem on G . They proved a Nordhaus Gaddum type theorem 5 leq gamma rm rik (G) gamma rm rik ( overline G ) leq n 3 for every graph G of order n geq 3 , where overline G is the complement of G . In this paper, we improve this result by showing that if G is not isomorphic to the 5 cycle, then 5 leq gamma rm rik (G) gamma rm rik ( overline G ) leq n 2 . Moreover, we show that the problem of deciding whether a graph has a k rainbow independent dominating function of a given weight is mathcal NP complete. Our results respond some open questions proposed by v S umenjak, et al.


A Variational Analysis Approach to Solving the Merton Problem. We address the Merton problem of maximizing the expected utility of terminal wealth using techniques from variational analysis. Under a general continuous semimartingale market model with stochastic parameters, we obtain a characterization of the optimal portfolio for general utility functions in terms of a forward backward stochastic differential equation (FBSDE) and derive solutions for a number of well known utility functions. Our results complement a previous studies conducted on optimal strategies in markets driven by Brownian noise with random drift and volatility parameters.


Intellectual Property Protection for Deep Learning Models Taxonomy, Methods, Attacks, and Evaluations. The training and creation of deep learning model is usually costly, thus it can be regarded as an intellectual property (IP) of the model creator. However, malicious users who obtain high performance models may illegally copy, redistribute, or abuse the models without permission. To deal with such security threats, a few deep neural networks (DNN) IP protection methods have been proposed in recent years. This paper attempts to provide a review of the existing DNN IP protection works and also an outlook. First, we propose the first taxonomy for DNN IP protection methods in terms of six attributes scenario, mechanism, capacity, type, function, and target models. Then, we present a survey on existing DNN IP protection works in terms of the above six attributes, especially focusing on the challenges these methods face, whether these methods can provide proactive protection, and their resistances to different levels of attacks. After that, we analyze the potential attacks on DNN IP protection methods from the aspects of model modifications, evasion attacks, and active attacks. Besides, a systematic evaluation method for DNN IP protection methods with respect to basic functional metrics, attack resistance metrics, and customized metrics for different application scenarios is given. Lastly, future research opportunities and challenges on DNN IP protection are presented.


Interpreting Neural Ranking Models using Grad CAM. Recently, applying deep neural networks in IR has become an important and timely topic. For instance, Neural Ranking Models(NRMs) have shown promising performance compared to the traditional ranking models. However, explaining the ranking results has become even more difficult with NRM due to the complex structure of neural networks. On the other hand, a great deal of research is under progress on Interpretable Machine Learning(IML), including Grad CAM. Grad CAM is an attribution method and it can visualize the input regions that contribute to the network's output. In this paper, we adopt Grad CAM for interpreting the ranking results of NRM. By adopting Grad CAM, we analyze how each query document term pair contributes to the matching score for a given pair of query and document. The visualization results provide insights on why a certain document is relevant to the given query. Also, the results show that neural ranking model captures the subtle notion of relevance. Our interpretation method and visualization results can be used for snippet generation and user query intent analysis.


Scaling in internally heated convection a unifying theory. We offer a unifying theory for turbulent purely internally heated convection, generalizing the unifying theories of Grossmann and Lohse (2000, 2001) for Rayleigh B 'enard turbulence and of Shishkina, Grossmann and Lohse (2016) for turbulent horizontal convection, which are both based on the splitting of the kinetic and thermal dissipation rates in respective boundary and bulk contributions. We obtain the mean temperature of the system and the Reynolds number (which are the response parameters) as function of the control parameters, namely the internal thermal driving strength (called, when nondimensionalized, the Rayleigh Roberts number) and the Prandtl number. The results of the theory are consistent with our direct numerical simulations.