Source author record

David Roundy

David Roundy 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

9works
5topics
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

9 published item(s)

preprint2020arXiv

Flat histogram method comparison on 2D Ising Model

We compare the convergence of several flat-histogram methods applied to the 2D Ising model, including the recently introduced stochastic approximation with a dynamic update factor (SAD) method. We compare this method with the Wang-Landau (WL) method, the $1/t$ variant of the WL method, and standard stochastic approximation Monte Carlo (SAMC). In addition, we consider a procedure WL followed by a "production run" with fixed weights that refines the estimation of the entropy. To our knowledge, this work is the first to test this approach against other methods. We find that WL followed by a production run \emph{does} converge to the true density of states, in contrast to pure WL. Three of the methods converge robustly: SAD, $1/t$-WL, and WL followed by a production run. Of these, SAD does not require \emph{a priori} knowledge of the energy range. This work also shows that WL followed by a production run performs superior to other forms of WL while ensuring both ergodicity and detailed balance.

preprint2019arXiv

Stochastic Approximation Monte Carlo with a Dynamic Update Factor

We present a new Monte Carlo algorithm based on the Stochastic Approximation Monte Carlo (SAMC) algorithm for directly calculating the density of states. The proposed method is Stochastic Approximation with a Dynamic update factor (SAD) which dynamically adjusts the update factor $γ_t$ during the course of the simulation. We test this method on the square-well fluid and the 31-atom Lennard-Jones cluster and compare the convergence behavior of several related Monte Carlo methods. We find that both the SAD and $1/t$-Wang-Landau ($1/t$-WL) methods rapidly converge to the correct density of states without the need for the user to specify an arbitrary tunable parameter $t_0$ as in the case of SAMC. SAD requires as input the temperature range of interest, in contrast to $1/t$-WL, which requires that the user identify the interesting range of energies. The convergence of the $1/t$-WL method is very sensitive to the energy range chosen for the low-temperature heat capacity of the Lennard-Jones cluster. Thus, SAD is more powerful in the common case in which the range of energies is not known in advance.

preprint2014arXiv

An efficient approach to approximating the pair distribution function of the inhomogeneous hard-sphere fluid

We introduce an approximation for the pair distribution function of the inhomogeneous hard sphere fluid. Our approximation makes use of our recently published averaged pair distribution function at contact which has been shown to accurately reproduce the averaged pair distribution function at contact for inhomogeneous density distributions. This approach achieves greater computational efficiency than previous approaches by enabling the use of exclusively fixed-kernel convolutions and thus allowing an implementation using fast Fourier transforms. We compare results for our pair distribution approximation with two previously published works and Monte-Carlo simulation, showing favorable results.

preprint2014arXiv

Experts' understanding of partial derivatives using the Partial Derivative Machine

Partial derivatives are used in a variety of different ways within physics. Most notably, thermodynamics uses partial derivatives in ways that students often find confusing. As part of a collaboration with mathematics faculty, we are at the beginning of a study of the teaching of partial derivatives, a goal of better aligning the teaching of multivariable calculus with the needs of students in STEM disciplines. As a part of this project, we have performed a pilot study of expert understanding of partial derivatives across three disciplines: physics, engineering and mathematics. Our interviews made use of the Partial Derivative Machine (PDM), which is a mechanical system featuring four observable and controllable properties, of which any two are independent. Using the PDM, we probed expert understanding of partial derivatives in an experimental context in which there is not a known functional form. Through these three interviews, we found that the mathematicians exhibited a striking difference in their understanding of derivatives relative to the other groups. The physicists and engineers were quick to use measurements to find a numeric approximation for a derivative. In contrast, the mathematicians repeatedly returned to speculation as to the functional form, and although they were comfortable drawing qualitative conclusions about the system from measurements, were reluctant to approximate the derivative through measurement. This pilot study led us to further questions. How do fields differ in their experts' concept image of partial derivatives? What representations of partial derivatives are preferred by experts? We plan to address these questions by means of further interviews with a wider range of disciplinary experts.

preprint2013arXiv

A Modified Histogram Method for Disordered Lattices

A new method is developed to extend the histogram method to lattices with any type of disorder. The Monte Carlo single- and multiple-histogram methods were developed to get the most out of only a few simulations, but are restricted to simulations with identical lattice configurations. The method introduced here expands on the histogram method to allow data from various disordered lattice configurations to be optimally combined in a single weighted average result. This method is applied to a simple Ising-like model of the relaxor ferroelectric perovskite solid solution $BaTiO_3$ - $Bi(Zn_{1/2}Ti_{1/2})O_3$ (BT-BZT), in which disorder plays a pivotal role.

preprint2013arXiv

Improved association in a classical density functional theory for water

We present a modification to our recently published SAFT-based classical density functional theory for water. We have recently developed and tested a functional for the averaged radial distribution function at contact of the hard-sphere fluid that is dramatically more accurate at interfaces than earlier approximations. We now incorporate this improved functional into the association term of our free energy functional for water, improving its description of hydrogen bonding. We examine the effect of this improvement by studying two hard solutes: a hard hydrophobic rod and a hard sphere. The improved functional leads to a moderate change in the density profile and a large decrease in the number of hydrogen bonds broken in the vicinity of the solutes.We present a modification to our recently published SAFT-based classical density functional theory for water. We have recently developed and tested a functional for the averaged radial distribution function at contact of the hard-sphere fluid that is dramatically more accurate at interfaces than earlier approximations. We now incorporate this improved functional into the association term of our free energy functional for water, improving its description of hydrogen bonding. We examine the effect of this improvement by studying two hard solutes: a hard hydrophobic rod and a hard sphere. The improved functional leads to a moderate change in the density profile and a large decrease in the number of hydrogen bonds broken in the vicinity of the solutes.

preprint2012arXiv

A Classical Density-Functional Theory for Describing Water Interfaces

We develop a classical density functional for water which combines the White Bear fundamental-measure theory (FMT) functional for the hard sphere fluid with attractive interactions based on the Statistical Associating Fluid Theory (SAFT-VR). This functional reproduces the properties of water at both long and short length scales over a wide range of temperatures, and is computationally efficient, comparable to the cost of FMT itself. We demonstrate our functional by applying it to systems composed of two hard rods, four hard rods arranged in a square and hard spheres in water.

preprint2012arXiv

Using Fundamental Measure Theory to Treat the Correlation Function of the Inhomogeneous Hard-Sphere Fluid

We investigate the value of the correlation function of an inhomogeneous hard-sphere fluid at contact. This quantity plays a critical role in Statistical Associating Fluid Theory (SAFT), which is the basis of a number of recently developed classical density functionals. We define two averaged values for the correlation function at contact, and derive formulas for each of them from the White Bear version of the Fundamental Measure Theory functional, using an assumption of thermodynamic consistency. We test these formulas, as well as two existing formulas against Monte Carlo simulations, and find excellent agreement between the Monte Carlo data and one of our averaged correlation functions.

preprint2010arXiv

Classical density-functional theory for water

We introduce a new computationally efficient and accurate classical density-functional theory for water and apply it to hydration of hard spheres and inert gas atoms. We find good agreement with molecular dynamics simulations for the hydration of hard spheres and promising agreement for the solvation of inert gas atoms in water. Finally, we explore the importance of the orientational ambiguity in state-of-the-art continuum theories of water, which are based on the molecular density only.