title: Simulation of fiber-reinforced viscoelastic structures subjected to finite strains: multiplicative approach

abstract: The study is devoted to the geometrically nonlinear simulation of fiber-reinforced composite structures. The applicability of the multiplicative approach to the simulation of viscoelastic properties of a composite material is assessed, certain improvements are suggested. For a greater accuracy in applications involving local compressive fiber buckling, a new family of hyperelastic potentials is introduced. This family allows us to account for the variable critical compressive stress, which depends on the fiber-matrix interaction. For the simulation of viscoelasticity, the well-established Sidoroff decomposition of the deformation gradient is implemented. To account for the viscosity of the matrix material, the model of Simo and Miehe (1992) is used; highly efficient iteration-free algorithms are implemented. The viscosity of the fiber is likewise described by the multiplicative decomposition of the deformation gradient, leading to a scalar differential equation; an efficient iteration-free algorithm is proposed for the implicit time stepping. The accuracy and convergence of the new iteration-free method is tested and compared to that of the standard scheme implementing the Newton iteration. To demonstrate the applicability of the approach, a pressurized multi-layer composite pipe is modelled; the so-called stretch inversion phenomenon is reproduced and explained. The stress distribution is obtained by a semi-analytical procedure; it may serve as a benchmark for FEM computations. Finally, the issue of the parameter identification is addressed.


title: An Aiding Tool for Building Design Generation, Thermal Assessment and Optimization -- EnergyPlus Interaction Overview

abstract: A building design aiding tool for space allocation and thermal performance optimization is being developed to help practitioners during the building space planning phase, predicting how it will behave regarding energy consumption and thermal comfort. The tool evaluates, ranks, and optimizes generated floor plans according to thermal performance criteria, using the dynamic simulation program EnergyPlus. The tool is currently able to use a wide variety of EnergyPlus objects, allowing for various template and detailed HVAC, DHW, and thermal and electrical energy production systems and components, as well as numerous internal gains types, construction elements and energy saving controls, to be accounted for and simulated in the generated buildings. This paper presents the tool overall concept as well as the main features regarding dynamic simulation. Some performance results are presented for distinct systems to illustrate the use and potential of the tool.


title: Hyper Space Exploration - A Multicriterial Quantitative Trade-Off Analysis for System Design in Complex Environment

abstract: Successful engineering requires environmentally adapted procedural and architectural approaches. While dealing with complicated issues has become an engineering standard mastering uncertainties in complex environment is still a major issue. Global trends, such as an increasing rate of disruptive technology changes or merging of technology fields, however, enforce the importance of this complex habitat. Missing experience in unknown technological territory faces engineers with two questions of paramount importance: 1) How can the best architectural solution within the space of potential alternatives be identified? 2) How can a proof-of-concept for considered solutions prior implementation be provided? Mastering lack of knowledge related risks and uncertainties states one of the most prominent tasks in according projects. The paper presents a novel methodology of system design in complex environment called Hyper Space Exploration (HSE). The HSE approach combines methods of virtual prototyping with those of design of virtual experiments based studies for statistical learning. Virtual prototyping allows an early feedback on system behavior with a proof-of-concept prior implementation. Statistical learning enables system architects to systematically build up the space of potential solution alternatives, model the effects of design and use-case variables on target indicators in complex territory, quantify target indicator trade-offs and finally identify Pareto-optimal system solutions. The first part of the paper characterizes engineering challenges in complex environment. Section two presents the HSE methodology with its two major constituents work flow and tool chain. Section three outlines first successful HSE applications that already proved its capabilities and universality. Final section four gives an outlook to further HSE applications as well as methodological future HSE extensions.


title: A goodness of fit test for the Pareto distribution

abstract: The Zenga (1984) inequality curve is constant in p for Type I Pareto distributions. This characterizing behavior will be exploited to obtain graphical and analytical tools for tail analysis and goodness of fit tests. A testing procedure for Pareto-type behavior based on a regression of technique will be introduced.


title: TTHRESH: Tensor Compression for Multidimensional Visual Data

abstract: Memory and network bandwidth are decisive bottlenecks when handling high-resolution multidimensional data sets in visualization applications, and they increasingly demand suitable data compression strategies. We introduce a novel lossy compression algorithm for multidimensional data over regular grids. It leverages the higher-order singular value decomposition (HOSVD), a generalization of the SVD to three dimensions and higher, together with bit-plane, run-length and arithmetic coding to compress the HOSVD transform coefficients. Our scheme degrades the data particularly smoothly and achieves lower mean squared error than other state-of-the-art algorithms at low-to-medium bit rates, as it is required in data archiving and management for visualization purposes. Further advantages of the proposed algorithm include very fine bit rate selection granularity and the ability to manipulate data at very small cost in the compression domain, for example to reconstruct filtered and/or subsampled versions of all (or selected parts) of the data set.


title: Controllable Semantic Image Inpainting

abstract: We develop a method for user-controllable semantic image inpainting: Given an arbitrary set of observed pixels, the unobserved pixels can be imputed in a user-controllable range of possibilities, each of which is semantically coherent and locally consistent with the observed pixels. We achieve this using a deep generative model bringing together: an encoder which can encode an arbitrary set of observed pixels, latent variables which are trained to represent disentangled factors of variations, and a bidirectional PixelCNN model. We experimentally demonstrate that our method can generate plausible inpainting results matching the user-specified semantics, but is still coherent with observed pixels. We justify our choices of architecture and training regime through more experiments.


title: Study of external electron injection and trapping in the blow-out regime

abstract: In this work we study electron side-injection and trapping in the blow-out regime in deep plasma channels. We analyze the maximum angle of injection, for which at least 90\% of the injected electrons are trapped. We discuss the dependence of this angle on the electrons' initial energy and their injection positions. In the scope of a semi-analytical blow-out model we show that the injection position is a less critical factor for trapping if electrons are injected into deep plasma channels. PIC simulations and analytical approximations support our results from the semi-analytical model.


title: Distribution Laws of Smooth Divisors

abstract: A classical result due to Deshouillers, Dress and Tenenbaum asserts that on average the distribution of the divisors of the integers follows the arcsine law. In this paper, we investigate the distribution of smooth divisors of the integers, that is, those divisors which are free of large prime factors. We show that on average these divisors are distributed according to a probability law that we will describe.


title: Flow Logic

abstract: Flow networks have attracted a lot of research in computer science. Indeed, many questions in numerous application areas can be reduced to questions about flow networks. Many of these applications would benefit from a framework in which one can formally reason about properties of flow networks that go beyond their maximal flow. We introduce Flow Logics: modal logics that treat flow functions as explicit first-order objects and enable the specification of rich properties of flow networks. The syntax of our logic BFL* (Branching Flow Logic) is similar to the syntax of the temporal logic CTL*, except that atomic assertions may be flow propositions, like $> \gamma$ or $\geq \gamma$, for $\gamma \in \mathbb{N}$, which refer to the value of the flow in a vertex, and that first-order quantification can be applied both to paths and to flow functions. We present an exhaustive study of the theoretical and practical aspects of BFL*, as well as extensions and fragments of it. Our extensions include flow quantifications that range over non-integral flow functions or over maximal flow functions, path quantification that ranges over paths along which non-zero flow travels, past operators, and first-order quantification of flow values. We focus on the model-checking problem and show that it is PSPACE-complete, as it is for CTL*. Handling of flow quantifiers, however, increases the complexity in terms of the network to ${\rm P}^{\rm NP}$, even for the LFL and BFL fragments, which are the flow-counterparts of LTL and CTL. We are still able to point to a useful fragment of BFL* for which the model-checking problem can be solved in polynomial time. Finally, we introduce and study the query-checking problem for BFL*, where under-specified BFL* formulas are used for network exploration.


title: Computing several eigenvalues of nonlinear eigenvalue problems by selection

abstract: Computing more than one eigenvalue for (large sparse) one-parameter polynomial and general nonlinear eigenproblems, as well as for multiparameter linear and nonlinear eigenproblems, is a much harder task than for standard eigenvalue problems. We present simple but efficient selection methods based on divided differences to do this. In contrast to locking techniques, it is not necessary to keep converged eigenvectors in the search space, so that the entire search space may be devoted to new information. The techniques are applicable to many types of matrix eigenvalue problems; standard deflation is possible only for linear one-parameter problems. The methods are easy to understand and implement. Although divided differences are well-known in the context of nonlinear eigenproblems, the proposed selection techniques are new for one-parameter problems. For multiparameter problems, we improve on and generalize our previous work. We also show how to use divided differences in the framework of homogeneous coordinates, which may be appropriate for generalized eigenvalue problems with infinite eigenvalues. While the approaches are valuable alternatives for one-parameter nonlinear eigenproblems, they seem the only option for multiparameter problems.


title: Higher order spacing ratios in random matrix theory and complex quantum systems

