# Publications

David M. Rogers, "Dual Characterization of the Ornstein-Zernike Equation in Moment Space." submitted, 2019.
I re-write Ornstein-Zernike theory in a convenient matrix basis so that future computational implementations can be made robust and that numerical and fitting error can be tightly controlled. It is hoped that the framework can be used in the near future to compute molecular solvation free energies and to realize its enormous potential economic and environmental benefits for formulating industrial fluids and consumer products. |

Phillip S. Hudson, Stefan Boresch, David M. Rogers, and H. Lee Woodcock., "Accelerating QM/MM Free Energy Computations via Intramolecular Force Matching" J. Chem. Theory Comput., 14 (12):6327–35, 2018.
My co-authors apply my algorithm and software implementation of Bayesian generalized linear model regression with linear inequality constraints to estimate molecular modeling parameters from quantum their mechanical calculations. They find the method gives robust results that greatly increase efficiency of additional simulations. |

David M. Rogers, "Extension of Kirkwood-Buff theory to the canonical ensemble." J. Chem. Phys., 148:054102, 2018.
We present a way to utilize Essmann et. al.'s smooth particle mesh Ewald algorithm (implemented here) to get the low-angle direct correlation function from canonical ensemble simulation data. The result avoids well-documented issues with long wavelength (small-angle scattering) modes that appear when using real-space methods. This result is combined with new estimates of finite-size effects and grand-canonical ensemble corrections to provide definitive results on the problem of extrapolating Kirkwood-Buff integrals. |

Juan M. Vanegas, Frank Heinrich, David M.Rogers, Bryan D. Carson, Sadie La Bauve, Briana C. Vernon, Bulent Akgun, Sushil Satija, Aihua Zheng, Margaret Kielian, Susan B. Rempe, and Michael S. Kent, "Insertion of Dengue E into lipid bilayers studied by neutron reflectivity and molecular dynamics simulations." BBA 1860(5):1216-1230, 2018.
We compare neutron reflectivity experiments to molecular dynamics calculations on the orientation and binding position of the Dengue envelope protein responsible for viral escape from the host cell's endosome. In addition to the hydrophobic fusion tip, we find important roles for three positively-charged residues in the viral protein that contribute to host membrane binding. These could potentially be targeted by new anti-viral medicines. |

David M. Rogers, "An information theory model for dissipation in open quantum systems." J. Phys., Conference Series 880(1):012039, 2017.
This paper presents a new, simple ansatz for adding dissipation to arbitrary stochastic forcing of a quantum dynamical system. For Gaussian random forces, it predicts a Lindblad equation identical to the Caldeira-Leggett model up to order <math>\beta^2</math>, where the same term is present with a different prefactor. The system-centric, phase space picture here shows that the <math>\beta^2</math> term represents a quantum confinement effect. |

David M. Rogers, "Einstein-Podolsky-Rosen paradox implies a minimum achievable temperature." Phys. Rev. E 95, 012149, 2017.
This paper provides measurement-based definitions of heat and work that can be realized in current laboratory setups. The first and second laws are proved despite the fact that temperature is treated completely as as a property of the interacting reservoir. Measurements of the work are subject to the famous EPR paradox because the work exchanged between two quantum systems is not defined until a measurement is performed. Based on this, we show that even an environment at absolute zero cannot lower a system's temperature below a minimum characteristic of the way the environment is coupled to the system. |

Guy W. Dayhoff II and David M. Rogers, "Driving forces in MD simulations of transition and ‘Free’ flows." Mol. Sim. 43(5-6), pp. 467-477, 2017. (special issue on Surface Chemistry)
We set out to test the Joule-Thomson analysis of thermodynamics of porous flow for gasses through a nanopore and found that while local equilibrium is established in the steady-state, finite-size effects cause heat flow opposite the flow direction that violates the assumption of an adiabatic porous plug. |

David M. Rogers, "Efficient Primitives for Standard Tensor Linear Algebra." Proc. XSEDE16 Conference on Diversity, Big Data, and Science at Scale, no. 14, 2016.
This paper introduces 3 basic functions that generalize BLAS to tensors and presents a code generation strategy for their efficient execution on GPUs that achieves peak performance on the same order of magnitude as for traditional, vendor-optimized matrix-multiplications. |

David M. Rogers, "Overcoming the Minimum Image Constraint Using the Closest Point Search." J. Mol. Graph. Model 68, pp. 197–205, 2016.
An elegant solution to the problem of finding periodic images in non-rectangular lattices is provided based on the closest vector problem. Related code is here: [1] |

Elisa La Bauve, Briana C. Vernon, Dongmei Ye, David M. Rogers, Cathryn M. Siegrist, Bryan Carson. Susan L. Rempe, Aihua Zheng, Margaret C. Kielian, Andrew P. Shreve, and Michael S. Kent. "Method for measuring the unbinding energy of strongly-bound membrane-associated proteins." BBA Biomembranes 1858(11): 2753–62, 2016.
This paper gives multiple experimental measurements of binding energy between the Dengue virus envelope protein and host membranes that largely confirm our computational predictions from 2015. I contributed all the theory for terminal velocity during sedimentation, along with a novel kinetic analysis providing the free energy and enthalpy of the dissociation barrier (all the details are at the end of the appendix). |

