Related papers: Minimum relaxation n-body simulations using orthog…
Common time-explicit numerical methods for kinetic simulations of plasmas in the low-collisions limit fall into two classes of algorithms: momentum conserving and energy conserving. Each has certain drawbacks. The PIC algorithm does not…
This paper presents computationally feasible rank-one relaxation algorithms for the efficient simulation of a time-incremental damage model with nonconvex incremental stress potentials in multiple spatial dimensions. While the standard…
We discuss a basis set developed to calculate perturbation coefficients in an expansion of the general N-body problem. This basis has two advantages. First, the basis is complete order-by-order for the perturbation series. Second, the…
We study dissipation in a small quantum system coupled to an environment held in thermodynamic equilibrium. The relaxation dynamics of a system subject to an abrupt quench in the parameters of the underlying Hamiltonian is investigated…
Sorting input objects is an important step in many machine learning pipelines. However, the sorting operator is non-differentiable with respect to its inputs, which prohibits end-to-end gradient-based optimization. In this work, we propose…
The relaxation function is the cornerstone to perform calculations in weakly driven processes. Properties that such a function should obey are already established, but the difficulty in its calculation is still an issue to be overcome. In…
We propose a dual decomposition and linear program relaxation of the NP -hard minimum cost multicut problem. Unlike other polyhedral relaxations of the multicut polytope, it is amenable to efficient optimization by message passing. Like…
The work explores a specific scenario for structural computational optimization based on the following elements: (a) a relaxed optimization setting considering the ersatz (bi-material) approximation, (b) a treatment based on a nonsmoothed…
Building mathematical models of brains is difficult because of the sheer complexity of the problem. One potential starting point is through basal cognition, which gives abstract representation of a range of organisms without central nervous…
We carry out overdamped simulations in a simple model of jamming - a collection of bi-disperse soft core frictionless disks in two dimensions - with the aim to explore the finite size dependence of different quantities, both the relaxation…
We present two modifications of the standard cell list algorithm for nonequilibrium molecular dynamics simulations of homogeneous, linear flows. When such a flow is modeled with periodic boundary conditions, the simulation box deforms with…
One of the fundamental aspects of statistical behaviour in many-body systems is exponential divergence of neighbouring orbits, which is often discussed in terms of Liapounov exponents. Here we study this topic for the classical…
In this paper, a new mixed integer nonlinear programming formulation is proposed for optimally placing and operating pressure reducing valves and chlorine booster stations in water distribution networks. The objective is the minimization of…
Generalized Fourier series with orthogonal polynomial bases have useful applications in several fields, including differential equations, pattern recognition, and image and signal processing. However, computing the generalized Fourier…
Symmetries (and their spontaneous rupturing) can be used to protect and engender novel quantum phases and lead to interesting collective phenomena. In Ref. 1, the authors described a general dynamical decoupling (polyfractal) protocol that…
It is common practice in molecular dynamics and Monte Carlo computer simulations to run multiple, separately-initialized simulations in order to improve the sampling of independent microstates. Here we examine the utility of an extreme case…
We consider a class of relaxation problems mixing slow and fast variations which can describe population dynamics models or hyperbolic systems, with varying stiffness (from non-stiff to strongly dissipative), and develop a multi-scale…
This paper is about partitioning in parallel and distributed simulation. That means decomposing the simulation model into a numberof components and to properly allocate them on the execution units. An adaptive solution based on…
We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the…
Cosmological N-Body simulations are used for a variety of applications. Indeed progress in the study of large scale structures and galaxy formation would have been very limited without this tool. For nearly twenty years the limitations…