Source author record

Thomas E. Markland

Thomas E. Markland appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

24 published item(s)

preprint2022arXiv

Optically induced anisotropy in time-resolved scattering: Imaging molecular scale structure and dynamics in disordered media with experiment and theory

Time-resolved scattering experiments enable imaging of materials at the molecular scale with femtosecond time resolution. However, in disordered media they provide access to just one radial dimension thus limiting the study of orientational structure and dynamics. Here we introduce a rigorous and practical theoretical framework for predicting and interpreting experiments combining optically induced anisotropy and time-resolved scattering. Using impulsive nuclear Raman and ultrafast X-ray scattering experiments of chloroform and simulations, we demonstrate that this framework can accurately predict and elucidate both the spatial and temporal features of these experiments.

preprint2022arXiv

Simulating nuclear and electronic quantum effects in enzymes

An accurate treatment of the structures and dynamics that lead to enhanced chemical reactivity in enzymes requires explicit treatment of both electronic and nuclear quantum effects. The former can be captured in ab initio molecular dynamics (AIMD) simulations while the latter can be included by performing ab initio path integral molecular dynamics (AI-PIMD) simulations. Both AIMD and AI-PIMD simulations have traditionally been computationally prohibitive for large enzymatic systems. Recent developments in streaming computer architectures and new algorithms to accelerate path integral simulations now make these simulations practical for biological systems, allowing elucidation of enzymatic reactions in unprecedented detail. In this chapter, we summarize these recent developments and discuss practical considerations for applying AIMD and AI-PIMD simulations to enzymes.

preprint2020arXiv

Excited state diabatization on the cheap using DFT: Photoinduced electron and hole transfer

Excited state electron and hole transfer underpin fundamental steps in processes such as exciton dissociation at photovoltaic heterojunctions, photoinduced charge transfer at electrodes, and electron transfer in photosynthetic reaction centers. Diabatic states corresponding to charge or excitation localized species, such as locally excited and charge transfer states, provide a physically intuitive framework to simulate and understand these processes. However, obtaining accurate diabatic states and their couplings from adiabatic electronic states generally leads to inaccurate results when combined with low-tier electronic structure methods, such as time dependent density functional theory (TDDFT), and exorbitant computational cost when combined with high-level wavefunction-based methods. Here we introduce a DFT-based diabatization scheme, Δ-ALMO(MSDFT2), which directly constructs the diabatic states using absolutely localized molecular orbitals (ALMOs). We demonstrate that our method, which combines ALMO calculations with the ΔSCF technique to construct electronically excited diabatic states and obtains their couplings with charge-transfer states using our MSDFT2 scheme, gives accurate results for excited state electron and hole transfer in both charged and uncharged systems that underlie DNA repair, charge separation in donor-acceptor dyads, chromophore-to-solvent electron transfer, and singlet fission. This framework for the accurate and efficient construction of excited state diabats and evaluation of their couplings directly from DFT thus offers a route to simulate and elucidate photoinduced electron and hole transfer in large disordered systems, such as those encountered in the condensed phase.

preprint2020arXiv

Exploiting machine learning to efficiently predict multidimensional optical spectra in complex environments

The excited state dynamics of chromophores in complex environments determine a range of vital biological and energy capture processes. Time-resolved, multidimensional optical spectroscopies provide a key tool to investigate these processes. Although theory has the potential to decode these spectra in terms of the electronic and atomistic dynamics, the need for large numbers of excited state electronic structure calculations severely limits first principles predictions of multidimensional optical spectra for chromophores in the condensed phase. Here, we leverage the locality of chromophore excitations to develop machine learning models to predict the excited state energy gap of chromophores in complex environments for efficiently constructing linear and multidimensional optical spectra. By analyzing the performance of these models, which span a hierarchy of physical approximations, across a range of chromophore-environment interaction strengths, we provide strategies for the construction of ML models that greatly accelerate the calculation of multidimensional optical spectra from first principles.

