BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tom Kennedy

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

ResearcherUnclaimed source record
math-phmath.MPmath.PRcond-mat.stat-mechcond-mat.str-elhep-th

Catalog footprint

What is connected

18works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

18 published item(s)

preprint2023arXiv

Tensor RG approach to high-temperature fixed point

We study a renormalization group (RG) map for tensor networks that include two-dimensional lattice spin systems such as the Ising model. Numerical studies of such RG maps have been quite successful at reproducing the known critical behavior. In those numerical studies the RG map must be truncated to keep the dimension of the legs of the tensors bounded. Our tensors act on an infinite-dimensional Hilbert space, and our RG map does not involve any truncations. Our RG map has a trivial fixed point which represents the high-temperature fixed point. We prove that if we start with a tensor that is close to this fixed point tensor, then the iterates of the RG map converge in the Hilbert-Schmidt norm to the fixed point tensor. It is important to emphasize that this statement is not true for the simplest tensor network RG map in which one simply contracts four copies of the tensor to define the renormalized tensor. The linearization of this simple RG map about the fixed point is not a contraction due to the presence of so-called CDL tensors. Our work provides a first step towards the important problem of the rigorous study of RG maps for tensor networks in a neighborhood of the critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Absence of renormalization group pathologies in some critical Dyson-Ising ferromagnets

The Dyson-Ising ferromagnet is a one-dimensional Ising model with a power law interaction. When the power is between -1 and -2, the model has a phase transition. Van Enter and Le Ny proved that at sufficiently low temperatures the decimation renormalization group transformation is not defined in the sense that the renormalized measure is not a Gibbs measure. We consider a modified model in which the nearest neighbor couplings are much larger than the other couplings. For a family of Hamiltonians which includes critical cases, we prove that the first step of the renormalization group transformation can be rigorously defined for majority rule and decimation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The difference between a discrete and continuous harmonic measure

We consider a discrete-time, continuous-state random walk with steps uniformly distributed in a disk of radius of $h$. For a simply connected domain $D$ in the plane, let $ω_h(0,\cdot;D)$ be the discrete harmonic measure at $0\in D$ associated with this random walk, and $ω(0,\cdot;D)$ be the (continuous) harmonic measure at $0$. For domains $D$ with analytic boundary, we prove there is a bounded continuous function $σ_D(z)$ on $\partial D$ such that for functions $g$ which are in $C^{2+α}(\partial D)$ for some $α>0$ $$ \lim_{h\downarrow 0} \frac{\int_{\partial D} g(ξ) ω_h(0,|dξ|;D) -\int_{\partial D} g(ξ)ω(0,|dξ|;D)}{h} = \int_{\partial D}g(z) σ_D(z) |dz|. $$ We give an explicit formula for $σ_D$ in terms of the conformal map from $D$ to the unit disc. The proof relies on some fine approximations of the potential kernel and Green's function of the random walk by their continuous counterparts, which may be of independent interest.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The first order correction to the exit distribution for some random walks

We study three different random walk models on several two-dimensional lattices by Monte Carlo simulations. One is the usual nearest neighbor random walk. Another is the nearest neighbor random walk which is not allowed to backtrack. The final model is the smart kinetic walk. For all three of these models the distribution of the point where the walk exits a simply connected domain $D$ in the plane converges weakly to harmonic measure on $\partial D$ as the lattice spacing $δ\rightarrow 0$. Let $ω(0,|dz|;D)$ be harmonic measure for $D$, and let $ω_δ(0,|dz|;D)$ be the discrete harmonic measure for one of the random walk models. Our definition of the random walk models is unusual in that we average over the orientation of the lattice with respect to the domain. We are interested in the limit of $(ω_δ(0,|dz|;D)- ω(0,|dz|;D))/δ$. Our Monte Carlo simulations of the three models lead to the conjecture that this limit equals $c_{M,L} \, ρ_D(z) |dz|$, where the function $ρ_D(z)$ depends on the domain, but not on the model or lattice, and the constant $c_{M,L}$ depends on the model and on the lattice, but not on the domain. So there is a form of universality for this first order correction. We also give an explicit formula for the conjectured density $ρ_D$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Smart Kinetic Self-Avoiding Walk and Schramm-Loewner Evolution