abstract: The distribution of the ratios of nearest neighbor level spacings has become a popular indicator of spectral fluctuations in complex quantum systems like interacting many-body localized and thermalization phases, quantum chaotic systems, and also in atomic and nuclear physics. In contrast to the level spacing distribution, which requires the cumbersome and at times ambiguous unfolding procedure, the ratios of spacings do not require unfolding and are easier to compute. In this work, for the class of Wigner-Dyson random matrices with nearest neighbor spacing ratios $r$ distributed as $P_{\beta}(r)$ for the three ensembles indexed by $\beta=1,2, 4$, their $k-$th order spacing ratio distributions are shown to be identical to $P_{\beta'}(r)$, where $\beta'$, an integer, is a function of $\beta$ and $k$. This result is shown for Gaussian and circular ensembles of random matrix theory and for several physical systems such as spin chains, chaotic billiards, Floquet systems and measured nuclear resonances.

categories: quant-ph cond-mat.dis-nn math-ph math.MP nlin.CD physics.data-an


title: The Bass diffusion model on finite Barabasi-Albert networks

abstract: Using a mean-field network formulation of the Bass innovation diffusion model and exact results by Fotouhi and Rabbat on the degree correlations of Barabasi-Albert networks, we compute the times of the diffusion peak and compare them with those on scale-free networks which have the same scale-free exponent but different assortativity properties. We compare our results with those obtained by Caldarelli et al. for the SIS epidemic model with the spectral method applied to adjacency matrices. It turns out that diffusion times on finite Barabasi-Albert networks are at a minimum. This may be due to a little-known property of these networks: although the value of the assortativity coefficient is close to zero, they look disassortative if one considers only a bounded range of degrees, including the smallest ones, and slightly assortative on the range of the higher degrees. We also find that if the trickle-down character of the diffusion process is enhanced by a larger initial stimulus on the hubs (via a inhomogeneous linear term in the Bass model), the relative difference between the diffusion times for BA networks and uncorrelated networks is even larger, reaching for instance the 34% in a typical case on a network with $10^4$ nodes.


title: Quantum computing with classical bits

abstract: A bit-quantum map relates probabilistic information for Ising spins or classical bits to quantum spins or qubits. Quantum systems are subsystems of classical statistical systems. The Ising spins can represent macroscopic two-level observables, and the quantum subsystem employs suitable expectation values and correlations. We discuss static memory materials based on Ising spins for which boundary information can be transported through the bulk in a generalized equilibrium state. They can realize quantum operations as the Hadamard or CNOT-gate for the quantum subsystem. Classical probabilistic systems can account for the entanglement of quantum spins. An arbitrary unitary evolution for an arbitrary number of quantum spins can be described by static memory materials for an infinite number of Ising spins which may, in turn, correspond to continuous variables or fields. We discuss discrete subsets of unitary operations realized by a finite number of Ising spins. They may be useful for new computational structures. We suggest that features of quantum computation or more general probabilistic computation may be realized by neural networks, neuromorphic computing or perhaps even the brain. This does neither require small isolated entities nor very low temperature. In these systems the processing of probabilistic information can be more general than for static memory materials. We propose a general formalism for probabilistic computing for which deterministic computing and quantum computing are special limiting cases.

categories: quant-ph cond-mat.stat-mech cs.IT hep-th math.IT


title: Field-free Magnetization Switching by Utilizing the Spin Hall Effect and Interlayer Exchange Coupling of Iridium

abstract: Magnetization switching by spin-orbit torque (SOT) via spin Hall effect represents as a competitive alternative to that by spin-transfer torque (STT) used for magnetoresistive random access memory (MRAM), as it does not require high-density current to go through the tunnel junction. For perpendicular MRAM, however, SOT driven switching of the free layer requires an external in-plane field, which poses limitation for viability in practical applications. Here we demonstrate field-free magnetization switching of a perpendicular magnet by utilizing an Iridium (Ir) layer. The Ir layer not only provides SOTs via spin Hall effect, but also induce interlayer exchange coupling with an in-plane magnetic layer that eliminates the need for the external field. Such dual functions of the Ir layer allows future build-up of magnetoresistive stacks for memory and logic applications. Experimental observations show that the SOT driven field-free magnetization reversal is characterized as domain nucleation and expansion. Micromagnetic modeling is carried out to provide in-depth understanding of the perpendicular magnetization reversal process in the presence of an in-plane exchange coupling field.


title: A characterization of linearized polynomials with maximum kernel

abstract: We provide sufficient and necessary conditions for the coefficients of a $q$-polynomial $f$ over $\mathbb{F}_{q^n}$ which ensure that the number of distinct roots of $f$ in $\mathbb{F}_{q^n}$ equals the degree of $f$. We say that these polynomials have maximum kernel. As an application we study in detail $q$-polynomials of degree $q^{n-2}$ over $\mathbb{F}_{q^n}$ which have maximum kernel and for $n\leq 6$ we list all $q$-polynomials with maximum kernel. We also obtain information on the splitting field of an arbitrary $q$-polynomial. Analogous results are proved for $q^s$-polynomials as well, where $\gcd(s,n)=1$.


title: Thin Reaction Zones in Highly Turbulent Medium

abstract: In a highly turbulent medium characterized by a low Damk\"ohler number Da, reactions are commonly considered to occur in distributed zones broadened by small-scale turbulent eddies. In the present communication, an alternative regime of propagation of reaction waves in a highly turbulent medium is introduced and studied theoretically and numerically. More specifically, propagation of an infinitely thin reaction sheet in a turbulent medium is analyzed, with molecular mixing of the reactant and product being allowed in wide layers. In this limiting case, an increase in the consumption velocity by turbulence is solely controlled by an increase in the reaction-sheet area. Based on physical reasoning and estimates, the area is hypothesized to be close to the mean area of an inert iso-scalar surface at the same turbulent Reynolds number. This hypothesis leads to a relation for the turbulent consumption velocity, which is similar to the well-known Damk\"ohler scaling associated commonly with distributed reaction zones at a low Da. The obtained theoretical results are validated by analyzing a big database (23 cases characterized by 0.01<Da<1) created recently in 3D direct numerical simulations of propagation of a statistically planar, one-dimensional, dynamically passive reaction wave in statistically stationary, homogeneous, isotropic turbulence. The DNS data well support the aforementioned relation. They also show that the reaction is localized to thin zones even at Da as low as 0.01, with a ratio of the turbulent and laminar consumption velocities being mainly controlled by the reaction-zone-surface area.


title: From probabilistic graphical models to generalized tensor networks for supervised learning

abstract: Tensor networks have found a wide use in a variety of applications in physics and computer science, recently leading to both theoretical insights as well as practical algorithms in machine learning. In this work we explore the connection between tensor networks and probabilistic graphical models, and show that it motivates the definition of generalized tensor networks where information from a tensor can be copied and reused in other parts of the network. We discuss the relationship between generalized tensor network architectures used in quantum physics, such as string-bond states, and architectures commonly used in machine learning. We provide an algorithm to train these networks in a supervised-learning context and show that they overcome the limitations of regular tensor networks in higher dimensions, while keeping the computation efficient. A method to combine neural networks and tensor networks as part of a common deep learning architecture is also introduced. We benchmark our algorithm for several generalized tensor network architectures on the task of classifying images and sounds, and show that they outperform previously introduced tensor-network algorithms. The models we consider also have a natural implementation on a quantum computer and may guide the development of near-term quantum machine learning architectures.


title: Transience and Recurrence of Markov Processes with Constrained Local Time

abstract: We study Markov processes conditioned so that their local time must grow slower than a prescribed function. Building upon recent work on Brownian motion with constrained local time in [5] and [33], we study transience and recurrence for a broad class of Markov processes. In order to understand the distribution of the local time, we determine the distribution of a non-decreasing L\'evy process (the inverse local time) conditioned to remain above a given level which varies in time. We study a time-dependent region, in contrast to previous works in which a process is conditioned to remain in a fixed region (e.g. [21,27]), so we must study boundary crossing probabilities for a family of curves, and thus obtain uniform asymptotics for such a family. Main results include necessary and sufficient conditions for transience or recurrence of the conditioned Markov process. We will explicitly determine the distribution of the inverse local time for the conditioned process, and in the transient case, we explicitly determine the law of the conditioned Markov process. In the recurrent case, we characterise the "entropic repulsion envelope" via necessary and sufficient conditions.


title: Cartan Geometry in Modal Homotopy Type Theory

abstract: This article constructs the moduli stack of torsionfree $G$-structures in homotopy type theory with one monadic modality. This yields a construction of this moduli stack for any $\infty$-topos equipped with any stable factorization systems. In the intended applications of this theory, the factorization systems are given by the deRham-Stack construction. Homotopy type theory allows a formulation of this abstract theory with surprising low complexity. This is witnessed by the accompanying formalization of large parts of this work.


title: Dark matter for $b\to s \mu^+ \mu^-$ anomaly in a gauged $U(1)_X$ model

abstract: We propose a new physics model which has a cold dark matter candidate and can explain the $b \to s \mu^+\mu^-$ anomaly at the same time. Our model includes a scalar quark $\widetilde{q}$ and a scalar lepton $\widetilde{\ell}$ which are $SU(2)_L$-doublet as well as a Dirac fermion $N$ which is $SU(2)_L$-singlet. The new particles are charged under a gauged $U(1)_X$ group which is spontaneously broken to a discrete $Z_2$ symmetry by a dark scalar $S$. The remnant $Z_2$ symmetry stabilizes the dark matter. Box diagrams with $\widetilde{q}$, $\widetilde{\ell}$, and $N$ running inside the loop can generate the correct Wilson coefficients $C_9^\mu = -C_{10}^\mu$ to accommodate the $b \to s \mu^+\mu^-$ anomaly while avoiding constraints such as $B_s-\overline{B}_s$ mixing. The dark matter annihilation into a second generation lepton pair via $t$-channel $\widetilde{\ell}$-exchanging process plays an important role in producing the current dark matter relic abundance of the universe, showing a strong interplay between the flavor and dark matter physics. We also discuss dark-gauge-interaction-dominated and Higgs-portal-dominated scenarios for dark matter physics.


title: A closed formula for the generating function of $p$-Bernoulli numbers: An elementary proof

abstract: For a two parameter family of Bernoulli numbers $B_{n, p}$ the exponential generating function is derived by elementary methods.


title: Magnetic tilting and emergent Majorana spin connection in topological superconductors

abstract: Due to the charge neutral and localized nature of surface Majorana modes, detection schemes usually rely on local spectroscopy or interference through the Josephson effect. Here, we theoretically study the magnetic response of a two-dimensional cone of Majorana fermions localized at the surface of class DIII Topological Superconductors. For a field parallel to the surface the Zeeman term vanishes and the orbital term induces a Doppler shift of the Andreev levels resulting in a tilting of the surface Majorana cone. For fields larger than a critical threshold field $H^*$ the system undergoes a transition from type I to type II Dirac-Majorana cone. In a spherical geometry the surface curvature leads to the emergence of the Majorana spin connection in the tilting term via an interplay between orbital and Zeeman, that generates a finite non-trivial coupling between negative and positive energy states. Majorana modes are thus expected to show a finite response to the applied field, that acquires a universal character in finite geometries and opens the way to detection of Majorana modes via time-dependent magnetic fields.


title: Elastic and diffractive scattering at the LHC

abstract: Inspired by the new TOTEM data on elastic $pp$ scattering at 13 TeV, we study the possibility to describe all the diffractive collider data ($\sigma_{\rm tot}, d\sigma_{\rm el}/dt$, $\rho\equiv$Re$A$/Im$A$, \sd) in a wide interval of energy (0.0625 to 13 TeV) in the framework of a two-channel eikonal model. We show that a satisfactory description can be achieved without an odd-signature (Odderon) exchange contribution. We consider the possible role of the QCD Odderon which may improve the description of $\rho$ and discuss the importance of the odd-signature term if the amplitude were to exceed the black disc limit.


title: An Enhanced Binary Particle-Swarm Optimization (E-BPSO) Algorithm for Service Placement in Hybrid Cloud Platforms

abstract: Nowadays, hybrid cloud platforms stand as an attractive solution for organizations intending to implement combined private and public cloud applications, in order to meet their profitability requirements. However, this can only be achieved through the utilization of available resources while speeding up execution processes. Accordingly, deploying new applications entails dedicating some of these processes to a private cloud solution, while allocating others to the public cloud. In this context, the present work is set to help minimize relevant costs and deliver effective choices for an optimal service placement solution within minimal execution time. Several evolutionary algorithms have been applied to solve the service placement problem and are used when dealing with complex solution spaces to provide an optimal placement and often produce a short execution time. The standard BPSO algorithm is found to display a significant disadvantage, namely, of easily trapping into local optima, in addition to demonstrating a noticeable lack of robustness in dealing with service placement problems. Hence, to overcome critical shortcomings associated with the standard BPSO, an Enhanced Binary Particle Swarm Optimization (E-BPSO) algorithm is proposed, consisting of a modification of the particle position updating equation, initially inspired from the continuous PSO. Our proposed E-BPSO algorithm is shown to outperform state-of-the-art approaches in terms of both cost and execution time, using a real benchmark.


title: Cancellable elements of the lattice of semigroup varieties: varieties satisfying a permutational identity of length 3

abstract: We completely determine all semigroup varieties satisfiyng a permutational identity of length 3 that are cancellable elements of the lattice of all semigroup varieties. Using this result, we provide a series of new examples of semigroup varieties that are a modular but not cancellable elements of this lattice.


title: Riesz bases associated with regular representations of semidirect product groups

abstract: This work is devoted to the study of Bessel and Riesz systems of the type $\big\{L_{\gamma}\mathsf{f}\big\}_{\gamma\in \Gamma}$ obtained from the action of the left regular representation $L_{\gamma}$ of a discrete non abelian group $\Gamma$ which is a semidirect product, on a function $\mathsf{f}\in \ell^2(\Gamma)$. The main features about these systems can be conveniently studied by means of a simple matrix-valued function $\mathbf{F}(\xi)$. These systems allow to derive sampling results in principal $\Gamma$-invariant spaces, i.e., spaces obtained from the action of the group $\Gamma$ on a element of a Hilbert space. Since the systems $\big\{L_{\gamma}\mathsf{f}\big\}_{\gamma\in \Gamma}$ are closely related to convolution operators, a connection with $C^*$-algebras is also established.


title: Boosted Training of Convolutional Neural Networks for Multi-Class Segmentation

abstract: Training deep neural networks on large and sparse datasets is still challenging and can require large amounts of computation and memory. In this work, we address the task of performing semantic segmentation on large volumetric data sets, such as CT scans. Our contribution is threefold: 1) We propose a boosted sampling scheme that uses a-posterior error maps, generated throughout training, to focus sampling on difficult regions, resulting in a more informative loss. This results in a significant training speed up and improves learning performance for image segmentation. 2) We propose a novel algorithm for boosting the SGD learning rate schedule by adaptively increasing and lowering the learning rate, avoiding the need for extensive hyperparameter tuning. 3) We show that our method is able to attain new state-of-the-art results on the VISCERAL Anatomy benchmark.