preprint2020arXiv

Quantum kinetic energy and isotope fractionation in aqueous ionic solutions

At room temperature, the quantum contribution to the kinetic energy of a water molecule exceeds the classical contribution by an order of magnitude. The quantum kinetic energy (QKE) of a water molecule is modulated by its local chemical environment and leads to uneven partitioning of isotopes between different phases in thermal equilibrium, which would not occur if the nuclei behaved classically. In this work, we use ab initio path integral simulations to show that QKEs of the water molecules and the equilibrium isotope fractionation ratios of the oxygen and hydrogen isotopes are sensitive probes of the hydrogen bonding structures in aqueous ionic solutions. In particular, we demonstrate how the QKE of water molecules in path integral simulations can be decomposed into translational, rotational and vibrational degrees of freedom, and use them to determine the impact of solvation on different molecular motions. By analyzing the QKEs and isotope fractionation ratios, we show how the addition of the Na$^+$, Cl$^-$ and HPO$_4^{2-}$ ions perturbs the competition between quantum effects in liquid water and impacts their local solvation structures.

preprint2019arXiv

Accurate and efficient DFT-based diabatization for hole and electron transfer using absolutely localized molecular orbitals

Diabatic states and the couplings between them are important for quantifying, elucidating, and predicting the rates and mechanisms of many chemical and biochemical processes. Here, we propose and investigate approaches to accurately compute diabatic couplings from density functional theory (DFT) using absolutely localized molecular orbitals (ALMOs). ALMOs provide an appealing approach to generate variationally optimized diabatic states and obtain their associated forces that allows for the relaxation of the donor and acceptor orbitals in a way that is internally consistent in how the method treats both the donor and acceptor states. Here, we show that one can obtain more accurate electronic couplings between ALMO-based diabats by employing the symmetrized transition density matrix to evaluate the exchange-correlation contribution. We demonstrate that this approach yields accurate results in comparison to other commonly used DFT-based diabatization methods across a wide array of electron and hole transfer processes occurring in systems ranging from conjugated organic molecules, such as thiophene and pentacene, to DNA base pairs. We also show that this approach yields accurate diabatic couplings even when combined with lower tiers of the DFT hierarchy, opening the door to combining it with quantum dynamics approaches to provide an ab initio treatment of nonadiabatic processes in the condensed phase.

preprint2016arXiv

Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory

Path integral molecular dynamics simulations, combined with an ab initio evaluation of interactions using electronic structure theory, incorporate the quantum mechanical nature of both the electrons and nuclei, which are essential to accurately describe systems containing light nuclei. However, path integral simulations have traditionally required a computational cost around two orders of magnitude greater than treating the nuclei classically, making them prohibitively costly for most applications. Here we show that the cost of path integral simulations can be dramatically reduced by extending our ring polymer contraction approach to ab initio molecular dynamics simulations. By using density functional tight binding as a reference system, we show that our ab initio ring polymer contraction (AI-RPC) scheme gives rapid and systematic convergence to the full path integral density functional theory result. We demonstrate the efficiency of this approach in ab initio simulations of liquid water and the reactive protonated and deprotonated water dimer systems. We find that the vast majority of the nuclear quantum effects are accurately captured using contraction to just the ring polymer centroid, which requires the same number of density functional theory calculations as a classical simulation. Combined with a multiple time step scheme using the same reference system, which allows the time step to be increased, our implementation is as fast as a typical classical ab initio molecular dynamics simulation and 35$\times$ faster than a full path integral calculation, while still exactly including the quantum sampling of nuclei. This development thus offers a route to routinely include nuclear quantum effects in ab initio molecular dynamics simulations at negligible computational cost.

preprint2016arXiv

Generalized Quantum Master Equations In and Out of Equilibrium: When Can One Win?

