BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. D. Doll

J. D. Doll contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.stat-mechphysics.chem-phphysics.comp-phcond-mat.mtrl-scicond-mat.softmath-phmath.MPmath.OCmath.PRquant-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

7 published item(s)

preprint2016arXiv

A large deviations analysis of certain qualitative properties of parallel tempering and infinite swapping algorithms

Parallel tempering, or replica exchange, is a popular method for simulating complex systems. The idea is to run parallel simulations at different temperatures, and at a given swap rate exchange configurations between the parallel simulations. From the perspective of large deviations it is optimal to let the swap rate tend to infinity and it is possible to construct a corresponding simulation scheme, known as infinite swapping. In this paper we propose a novel use of large deviations for empirical measures for a more detailed analysis of the infinite swapping limit in the setting of continuous time jump Markov processes. Using the large deviations rate function and associated stochastic control problems we consider a diagnostic based on temperature assignments, which can be easily computed during a simulation. We show that the convergence of this diagnostic to its a priori known limit is a necessary condition for the convergence of infinite swapping. The rate function is also used to investigate the impact of asymmetries in the underlying potential landscape, and where in the state space poor sampling is most likely to occur.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On Performance Measures for Infinite Swapping Monte Carlo Methods

We introduce and illustrate a number of performance measures for rare-event sampling methods. These measures are designed to be of use in a variety of expanded ensemble techniques including parallel tempering as well as infinite and partial infinite swapping approaches. Using a variety of selected applications we address questions concerning the variation of sampling performance with respect to key computational ensemble parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rare-Event Sampling: Occupation-Based Performance Measures for Parallel Tempering and Infinite Swapping Monte Carlo Methods

In the present paper we identify a rigorous property of a number of tempering-based Monte Carlo sampling methods, including parallel tempering as well as partial and infinite swapping. Based on this property we develop a variety of performance measures for such rare-event sampling methods that are broadly applicable, informative, and straightforward to implement. We illustrate the use of these performance measures with a series of applications involving the equilibrium properties of simple Lennard-Jones clusters, applications for which the performance levels of partial and infinite swapping approaches are found to be higher than those of conventional parallel tempering.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

An Infinite Swapping Approach to the Rare-Event Sampling Problem

We describe a new approach to the rare-event Monte Carlo sampling problem. This technique utilizes a symmetrization strategy to create probability distributions that are more highly connected and thus more easily sampled than their original, potentially sparse counterparts. After discussing the formal outline of the approach and devising techniques for its practical implementation, we illustrate the utility of the technique with a series of numerical applications to Lennard-Jones clusters of varying complexity and rare-event character.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Comparative Monte Carlo Efficiency by Monte Carlo Analysis

We propose a modified power method for computing the subdominant eigenvalue $λ_2$ of a matrix or continuous operator. Here we focus on defining simple Monte Carlo methods for its application. The methods presented use random walkers of mixed signs to represent the subdominant eigenfuction. Accordingly, the methods must cancel these signs properly in order to sample this eigenfunction faithfully. We present a simple procedure to solve this sign problem and then test our Monte Carlo methods by computing the $λ_2$ of various Markov chain transition matrices. We first computed ${λ_2}$ for several one and two dimensional Ising models, which have a discrete phase space, and compared the relative efficiencies of the Metropolis and heat-bath algorithms as a function of temperature and applied magnetic field. Next, we computed $λ_2$ for a model of an interacting gas trapped by a harmonic potential, which has a mutidimensional continuous phase space, and studied the efficiency of the Metropolis algorithm as a function of temperature and the maximum allowable step size $Δ$. Based on the $λ_2$ criterion, we found for the Ising models that small lattices appear to give an adequate picture of comparative efficiency and that the heat-bath algorithm is more efficient than the Metropolis algorithm only at low temperatures where both algorithms are inefficient. For the harmonic trap problem, we found that the traditional rule-of-thumb of adjusting $Δ$ so the Metropolis acceptance rate is around 50% range is often sub-optimal. In general, as a function of temperature or $Δ$, $λ_2$ for this model displayed trends defining optimal efficiency that the acceptance ratio does not. The cases studied also suggested that Monte Carlo simulations for a continuum model are likely more efficient than those for a discretized version of the model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convergence Characteristics of the Cumulant Expansion for Fourier Path Integrals

The cumulant representation of the Fourier path integral method is examined to determine the asymptotic convergence characteristics of the imaginary-time density matrix with respect to the number of path variables $N$ included. It is proved that when the cumulant expansion is truncated at order $p$, the asymptotic convergence rate of the density matrix behaves like $N^{-(2p+1)}$. The complex algebra associated with the proof is simplified by introducing a diagrammatic representation of the contributing terms along with an associated linked-cluster theorem. The cumulant terms at each order are expanded in a series such that the the asymptotic convergence rate is maintained without the need to calculate the full cumulant at order $p$. Using this truncated expansion of each cumulant at order $p$, the numerical cost in developing Fourier path integral expressions having convergence order $N^{-(2p+1)}$ is shown to be approximately linear in the number of required potential energy evaluations making the method promising for actual numerical implementation.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A Simple Method for Calculating Quantum Effects on the Temperature Dependence of Bimolecular Reaction Rates: An Application to $CH_4 + H \to CH_3 + H_2$

The temperature dependence of the rate of the reaction CH_4+H \to CH_3+H_2 is studied using classical collision theory with a temperature-dependent effective potential derived from a path integral analysis. Analytical expressions are obtained for the effective potential and for the rate constant. The rate constant expressions use a temperature-dependent activation energy. They give better agreement with the available experimental results than do previous empirical fits. Since all but one of the parameters in the present expressions are obtained from theory, rather than by fitting to experimental reaction rates, the expressions can be expected to be more dependable than purely empirical expressions at temperatures above 2000 K or below 350 K, where experimental results are not available.

0 citations0 reviewsNo code linkedOpen access