title: Structured Variational Learning of Bayesian Neural Networks with Horseshoe Priors

abstract: Bayesian Neural Networks (BNNs) have recently received increasing attention for their ability to provide well-calibrated posterior uncertainties. However, model selection---even choosing the number of nodes---remains an open question. Recent work has proposed the use of a horseshoe prior over node pre-activations of a Bayesian neural network, which effectively turns off nodes that do not help explain the data. In this work, we propose several modeling and inference advances that consistently improve the compactness of the model learned while maintaining predictive performance, especially in smaller-sample settings including reinforcement learning.


title: Holographic thermalization in AdS-Gauss-Bonnet gravity for small entangled regions

abstract: In this paper we study the propagation of entanglement entropy after a global instantaneous quench on the CFT boundary of AdS bulk. We consider the Gauss-Bonnet model as a higher curvature gravity model for which we correct the RT(HRT) proposal to compute the holographic entanglement entropy(HEE). To obtain an analytical solution we perform an approximation approach which bounds our computations to the small subregions and we compare its thermalization regimes to the result of large subsystem case. We can see tsunami picture where the evolution of entanglement breaks down for the large systems and so its details depends just on the shape and size of entangled region and also the used gravity model. We can see the phase transition in this regime is always continuous regardless the shape and size, in contrary with large subregions.


title: Understanding Complex Systems: From Networks to Optimal Higher-Order Models

abstract: To better understand the structure and function of complex systems, researchers often represent direct interactions between components in complex systems with networks, assuming that indirect influence between distant components can be modelled by paths. Such network models assume that actual paths are memoryless. That is, the way a path continues as it passes through a node does not depend on where it came from. Recent studies of data on actual paths in complex systems question this assumption and instead indicate that memory in paths does have considerable impact on central methods in network science. A growing research community working with so-called higher-order network models addresses this issue, seeking to take advantage of information that conventional network representations disregard. Here we summarise the progress in this area and outline remaining challenges calling for more research.

categories: physics.soc-ph cond-mat.dis-nn cs.SI physics.data-an


title: Uncertainty Estimations by Softplus normalization in Bayesian Convolutional Neural Networks with Variational Inference

abstract: We introduce a novel uncertainty estimation for classification tasks for Bayesian convolutional neural networks with variational inference. By normalizing the output of a Softplus function in the final layer, we estimate aleatoric and epistemic uncertainty in a coherent manner. The intractable posterior probability distributions over weights are inferred by Bayes by Backprop. Firstly, we demonstrate how this reliable variational inference method can serve as a fundamental construct for various network architectures. On multiple datasets in supervised learning settings (MNIST, CIFAR-10, CIFAR-100), this variational inference method achieves performances equivalent to frequentist inference in identical architectures, while the two desiderata, a measure for uncertainty and regularization are incorporated naturally. Secondly, we examine how our proposed measure for aleatoric and epistemic uncertainties is derived and validate it on the aforementioned datasets.


title: Substrate-induced magnetism and topological phase transition in silicene

abstract: Silicene has shown great application potential as a versatile material for nanoelectronics, particularly promising as building block for spintronic applications. Unfortunately, despite its intriguing properties, such as relatively large spin-orbit interactions, one of the biggest hurdles for silicene to be useful as a host spintronic material is the lack of magnetism or the topological phase transition owing to the silicene-substrate interactions, which influence its fundamental properties and has yet to be fully explored. Here, we show that when silicene is grown on CeO2 substrate, an appreciable robust magnetic moment appears in silicene covalently bonded to CeO2 (111), while a topological phase transition to a band insulator occurs regardless of van der Waals (vdWs) interaction or covalent bonding interaction at interface. The induced magnetism of silicene is due to the breaking of Si-Si {\pi}-bonding, also resulting in trivial topological phase. The silicene-substrate interaction, even weak vdWs force (equivalent to an electric field), can destroy quantum spin Hall effect (QSHE) of silicene. We propose a viable strategy --- constructing inverse symmetrical sandwich structure (protective layer/silicene/substrate) --- to preserve quantum spin Hall (QSH) state of silicene in weak vdWs interaction system. This work takes a critical step towards fundamental physics and realistic applications of silicene-based spintronic devices.


title: Isoscalar Giant Monopole, Dipole, and Quadrupole Resonances in $^{90,92}$Zr and $^{92}$Mo

abstract: The isoscalar giant monopole, dipole, and quadrupole strength distributions have been deduced in $^{90, 92}$Zr, and $^{92}$Mo from "background-free" spectra of inelastic $\alpha$-particle scattering at a beam energy of 385 MeV at extremely forward angles, including 0$^{\circ}$. These strength distributions were extracted by a multipole-decomposition analysis based on the expected angular distributions of the respective multipoles. All these strength distributions for the three nuclei practically coincide with each other, affirming that giant resonances, being collective phenomena, are not influenced by nuclear shell structure near $A\sim$90, contrary to the claim in a recent measurement.


title: Dynamical phase diagrams of a love capacity constrained prey-predator model

abstract: One interesting question in love relationships is: finally, what and when is the end of this love relationship? Using a prey-predator Verhulst-Lotka-Volterra (VLV) model we imply cooperation and competition tendency between people in order to describe a "love dilemma game". We select the most simple but immediately most complex case for studying the set of nonlinear differential equations, i.e. that implying three persons, being at the same time prey and predator. We describe four different scenarios in such a love game containing either a one-way love or a love triangle. Our results show that it is hard to love more than one person simultaneously. Moreover, to love several people simultaneously is an unstable state. We find some condition in which persons tend to have a friendly relationship and love someone in spite of their antagonistic interaction. We demonstrate the dynamics by displaying flow diagrams.


title: Accelerating delayed-acceptance Markov chain Monte Carlo algorithms