Generalized quantum master equations (GQMEs) are an important tool in modeling chemical and physical processes. For a large number of problems it has been shown that exact and approximate quantum dynamics methods can be made dramatically more efficient, and in the latter case more accurate, by proceeding via the GQME formalism. However, there are many situations where utilizing the GQME approach seems to offer no advantage over a direct evaluation of the property of interest. Here we provide a more detailed understanding of the conditions under which these methods will offer benefits. In particular, we derive exact expressions for the memory kernel for systems both in and out of equilibrium, and show the conditions under which these expressions will be guaranteed to return a result identical to that obtained from direct simulation. We also show the conditions which approximate methods must satisfy if they are to offer different results when used in conjunction with the GQME formalism. These exact analytical results thus provide new insights as to when proceeding via the GQME approach can improve the accuracy or efficiency of simulations.

preprint2015arXiv

Accurate nonadiabatic quantum dynamics on the cheap: making the most of mean field theory with master equations

In this article we show how Ehrenfest mean field theory can be made both a more accurate and efficient method to treat nonadiabatic quantum dynamics by combining it with the generalized quantum master equation framework. The resulting mean field generalized quantum master equation (MF-GQME) approach is a non-perturbative and non-Markovian theory to treat open quantum systems without any restrictions on the form of the Hamiltonian that it can be applied to. By studying relaxation dynamics in a wide range of dynamical regimes, typical of charge and energy transfer, we show that MF-GQME provides a much higher accuracy than a direct application of mean field theory. In addition, these increases in accuracy are accompanied by computational speed-ups of between one and two orders of magnitude that become larger as the system becomes more nonadiabatic. This combination of quantum-classical theory and master equation techniques thus makes it possible to obtain the accuracy of much more computationally expensive approaches at a cost lower than even mean field dynamics, providing the ability to treat the quantum dynamics of atomistic condensed phase systems for long times.

preprint2015arXiv

Efficient and accurate surface hopping for long time nonadiabatic quantum dynamics

The quantum-classical Liouville equation offers a rigorous approach to nonadiabatic quantum dynamics based on surface hopping type trajectories. However, in practice the applicability of this approach has been limited to short times owing to unfavorable numerical scaling. In this paper we show that this problem can be alleviated by combining it with a formally exact generalized quantum master equation treatment. This allows dramatic improvements in the efficiency of the approach in nonadiabatic regimes, making it computationally tractable to treat the quantum dynamics of complex systems for long times. We demonstrate our approach by applying it to a model of condensed phase charge transfer where our method is shown to be numerically exact in regimes where fewest-switches surface hopping and mean field approaches fail to obtain the either the correct rates or long-time populations.

preprint2015arXiv

Quantum delocalization of protons in the hydrogen bond network of an enzyme active site

Enzymes utilize protein architectures to create highly specialized structural motifs that can greatly enhance the rates of complex chemical transformations. Here we use experiments, combined with ab initio simulations that exactly include nuclear quantum effects, to show that a triad of strongly hydrogen bonded tyrosine residues within the active site of the enzyme ketosteroid isomerase (KSI) facilitates quantum proton delocalization. This delocalization dramatically stabilizes the deprotonation of an active site tyrosine residue, resulting in a very large isotope effect on its acidity. When an intermediate analog is docked, it is incorporated into the hydrogen bond network, giving rise to extended quantum proton delocalization in the active site. These results shed light on the role of nuclear quantum effects in the hydrogen bond network that stabilizes the reactive intermediate of KSI, and the behavior of protons in biological systems containing strong hydrogen bonds.

preprint2014arXiv

Quantum fluctuations and isotope effects in ab initio descriptions of water