David M. Rogers. "Thermodynamics of Maximum Transition Entropy for Quantum Assemblies." arXiv:1503.01232 submitted, 2016.
The work presents a new, general, theoretical foundation for the dynamics of open quantum systems modeled on the maximum entropy derivation of equilibrium statistical mechanics. Computational results are presented for three detailed systems to validate and reinforce the theory. It represents a significant advancement for the field, as it lucidly connects the dynamics of a single wavefunction plus environmental noise to the Caldeira-Leggett model for density matrices. |

Andriy Anishkin, Juan M. Vanegas, David M. Rogers, Philip L. Lorenzi, Wai Kin Chan, Preeti Purwaha, John N. Weinstein, Sergei Sukharev, and Susan B. Rempe. "Catalytic Role of the Substrate Defines Specificity of Therapeutic L-Asparaginase."
J. Mol. Biol. 427:2867-2885, 2015. We present an explanation for the (until now controversial) catalytic mechanism of type 2 bacterial L-asparaginase enzymes. By using the carboxylic acid of the substrate (asparagine) as the proton acceptor, this enzyme is able to preferentially carry out deamidation on asparagine more quickly than for the competing substrate, glutamine. The hypothesis, re-discovered from our MD simulations, was first put forward years ago in contested experimental studies, and now finds additional support from our MD and QM calculations. |

David M. Rogers. "Towards a Direct, By-Need Evaluator for Dependently Typed Languages."
arXiv:1509.07036 submitted, 2015. This paper describes the implementation of a new interpreted language for distributed parallel computing. It achieves its goal by maintaining pure functional semantics, allowing all terms in the language to be partially evaluated and serialized to network storage at any point during computation. |

Marielle Soniat, David M. Rogers, and Susan Rempe. "Dispersion- and Exchange-Corrected Density Functional Theory for Sodium Ion Hydration." J. Chem. Theory. Comput. 142:074101, 2015.
We studied the influence of dispersion energy corrections on the free energy of formation for sodium-water clusters computed with DFT and wound up discovering that dispersion and split-range exchange functionals can somewhat counter-balance each other. The charged sodium ion pulls on the water's electrons, clearly showing which density functionals over-polarize compared to CCSD. Split-range exchange can reduce this over-polarization, but results in reduced electrostatic interaction. Dispersion can lower the binding energy again to counter-balance. So, functionals fit to experimental formation energies need both effects to avoid lowering energies by over-polarizing. |

David M. Rogers. "Real-space quadrature: a convenient, efficient representation for multipole expansions." J. Chem. Phys. 142:074101, 2015. (Presentation)
I introduce sets of point charges that are able to simultaneously reproduce all multipole (spherical harmonic) expansions up to arbitrary order. The number of points is space-optimal. Translations are described from the usual harmonics and from Cartesian moments (dipole, quadrupole, etc.) on supersymmetric tensors to directional moments using the point weight distribution. Many applications are possible, including trivial implementation of multipoles in molecular mechanics and representing probability distributions over rotation space. |

David M. Rogers, Michael S. Kent, and Susan B. Rempe, "Molecular basis of endosomal-membrane association for the dengue virus envelope protein." BBA Biomembranes 1848(4):1041-52, 2015.
A fully atomistic potential of mean force for association of the viral envelope protein from Dengue virus was compared to a Poisson-Boltzmann electrostatic plus dispersion model. The results are consistent, showing hope for this type of combined scale simulation. |

Yaqin Fu, Binsong Li, Ying-Bing Jiang, Darren R. Dunphy, Andy Tsai, Siu-Yue Tam, Hongyou Fan, Hongxia Zhan, David Rogers, Susan Rempe, Plamen Atanassov, Joseph L. Cecchi, and C. Jeffrey Brinker "Atomic Layer Deposition of L-Alanine Polypeptide." JACS 136(45):15821–4, 2014.
This paper with our experimental collaborators carried out blocked peptide synthesis by vapor-depositing Boc-L-alanine to create a uniform thin film of polypeptides grown on a silica substrate activated by aminopropyltrimethoxysilane. |

David M. Rogers. "Silmaril, A Functional Language for Distributed Parallel Evaluation." Submitted version |

Mathias B. Andersen, David M. Rogers, Junyu Mai, Benjamin Schudel, Anson V. Hatch, Susan B. Rempe and Ali Mani. "Spatiotemporal pH dynamics in concentration polarization near ion-selective membranes." Langmuir, 30(26):7902–7912, 2014 |

W. K. Chan, P. L. Lorenzi, A. Anishkin, P. Purwaha, D. M. Rogers, S. Sukharev, S.B. Rempe, and J. N. Weinstein. "The glutaminase activity of l-asparaginase is not required for anticancer activity against ASNS-negative cells." Blood. 123(23):3596-606, 2014. |