abstract: Delayed-acceptance Markov chain Monte Carlo (DA-MCMC) samples from a probability distribution via a two-stages version of the Metropolis-Hastings algorithm, by combining the target distribution with a "surrogate" (i.e. an approximate and computationally cheaper version) of said distribution. DA-MCMC accelerates MCMC sampling in complex applications, while still targeting the exact distribution. We design a computationally faster, albeit approximate, DA-MCMC algorithm. We consider parameter inference in a Bayesian setting where a surrogate likelihood function is introduced in the delayed-acceptance scheme. When the evaluation of the likelihood function is computationally intensive, our scheme produces a 2-4 times speed-up, compared to standard DA-MCMC. However, the acceleration is highly problem dependent. Inference results for the standard delayed-acceptance algorithm and our approximated version are similar, indicating that our algorithm can return reliable Bayesian inference. As a computationally intensive case study, we introduce a novel stochastic differential equation model for protein folding data.


title: Online Variant of Parcel Allocation in Last-mile Delivery

abstract: We investigate the problem of last-mile delivery, where a large pool of citizen crowd-workers are hired to perform a variety of location-specific urban logistics parcel delivering tasks. Current approaches focus on offline scenarios, where all the spatio temporal information of parcels and workers are given. However, the offline scenarios can be impractical since parcels and workers appear dynamically in real applications, and their information is unknown in advance. In this paper, in order to solve the shortcomings of the offline setting, we first formalize the online parcel allocation in last-mile delivery problem, where all parcels were put in pop-stations in advance, while workers arrive dynamically. Then we propose an algorithm which provides theoretical guarantee for the parcel allocation in last-mile delivery. Finally, we verify the effectiveness and efficiency of the proposed method through extensive experiments on real and synthetic datasets.


title: Ego-Lane Analysis System (ELAS): Dataset and Algorithms

abstract: Decreasing costs of vision sensors and advances in embedded hardware boosted lane related research detection, estimation, and tracking in the past two decades. The interest in this topic has increased even more with the demand for advanced driver assistance systems (ADAS) and self-driving cars. Although extensively studied independently, there is still need for studies that propose a combined solution for the multiple problems related to the ego-lane, such as lane departure warning (LDW), lane change detection, lane marking type (LMT) classification, road markings detection and classification, and detection of adjacent lanes (i.e., immediate left and right lanes) presence. In this paper, we propose a real-time Ego-Lane Analysis System (ELAS) capable of estimating ego-lane position, classifying LMTs and road markings, performing LDW and detecting lane change events. The proposed vision-based system works on a temporal sequence of images. Lane marking features are extracted in perspective and Inverse Perspective Mapping (IPM) images that are combined to increase robustness. The final estimated lane is modeled as a spline using a combination of methods (Hough lines with Kalman filter and spline with particle filter). Based on the estimated lane, all other events are detected. To validate ELAS and cover the lack of lane datasets in the literature, a new dataset with more than 20 different scenes (in more than 15,000 frames) and considering a variety of scenarios (urban road, highways, traffic, shadows, etc.) was created. The dataset was manually annotated and made publicly available to enable evaluation of several events that are of interest for the research community (i.e., lane estimation, change, and centering; road markings; intersections; LMTs; crosswalks and adjacent lanes). ELAS achieved high detection rates in all real-world events and proved to be ready for real-time applications.


title: Simulating Coulomb gases and log-gases with hybrid Monte Carlo algorithms

abstract: Coulomb and log-gases are exchangeable singular Boltzmann-Gibbs measures appearing in mathematical physics at many places, in particular in random matrix theory. We explore experimentally an efficient numerical method for simulating such gases. It is an instance of the Hybrid or Hamiltonian Monte Carlo algorithm, in other words a Metropolis-Hastings algorithm with proposals produced by a kinetic or underdamped Langevin dynamics. This algorithm has excellent numerical behavior despite the singular interaction, in particular when the number of particles gets large. It is more efficient than the well known overdamped version previously used for such problems, and allows new numerical explorations. It suggests for instance to conjecture a universality of the Gumbel fluctuation at the edge of beta Ginibre ensembles for all beta.


title: Dynamical resonance quench and Fano interference in spontaneous Raman scattering from quasiparticle and collective excitations

abstract: Time-resolved spontaneous Raman spectroscopy serves as a probe for incoherent quasiparticle and collective excitation dynamics, and allows to distinguish symmetry changes across a photoinduced phase transition through the inelastic light scattering selection rules. Largely unexplored is the role of the Raman resonance enhancement in the time-domain, and the transient interaction between scattering from quasiparticles and collective excitations, with the latter interaction leading to a Fano interference. In this work, we report on the observation of a phonon Raman tensor quench and Fano interference after strong photoexcitation of an intrinsic semiconductor. We observed a dynamic phonon scattering rate asymmetry and spectral asymmetry through simultaneous detection of both the anti-Stokes and Stokes response. The asymmetric phonon scattering rate is ascribed to the combined effect of the transient phonon population and the reduction of the phonon Raman tensor resulting from the photoexcited hole population. This same hole population results in a strong enhancement of the Fano spectral asymmetry.


title: Efficient adaptive multilevel stochastic Galerkin approximation using implicit a posteriori error estimation

abstract: Partial differential equations (PDEs) with inputs that depend on infinitely many parameters pose serious theoretical and computational challenges. Sophisticated numerical algorithms that automatically determine which parameters need to be activated in the approximation space in order to estimate a quantity of interest to a prescribed error tolerance are needed. For elliptic PDEs with parameter-dependent coefficients, stochastic Galerkin finite element methods (SGFEMs) have been well studied. Under certain assumptions, it can be shown that there exists a sequence of SGFEM approximation spaces for which the energy norm of the error decays to zero at a rate that is independent of the number of input parameters. However, it is not clear how to adaptively construct these spaces in a practical and computationally efficient way. We present a new adaptive SGFEM algorithm that tackles elliptic PDEs with parameter-dependent coefficients quickly and efficiently. We consider approximation spaces with a multilevel structure---where each solution mode is associated with a finite element space on a potentially different mesh---and use an implicit a posteriori error estimation strategy to steer the adaptive enrichment of the space. At each step, the components of the error estimator are used to assess the potential benefits of a variety of enrichment strategies, including whether or not to activate more parameters. No marking or tuning parameters are required. Numerical experiments for a selection of test problems demonstrate that the new method performs optimally in that it generates a sequence of approximations for which the estimated energy error decays to zero at the same rate as the error for the underlying finite element method applied to the associated parameter-free problem.


title: A single parameter Hermite-Pad\'e series representations for Ap\'ery's constant

abstract: Inspired by the results of Rhin and Viola (2001), the purpose of this work is to elaborate on a series representation for $\zeta \left( 3\right)$ which only depends on one single integer parameter. This is accomplished by deducing a Hermite-Pad\'e approximation problem using ideas of Sorokin (1998). As a consequence we get a new recurrence relation for the approximation of $\zeta(3)$ as well as a corresponding new continued fraction expansion for $\zeta(3)$, which do no reproduce Ap\'ery's phenomenon, i.e., though the approaches are different, they lead to the same sequence of diophantine approximations to $\zeta \left( 3\right) $. Finally, the convergence rates of several series representations of $\zeta(3)$ are compared.


title: \"Uber die Funktionen des Irrationalit\"atsma\ss es f\"{u}r zwei irrationalen Zahlen

abstract: For real $\xi$ we consider irrationality measure function $\psi_\xi (t) = \min_{1\le q \le t, \, q\in \mathbb{Z}} ||q\xi||$. We prove that in the case $\alpha \pm \beta \not\in \mathbb{Z}$ there exist arbitrary large values of $t$ with $|\psi_\alpha (t) -\psi_\beta(t)| \ge\left(\sqrt{\frac{\sqrt{5}+1}{2}}-1\right) \min (\psi_\alpha(t), \psi_\beta (t))$. This result is optimal.


title: Electrostatic effects, band distortions and superconductivity in twisted graphene bilayers

abstract: Bilayer graphene twisted by a small angle shows a significant charge modulation away from neutrality, as the charge in the narrow bands near the Dirac point is mostly localized in the regions of the Moir\'e pattern with $AA$ stacking. The resulting electrostatic potential gives rise to the dominant contribution to the electron-electron interaction within this low energy band, which becomes significantly distorted. The changes in the band structure are described by new local interactions, and lead to an assisted electron-hopping process. These couplings favor superconductivity at certain fillings.


title: Searching for Dark Matter with Paleo-Detectors

abstract: A large experimental program is underway to extend the sensitivity of direct detection experiments, searching for interaction of Dark Matter with nuclei, down to the neutrino floor. However, such experiments are becoming increasingly difficult and costly due to the large target masses and exquisite background rejection needed for the necessary improvements in sensitivity. We investigate an alternative approach to the detection of Dark Matter-nucleon interactions: Searching for the persistent traces left by Dark Matter scattering in ancient minerals obtained from much deeper than current underground laboratories. We estimate the sensitivity of paleo-detectors, which extends far beyond current upper limits for a wide range of Dark Matter masses. The sensitivity of our proposal also far exceeds the upper limits set by Snowden-Ifft et al. more than three decades ago using ancient Mica in an approach similar to paleo-detectors.


title: A Unique Connection for Born Geometry

abstract: It has been known for a while that the effective geometrical description of compactified strings on $d$-dimensional target spaces implies a generalization of geometry with a doubling of the sets of tangent space directions. This generalized geometry involves an $O(d,d)$ pairing $\eta$ and an $O(2d)$ generalized metric $\mathcal{H}$. More recently it has been shown that in order to include T-duality as an effective symmetry, the generalized geometry also needs to carry a phase space structure or more generally a para-Hermitian structure encoded into a skew-symmetric pairing $\omega$. The consistency of string dynamics requires this geometry to satisfy a set of compatibility relations that form what we call a Born geometry. In this work we prove an analogue of the fundamental theorem of Riemannian geometry for Born geometry. We show that there exists a unique connection which preserves the Born structure $(\eta,\omega,\mathcal{H})$ and which is torsionless in a generalized sense. This resolves a fundamental ambiguity that is present in the double field theory formulation of effective string dynamics.


