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

David J. Wales

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechcond-mat.mtrl-sciphysics.chem-phcond-mat.dis-nnhep-thphysics.comp-phastro-ph.IMastro-ph.COBiomoleculesComputationhep-lathep-phmath.AGmath.NAnlin.CDphysics.atm-clusphysics.data-anquant-ph

Catalog footprint

What is connected

18works
19topics
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)

preprint2022arXiv

Exact electronic states with shallow quantum circuits through global optimisation

Quantum computers promise to revolutionise electronic simulations by overcoming the exponential scaling of many-electron problems. While electronic wave functions can be represented using a product of fermionic unitary operators, shallow quantum circuits for exact states have not yet been achieved. We construct universal wave functions from gate-efficient, symmetry-preserving fermionic operators by introducing an algorithm that globally optimises the wave function in the discrete ansatz design and the continuous parameter spaces. Our approach maximises the accuracy that can be obtained with near-term quantum circuits. Highly accurate numerical simulations on strongly correlated molecules, including water and molecular nitrogen, and the condensed-matter Hubbard model, demonstrate that our algorithm reliably advances the state-of-the-art, defining a new paradigm for quantum simulations featuring strong electron correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Nested sampling for physical scientists

We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Kinetic Transition Networks for the Thomson Problem and Smale's 7th Problem

The Thomson Problem, arrangement of identical charges on the surface of a sphere, has found many applications in physics, chemistry and biology. Here we show that the energy landscape of the Thomson Problem for $N$ particles with $N=132, 135, 138, 141, 144, 147$ and $150$ is single funnelled, characteristic of a structure-seeking organisation where the global minimum is easily accessible. Algorithmically constructing starting points close to the global minimum of such a potential with spherical constraints is one of Smale's 18 unsolved problems in mathematics for the 21st century because it is important in the solution of univariate and bivariate random polynomial equations. By analysing the kinetic transition networks, we show that a randomly chosen minimum is in fact always `close' to the global minimum in terms of the number of transition states that separate them, a characteristic of small world networks.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Structural analysis of high-dimensional basins of attraction

We propose an efficient Monte Carlo method for the computation of the volumes of high-dimensional bodies with arbitrary shape. We start with a region of known volume within the interior of the manifold and then use the multistate Bennett acceptance-ratio method to compute the dimensionless free-energy difference between a series of equilibrium simulations performed within this object. The method produces results that are in excellent agreement with thermodynamic integration, as well as a direct estimate of the associated statistical uncertainties. The histogram method also allows us to directly obtain an estimate of the interior radial probability density profile, thus yielding useful insight into the structural properties of such a high-dimensional body. We illustrate the method by analyzing the effect of structural disorder on the basins of attraction of mechanically stable packings of soft repulsive spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turning intractable counting into sampling: computing the configurational entropy of three-dimensional jammed packings

We report a numerical calculation of the total number of disordered jammed configurations $Ω$ of $N$ repulsive, three-dimensional spheres in a fixed volume $V$. To make these calculations tractable, we increase the computational efficiency of the approach of Xu et al. (Phys. Rev. Lett. 106, 245502 (2011)) and Asenjo et al. (Phys. Rev. Lett. 112, 098002 (2014)) and we extend the method to allow computation of the configurational entropy as a function of pressure. The approach that we use computes the configurational entropy by sampling the absolute volume of basins of attraction of the stable packings in the potential energy landscape. We find a surprisingly strong correlation between the pressure of a configuration and the volume of its basin of attraction in the potential energy landscape. This relation is well described by a power law. Our methodology to compute the number of minima in the potential energy landscape should be applicable to a wide range of other enumeration problems in statistical physics, string theory, cosmology and machine learning, that aim to find the distribution of the extrema of a scalar cost function that depends on many degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Analysing kinetic transition networks for rare events

The graph transformation approach is a recently proposed method for computing mean first passage times, rates, and committor probabilities for kinetic transition networks. Here we compare the performance to existing linear algebra methods, focusing on large, sparse networks. We show that graph transformation provides a much more robust framework, succeeding when numerical precision issues cause the other methods to fail completely. These are precisely the situations that correspond to rare event dynamics for which the graph transformation was introduced.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Certification and the Potential Energy Landscape

