Researcher profile

Robert A. DiStasio Jr.

Robert A. DiStasio Jr. contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

8 published item(s)

preprint2020arXiv

QM7-X: A comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules

We introduce QM7-X, a comprehensive dataset of 42 physicochemical properties for $\approx$ 4.2 M equilibrium and non-equilibrium structures of small organic molecules with up to seven non-hydrogen (C, N, O, S, Cl) atoms. To span this fundamentally important region of chemical compound space (CCS), QM7-X includes an exhaustive sampling of (meta-)stable equilibrium structures - comprised of constitutional/structural isomers and stereoisomers, e.g., enantiomers and diastereomers (including cis-/trans- and conformational isomers) - as well as 100 non-equilibrium structural variations thereof to reach a total of $\approx$ 4.2 M molecular structures. Computed at the tightly converged quantum-mechanical PBE0+MBD level of theory, QM7-X contains global (molecular) and local (atom-in-a-molecule) properties ranging from ground state quantities (such as atomization energies and dipole moments) to response quantities (such as polarizability tensors and dispersion coefficients). By providing a systematic, extensive, and tightly-converged dataset of quantum-mechanically computed physicochemical properties, we expect that QM7-X will play a critical role in the development of next-generation machine-learning based models for exploring greater swaths of CCS and performing in silico design of molecules with targeted properties.

preprint2015arXiv

A Systematic Study of Chloride Ion Solvation in Water using van der Waals Inclusive Hybrid Density Functional Theory

In this work, the solvation and electronic structure of the aqueous chloride ion solution was investigated using Density Functional Theory (DFT) based \textit{ab initio} molecular dynamics (AIMD). From an analysis of radial distribution functions, coordination numbers, and solvation structures, we found that exact exchange ($E_{\rm xx}$) and non-local van der Waals (vdW) interactions effectively \textit{weaken} the interactions between the Cl$^-$ ion and the first solvation shell. With a Cl-O coordination number in excellent agreement with experiment, we found that most configurations generated with vdW-inclusive hybrid DFT exhibit 6-fold coordinated distorted trigonal prism structures, which is indicative of a significantly disordered first solvation shell. By performing a series of band structure calculations on configurations generated from AIMD simulations with varying DFT potentials, we found that the solvated ion orbital energy levels (unlike the band structure of liquid water) strongly depend on the underlying molecular structures. In addition, these orbital energy levels were also significantly affected by the DFT functional employed for the electronic structure; as the fraction of $E_{\rm xx}$ was increased, the gap between the highest occupied molecular orbital of Cl$^-$ and the valence band maximum of liquid water steadily increased towards the experimental value.

preprint2015arXiv

Analytical nuclear gradients for the range-separated many-body dispersion model of noncovalent interactions

Accurate treatment of the long-range electron correlation energy, including van der Waals (vdW) or dispersion interactions, is essential for describing the structure, dynamics, and function of a wide variety of systems. Among the most accurate models for including dispersion into density functional theory (DFT) is the range-separated many-body dispersion (MBD) method [A. Ambrossetti et al., J. Chem. Phys. 140, 18A508 (2014)], in which the correlation energy is modeled at short-range by a semi-local density functional and at long-range by a model system of coupled quantum harmonic oscillators. In this work, we develop analytical gradients of the MBD energy with respect to nuclear coordinates, including all implicit coordinate dependencies arising from the partitioning of the charge density into Hirshfeld effective volumes. To demonstrate the efficiency and accuracy of these MBD gradients for geometry optimizations of systems with intermolecular and intramolecular interactions, we optimized conformers of the benzene dimer and isolated small peptides with aromatic side-chains. We find excellent agreement with the wavefunction theory reference geometries of these systems (at a fraction of the computational cost) and find that MBD consistently outperforms the popular TS and D3(BJ) dispersion corrections. To demonstrate the performance of the MBD model on a larger system with supramolecular interactions, we optimized the C$_{60}$@C$_{60}$H$_{28}$ buckyball catcher host-guest complex. Finally, we find that neglecting the implicit nuclear coordinate dependence arising from the charge density partitioning, as has been done in prior numerical treatments, leads to an unacceptable error in the MBD forces, with relative errors of ~20% (on average) that can extend well beyond 100%.

preprint2015arXiv

Inverse design of disordered stealthy hyperuniform spin chains