title: Torsion groups of elliptic curves over quadratic fields $\mathbb{Q}(\sqrt{d}),$ $0<d<100$

abstract: We prove results towards classifying the possible torsion subgroups of elliptic curves over quadratic fields $\mathbb{Q}(\sqrt{d})$, where $0<d<100$ is a square-free integer, and obtain a complete classification for 49 out of 60 such fields. Over the remaining 11 quadratic fields, we cannot rule out the possibility of the group $\mathbb{Z}/16\mathbb{Z}$ appearing as a torsion group of an elliptic curve.


title: Energy conservation in self-phase modulation

abstract: Spectral broadening of ultrashort laser pulses is simultaneously described by either self-phase modulation (SPM) or four-wave mixing (FWM). The latter implies the instantaneous conservation of both the photon number and energy, while the former describes a time-dependent frequency shift, implying a violation of the energy conservation if the number of photons is to be conserved in each time slice. We resolve this paradox by considering the transient energy storage in the propagation medium, that can be calculated in the SPM formalism via the dephasing between the incident pulse and the medium polarization leading to an effective imaginary part in the third-order susceptibility. In parallel, considering the temporal variation of the incident intensity in FWM offsets the instantaneous frequency.


title: An improved cosmological parameter inference scheme motivated by deep learning

abstract: Dark matter cannot be observed directly, but its weak gravitational lensing slightly distorts the apparent shapes of background galaxies, making weak lensing one of the most promising probes of cosmology. Several observational studies have measured the effect, and there are currently running, and planned efforts to provide even larger, and higher resolution weak lensing maps. Due to nonlinearities on small scales, the traditional analysis with two-point statistics does not fully capture all the underlying information. Multiple inference methods were proposed to extract more details based on higher order statistics, peak statistics, Minkowski functionals and recently convolutional neural networks (CNN). Here we present an improved convolutional neural network that gives significantly better estimates of $\Omega_m$ and $\sigma_8$ cosmological parameters from simulated convergence maps than the state of art methods and also is free of systematic bias. We show that the network exploits information in the gradients around peaks, and with this insight, we construct a new, easy-to-understand, and robust peak counting algorithm based on the 'steepness' of peaks, instead of their heights. The proposed scheme is even more accurate than the neural network on high-resolution noiseless maps. With shape noise and lower resolution its relative advantage deteriorates, but it remains more accurate than peak counting.


title: Search for the decay of a Higgs boson in the $\ell\ell\gamma$ channel in proton-proton collisions at $\sqrt{s} =$ 13 TeV

abstract: A search for a Higgs boson decaying into a pair of electrons or muons and a photon is described. Higgs boson decays to a Z boson and a photon (H $\to$ Z$\gamma\to\ell\ell\gamma$, $\ell =$ e or $\mu$), or to two photons, one of which has an internal conversion into a muon pair (H $\to\gamma^{*}\gamma\to\mu\mu\gamma$) were considered. The analysis is performed using a data set recorded by the CMS experiment at the LHC from proton-proton collisions at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 35.9 fb$^{-1}$. No significant excess above the background prediction has been found. Limits are set on the cross section for a standard model Higgs boson decaying to opposite-sign electron or muon pairs and a photon. The observed limits on cross section times the corresponding branching fractions vary between 1.4 and 4.0 (6.1 and 11.4) times the standard model cross section for H $\to\gamma^{*}\gamma\to\mu\mu\gamma$ (H $\to$ Z$\gamma\to\ell\ell\gamma$) in the 120-130 GeV mass range of the $\ell\ell\gamma$ system. The H $\to\gamma^*\gamma\to\mu\mu\gamma$ and H $\to$ Z$\gamma\to\ell\ell\gamma$ analyses are combined for $m_\mathrm{H} =$ 125 GeV, obtaining an observed (expected) 95% confidence level upper limit of 3.9 (2.0) times the standard model cross section.


title: A Dataset for Building Code-Mixed Goal Oriented Conversation Systems

abstract: There is an increasing demand for goal-oriented conversation systems which can assist users in various day-to-day activities such as booking tickets, restaurant reservations, shopping, etc. Most of the existing datasets for building such conversation systems focus on monolingual conversations and there is hardly any work on multilingual and/or code-mixed conversations. Such datasets and systems thus do not cater to the multilingual regions of the world, such as India, where it is very common for people to speak more than one language and seamlessly switch between them resulting in code-mixed conversations. For example, a Hindi speaking user looking to book a restaurant would typically ask, "Kya tum is restaurant mein ek table book karne mein meri help karoge?" ("Can you help me in booking a table at this restaurant?"). To facilitate the development of such code-mixed conversation models, we build a goal-oriented dialog dataset containing code-mixed conversations. Specifically, we take the text from the DSTC2 restaurant reservation dataset and create code-mixed versions of it in Hindi-English, Bengali-English, Gujarati-English and Tamil-English. We also establish initial baselines on this dataset using existing state of the art models. This dataset along with our baseline implementations is made publicly available for research purposes.


title: Asymptotics for the number of zero drift reflectable walks in a Weyl chamber of type A

abstract: We study lattice walks in a Weyl chamber of type A with fixed or free end points. For lattice walk models with zero drift that may be counted by means of a reflection argument, we determine asymptotics for the number of such walks as their length tends to infinity. These models are equivalent to the lock step model and the random turns model of vicious walkers. As special cases, our main results include various asymptotic formulas found in the literature.


title: Mumford-Shah Mesh Processing using the Ambrosio-Tortorelli Functional

abstract: The Mumford-Shah functional approximates a function by a piecewise smooth function. Its versatility makes it ideal for tasks such as image segmentation or restoration, and it is now a widespread tool of image processing. Recent work has started to investigate its use for mesh segmentation and feature lines detection, but we take the stance that the power of this functional could reach far beyond these tasks and integrate the everyday mesh processing toolbox. In this paper, we discretize an Ambrosio-Tortorelli approximation via a Discrete Exterior Calculus formulation. We show that, combined with a new shape optimization routine, several mesh processing problems can be readily tackled within the same framework. In particular, we illustrate applications in mesh denoising, normal map embossing, mesh inpainting and mesh segmentation.


title: Fluctuations for block spin Ising models

abstract: We analyze the high temperature fluctuations of the magnetization of the so-called Ising block model. This model was recently introduced by Berthet, Rigollet and Srivastava. We prove a Central Limit Theorems (CLT) for the magnetization in the high temperature regime. At the same time we show that this CLT breaks down at a line of critical temperatures. At this line we show the validity of a non-standard Central Limit Theorems for the magnetization.


title: Limited Information Strategies and Discrete Selectivity

abstract: We relate the property of discrete selectivity and its corresponding game, both recently introduced by V.V. Tkachuck, to a variety of selection principles and point picking games. In particular we show that player II can win the discrete selection game on \(C_p(X)\) if and only if player II can win a variant of the point open game on \(X\). We also show that the existence of limited information strategies in the discrete selection game on \(C_p(X)\) for either player are equivalent to other well-known topological properties.


title: A multi-band analogue frequency sensor with sub-MHz resolution based on a Vortex Nano-Oscillator

abstract: As the Internet of Things (IoT) becomes more integral to how we live our lives and society becomes ever increasingly connected, the need for efficient nanoscale radio frequency detection is becoming more of a pressing issue. The commonly used technique of a superheterodyne conversion of the incoming signal to an intermediate frequency (IF), requires a relatively large electronic circuit, including a local oscillator, mixer, various filters and a rectifier (i.e. diode). In order to radically simplify this process, we propose a spintronics-based nanometric single device whose resistance varies linearly as a function of the incoming frequency. This analogue frequency sensor allows a direct measurement of the incoming frequency and can work over many frequency bands and could replace a significantly larger conventional electronics circuit, thus reducing the footprint and potentially the energy cost. Three types of operation of the frequency sensor are demonstrated: a HF/VHF sensor based on a subharmonic modulation working across 9 different bands with 1-10 MHz bandwidth and 200 kHz resolution, a VHF/UHF sensor based on direct gyrotropic excitation operating between 150-225MHz with 250 kHz resolution and a UHF/SHF sensor based on indirect spin wave excitation operating between 2-4GHz with a 2.5 MHz resolution.


title: On Machine Learning and Structure for Mobile Robots

abstract: Due to recent advances - compute, data, models - the role of learning in autonomous systems has expanded significantly, rendering new applications possible for the first time. While some of the most significant benefits are obtained in the perception modules of the software stack, other aspects continue to rely on known manual procedures based on prior knowledge on geometry, dynamics, kinematics etc. Nonetheless, learning gains relevance in these modules when data collection and curation become easier than manual rule design. Building on this coarse and broad survey of current research, the final sections aim to provide insights into future potentials and challenges as well as the necessity of structure in current practical applications.


title: Partially-Supervised Image Captioning

abstract: Image captioning models are becoming increasingly successful at describing the content of images in restricted domains. However, if these models are to function in the wild - for example, as assistants for people with impaired vision - a much larger number and variety of visual concepts must be understood. To address this problem, we teach image captioning models new visual concepts from labeled images and object detection datasets. Since image labels and object classes can be interpreted as partial captions, we formulate this problem as learning from partially-specified sequence data. We then propose a novel algorithm for training sequence models, such as recurrent neural networks, on partially-specified sequences which we represent using finite state automata. In the context of image captioning, our method lifts the restriction that previously required image captioning models to be trained on paired image-sentence corpora only, or otherwise required specialized model architectures to take advantage of alternative data modalities. Applying our approach to an existing neural captioning model, we achieve state of the art results on the novel object captioning task using the COCO dataset. We further show that we can train a captioning model to describe new visual concepts from the Open Images dataset while maintaining competitive COCO evaluation scores.


