Researcher profile

B. Leimkuhler

B. Leimkuhler contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
3topics
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

4 published item(s)

preprint2020arXiv

2D mobile breather scattering in a hexagonal crystal lattice

We describe, for the first time, the full 2D scattering of long-lived breathers in a model hexagonal lattice of atoms. The chosen system, representing an idealized model of mica, combines a Lennard-Jones interatomic potential with an "egg-box" harmonic potential well surface. We investigate the dependence of breather properties on the ratio of the well depths associated to the interaction and on-site potentials. High values of this ratio lead to large spatial displacements in adjacent chains of atoms and thus enhance the two dimensional character of the quasi-one-dimensional breather solutions. This effect is further investigated during breather-breather collisions by following the constrained energy density function in time for a set of randomly excited mobile breather solutions. Certain collisions lead to 60$^\circ$ scattering, and collisions of mobile and stationary breathers can generate a rich variety of states.

preprint2014arXiv

Nonlinear propagating localized modes in a 2D hexagonal crystal lattice

In this paper we consider a 2D hexagonal crystal lattice model first proposed by Marin, Eilbeck and Russell in 1998. We perform a detailed numerical study of nonlinear propagating localized modes, that is, propagating discrete breathers and kinks. The original model is extended to allow for arbitrary atomic interactions, and to allow atoms to travel out of the unit cell. A new on-site potential is considered with a periodic smooth function with hexagonal symmetry. We are able to confirm the existence of long-lived propagating discrete breathers. Our simulations show that, as they evolve, breathers appear to localize in frequency space, i.e. the energy moves from sidebands to a main frequency band. Our numerical findings contribute to the open question of whether exact moving breather solutions exist in 2D hexagonal layers in physical crystal lattices.

preprint2014arXiv

Numerical simulations of nonlinear modes in mica: past, present and future

We review research on the role of nonlinear coherent phenomena (e.g breathers and kinks) in the formation linear decorations in mica crystal. The work is based on a new model for the motion of the mica hexagonal K layer, which allows displacement of the atoms from the unit cell. With a simple piece-wise polynomial inter-particle potential, we verify the existence of localized long-lived breathers in an idealized lattice at 0K. Moreover, our model allows us to observe long-lived localized kinks. We study the interactions of such localized modes along a lattice direction, and in addition demonstrate fully two dimensional scattering of such pulses for the first time. For large interatomic forces we observe a spreading horseshoe-shaped wave, a type of shock wave but with a breather profile.

preprint2014arXiv

On the long-time integration of stochastic gradient systems

This article addresses the weak convergence of numerical methods for Brownian dynamics. Typical analyses of numerical methods for stochastic differential equations focus on properties such as the weak order which estimates the asymptotic (stepsize $h \to 0$) convergence behavior of the error of finite time averages. Recently it has been demonstrated, by study of Fokker-Planck operators, that a non-Markovian numerical method [Leimkuhler and Matthews, 2013; Leimkuhler et al., 2013] generates approximations in the long time limit with higher accuracy order (2nd order) than would be expected from its weak convergence analysis (finite-time averages are 1st order accurate). In this article we describe the transition from the transient to the steady-state regime of this numerical method by estimating the time-dependency of the coefficients in an asymptotic expansion for the weak error, demonstrating that the convergence to 2nd order is exponentially rapid in time. Moreover, we provide numerical tests of the theory, including comparisons of the efficiencies of the Euler-Maruyama method, the popular 2nd order Heun method, and the non-Markovian method.