The smart kinetic self-avoiding walk (SKSAW) is a random walk which never intersects itself and grows forever when run in the full-plane. At each time step the walk chooses the next step uniformly from among the allowable nearest neighbors of the current endpoint of the walk. In the full-plane a nearest neighbor is allowable if it has not been visited before and there is a path from the nearest neighbor to infinity through sites that have not been visited before. It is well known that on the hexagonal lattice the SKSAW in a bounded domain between two boundary points is equivalent to an interface in critical percolation, and hence its scaling limit is the chordal Schramm-Loewner evolution with kappa=6 (SLE_6). Like SLE there are variants of the SKSAW depending on the domain and the initial and terminal points. On the hexagonal lattice these variants have been shown to converge to the corresponding version of SLE_6. We conjecture that the scaling limit of all these variants on any regular lattice is the corresponding version of SLE_6. We test this conjecture for the square lattice by simulating the SKSAW in the full-plane and find excellent agreement with the predictions of full-plane SLE_6.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Conformal invariance predictions for the three-dimensional self-avoiding walk

If the three dimensional self-avoiding walk (SAW) is conformally invariant, then one can compute the hitting densities for the SAW in a half-space and in a sphere. The ensembles of SAW's used to define these hitting densities involve walks of arbitrary lengths, and so these ensembles cannot be directly studied by the pivot Monte Carlo algorithm for the SAW. We show that these mixed length ensembles should have the same scaling limit as certain weighted ensembles that only involve walks with a single length, thus providing a fast method for simulating these ensembles. Preliminary simulations which found good agreement between the predictions and Monte Carlo simulations for the SAW were reported in [14]. In this paper we present more accurate simulations testing the predictions and find even stronger support for the prediction that the SAW is conformally invariant in three dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Self-avoiding walks in a rectangle

A celebrated problem in numerical analysis is to consider Brownian motion originating at the centre of a $10 \times 1$ rectangle, and to evaluate the ratio of probabilities of a Brownian path hitting the short ends of the rectangle before hitting one of the long sides. For Brownian motion this probability can be calculated exactly \cite{BLWW04}. Here we consider instead the more difficult problem of a self-avoiding walk in the scaling limit, and pose the same question. Assuming that the scaling limit of SAW is conformally invariant, we evaluate, asymptotically, the same ratio of probabilities. For the SAW case we find the probability ratio is approximately 200 times greater than for Brownian motion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Lattice effects in the scaling limit of the two-dimensional self-avoiding walk

We consider the two-dimensional self-avoiding walk (SAW) in a simply connected domain that contains the origin. The SAW starts at the origin and ends somewhere on the boundary. The distribution of the endpoint along the boundary is expected to differ from the SLE partition function prediction for this distribution because of lattice effects that persist in the scaling limit. We give a precise conjecture for how to compute this lattice effect correction and support our conjecture with simulations. We also give a precise conjecture for the lattice corrections that persist in the scaling limit of the lambda-SAW walk.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Simulating self-avoiding walks in bounded domains