title: Co-orbital resonance with a migrating proto-giant planet

abstract: In this work we pose the possibility that, at an early stage, the migration of a proto--giant planet caused by the presence of a gaseous circumstellar disk could explain the continuous feeding of small bodies into its orbit. Particularly, we study the probability of capture and permanence in co--orbital resonance of these small bodies, as planets of diverse masses migrate by interaction with the gaseous disk, and the drag induced by this disk dissipates energy from these small objects, making capture more likely. Also, we study the relevance of the circumplanetary disk, a structure formed closely around the planet where the gas density is enhanced, in the process of capture. It is of great interest for us to study the capture of small bodies in 1:1 resonance because it could account for the origin of the Trojan population, which has been proposed \citep{2011Icar..215..669K} as a mechanism of quasi-satellites and irregular satellites capture.


title: Predicting polymorphism in molecular crystals using orientational entropy

abstract: We introduce a computational method to discover polymorphs in molecular crystals at finite temperature. The method is based on reproducing the crystallization process starting from the liquid and letting the system discover the relevant polymorphs. This idea, however, conflicts with the fact that crystallization has a time scale much longer than that of molecular simulations. In order to bring the process within affordable simulation time, we enhance the fluctuations of a collective variable by constructing a bias potential with well tempered metadynamics. We use as collective variable an entropy surrogate based on an extended pair correlation function that includes the correlation between the orientation of pairs of molecules. We also propose a similarity metric between configurations based on the extended pair correlation function and a generalized Kullback-Leibler divergence. In this way, we automatically classify the configurations as belonging to a given polymorph using our metric and a hierarchical clustering algorithm. We find all relevant polymorphs for both substances and we predict new polymorphs. One of them is stabilized at finite temperature by entropic effects.


title: Golden ratios, Lucas Sequences and the Quadratic Family

abstract: It is conjectured that there is a converging sequence of some generalized Fibonacci ratios, given the difference between consecutive ratios, such as the Golden Ratio, $\varphi^1$, and the next golden ratio $\varphi^2$. Moreover, the graphic depiction of those ratios show some overlap with the quadratic family, and some numerical evidence suggest that everyone of those ratios in the finite set obtained, belong to at least one quadratic family, and finally a proof is presented that the converging sequence of some generalized Fibonacci ratios belong to at least one quadratic family.


title: An Optomechanical Quantum Cavendish Experiment

abstract: An open question in experimental physics is the characterization of gravitational effects in quantum regimes. We propose an experimental set-up that uses well-tested techniques in cavity optomechanics to observe the effects of the gravitational interaction between two micro-mechanical oscillators on the interference of the cavity photons through the shifts in the visibility of interfering photons. The gravitational coupling leads to a shift in the period and magnitude of the visibility whose observability is within reach of current technology. We discuss the feasibility of the set-up as well as the effects on entanglement due to gravitational interaction.


title: A posteriori error estimates for the Stokes problem with singular sources

abstract: We propose a posteriori error estimators for classical low-order inf-sup stable and stabilized finite element approximations of the Stokes problem with singular sources in two and three dimensional Lipschitz, but not necessarily convex, polytopal domains. The designed error estimators are proven to be reliable and locally efficient. On the basis of these estimators we design a simple adaptive strategy that yields optimal rates of convergence for the numerical examples that we perform.


title: A Self-Replication Basis for Designing Complex Agents

abstract: In this work, we describe a self-replication-based mechanism for designing agents of increasing complexity. We demonstrate the validity of this approach by solving simple, standard evolutionary computation problems in simulation. In the context of these simulation results, we describe the fundamental differences of this approach when compared to traditional approaches. Further, we highlight the possible advantages of applying this approach to the problem of designing complex artificial agents, along with the potential drawbacks and issues to be addressed in the future.


title: Bounds on the number of 2-level polytopes, cones and configurations

abstract: We prove an upper bound of the form $2^{O(d^2 \mathrm{polylog}\,d)}$ on the number of affine (resp. linear) equivalence classes of, by increasing order of generality, 2-level d-polytopes, d-cones and d-configurations. This in particular answers positively a conjecture of Bohn et al. on 2-level polytopes. We obtain our upper bound by relating affine (resp. linear) equivalence classes of 2-level d-polytopes, d-cones and d-configurations to faces of the correlation cone. We complement this with a $2^{\Omega(d^2)}$ lower bound, by estimating the number of nonequivalent stable set polytopes of bipartite graphs.


title: A New Decision Theoretic Sampling Plan for Exponential Distribution under Type-I Censoring

abstract: In this paper a new decision theoretic sampling plan (DSP) is proposed for Type-I censored exponential distribution. The proposed DSP is based on a new estimator of the expected lifetime of an exponential distribution which always exists, unlike the usual maximum likelihood estimator. The DSP is a modification of the Bayesian variable sampling plan of \cite{Lam:1994}. An optimum DSP is derived in the sense that it minimizes the Bayes risk. In terms of the Bayes risks, it performs better than Lam's sampling plan and its performance is as good as the Bayesian sampling plan of \cite{LLH:2002}, although implementation of the DSP is very simple. Analytically it is more tractable than the Bayesian sampling plan of \cite{LLH:2002}, and it can be easily generalized for any other loss functions also. A finite algorithm is provided to obtain the optimal plan and the corresponding minimum Bayes risk is calculated. Extensive numerical comparisons with the optimal Bayesian sampling plan proposed by \cite{LLH:2002} are made. The results have been extended for three degree polynomial loss function and for Type-I hybrid censoring scheme.


title: Markov Process Simulation with Quantum Computer

abstract: This paper shows a novel way of simulating a Markov process by a quantum computer. The main purpose of the paper is to show a particular application of quantum computing in the field of stochastic processes analysis. Using a Quantum computer, the process could be superposed, where the random variables of the Markov chain are represented by entangled qubit states, which gives the great opportunity of having all the possible scenarios simultaneously.


title: Basics of Laser-Plasma Interaction: a Selection of Topics

abstract: A short, tutorial introduction to some basic concepts of laser-plasma interactions at ultra-high intensities is given. The selected topics include a) elements of the relativistic dynamics of an electron in electromagnetic fields, including the ponderomotive force and classical radiation friction; b) the "relativistic" nonlinear optical transparency and self-focusing; c) the moving mirror concept and its application to light sail acceleration and high harmonic generation, with a note on related instabilities; d) some specific phenomena related to the absorption of energy, kinetic momentum and angular momentum from the laser light.


title: Probing the scale of non-commutativity of space

abstract: Examining quantum electrodynamics in non-commutative (NC) spaces along with composite operators in these spaces, we show that i) any charge g for a fermion matter field is allowed provided the basic NC photon-photon coupling is g, however no other multiples of g are permitted and ii) composite operators do not have a simple transformation which can be attributed to the effective total charge of the composite particle. Taken together these results place a limit on the scale of non-commutativity to be at most smaller that current LHC limits for compositeness. Furthermore, they also suggest that a substructure at still smaller scales is needed if such spaces are to be a physical reality.


title: Accelerating molecular dynamics simulations with population annealing

abstract: Population annealing is a powerful tool for large-scale Monte Carlo simulations. We adapt this method to molecular dynamics simulations and demonstrate its excellent accelerating effect by simulating the folding of a short peptide commonly used to gauge the performance of algorithms. The method is compared to the well established parallel tempering approach and is found to yield similar performance for the same computational resources. In contrast to other methods, however, population annealing scales to a nearly arbitrary number of parallel processors and it is thus a unique tool that enables molecular dynamics to tap into the massively parallel computing power available in supercomputers that is so much needed for a range of difficult computational problems.


title: Homonym Detection in Curated Bibliographies: Learning from dblp's Experience (full version)

abstract: Identifying (and fixing) homonymous and synonymous author profiles is one of the major tasks of curating personalized bibliographic metadata repositories like the dblp computer science bibliography. In this paper, we present and evaluate a machine learning approach to identify homonymous author bibliographies using a simple multilayer perceptron setup. We train our model on a novel gold-standard data set derived from the past years of active, manual curation at the dblp computer science bibliography.


title: Probing for Extra Top Yukawa Couplings in Light of $t\bar th(125)$ Observation

abstract: The observation of $t\bar th(125)$ production at the Large Hadron Collider (LHC) is the first direct measurement of the top Yukawa coupling. It opens the window on an extra top Yukawa coupling, $\rho_{tt}$, from a second Higgs doublet, without a $Z_2$ symmetry to forbid flavor changing neutral Higgs couplings. We show that $t\bar th$ and Higgs property measurements at the High Luminosity LHC can constrain the ${\rm Re}\,\rho_{tt}$--${\rm Im}\,\rho_{tt}$ parameter space that could drive electroweak baryogenesis, but the $\Gamma_h$ width measurement must be considerably improved beyond current projections.


title: Shifted-antimagic Labelings for Graphs

abstract: The concept of antimagic labelings of a graph is to produce distinct vertex sums by labeling edges through consecutive numbers starting from one. A long-standing conjecture is that every connected graph, except a single edge, is antimagic. Some graphs are known to be antimagic, but little has been known about sparse graphs, not even trees. This paper studies a weak version called $k$-shifted-antimagic labelings which allow the consecutive numbers starting from $k+1$, instead of starting from 1, where $k$ can be any integer. This paper establishes connections among various concepts proposed in the literature of antimagic labelings and extends previous results in three aspects: $\bullet$ Some classes of graphs, including trees and graphs whose vertices are of odd degrees, which have not been verified to be antimagic are shown to be $k$-shifted-antimagic for sufficiently large $k$. $\bullet$ Some graphs are proved $k$-shifted-antimagic for all $k$, while some are proved not for some particular $k$. $\bullet$ Disconnected graphs are also considered.


title: Imbalanced Randomization in Clinical Trials