Typically, there is no guarantee that a numerical approximation obtained using standard nonlinear equation solvers is indeed an actual solution, meaning that it lies in the quadratic convergence basin. Instead, it may lie only in the linear convergence basin, or even in a chaotic region, and hence not converge to the corresponding stationary point when further optimization is attempted. In some cases, these non-solutions could be misleading. Proving that a numerical approximation will quadratically converge to a stationary point is termed \textit{certification}. In this report, we provide details of how Smale's $α$-theory can be used to certify numerically obtained stationary points of a potential energy landscape, providing a \textit{mathematical proof} that the numerical approximation does indeed correspond to an actual stationary point, independent of the precision employed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Design principles for Bernal spirals and helices with tunable pitch

Using the framework of potential energy landscape theory, we describe two in silico designs for self-assembling helical colloidal superstructures based upon dipolar dumbbells and Janus-type building blocks, respectively. Helical superstructures with controllable pitch length are obtained using external magnetic field driven assembly of asymmetric dumbbells involving screened electrostatic as well as magnetic dipolar interactions. The pitch of the helix is tuned by modulating the Debye screening length over an experimentally accessible range. The second design is based on building blocks composed of rigidly linked spheres with short-range anisotropic interactions, which are predicted to self-assemble into Bernal spirals. These spirals are quite flexible, and longer helices undergo rearrangements via cooperative, hinge-like moves, in agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Superposition Enhanced Nested Sampling

The theoretical analysis of many problems in physics, astronomy and applied mathematics requires an efficient numerical exploration of multimodal parameter spaces that exhibit broken ergodicity. Monte Carlo methods are widely used to deal with these classes of problems, but such simulations suffer from a ubiquitous sampling problem: the probability of sampling a particular state is proportional to its entropic weight. Devising an algorithm capable of sampling efficiently the full phase space is a long-standing problem. Here we report a new hybrid method for the exploration of multimodal parameter spaces exhibiting broken ergodicity. Superposition enhanced nested sampling (SENS) combines the strengths of global optimization with the unbiased/athermal sampling of nested sampling, greatly enhancing its efficiency with no additional parameters. We report extensive tests of this new approach for atomic clusters that are known to have energy landscapes for which conventional sampling schemes suffer from broken ergodicity. We also introduce a novel parallelization algorithm for nested sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Certifying the Potential Energy Landscape

It is highly desirable for a numerical approximation of a stationary point for a potential energy landscape to lie in the quadratic convergence basin of that stationary point. However, it is possible that an approximation may lie only in the linear convergence basin, or even in a chaotic region, and hence not converge to the actual stationary point when further optimization is attempted. Proving that a numerical approximation will quadratically converge to the associated stationary point is termed certifying the numerical approximation. We employ Smale's α-theory to stationary points, providing a certification that serves as a mathematical proof that the numerical approximation does indeed correspond to an actual stationary point, independent of the precision employed. As a practical example, employing recently developed certification algorithms, we show how the α-theory can be used to certify all the known minima and transition states of Lennard-Jones LJ$_{N}$ atomic clusters for N = 7, ...,14.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Defect Motifs for Constant Mean Curvature Surfaces

The energy landscapes of electrostatically charged particles embedded on constant mean curvature surfaces are analysed for a wide range of system size, curvature, and interaction potentials. The surfaces are taken to be rigid, and the basin-hopping method is used to locate the putative global minimum structures. The defect motifs favoured by potential energy agree with experimental observations for colloidal systems: extended defects (scars and pleats) for weakly positive and negative Gaussian curvatures, and isolated defects for strongly negative Gaussian curvatures. Near the phase boundary between these regimes the two motifs are in strong competition, as evidenced from the appearance of distinct funnels in the potential energy landscape. We also report a novel defect motif consisting of pentagon pairs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Energy landscapes, structural topologies and rearrangement mechanisms in clusters of dipolar particles

Clusters of spherical particles with isotropic attraction favour compact structures that maximise the number of energetically optimal nearest-neighbour interactions. In contrast, dipolar interactions lead to the formation of chains with a low coordination number. When both isotropic and dipolar interactions are present, the competition between them can lead to intricate knot, link and coil structures. Here, we investigate how these structures may self-organise and interconvert in clusters bound by the Stockmayer potential (Lennard-Jones plus point dipole). We map out the low-lying region of the energy landscape using disconnectivity graphs to follow how it evolves as the strength of the dipolar interactions increases. From comprehensive surveys of isomerisation pathways, we identify a number of rearrangement mechanisms that allow the topology of chain-like structures to interconvert.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Potential Energy Landscapes for the 2D XY Model: Minima, Transition States and Pathways

