Source author record

S. J. Thwaite

S. J. Thwaite appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math-ph math.MP cond-mat.other cond-mat.stat-mech Discrete Mathematics math.CA math.QA math.RA math.RT math.ST Statistics Theory

Catalog footprint

What is connected

5works

11topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2015arXiv

Walk-Sums, Continued Fractions and Unique Factorisation on Digraphs

We show that the series of all walks between any two vertices of any (possibly weighted) directed graph $\mathcal{G}$ is given by a universal continued fraction of finite depth and breadth involving the simple paths and simple cycles of $\mathcal{G}$. A simple path is a walk forbidden to visit any vertex more than once. We obtain an explicit formula giving this continued fraction. Our results are based on an equivalent to the fundamental theorem of arithmetic: we demonstrate that arbitrary walks on $\mathcal{G}$ factorize uniquely into nesting products of simple paths and simple cycles, where nesting is a product operation between walks that we define. We show that the simple paths and simple cycles are the prime elements of the set of all walks on $\mathcal{G}$ equipped with the nesting product. We give an algorithm producing the prime factorization of individual walks, and obtain a recursive formula producing the prime factorization of sets of walks. Our results have already found applications in machine learning, matrix computations and quantum mechanics.

preprint2014arXiv

An Exact Formulation of the Time-Ordered Exponential using Path-Sums

We present the path-sum formulation for $\mathsf{OE}[\mathsf{H}](t',t)=\mathcal{T}\,\text{exp}\big(\int_{t}^{t'}\!\mathsf{H}(τ)\,dτ\big)$, the time-ordered exponential of a time-dependent matrix $\mathsf{H}(t)$. The path-sum formulation gives $\mathsf{OE}[\mathsf{H}]$ as a branched continued fraction of finite depth and breadth. The terms of the path-sum have an elementary interpretation as self-avoiding walks and self-avoiding polygons on a graph. Our result is based on a representation of the time-ordered exponential as the inverse of an operator, the mapping of this inverse to sums of walks on graphs and the algebraic structure of sets of walks. We give examples demonstrating our approach. We establish a super-exponential decay bound for the magnitude of the entries of the time-ordered exponential of sparse matrices. We give explicit results for matrices with commonly encountered sparse structures.

preprint2014arXiv

Exact Inference on Gaussian Graphical Models of Arbitrary Topology using Path-Sums

We present the path-sum formulation for exact statistical inference of marginals on Gaussian graphical models of arbitrary topology. The path-sum formulation gives the covariance between each pair of variables as a branched continued fraction of finite depth and breadth. Our method originates from the closed-form resummation of infinite families of terms of the walk-sum representation of the covariance matrix. We prove that the path-sum formulation always exists for models whose covariance matrix is positive definite: i.e.~it is valid for both walk-summable and non-walk-summable graphical models of arbitrary topology. We show that for graphical models on trees the path-sum formulation is equivalent to Gaussian belief propagation. We also recover, as a corollary, an existing result that uses determinants to calculate the covariance matrix. We show that the path-sum formulation formulation is valid for arbitrary partitions of the inverse covariance matrix. We give detailed examples demonstrating our results.

preprint2012arXiv

Evaluating Matrix Functions by Resummations on Graphs: the Method of Path-Sums

We introduce the method of path-sums which is a tool for exactly evaluating a function of a discrete matrix with possibly non-commuting entries, based on the closed-form resummation of infinite families of terms in the corresponding Taylor series. If the matrix is finite, our approach yields the exact result in a finite number of steps. We achieve this by combining a mapping between matrix powers and walks on a weighted directed graph with a universal graph-theoretic result on the structure of such walks. We present path-sum expressions for a matrix raised to a complex power, the matrix exponential, matrix inverse, and matrix logarithm. We show that the quasideterminants of a matrix can be naturally formulated in terms of a path-sum, and present examples of the application of the path-sum method. We show that obtaining the inversion height of a matrix inverse and of quasideterminants is an NP-complete problem.

preprint2007arXiv

Pairing mean-field theory for the dynamics of dissociation of molecular Bose-Einstein condensates

We develop a pairing mean-field theory to describe the quantum dynamics of the dissociation of molecular Bose-Einstein condensates into their constituent bosonic or fermionic atoms. We apply the theory to one, two, and three-dimensional geometries and analyze the role of dimensionality on the atom production rate as a function of the dissociation energy. As well as determining the populations and coherences of the atoms, we calculate the correlations that exist between atoms of opposite momenta, including the column density correlations in 3D systems. We compare the results with those of the undepleted molecular field approximation and argue that the latter is most reliable in fermionic systems and in lower dimensions. In the bosonic case we compare the pairing mean-field results with exact calculations using the positive-$P$ stochastic method and estimate the range of validity of the pairing mean-field theory. Comparisons with similar first-principle simulations in the fermionic case are currently not available, however, we argue that the range of validity of the present approach should be broader for fermions than for bosons in the regime where Pauli blocking prevents complete depletion of the molecular condensate.