abstract: Randomization is a common technique used in clinical trials to eliminate potential bias and confounders in a patient population. Equal allocation to treatment groups is the standard due to its optimal efficiency in many cases. However, in certain scenarios, unequal allocation can improve efficiency. In superiority trials with more than two groups, the optimal randomization is not always a balanced randomization. In non-inferiority trials, additive margin with equal variance is the only instance with balanced randomization. Optimal randomization for non-inferiority trials can be far from 1:1 and can greatly improve efficiency, a fact which is commonly overlooked. A tool for sample size calculation for non-inferiority trials with additive or multiplicative margin with normal, binomial or Poisson distribution is available at http://www.statlab.wisc.edu/shiny/SSNI/.


title: Yamabe flow on non-compact manifolds with unbounded initial curvature

abstract: We prove global existence of Yamabe flows on non-compact manifolds $M$ of dimension $m\geq3$ under the assumption that the initial metric $g_0=u_0g_M$ is conformally equivalent to a complete background metric $g_M$ of bounded, non-positive scalar curvature and positive Yamabe invariant with conformal factor $u_0$ bounded from above and below. We do not require initial curvature bounds. In particular, the scalar curvature of $(M,g_0)$ can be unbounded from above and below without growth condition.


title: On finite sets of small tripling or small alternation in arbitrary groups

abstract: We prove Bogolyubov-Ruzsa-type results for finite subsets of groups with small tripling, $|A^3|\leq O(|A|)$, or small alternation, $|AA^{\text{-}1} A|\leq O(|A|)$. As applications, we obtain a qualitative analog of Bogolyubov's Lemma for dense sets in arbitrary finite groups, as well as a quantitative arithmetic regularity lemma for sets of bounded VC-dimension in finite groups of bounded exponent. The latter result generalizes the abelian case, due to Alon, Fox, and Zhao, and gives a quantitative version of previous work of the author, Pillay, and Terry.


title: Explicit expressions for the related numbers of higher order Appell polynomials

abstract: In this note, by using the Hasse-Teichm\"uller derivatives, we obtain two explicit expressions for the related numbers of higher order Appell polynomials. One of them presents a determinant expression for the related numbers of higher order Appell polynomials, which involves several determinant expressions of special numbers, such as the higher order generalized hypergeometric Bernoulli and Cauchy numbers, thus recovers the classical determinant expressions of Bernoulli and Cauchy numbers stated in an article by Glaisher in 1875.


title: Linear Single- and Three-Phase Voltage Forecasting and Bayesian State Estimation with Limited Sensing

abstract: Implementing state estimation in low and medium voltage power distribution is still challenging given the scale of many networks and the reliance of traditional methods on a large number of measurements. This paper proposes a method to improve voltage predictions in real-time by leveraging a limited set of real-time measurements. The method relies on Bayesian estimation formulated as a linear least squares estimation problem, which resembles the classical weighted least-squares (WLS) approach for scenarios where full network observability is not available. We build on recently developed linear approximations for unbalanced three-phase power flow to construct voltage predictions as a linear mapping of load predictions constructed with Gaussian processes. The estimation step to update the voltage forecasts in real-time is a linear computation allowing fast high-resolution state estimate updates. The uncertainty in forecasts can be determined a priori and smoothed a posteriori, making the method useful for both planning, operation and post-hoc analysis. The method outperforms conventional WLS and is applied to different test feeders and validated on a real test feeder with the utility Alliander in The Netherlands.


title: Low-Frequency Quasi-Periodic Oscillations in the X-ray Nova MAXI J1535-571 at the Initial Stage of Its 2017 Outburst

abstract: We report the discovery of low-frequency quasi-periodic oscillations (QPOs) in the power spectrum of the X-ray nova MAXI J1535-571 at the initial stage of its outburst in September 2017. Based on data from the SWIFT and INTEGRAL instruments, we have traced the evolution of the QPO parameters (primarily their frequency) with time and their correlation with changes in the X-ray spectrum of the source (changes in the emission flux and hardness). We place constraints on the theoretical QPO generation models.


title: Fingerprints of stacking order in graphene layers from ab initio second-order Raman spectra

abstract: We present an \textit{ab initio} study based on density-functional theory of first- and second-order Raman spectra of graphene-based materials with different stacking arrangements and numbers of layers. Going from monolayer and bilayer graphene to periodic graphitic structures, we investigate the behavior of the first-order G-band and of the second-order 2D-band excited by the same set of photon energies. The former turns out to be very similar in all considered graphene-based materials, while in the latter we find the signatures of individual structures. With a systematic analysis of the second-order Raman spectra at varying frenquency of the incident radiation, we monitor the Raman signal and identify the contributions from different phonon modes that are characteristic of each specific arrangement. Supported by good agreement with experimental findings and with previous theoretical studies based on alternative approaches, our results propose an effective tool to probe and analyze the fingerprints of graphene-based and other low-dimensional materials.


title: A Leslie-Gower predator-prey model with a free boundary

abstract: In this paper, we consider a Leslie-Gower predator-prey model in one-dimensional environment. We study the asymptotic behavior of two species evolving in a domain with a free boundary. Sufficient conditions for spreading success and spreading failure are obtained. We also derive sharp criteria for spreading and vanishing of the two species. Finally, when spreading is successful, we show that the spreading speed is between the minimal speed of traveling wavefront solutions for the predator-prey model on the whole real line (without a free boundary) and an elliptic problem that follows from the original model.


title: A new characterization of the Gamma distribution and associated goodness of fit tests

abstract: We propose a class of weighted $L_2$-type tests of fit to the Gamma distribution. Our novel procedure is based on a fixed point property of a new transformation connected to a Steinian characterization of the family of Gamma distributions. We derive the weak limits of the statistic under the null hypothesis and under contiguous alternatives. Further, we establish the global consistency of the tests and apply a parametric bootstrap technique in a Monte Carlo simulation study to show the competitiveness to existing procedures.


title: One-Shot Unsupervised Cross Domain Translation

abstract: Given a single image x from domain A and a set of images from domain B, our task is to generate the analogous of x in B. We argue that this task could be a key AI capability that underlines the ability of cognitive agents to act in the world and present empirical evidence that the existing unsupervised domain translation methods fail on this task. Our method follows a two step process. First, a variational autoencoder for domain B is trained. Then, given the new sample x, we create a variational autoencoder for domain A by adapting the layers that are close to the image in order to directly fit x, and only indirectly adapt the other layers. Our experiments indicate that the new method does as well, when trained on one sample x, as the existing domain transfer methods, when these enjoy a multitude of training samples from domain A. Our code is made publicly available at https://github.com/sagiebenaim/OneShotTranslation


title: A PDE approach to fractional diffusion: a space-fractional wave equation

abstract: We study solution techniques for an evolution equation involving second order derivative in time and the spectral fractional powers, of order $s \in (0,1)$, of symmetric, coercive, linear, elliptic, second-order operators in bounded domains $\Omega$. We realize fractional diffusion as the Dirichlet-to-Neumann map for a nonuniformly elliptic problem posed on the semi-infinite cylinder $\mathcal{C} = \Omega \times (0,\infty)$. We thus rewrite our evolution problem as a quasi-stationary elliptic problem with a dynamic boundary condition and derive space, time, and space-time regularity estimates for its solution. The latter problem exhibits an exponential decay in the extended dimension and thus suggests a truncation that is suitable for numerical approximation. We propose and analyze two fully discrete schemes. The discretization in time is based on finite difference discretization techniques: trapezoidal and leapfrog schemes. The discretization in space relies on the tensorization of a first-degree FEM in $\Omega$ with a suitable $hp$-FEM in the extended variable. For both schemes we derive stability and error estimates.


title: Two Algorithms in Group Theory

abstract: We present a new algorithm deciding if the intersection of a quasiconvex subgroup of a negatively curved group with a conjugate is finite. We also give a short proof of decidability of the membership problem for quasiconvex subgroups of finitely generated groups with decidable word problem.


title: VIoLET: A Large-scale Virtual Environment for Internet of Things

abstract: IoT deployments have been growing manifold, encompassing sensors, networks, edge, fog and cloud resources. Despite the intense interest from researchers and practitioners, most do not have access to large-scale IoT testbeds for validation. Simulation environments that allow analytical modeling are a poor substitute for evaluating software platforms or application workloads in realistic computing environments. Here, we propose VIoLET, a virtual environment for defining and launching large-scale IoT deployments within cloud VMs. It offers a declarative model to specify container-based compute resources that match the performance of the native edge, fog and cloud devices using Docker. These can be inter-connected by complex topologies on which private/public networks, and bandwidth and latency rules are enforced. Users can configure synthetic sensors for data generation on these devices as well. We validate VIoLET for deployments with > 400 devices and > 1500 device-cores, and show that the virtual IoT environment closely matches the expected compute and network performance at modest costs. This fills an important gap between IoT simulators and real deployments.


title: The Minimum Principle for Convex Subequations

abstract: A subequation on an open subset $X\subset \mathbb R^n$ is a subset $F$ of the space of $2$-jets on $X$ with certain properties. A smooth function is said to be $F$-subharmonic if all of its $2$-jets lie in $F$, and using the viscosity technique one can extend the notion of $F$-subharmonicity to any upper-semicontinuous function. Let $\mathcal P$ denote the subequation consisting of those $2$-jets whose Hessian part is semipositive. We introduce a notion of product subequation $F\#\mathcal P$ on $X\times \mathbb R^{m}$ and prove, under suitable hypotheses, that if $F$ is convex and $f(x,y)$ is $F\#\mathcal P$-subharmonic then the marginal function $$ g(x):= \inf_y f(x,y)$$ is $F$-subharmonic. This generalises the classical statement that the marginal function of a convex function is again convex. We also prove a complex version of this result that generalises the Kiselman minimum principle for the marginal function of a plurisubharmonic function.


title: The Toybox Dataset of Egocentric Visual Object Transformations

