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 M. Wilkins

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

ResearcherUnclaimed source record
physics.chem-phquant-phmath-phmath.MPphysics.comp-phBiological Physicscond-mat.stat-mechQuantitative Methods

Catalog footprint

What is connected

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

8 published item(s)

preprint2022arXiv

Charge gradients around dendritic voids cause nanoscale inhomogeneities in liquid water

Water is the matrix of life and is generally considered a homogeneous, uniform, liquid. Recent experiments report that water is instead a two-state liquid. However, subsequently, these findings were contested. The structure of water and whether we should think of it as uniform therefore remains an open question. Here, we report femtosecond elastic second harmonic scattering (fs-ESHS) of liquid water in comparison to an isotropic liquid (CCl4) and show that water is indeed a non-uniform liquid, The coherent fs-ESHS intensity was interpreted, using molecular dynamics simulations, as arising from charge density fluctuations and consequentially enhanced nanoscale polarizabilities around transient voids having an average lifetime of 300 fs. Although voids were also present in CCl4, they were not characterized by hydrogen bond defects and did not show strong polarizability fluctuations, leading to fs-ESHS of an isotropic liquid. The voids increased in number at higher temperatures above room temperature, in agreement with the fs-ESHS results. The measured nanoscale-femtosecond inhomogeneities therefore do not necessarily relate to the proposed two state model of water, but instead underscore the elusive nature of liquid water, and undoubtedly have an impact on any type of transition that occurs in liquid water.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Hydrogen Bonding and Nuclear Quantum Effects in Clays

Hydrogen bonds are of paramount importance in the chemistry of clays, mediating the interaction between the clay surface and water, and for some materials between separate layers. It is well-established that the accuracy of a computational model for clays depends on the level of theory at which the electronic structure is treated. However, for hydrogen-bonded systems the motion of light H nuclei on the electronic potential energy surface is often affected by quantum delocalization. Using path integral molecular dynamics, we show that nuclear quantum effects lead to a relatively small change in the structure of clays, but one that is comparable to the variation incurred by treating the clay at different levels of electronic structure theory. Accounting for quantum effects weakens the hydrogen bonds in clays, with H-bonds between different layers of the clay affected more than those within the same layer; this is ascribed to the fact that the confinement of an H atom inside a layer is independent of its participation in hydrogen-bonding. More importantly, the weakening of hydrogen bonds by nuclear quantum effects causes changes in the vibrational spectra of these systems, significantly shifting the O--H stretching peaks and meaning that in order to fully understand these spectra by computational modelling, both electronic and nuclear quantum effects must be included. We show that after reparametrization of the popular CLAYFF model for clays, the O--H stretching region of their vibrational spectra better matches the experimental one, with no detriment to the model's agreement with other experimental properties.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Inexpensive modelling of quantum dynamics using path integral generalized Langevin equation thermostats

The properties of molecules and materials containing light nuclei are affected by their quantum mechanical nature. Modelling these quantum nuclear effects accurately requires computationally demanding path integral techniques. Considerable success has been achieved in reducing the cost of such simulations by using generalized Langevin dynamics to induce frequency-dependent fluctuations. Path integral generalized Langevin equation methods, however, have this far been limited to the study of static, thermodynamic properties due to the large perturbation to the system's dynamics induced by the aggressive thermostatting. Here we introduce a post-processing scheme, based on analytical estimates of the dynamical perturbation induced by the generalized Langevin dynamics, that makes it possible to recover meaningful time correlation properties from a thermostatted trajectory. We show that this approach yields spectroscopic observables for model and realistic systems which have an accuracy comparable to much more demanding approximate quantum dynamics techniques based on full path integral simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Theoretical Investigation Into Energy Transfer In Photosynthetic Open Quantum Systems

This thesis looks at the electronic energy transfer in the Fenna-Matthews-Olson complex, in which evidence of long-lived coherence has been observed in 2-dimensional infrared experiments. I use three techniques: the numerically exact Hierarchical Equations of Motion, and the perturbative Redfield and Foerster theories, the latter of which ignores quantum coherence in the transfer. Both of the approximate methods perform very well - and while oscillations in site populations (a hallmark of coherence) are present in the exact transfer dynamics and absent in the dynamics of Foerster theory, the latter gives a reasonable prediction of transfer rates and steady-state populations, despite being incoherent - suggesting that coherence is not vital for the dynamics of transfer. Since Foerster theory is very inexpensive to run and performs so well, I then apply it to calculate the effects of static disorder in bacteriochlorophyll site energies and of a more structured spectral density. Ultimately, the energy transfer in the complex is found to be very robust to changes in its environment, which is advantageous for its biological function.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nuclear quantum effects in water exchange around lithium and fluoride ions