Let D be a domain in the plane containing the origin. We are interested in the ensemble of self-avoiding walks (SAW's) in D which start at the origin and end on the boundary of the domain. We introduce an ensemble of SAW's that we expect to have the same scaling limit. The advantage of our ensemble is that it can be simulated using the pivot algorithm. Our ensemble makes it possible to accurately study SLE predictions for the SAW in bounded simply connected domains. One such prediction is the distribution along the boundary of the endpoint of the SAW. We use the pivot algorithm to simulate our ensemble and study this density. In particular the lattice effects in this density that persist in the scaling limit are seen to be given by a purely local function.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The self-avoiding walk in a strip

We review the existence of the infinite length self-avoiding walk in the half plane and its relationship to bridges. We prove that this probability measure is also given by the limit as $β\rightarrow β_c-$ of the probability measure on all finite length walks $ω$ with the probability of $ω$ proportional to $β_c^{|ω|}$ where $|ω|$ is the number of steps in $ω$. The self-avoiding walk in a strip $\{z : 0<\Im(z)<y\}$ is defined by considering all self-avoiding walks $ω$ in the strip which start at the origin and end somewhere on the top boundary with probability proportional to $β_c^{|ω|}$ We prove that this probability measure may be obtained by conditioning the SAW in the half plane to have a bridge at height $y$. This observation is the basis for simulations to test conjectures on the distribution of the endpoint of the SAW in a strip and the relationship between the distribution of this strip SAW and SLE$_{8/3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transforming fixed-length self-avoiding walks into radial SLE_8/3

We conjecture a relationship between the scaling limit of the fixed-length ensemble of self-avoiding walks in the upper half plane and radial SLE with kappa=8/3 in this half plane from 0 to i. The relationship is that if we take a curve from the fixed-length scaling limit of the SAW, weight it by a suitable power of the distance to the endpoint of the curve and then apply the conformal map of the half plane that takes the endpoint to i, then we get the same probability measure on curves as radial SLE. In addition to a non-rigorous derivation of this conjecture, we support it with Monte Carlo simulations of the SAW. Using the conjectured relationship between the SAW and radial SLE, our simulations give estimates for both the interior and boundary scaling exponents. The values we obtain are within a few hundredths of a percent of the conjectured values.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Renormalization group maps for Ising models in lattice gas variables

Real space renormalization group maps, e.g., the majority rule transformation, map Ising type models to Ising type models on a coarser lattice. We show that each coefficient of the renormalized Hamiltonian in the lattice gas variables depends on only a finite number of values of the renormalized Hamiltonian. We introduce a method which computes the values of the renormalized Hamiltonian with high accuracy and so computes the coefficients in the lattice gas variables with high accuracy. For the critical nearest neighbor Ising model on the square lattice with the majority rule transformation, we compute over 1,000 different coefficients in the lattice gas variable representation of the renormalized Hamiltonian and study the decay of these coefficients. We find that they decay exponentially in some sense but with a slow decay rate. We also show that the coefficients in the spin variables are sensitive to the truncation method used to compute them.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical computations for the Schramm-Loewner Evolution

We review two numerical methods related to the Schramm-Loewner evolution (SLE). The first simulates SLE itself. More generally, it finds the curve in the half-plane that results from the Loewner equation for a given driving function. The second method can be thought of as the inverse problem. Given a simple curve in the half-plane it computes the driving function in the Loewner equation. This algorithm can be used to test if a given random family of curves in the half-plane is SLE by computing the driving process for the curves and testing if it is Brownian motion. More generally, this algorithm can be used to compute the driving process for random curves that may not be SLE. Most of the material presented here has appeared before. Our goal is to give a pedagogic review, illustrate some of the practical issues that arise in these computations and discuss some open problems.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Computing the Loewner driving process of random curves in the half plane

We simulate several models of random curves in the half plane and numerically compute their stochastic driving process (as given by the Loewner equation). Our models include models whose scaling limit is the Schramm-Loewner evolution (SLE) and models for which it is not. We study several tests of whether the driving process is Brownian motion. We find that just testing the normality of the process at a fixed time is not effective at determining if the process is Brownian motion. Tests that involve the independence of the increments of Brownian motion are much more effective. We also study the zipper algorithm for numerically computing the driving function of a simple curve. We give an implementation of this algorithm which runs in a time O(N^1.35) rather than the usual O(N^2), where N is the number of points on the curve.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A Fast Algorithm for Simulating the Chordal Schramm-Loewner Evolution

The Schramm-Loewner evolution (SLE) can be simulated by dividing the time interval into N subintervals and approximating the random conformal map of the SLE by the composition of N random, but relatively simple, conformal maps. In the usual implementation the time required to compute a single point on the SLE curve is O(N). We give an algorithm for which the time to compute a single point is O(N^p) with p<1. Simulations with kappa=8/3 and kappa=6 both give a value of p of approximately 0.4.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Conformal Invariance and Stochastic Loewner Evolution Predictions for the 2D Self-Avoiding Walk - Monte Carlo Tests

Simulations of the self-avoiding walk (SAW) are performed in a half-plane and a cut-plane (the complex plane with the positive real axis removed) using the pivot algorithm. We test the conjecture of Lawler, Schramm and Werner that the scaling limit of the two-dimensional SAW is given by Schramm's Stochastic Loewner Evolution (SLE). The agreement is found to be excellent. The simulations also test the conformal invariance of the SAW since conformal invariance would imply that if we map the walks in the cut-plane into the half plane using the conformal map z -> sqrt(z), then the resulting walks will have the same distribution as the SAW in the half plane. The simulations show excellent agreement between the distributions.

0 citations0 reviewsNo code linkedOpen access