Nuclear quantum effects, such as zero-point energy and tunneling, cause significant changes to the structure and dynamics of hydrogen bonded systems such as liquid water. However, due to the current inability to simulate liquid water using an exact description of its electronic structure, the interplay between nuclear and electronic quantum effects remains unclear. Here we use simulations that incorporate the quantum mechanical nature of both the nuclei and electrons to provide a fully ab initio determination of the particle quantum kinetic energies, free energy change upon exchanging hydrogen for deuterium and the isotope fractionation ratio in water. These properties, which selectively probe the quantum nature of the nuclear degrees of freedom, allow us to make direct comparison to recent experiments and elucidate how electronic exchange and correlation and nuclear quantum fluctuations determine the structure of the hydrogen bond in water.

preprint2013arXiv

Interface limited growth of heterogeneously nucleated ice in supercooled water

Heterogeneous ice growth exhibits a maximum in freezing rate arising from the competition between kinetics and the thermodynamic driving force between the solid and liquid states. Here, we use molecular dynamics simulations to elucidate the atomistic details of this competition, focusing on water properties in the interfacial region along the secondary prismatic direction. The crystal growth velocity is maximized when the efficiency of converting interfacial water molecules to ice, collectively known as the attachment kinetics, is greatest. We find water molecules that contact the intermediate ice layer in concave regions along the atomistically roughened surface are more likely to freeze directly. The increased roughening of the solid surface at large undercoolings consequently plays an important limiting role on the rate of ice growth, as water molecules are unable to integrate into increasingly deeper surface pockets. These results provide insights into the molecular mechanisms for self-assembly of solid phases that are important in many biological and atmospheric processes.

preprint2013arXiv

Multiple Time Step Integrators in Ab Initio Molecular Dynamics

Multiple time-scale algorithms exploit the natural separation of time-scales in chemical systems to greatly accelerate the efficiency of molecular dynamics simulations. Although the utility of these methods in systems where the interactions are described by empirical potentials is now well established, their application to ab initio molecular dynamics calculations has been limited by difficulties associated with splitting the ab initio potential into fast and slowly varying components. Here we show that such a timescale separation is possible using two different schemes: one based on fragment decomposition and the other on range separation of the Coulomb operator in the electronic Hamiltonian. We demonstrate for both water clusters and a solvated hydroxide ion that multiple time-scale molecular dynamics allows for outer time steps of 2.5 fs, which are as large as those obtained when such schemes are applied to empirical potentials, while still allowing for bonds to be broken and reformed throughout the dynamics. This permits computational speedups of up to 4.4x, compared to standard Born-Oppenheimer ab initio molecular dynamics with a 0.5 fs time step, while maintaining the same energy conservation and accuracy.

preprint2013arXiv

Unraveling quantum mechanical effects in water using isotopic fractionation

When two phases of water are at equilibrium, the ratio of hydrogen isotopes in each is slightly altered due to their different phase affinities. This isotopic fractionation process can be utilized to analyze water's movement in the world's climate. Here we show that equilibrium fractionation ratios, an entirely quantum mechanical property, also provide a sensitive probe to assess the magnitude of nuclear quantum fluctuations in water. By comparing the predictions of a series of water models, we show that those describing the OH chemical bond as rigid or harmonic greatly over-predict the magnitude of isotope fractionation. Models that account for anharmonicity in this coordinate are shown to provide much more accurate results due to their ability to give partial cancellation between inter and intra-molecular quantum effects. These results give evidence of the existence of competing quantum effects in water and allow us to identify how this cancellation varies across a wide range of temperatures. In addition, this work demonstrates that simulation can provide accurate predictions and insights into hydrogen fractionation.

preprint2012arXiv

Efficient methods and practical guidelines for simulating isotope effects

