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

Brian Van Koten

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

ResearcherUnclaimed source record
math.NAcond-mat.mtrl-sciphysics.comp-phBiomoleculescond-mat.stat-mechmath.PRMethodologyNumerical Analysis

Catalog footprint

What is connected

7works
8topics
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

7 published item(s)

preprint2020arXiv

Stratification as a general variance reduction method for Markov chain Monte Carlo

The Eigenvector Method for Umbrella Sampling (EMUS) belongs to a popular class of methods in statistical mechanics which adapt the principle of stratified survey sampling to the computation of free energies. We develop a detailed theoretical analysis of EMUS. Based on this analysis, we show that EMUS is an efficient general method for computing averages over arbitrary target distributions. In particular, we show that EMUS can be dramatically more efficient than direct MCMC when the target distribution is multimodal or when the goal is to compute tail probabilities. To illustrate these theoretical results, we present a tutorial application of the method to a problem from Bayesian statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Eigenvector method for umbrella sampling enables error analysis

Umbrella sampling efficiently yields equilibrium averages that depend on exploring rare states of a model by biasing simulations to windows of coordinate values and then combining the resulting data with physical weighting. Here, we introduce a mathematical framework that casts the step of combining the data as an eigenproblem. The advantage to this approach is that it facilitates error analysis. We discuss how the error scales with the number of windows. Then, we derive a central limit theorem for averages that are obtained from umbrella sampling. The central limit theorem suggests an estimator of the error contributions from individual windows, and we develop a simple and computationally inexpensive procedure for implementing it. We demonstrate this estimator for simulations of the alanine dipeptide and show that it emphasizes low free energy pathways between stable states in comparison to existing approaches for assessing error contributions. We discuss the possibility of using the estimator and, more generally, the eigenvector method for umbrella sampling to guide adaptation of the simulation parameters to accelerate convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sharp entrywise perturbation bounds for Markov chains

For many Markov chains of practical interest, the invariant distribution is extremely sensitive to perturbations of some entries of the transition matrix, but insensitive to others; we give an example of such a chain, motivated by a problem in computational statistical physics. We have derived perturbation bounds on the relative error of the invariant distribution that reveal these variations in sensitivity. Our bounds are sharp, we do not impose any structural assumptions on the transition matrix or on the perturbation, and computing the bounds has the same complexity as computing the invariant distribution or computing other bounds in the literature. Moreover, our bounds have a simple interpretation in terms of hitting times, which can be used to draw intuitive but rigorous conclusions about the sensitivity of a chain to various types of perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysis of Blended Atomistic/Continuum Hybrid Methods

We present a comprehensive error analysis of two prototypical atomistic-to-continuum coupling methods of blending type: the energy-based and the force-based quasicontinuum methods. Our results are valid in two and three dimensions, for finite range many-body interactions (e.g., EAM type), and in the presence of lattice defects (we consider point defects and dislocations). The two key ingredients in the analysis are (i) new force and energy consistency error estimates; and (ii) a new technique for proving energy norm stability of a/c couplings that requires only the assumption that the exact atomistic solution is a stable equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Symmetries of 2-lattices and second order accuracy of the Cauchy--Born Model

We show that the Cauchy--Born model of a single-species 2-lattice is second order if the atomistic and continuum kinematics are connected in a novel way. Our proof uses a generalization to 2-lattices of the point symmetry of Bravais lattices. Moreover, by identifying similar symmetries in multi-species pair interaction models, we construct a new stored energy density, using shift-gradients but not strain gradients, that is also second order accurate. These results can be used to develop highly accurate continuum models and atomistic/continuum coupling methods for materials such as graphene, hcp metals, and shape memory alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Computational and Theoretical Investigation of the Accuracy of Quasicontinuum Methods

We give computational results to study the accuracy of several quasicontinuum methods for two benchmark problems - the stability of a Lomer dislocation pair under shear and the stability of a lattice to plastic slip under tensile loading. We find that our theoretical analysis of the accuracy near instabilities for one-dimensional model problems can successfully explain most of the computational results for these multi-dimensional benchmark problems. However, we also observe some clear discrepancies, which suggest the need for additional theoretical analysis and benchmark problems to more thoroughly understand the accuracy of quasicontinuum methods.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of Energy-Based Blended Quasicontinuum Approximations

The development of patch test consistent quasicontinuum energies for multi-dimensional crystalline solids modeled by many-body potentials remains a challenge. The original quasicontinuum energy (QCE) has been implemented for many-body potentials in two and three space dimensions, but it is not patch test consistent. We propose that by blending the atomistic and corresponding Cauchy-Born continuum models of QCE in an interfacial region with thickness of a small number $k$ of blended atoms, a general quasicontinuum energy (BQCE) can be developed with the potential to significantly improve the accuracy of QCE near lattice instabilities such as dislocation formation and motion. In this paper, we give an error analysis of the blended quasicontinuum energy (BQCE) for a periodic one-dimensional chain of atoms with next-nearest neighbor interactions. Our analysis includes the optimization of the blending function for an improved convergence rate. We show that the $\ell^2$ strain error for the non-blended QCE energy (QCE), which has low order $\text{O}(ε^{1/2})$ where $ε$ is the atomistic length scale, can be reduced by a factor of $k^{3/2}$ for an optimized blending function where $k$ is the number of atoms in the blending region. The QCE energy has been further shown to suffer from a O$(1)$ error in the critical strain at which the lattice loses stability. We prove that the error in the critical strain of BQCE can be reduced by a factor of $k^2$ for an optimized blending function, thus demonstrating that the BQCE energy for an optimized blending function has the potential to give an accurate approximation of the deformation near lattice instabilities such as crack growth.

0 citations0 reviewsNo code linkedOpen access