Positioned between crystalline solids and liquids, disordered many-particle systems which are stealthy and hyperuniform represent new states of matter that are endowed with novel physical and thermodynamic properties. Such stealthy and hyperuniform states are unique in that they are transparent to radiation for a range of wavenumbers around the origin. In this work, we employ recently developed inverse statistical-mechanical methods, which seek to obtain the optimal set of interactions that will spontaneously produce a targeted structure or configuration as a unique ground state, to investigate the spin-spin interaction potentials required to stabilize disordered stealthy hyperuniform one-dimensional (1D) Ising-like spin chains. By performing an exhaustive search over the spin configurations that can be enumerated on periodic 1D integer lattices containing $N=2,3,\ldots,36$ sites, we were able to identify and structurally characterize \textit{all} stealthy hyperuniform spin chains in this range of system sizes. Within this pool of stealthy hyperuniform spin configurations, we then utilized such inverse optimization techniques to demonstrate that stealthy hyperuniform spin chains can be realized as either unique or degenerate disordered ground states of radial long-ranged (relative to the spin chain length) spin-spin interactions. Such exotic ground states are distinctly different from spin glasses in both their inherent structural properties and the nature of the spin-spin interactions required to stabilize them. As such, the implications and significance of the existence of such disordered stealthy hyperuniform ground state spin systems warrants further study, including whether their bulk physical properties and excited states, like their many-particle system counterparts, are singularly remarkable, and can be experimentally realized.

preprint2015arXiv

Local Structure Analysis in $Ab$ $Initio$ Liquid Water

Within the framework of density functional theory, the inclusion of exact exchange and non-local van der Waals/dispersion (vdW) interactions is crucial for predicting a microscopic structure of ambient liquid water that quantitatively agrees with experiment. In this work, we have used the local structure index (LSI) order parameter to analyze the local structure in such highly accurate $ab$ $initio$ liquid water. At ambient conditions, the LSI probability distribution, P($I$), was unimodal with most water molecules characterized by more disordered high-density-like local environments. With thermal excitations removed, the resultant bimodal P($I$) in the inherent potential energy surface (IPES) exhibited a 3:1 ratio between high- and low-density-like molecules, with the latter forming small connected clusters amid the predominant population. By considering the spatial correlations and hydrogen bond network topologies $among$ water molecules with the same LSI identities, we demonstrate that the signatures of the experimentally observed low- (LDA) and high-density (HDA) amorphous phases of ice are present in the IPES of ambient liquid water. Analysis of the LSI autocorrelation function uncovered a persistence time of $\sim$ 4 ps---a finding consistent with the fact that natural thermal fluctuations are responsible for transitions between these distinct yet transient local aqueous environments in ambient liquid water.

preprint2014arXiv

The Individual and Collective Effects of Exact Exchange and Dispersion Interactions on the Ab Initio Structure of Liquid Water

In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local vdW/dispersion interactions, via a fully self-consistent density-dependent dispersion correction, and approximate nuclear quantum effects (aNQE), via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx, vdW, and aNQE as resulting from a large-scale AIMD simulation of (H$_2$O)$_{128}$ at the PBE0+vdW level of theory, significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, $S_{\rm OO}(Q)$, and corresponding oxygen-oxygen radial distribution function, $g_{\rm OO}(r)$, that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment as demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of the PBE0 hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions.

preprint2013arXiv

Long-range correlation energy calculated from coupled atomic response functions

An accurate determination of the electron correlation energy is essential for describing the structure, stability, and function in a wide variety of systems, ranging from gas-phase molecular assemblies to condensed matter and organic/inorganic interfaces. Even small errors in the correlation energy can have a large impact on the description of chemical and physical properties in the systems of interest. In this context, the development of efficient approaches for the accurate calculation of the long-range correlation energy (and hence dispersion) is the main challenge. In the last years a number of methods have been developed to augment density functional approximations via dispersion energy corrections, but most of these approaches ignore the intrinsic many-body nature of correlation effects, leading to inconsistent and sometimes even qualitatively incorrect predictions. Here we build upon the recent many-body dispersion (MBD) framework, which is intimately linked to the random-phase approximation for the correlation energy. We separate the correlation energy into short-range contributions that are modeled by semi-local functionals and long-range contributions that are calculated by mapping the complex all-electron problem onto a set of atomic response functions coupled in the dipole approximation. We propose an effective range-separation of the coupling between the atomic response functions that extends the already broad applicability of the MBD method to non-metallic materials with highly anisotropic responses, such as layered nanostructures. Application to a variety of high-quality benchmark datasets illustrates the accuracy and applicability of the improved MBD approach, which offers the prospect of first-principles modeling of large structurally complex systems with an accurate description of the long-range correlation energy.

preprint2012arXiv

Many-Body Dispersion Interactions in Molecular Crystal Polymorphism

Polymorphs in molecular crystals are often very close in energy, yet they may possess markedly different physical and chemical properties. The understanding and prediction of polymorphism is of paramount importance for a variety of applications, including pharmaceuticals, non-linear optics, and hydrogen storage. Here, we show that the non-additive many-body dispersion (MBD) energy beyond the standard pairwise approximation is crucial for the correct qualitative and quantitative description of polymorphism in molecular crystals. This is rationalized by the sensitive dependence of the MBD energy on the polymorph geometry and the ensuing dynamic electric fields inside molecular crystals. We use the glycine crystal as a fundamental and stringent benchmark case to demonstrate the accuracy of the DFT+MBD method.