The shift in chemical equilibria due to isotope substitution is often exploited to gain insight into a wide variety of chemical and physical processes. It is a purely quantum mechanical effect, which can be computed exactly using simulations based on the path integral formalism. Here we discuss how these techniques can be made dramatically more efficient, and how they ultimately outperform quasi-harmonic approximations to treat quantum liquids not only in terms of accuracy, but also in terms of computational efficiency. To achieve this goal we introduce path integral quantum mechanics estimators based on free energy perturbation, which enable the evaluation of isotope effects using only a single path integral molecular dynamics trajectory of the naturally abundant isotope. We use as an example the calculation of the free energy change associated with H/D and 16O/18O substitutions in liquid water, and of the fractionation of those isotopes between the liquid and the vapor phase. In doing so, we demonstrate and discuss quantitatively the relative benefits of each approach, thereby providing a set of guidelines that should facilitate the choice of the most appropriate method in different, commonly encountered scenarios. The efficiency of the estimators we introduce and the analysis that we perform should in particular facilitate accurate ab initio calculation of isotope effects in condensed phase systems.

preprint2012arXiv

Growing point-to-set length scale correlates with growing relaxation times in model supercooled liquids

It has been demonstrated recently that supercooled liquids sharing simple structural features (e.g. pair distribution functions) may exhibit strikingly distinct dynamical behavior. Here we show that a more subtle structural feature correlates with relaxation times in three simulated systems that have nearly identical radial distribution functions but starkly different dynamical behavior. In particular, for the first time we determine the thermodynamic "point-to-set" length scale in several canonical model systems and demonstrate the quantitative connection between this length scale and the growth of relaxation times. Our results provide clues necessary for distinguishing competing theories of the glass transition.

preprint2011arXiv

Reduced density matrix hybrid approach: An efficient and accurate method for adiabatic and non-adiabatic quantum dynamics

We present a new approach to calculate real-time quantum dynamics in complex systems. The formalism is based on the partitioning of a system's environment into "core" and "reservoir" modes, with the former to be treated quantum mechanically and the latter classically. The presented method only requires the calculation of the system's reduced density matrix averaged over the quantum core degrees of freedom which is then coupled to a classically evolved reservoir to treat the remaining modes. We demonstrate our approach by applying it to the spin-boson problem using the noninteracting blip approximation to treat the system and core, and Ehrenfest dynamics to treat the reservoir. The resulting hybrid methodology is accurate for both fast and slow baths, since it naturally reduces to its composite methods in their respective regimes of validity. In addition, our combined method is shown to yield good results in intermediate regimes where neither approximation alone is accurate and to perform equally well for both strong and weak system-bath coupling. Our approach therefore provides an accurate and efficient methodology for calculating quantum dynamics in complex systems.

preprint2011arXiv

Reduced density matrix hybrid approach: Application to electronic energy transfer

Electronic energy transfer in the condensed phase, such as that occurring in photosynthetic complexes, frequently occurs in regimes where the energy scales of the system and environment are similar. This situation provides a challenge to theoretical investigation since most approaches are accurate only when a certain energetic parameter is small compared to others in the problem. Here we show that in these difficult regimes, the Ehrenfest approach provides a good starting point for a dynamical description of the energy transfer process due to its ability to accurately treat coupling to slow environmental modes. To further improve on the accuracy of the Ehrenfest approach, we use our reduced density matrix hybrid framework to treat the faster environmental modes quantum mechanically, at the level of a perturbative master equation. This combined approach is shown to provide an efficient and quantitative description of electronic energy transfer in a model dimer and the Fenna-Matthews-Olson complex and is used to investigate the effect of environmental preparation on the resulting dynamics.

preprint2011arXiv

Theory and simulations of quantum glass forming liquids

A comprehensive microscopic dynamical theory is presented for the description of quantum fluids as they transform into glasses. The theory is based on a quantum extension of mode-coupling theory. Novel effects are predicted, such as reentrant behavior of dynamical relaxation times. These predictions are supported by path integral ring polymer molecular dynamics simulations. The simulations provide detailed insight into the factors that govern slow dynamics in glassy quantum fluids. Connection to other recent work on both quantum glasses as well as quantum optimization problems is presented.

preprint2010arXiv

Competing quantum effects in the dynamics of a flexible water model