David M. Rogers, Dian Jiao, Lawrence Pratt, and Susan B. Rempe. "Structural Models and Molecular Thermodynamics of Hydration of Ions and Small Molecules" Annu. Rep. Comp. Chem. 8:71–127, 2012. |

David M. Rogers and Susan B. Rempe. "Irreversible Thermodynamics." J. Phys.: Conf. Ser. 402:012014, 2012. |

David M. Rogers, Thomas L. Beck, and Susan B. Rempe. "An Information Theory Approach to Nonlinear, Nonequilibrium Thermodynamics." J. Stat. Phys. 145(2):385-409, 2011
We show how the interpretation of thermodynamic states as representing system information leads naturally to thermodynamic cycles and the first and second laws of thermodynamics as well as similar formulations for nontrivial nonequilibrium problems. The logical development of the theory also leads naturally to correct indistinguishability factors in the partition function. |

Sameer Varma, David M. Rogers, Lawrence R. Pratt, and Susan B. Rempe. "Perspectives on Ion Selectivity: Design Principles for K+ Selectivity in Membrane Transport." J. Gen. Physiol., 137(6):479-488, 2011.
We review the development of models for understanding the physical basis of selectivity for K+ ions over Na+, its sibling only one row behind, in membrane channels and transporters. Although the problem is subtle because of the morass of competing effects, we emphasize work analyzing the systematic influence of the environment on tipping local binding site structure toward selective configurations. |

David M. Rogers and Susan B. Rempe. “Probing the Thermodynamics of Competitive Ion Binding Using Minimum Energy Structures.” J. Phys. Chem. B, 115(29):9116-29, 2011.
We presented an extension of the Quasi-Chemical theory for quantifying the impact of local structure on ion complexation thermodynamics. The theory can be simply represented using a set of thermodynamic cycles involving binding site structural and compositional states as reaction intermediates. |

Susan B. Rempe and David M. Rogers; et. al. “Computational and experimental platform for understanding and optimizing water flux and salt rejection in nanoporous membranes.” Sandia Technical Report, SAND2010-6735, 2010.
We summarize work on designing polymer coatings for salt exclusion in water transporting nanopores. In this work, I collected available molecular dynamics results for these systems and performed a novel energy efficiency analysis able to relate atomistic and experimental scales as well as identify important design goals and chemical principles for material performance. |

David M. Rogers, 'Using Bayes' Theorem for Free Energy Calculations', 2009.
We investigated the central quantity of free energies in a Bayesian context and provide estimators for solvation free energies as well as optimal potential of mean force approximations to model polymer coarse-grained dynamics from atomistic simulations. |

Zhen Zhao, David M. Rogers and Thomas L. Beck. "Polarization and Charge Transfer in the Hydration of Chloride Ions." J. Chem. Phys., 132:014502, 2010.
Dr. Zhao's ab-initio analysis of the charge distribution in water-ion clusters highlighted the importance of many-body water-water interactions and charge transfer effects in determining cluster structural and energetic properties. These are still challenging to represent in modern polarizable forcefields and have implications for anion properties at interfaces. |

David M. Rogers and Thomas L. Beck. "Quasi-Chemical and Structural Analysis of Polarizable Anion Hydration." J. Chem. Phys., 132:014505, 2010.
The role of polarizability in forcefield-based models of ions and water was examined. Utilizing some of our recent developments on quasi-chemical theory, we have been able to quantify the tightened, asymmetric nature of the ion's local solvation waters induced by increased polarizability as well as the exact effects of polarization on the solvation free energy. The results suggest some potential problems and diagnostics for such models. |

David M. Rogers and Thomas L. Beck. Force Solve (Sourceforge, Chicago IL, 2008).
This force matching software implements and tests coarse-graining for general molecular systems in a mere 4000 lines of code. It is able to parametrize coarse Hamiltonians from atomic trajectory data given arbitrary definitions of coarse united-atom type models as well as carry out short Langevin Dynamics simulations on the coarse scale. The program's main drawbacks are its slow speed and high memory usage due to its simplistic design, attributable to the interpreted nature of python. |

David M. Rogers and Thomas L. Beck. "Resolution and Scale Independent Nonparametric Function Matching Using a String Energy Penalized Spline Prior." 2008. arXiv:1003.4741v1 (stat.ML).
Fresh insight is provided into long-standing mathematical issues surrounding computational modeling of continuous functions from a few sampled data points. The present research lays the groundwork for predicting the behavior of complicated many-body systems using advanced regression techniques. |

David M. Rogers and Thomas L. Beck. "Modeling molecular and ionic absolute solvation free energies with quasichemical theory bounds." J. Chem. Phys., 129:134505, 2008.
We develop a Bayesian method for computing (with error bars) the free energy for forming a nano-bubble in an arbitrary solvent system. This forms the first step of a thermodynamic cycle for dissolving a real solute. We prove that upper and lower bounds for that solvation free energy can be obtained from two simulations (with and without the solute present). The method is excellent for dissolving gas in water, while the upper/lower bounds are larger for dissolving water or ions. |