We employ classical and ring polymer molecular dynamics simulations to study the effect of nuclear quantum fluctuations on the structure and the water exchange dynamics of aqueous solutions of lithium and fluoride ions. While we obtain reasonably good agreement with experimental data for solutions of lithium by augmenting the Coulombic interactions between the ion and the water molecules with a standard Lennard-Jones ion-oxygen potential, the same is not true for solutions of fluoride, for which we find that a potential with a softer repulsive wall gives much better agreement. A small degree of destabilization of the first hydration shell is found in quantum simulations of both ions when compared with classical simulations, with the shell becoming less sharply defined and the mean residence time of the water molecules in the shell decreasing. In line with these modest differences, we find that the mechanisms of the exchange processes are unaffected by quantization, so a classical description of these reactions gives qualitatively correct and quantitatively reasonable results. We also find that the quantum effects in solutions of lithium are larger than in solutions of fluoride. This is partly due to the stronger interaction of lithium with water molecules, partly due to the lighter mass of lithium, and partly due to competing quantum effects in the hydration of fluoride, which are absent in the hydration of lithium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Why quantum coherence is not important in the Fenna-Matthews-Olsen Complex

We develop and present an improvement to the conventional technique for solving the Hierarchical Equations of Motion which reduces the memory cost by more than 75% while retaining the same convergence rate and accuracy. This allows for a full calculation of the population dynamics of the 24-site FMO trimer for long timescales with very little effort, and we present the first fully converged, exact results for the 7-site subsystem of the monomer, and for the full 24-site trimer. Owing to this new approach, our numerically exact 24-site, 2-exponential results are the most demanding HEOM calculations performed to date, surpassing the 50-site, 1-exponential results of Strumpfer and Schulten [2012, J. Chem. Thy. & Comp., 8, 2808]. We then show where our exact 7-site results deviate from the approximation of Ishizaki and Fleming [2009, Proc. Natl. Acad. Sci. USA, 106, 17255]. Our exact results are then compared to calculations using the incoherent Forster theory, and it is found that the energy transfer from the antenna to the reaction center occurs more than 50 times faster than the fluorescence lifetime of the excitation, whether or not coherence is considered. This means that coherence is not likely to improve the efficiency of the photosynthesis. In fact, the incoherent theory often tends to over-predict the rates of energy transfer, suggesting that in some cases electronic coherence may actually slow down the photosynthetic process.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analytic influence functionals for numerical Feynman integrals in most open quantum systems

Fully analytic formulas, which do not involve any numerical integration, are derived for the discretized influence functionals of a very extensive assortment of spectral distributions. For Feynman integrals derived using the Trotter splitting and Strang splitting, we present general formulas for the discretized influence functionals in terms of proper integrals of the bath response function. When an analytic expression exists for the bath response function, these integrals can almost always be evaluated analytically. In cases where these proper integrals cannot be integrated analytically, numerically computing them is much faster and less error-prone than calculating the discretized influence functionals in the traditional way, which involves numerically calculating integrals whose bounds are both infinite. As an example, we present the analytic discretized influence functional for a bath response function of the form α(t)=\sum_{j}^{K}p_{j}e^{Ω_{j}t}, which is a natural form for many spectral distribution functions (including the very popular Lorentz-Drude/Debye function), and for other spectral distribution functions it is a form that is easily obtainable by a least-squares fit . Evaluating our analytic formulas for this example case is much faster and easier to implement than numerically calculating the discretized influence funcitonal in the traditional way. In the appendix we provide analytic expressions for p_{j} and Ω_{j} for a variety of spectral distribution forms, and as a second example we provide the analytic bath response function and analytic influence functionals for spectral distributions of the form J(ω)\proptoω^{s}e^{-(ω/ω_{c})^{q}}. The value of the analytic expression for this bath response function extends beyond its use for calculating Feynman integrals.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Optimal representation of the bath response function & fast calculation of influence functional coefficients in open quantum systems with BATHFIT 1

Today's most popular techniques for accurately calculating the dynamics of the reduced density operator in an open quantum system, either require, or gain great computational benefits, from representing the bath response function a(t) in the form a(t)=Σ_k^K p_k e^{O_k t} . For some of these techniques, the number of terms in the series K plays the lead role in the computational cost of the calculation, and is therefore often a limiting factor in simulating open quantum system dynamics. We present an open source MATLAB program called BATHFIT 1, whose input is any spectral distribution functions J(w) or bath response function, and whose output attempts to be the set of parameters {p_k,w_k}_k=1^K such that for a given value of K, the series Σ_k^k p_k e^{O_k t} is as close as possible to a(t). This should allow the user to represent a(t) as accurately as possible with as few parameters as possible. The program executes non-linear least squares fitting, and for a very wide variety of spectral distribution functions, competent starting parameters are used for these fits. For most forms of J(w), these starting parameters, and the exact a(t) corresponding to the given J(w), are calculated using the recent Pade decomposition technique - therefore this program can also be used to merely implement the Pade decomposition for these spectral distribution functions; and it can also be used just to efficiently and accurately calculate a(t) for any given J(w) . The program also gives the J(w) corresponding to a given a(t), which may allow one to assess the quality (in the w-domain) of a representation of a(t) being used. Finally, the program can calculate the discretized influence functional coefficients for any J(w), and this is computed very efficiently for most forms of J(w) by implementing the recent technique published in [Quantum Physics Letters (2012) 1 (1) pg. 35].

0 citations0 reviewsNo code linkedOpen access