Numerous studies have identified large quantum mechanical effects in the dynamics of liquid water. In this paper, we suggest that these effects may have been overestimated due to the use of rigid water models and flexible models in which the intramolecular interactions were described using simple harmonic functions. To demonstrate this, we introduce a new simple point charge model for liquid water, q-TIP4P/F, in which the O--H stretches are described by Morse-type functions. We have parameterized this model to give the correct liquid structure, diffusion coefficient, and infra-red absorption frequencies in quantum (path integral-based) simulations. By comparing classical and quantum simulations of the liquid, we find that quantum mechanical fluctuations increase the rates of translational diffusion and orientational relaxation in our model by a factor of around 1.15. This effect is much smaller than that observed in all previous simulations of simple empirical water models, which have found a quantum effect of at least 1.4 regardless of the quantum simulation method or the water model employed. The small quantum effect in our model is a result of two competing phenomena. Intermolecular zero point energy and tunneling effects destabilize the hydrogen bonding network, leading to a less viscous liquid with a larger diffusion coefficient. However this is offset by intramolecular zero point motion, which changes the average water monomer geometry resulting in a larger dipole moment, stronger intermolecular interactions, and slower diffusion. We end by suggesting, on the basis of simulations of other potential energy models, that the small quantum effect we find in the diffusion coefficient is associated with the ability of our model to produce a single broad O-H stretching band in the infra-red absorption spectrum.

preprint2010arXiv

Efficient multiple time scale molecular dynamics: using colored noise thermostats to stabilize resonances

Multiple time scale molecular dynamics enhances computational efficiency by updating slow motions less frequently than fast motions. However, in practice the largest outer time step possible is limited not by the physical forces but by resonances between the fast and slow modes. In this paper we show that this problem can be alleviated by using a simple colored noise thermostatting scheme which selectively targets the high frequency modes in the system. For two sample problems, flexible water and solvated alanine dipeptide, we demonstrate that this allows the use of large outer time steps while still obtaining accurate sampling and minimizing the perturbation of the dynamics. Furthermore, this approach is shown to be comparable to constraining fast motions, thus providing an alternative to molecular dynamics with constraints.

preprint2010arXiv

Efficient stochastic thermostatting of path integral molecular dynamics

The path integral molecular dynamics (PIMD) method provides a convenient way to compute the quantum mechanical structural and thermodynamic properties of condensed phase systems at the expense of introducing an additional set of high-frequency normal modes on top of the physical vibrations of the system. Efficiently sampling such a wide range of frequencies provides a considerable thermostatting challenge. Here we introduce a simple stochastic path integral Langevin equation (PILE) thermostat which exploits an analytic knowledge of the free path integral normal mode frequencies. We also apply a recently-developed colored-noise thermostat based on a generalized Langevin equation (GLE), which automatically achieves a similar, frequency-optimized sampling. The sampling efficiencies of these thermostats are compared with that of the more conventional Nosé-Hoover chain (NHC) thermostat for a number of physically relevant properties of the liquid water and hydrogen-in-palladium systems. In nearly every case, the new PILE thermostat is found to perform just as well as the NHC thermostat while allowing for a computationally more efficient implementation. The GLE thermostat also proves to be very robust delivering a near-optimum sampling efficiency in all of the cases considered. We suspect that these simple stochastic thermostats will therefore find useful application in many future PIMD simulations.

preprint2010arXiv

Quantum fluctuations can promote or inhibit glass formation

The very nature of glass is somewhat mysterious: while relaxation times in glasses are of sufficient magnitude that large-scale motion on the atomic level is essentially as slow as it is in the crystalline state, the structure of glass appears barely different than that of the liquid that produced it. Quantum mechanical systems ranging from electron liquids to superfluid helium appear to form glasses, but as yet no unifying framework exists connecting classical and quantum regimes of vitrification. Here we develop new insights from theory and simulation into the quantum glass transition that surprisingly reveal distinct regions where quantum fluctuations can either promote or inhibit glass formation.