21,039 research outputs found
Simulating Brownian suspensions with fluctuating hydrodynamics
Fluctuating hydrodynamics has been successfully combined with several
computational methods to rapidly compute the correlated random velocities of
Brownian particles. In the overdamped limit where both particle and fluid
inertia are ignored, one must also account for a Brownian drift term in order
to successfully update the particle positions. In this paper, we present an
efficient computational method for the dynamic simulation of Brownian
suspensions with fluctuating hydrodynamics that handles both computations and
provides a similar approximation as Stokesian Dynamics for dilute and
semidilute suspensions. This advancement relies on combining the fluctuating
force-coupling method (FCM) with a new midpoint time-integration scheme we
refer to as the drifter-corrector (DC). The DC resolves the drift term for
fluctuating hydrodynamics-based methods at a minimal computational cost when
constraints are imposed on the fluid flow to obtain the stresslet corrections
to the particle hydrodynamic interactions. With the DC, this constraint need
only be imposed once per time step, reducing the simulation cost to nearly that
of a completely deterministic simulation. By performing a series of
simulations, we show that the DC with fluctuating FCM is an effective and
versatile approach as it reproduces both the equilibrium distribution and the
evolution of particulate suspensions in periodic as well as bounded domains. In
addition, we demonstrate that fluctuating FCM coupled with the DC provides an
efficient and accurate method for large-scale dynamic simulation of colloidal
dispersions and the study of processes such as colloidal gelation
One-dimensional Chern-Simons theory
We study a one-dimensional toy version of the Chern-Simons theory. We
construct its simplicial version which comprises features of a low-energy
effective gauge theory and of a topological quantum field theory in the sense
of Atiyah.Comment: 37 page
Rank Centrality: Ranking from Pair-wise Comparisons
The question of aggregating pair-wise comparisons to obtain a global ranking
over a collection of objects has been of interest for a very long time: be it
ranking of online gamers (e.g. MSR's TrueSkill system) and chess players,
aggregating social opinions, or deciding which product to sell based on
transactions. In most settings, in addition to obtaining a ranking, finding
`scores' for each object (e.g. player's rating) is of interest for
understanding the intensity of the preferences.
In this paper, we propose Rank Centrality, an iterative rank aggregation
algorithm for discovering scores for objects (or items) from pair-wise
comparisons. The algorithm has a natural random walk interpretation over the
graph of objects with an edge present between a pair of objects if they are
compared; the score, which we call Rank Centrality, of an object turns out to
be its stationary probability under this random walk. To study the efficacy of
the algorithm, we consider the popular Bradley-Terry-Luce (BTL) model
(equivalent to the Multinomial Logit (MNL) for pair-wise comparisons) in which
each object has an associated score which determines the probabilistic outcomes
of pair-wise comparisons between objects. In terms of the pair-wise marginal
probabilities, which is the main subject of this paper, the MNL model and the
BTL model are identical. We bound the finite sample error rates between the
scores assumed by the BTL model and those estimated by our algorithm. In
particular, the number of samples required to learn the score well with high
probability depends on the structure of the comparison graph. When the
Laplacian of the comparison graph has a strictly positive spectral gap, e.g.
each item is compared to a subset of randomly chosen items, this leads to
dependence on the number of samples that is nearly order-optimal.Comment: 45 pages, 3 figure
Least Squares Ranking on Graphs
Given a set of alternatives to be ranked, and some pairwise comparison data,
ranking is a least squares computation on a graph. The vertices are the
alternatives, and the edge values comprise the comparison data. The basic idea
is very simple and old: come up with values on vertices such that their
differences match the given edge data. Since an exact match will usually be
impossible, one settles for matching in a least squares sense. This formulation
was first described by Leake in 1976 for rankingfootball teams and appears as
an example in Professor Gilbert Strang's classic linear algebra textbook. If
one is willing to look into the residual a little further, then the problem
really comes alive, as shown effectively by the remarkable recent paper of
Jiang et al. With or without this twist, the humble least squares problem on
graphs has far-reaching connections with many current areas ofresearch. These
connections are to theoretical computer science (spectral graph theory, and
multilevel methods for graph Laplacian systems); numerical analysis (algebraic
multigrid, and finite element exterior calculus); other mathematics (Hodge
decomposition, and random clique complexes); and applications (arbitrage, and
ranking of sports teams). Not all of these connections are explored in this
paper, but many are. The underlying ideas are easy to explain, requiring only
the four fundamental subspaces from elementary linear algebra. One of our aims
is to explain these basic ideas and connections, to get researchers in many
fields interested in this topic. Another aim is to use our numerical
experiments for guidance on selecting methods and exposing the need for further
development.Comment: Added missing references, comparison of linear solvers overhauled,
conclusion section added, some new figures adde
Ground-state Stabilization of Open Quantum Systems by Dissipation
Control by dissipation, or environment engineering, constitutes an important
methodology within quantum coherent control which was proposed to improve the
robustness and scalability of quantum control systems. The system-environment
coupling, often considered to be detrimental to quantum coherence, also
provides the means to steer the system to desired states. This paper aims to
develop the theory for engineering of the dissipation, based on a ground-state
Lyapunov stability analysis of open quantum systems via a Heisenberg-picture
approach. Algebraic conditions concerning the ground-state stability and
scalability of quantum systems are obtained. In particular, Lyapunov stability
conditions expressed as operator inequalities allow a purely algebraic
treatment of the environment engineering problem, which facilitates the
integration of quantum components into a large-scale quantum system and draws
an explicit connection to the classical theory of vector Lyapunov functions and
decomposition-aggregation methods for control of complex systems. The
implications of the results in relation to dissipative quantum computing and
state engineering are also discussed in this paper.Comment: 18 pages, to appear in Automatic
- …