We describe a numerical study of the potential energy landscape for the two-dimensional XY model (with no disorder), considering up to 100 spins and CPU and GPU implementations of local optimization, focusing on minima and saddles of index one (transition states). We examine both periodic and anti-periodic boundary conditions, and show that the number of stationary points located increases exponentially with increasing lattice size. The corresponding disconnectivity graphs exhibit funneled landscapes; the global minima are readily located because they exhibit relatively large basins of attraction compared to the higher energy minima as the lattice size increases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Visualizing Basins of Attraction for Different Minimization Algorithms

We report a study of the basins of attraction for potential energy minima defined by different minimisation algorithms for an atomic system. We find that whereas some minimisation algorithms produce compact basins, others produce basins with complex boundaries or basins consisting of disconnected parts. Such basins deviate from the "correct" basin of attraction defined by steepest-descent pathways, and the differences can be controlled to some extent by adjustment of the maximum step size. The choice of the most convenient minimisation algorithm depends on the problem in hand. We show that while L-BFGS is the fastest minimiser, the FIRE algorithm is also quite fast, and can lead to less fragmented basins of attraction.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydrogen on graphene under stress: Molecular dissociation and gap opening

Density functional calculations are employed to study the molecular dissociation of hydrogen on graphene, the diffusion of chemisorbed atomic species, and the electronic properties of the resulting hydrogen on graphene system. Our results show that applying stress to the graphene substrate can lower the barrier to dissociation of molecular hydrogen by a factor of six, and change the process from endothermic to exothermic. These values for the barrier and the heat of reaction, unlike the zero stress values, are compatible with the time scales observed in experiments. Diffusion, on the other hand, is not greatly modified by stress. We analyse the electronic structure for configurations relevant to molecular dissociation and adsorption of atomic hydrogen on a graphene single layer. An absolute band gap of 0.5 eV is found for the equilibrium optimum configuration for a narrow range of coverages ($θ\approx 0.25$). This value is in good agreement with experiment [Elias et al., Science {\bf 323}, 610 (2009)].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Energy Landscape, Folding Pathways and the Kinetics of a Knotted Protein

The folding pathway and rate coefficients of the folding of a knotted protein are calculated for a potential energy function with minimal energetic frustration. A kinetic transition network is constructed using the discrete path sampling approach, and the resulting potential energy surface is visualized by constructing disconnectivity graphs. Owing to topological constraints, the low-lying portion of the landscape consists of three distinct regions, corresponding to the native knotted state and to configurations where either the N- or C-terminus is not yet folded into the knot. The fastest folding pathways from denatured states exhibit early formation of the N-terminus portion of the knot and a rate-determining step where the C-terminus is incorporated. The low-lying minima with the N-terminus knotted and the C-terminus free therefore constitute an off-pathway intermediate for this model. The insertion of both the N- and C-termini into the knot occur late in the folding process, creating large energy barriers that are the rate limiting steps in the folding process. When compared to other protein folding proteins of a similar length, this system folds over six orders of magnitude more slowly.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Supercooled Lennard-Jones Liquids and Glasses: a Kinetic Monte Carlo Approach

A kinetic Monte Carlo (KMC) method is used to study the structural properties and dynamics of a supercooled binary Lennard-Jones liquid around the glass transition temperature. This technique permits us to explore the potential energy surface and barrier distributions without suffering the exponential slowing down at low temperature that affects molecular dynamics simulations. In agreement with previous studies we observe a distinct change in behaviour around $T=0.45$, close to the dynamical transition temperature $T_c$ of mode coupling theory (MCT). Below this temperature the number of different local minima visited by the system for the same number of KMC steps decreases by more than an order of magnitude. The mean number of atoms involved in each jump between local minima and the average distance they move also decreases significantly, and new features appear in the partial structure factor. Above $T~0.45$ the probability distribution for the magnitude of the atomic displacement per KMC step exhibits an exponential decay, which is only weakly temperature dependent.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Rearrangements of Water Dimer and Hexamer

Rearrangement mechanisms of the water dimer and the cage form of the water hexamer are examined theoretically with particular reference to tunneling splittings and spectroscopy. The three lowest barrier rearrangements of the water dimer are characterized by {\it ab initio} methods and compared with the results of previous constrained calculations. The acceptor-tunneling pathway does not proceed via a direct rotation around the $C_2$ axis of the acceptor, but rather via relatively asynchronous rotation of the donor about the hydrogen bond and an associated `wag' of the acceptor. Rearrangements between different cage isomers of the water hexamer are studied for two empirical potentials. The experimentally observed triplet splittings may be the result of flip and bifurcation rearrangements of the two single-donor, single-acceptor monomers. Two-dimensional quantum calculations of the nuclear dynamics suggest that delocalization over more than one cage isomer may occur, especially in excited states.

0 citations0 reviewsNo code linkedOpen access