abstract: In object recognition research, many commonly used datasets (e.g., ImageNet and similar) contain relatively sparse distributions of object instances and views, e.g., one might see a thousand different pictures of a thousand different giraffes, mostly taken from a few conventionally photographed angles. These distributional properties constrain the types of computational experiments that are able to be conducted with such datasets, and also do not reflect naturalistic patterns of embodied visual experience. As a contribution to the small (but growing) number of multi-view object datasets that have been created to bridge this gap, we introduce a new video dataset called Toybox that contains egocentric (i.e., first-person perspective) videos of common household objects and toys being manually manipulated to undergo structured transformations, such as rotation, translation, and zooming. To illustrate potential uses of Toybox, we also present initial neural network experiments that examine 1) how training on different distributions of object instances and views affects recognition performance, and 2) how viewpoint-dependent object concepts are represented within the hidden layers of a trained network.


title: Customized Local Differential Privacy for Multi-Agent Distributed Optimization

abstract: Real-time data-driven optimization and control problems over networks may require sensitive information of participating users to calculate solutions and decision variables, such as in traffic or energy systems. Adversaries with access to coordination signals may potentially decode information on individual users and put user privacy at risk. We develop local differential privacy, which is a strong notion that guarantees user privacy regardless of any auxiliary information an adversary may have, for a larger family of convex distributed optimization problems. The mechanism allows agent to customize their own privacy level based on local needs and parameter sensitivities. We propose a general sampling based approach for determining sensitivity and derive analytical bounds for specific quadratic problems. We analyze inherent trade-offs between privacy and suboptimality and propose allocation schemes to divide the maximum allowable noise, a privacy budget, among all participating agents. Our algorithm is implemented to enable privacy in distributed optimal power flow for electric grids.


title: Flavor anomalies and dark matter in SUSY with an extra U(1)

abstract: Motivated by the recent anomalies in $b\rightarrow s$ transitions that emerged at LHCb, we consider a model with an $L_\mu - L_\tau$ gauge symmetry and additional vector-like fermions. We find that by introducing supersymmetry the model can be made consistent with the long-standing deviation in the measured value of the anomalous magnetic moment of the muon, $g-2$, and neutralino dark matter of broad mass ranges and properties. In particular, dark matter candidates include the well-known 1 TeV higgsino, which in the MSSM is typically not compatible with solutions to the $g-2$ puzzle. Moreover, its spin-independent cross section could be at the origin of the recent small excess in XENON-1T data. We apply to the model constraints arising from flavor precision measurements and direct searches at the Large Hadron Collider and show that they do not currently exclude the relevant parameter space regions.


title: Differential Evolution Algorithm Aided Turbo Channel Estimation and Multi-User Detection for G.Fast Systems in the Presence of FEXT

abstract: The ever-increasing demand for broadband Internet access has motivated the further development of the digital subscriber line to the G.fast standard in order to expand its operational band from 106 MHz to 212 MHz. Conventional far-end crosstalk (FEXT) based cancellers falter in the upstream transmission of this emerging G.fast system. In this paper, we propose a novel differential evolution algorithm (DEA) aided turbo channel estimation (CE) and multi-user detection (MUD) scheme for the G.fast upstream including the frequency band up to 212 MHz, which is capable of approaching the optimal Cramer-Rao lower bound of the channel estimate, whilst approaching the optimal maximum likelihood (ML) MUD's performance associated with perfect channel state information, and yet only imposing about 5% of its computational complexity. Explicitly, the turbo concept is exploited by iteratively exchanging information between the continuous value-based DEA assisted channel estimator and the discrete value-based DEA MUD. Our extensive simulations show that 18 dB normalized mean square error gain is attained by the channel estimator and 10 dB signal-to-noise ratio gain can be achieved by the MUD upon exploiting this iteration gain. We also quantify the influence of the CE error, of the copper length and of the impulse noise. Our study demonstrates that the proposed DEA aided turbo CE and MUD scheme is capable of offering near-capacity performance at an affordable complexity for the emerging G.fast systems.


title: The merger that led to the formation of the Milky Way's inner stellar halo and thick disk

abstract: The assembly process of our Galaxy can be retrieved using the motions and chemistry of individual stars. Chemo-dynamical studies of the nearby halo have long hinted at the presence of multiple components such as streams, clumps, duality and correlations between the stars' chemical abundances and orbital parameters. More recently, the analysis of two large stellar surveys have revealed the presence of a well-populated chemical elemental abundance sequence, of two distinct sequences in the colour-magnitude diagram, and of a prominent slightly retrograde kinematic structure all in the nearby halo, which may trace an important accretion event experienced by the Galaxy. Here report an analysis of the kinematics, chemistry, age and spatial distribution of stars in a relatively large volume around the Sun that are mainly linked to two major Galactic components, the thick disk and the stellar halo. We demonstrate that the inner halo is dominated by debris from an object which at infall was slightly more massive than the Small Magellanic Cloud, and which we refer to as Gaia-Enceladus. The stars originating in Gaia-Enceladus cover nearly the full sky, their motions reveal the presence of streams and slightly retrograde and elongated trajectories. Hundreds of RR Lyrae stars and thirteen globular clusters following a consistent age-metallicity relation can be associated to Gaia-Enceladus on the basis of their orbits. With an estimated 4:1 mass-ratio, the merger with Gaia-Enceladus must have led to the dynamical heating of the precursor of the Galactic thick disk and therefore contributed to the formation of this component approximately 10 Gyr ago. These findings are in line with simulations of galaxy formation, which predict that the inner stellar halo should be dominated by debris from just a few massive progenitors.


title: (K,L)-eigenvectors in max-min algebra

abstract: Using the concept of (K,L)-eigenvector, we investigate the structure of the max-min eigenspace associated with a given eigenvalue of a matrix in the max-min algebra (also known as fuzzy algebra). In our approach, the max-min eigenspace is split into several regions according to the order relations between the eigenvalue and the components of x. The resulting theory of (K,L)-eigenvectors, being based on the fundamental results of Gondran and Minoux, allows to describe the whole max-min eigenspace explicitly and in more detail.


title: Materials Informatics for Dark Matter Detection

abstract: Dark Matter particles are commonly assumed to be weakly interacting massive particles (WIMPs) with a mass in the GeV to TeV range. However, recent interest has shifted towards lighter WIMPs, which are more difficult to probe experimentally. A detection of sub-GeV WIMPs would require the use of small gap materials in sensors. Using recent estimates of the WIMP mass, we identify the relevant target space towards small gap materials (100-10 meV). Dirac Materials, a class of small- or zero-gap materials, emerge as natural candidates for sensors for Dark Matter detection. We propose the use of informatics tools to rapidly assay materials band structures to search for small gap semiconductors and semimetals, rather than focusing on a few preselected compounds. As a specific example of the proposed strategy, we use the organic materials database (omdb.diracmaterials.org) to identify organic candidates for sensors: the narrow band gap semiconductors BNQ-TTF and DEBTTT with gaps of 40 and 38 meV, and the Dirac-line semimetal (BEDT-TTF)$\cdot$Br which exhibits a tiny gap of $\approx$ 50 meV when spin-orbit coupling is included. We outline a novel and powerful approach to search for dark matter detection sensor materials by means of a rapid assay of materials using informatics tools.

categories: cond-mat.mtrl-sci astro-ph.HE astro-ph.IM hep-ex hep-ph


title: Self-contained two-layer shallow water theory of strong internal bores

abstract: We show that interfacial gravity waves comprising strong hydraulic jumps (bores) can be described by a two-layer hydrostatic shallow-water (SW) approximation without invoking additional front conditions. The theory is based on a new SW momentum equation which is derived in locally conservative form containing a free parameter $\alpha.$ This parameter, which defines the relative contribution of each layer to the pressure at the interface, affects only hydraulic jumps but not continuous waves. The Rankine-Hugoniot jump conditions for the momentum and mass conservation equations are found to be mathematically equivalent to the classical front conditions, which were previously thought to be outside the scope of SW approximation. Dimensional arguments suggest that $\alpha$ depends on the density ratio. For nearly equal densities, both layers are expected to affect interfacial pressure with approximately equal weight coefficients, which corresponds to $\alpha\approx0.$ The front propagation velocity for $\alpha=0$ agrees well with experimental and numerical results in a wide range of bore strengths. A remarkably better agreement with high-accuracy numerical results is achieved by $\alpha=\sqrt{5}-2,$ which yields the largest height that a stable gravity current can have.


title: Pathways to proposing causes for unexpected experimental results

abstract: Models of physical systems are used to explain and predict experimental results and observations. When students encounter discrepancies between the actual and expected behavior of a system, they revise their models to include the newly acquired observations, or change their apparatus to better represent their models. The Modeling Framework for Experimental Physics describes the process of matching measurements and observations to models by making revisions to resolve discrepancies. As part of a larger effort to create assessments of students' modeling abilities in the context of upper-division electronics courses, we used the Modeling Framework to develop and code think-aloud problem-solving activities centered on troubleshooting an inverting amplifier circuit. We observed that some participants iteratively and continuously made measurements and revisions if they could not immediately propose a cause for an observed discrepancy. This pathway has not been previously discussed in the Modeling Framework. In this paper, we discuss two episodes where students undergo this process to converge on a proposed cause post hoc. We conclude by discussing implications for a modeling assessment based on the observed modeling behavior.


title: AVX-512 extension to OpenQCD 1.6

abstract: We publish an extension of openQCD-1.6 with AVX-512 vector instructions using Intel intrinsics. Recent Intel processors support extended instruction sets with operations on 512-bit wide vectors, increasing both the capacity for floating point operations and register memory. Optimal use of the new capabilities requires reorganising data and floating point operations into these wider vector units. We report on the implementation and performance of the AVX-512 OpenQCD extension on clusters using Intel Knights Landing and Xeon Scalable (Skylake) CPUs. In complete HMC trajectories with physically relevant parameters we observe a performance increase of 5% to 10%.