|
Specific application areas targeted by our group are as follows.
Understanding free-energy calculationsThe free-energy perturbation (FEP) method is a core technique in computational chemistry. Despite its importance and wide use, it is frequently applied in ways that are inefficient at best, and incorrect at worst. FEP calculations can be performed in either of two directions, depending on which system is taken as the reference. It is well known that these two calculations lead to results that differ systematically. One mistake commonly made is to assume that both calculations are equally wrong, and the best result is obtained by splitting the difference. By applying modeling concepts to the simulation process itself, we have shown that this heuristic is incorrect, and we have formulated one that is much more reliable. Staged FEP calculations are usually optimized by equalizing the free-energy difference between stages. We have found that this heuristic too is incorrect, and can lead to imprecisions that are several orders of magnitude worse than an approach based on consideration of entropy differences.Efficient evaluation of free-energies and phase equilibria in solidsCalculating the free energy of solid phases is more difficult than in fluid phases. Fluid free-energies are simpler because they can be evaluated by a process of inserting a single molecule. Because this is not such a big change (it describes a thermodynamically intensive property, the chemical potential), free-energy methods can handle it without great difficulty. In solids, molecule insertion cannot be completed without introducing a defect in the lattice, and this complicates the use of such methods. Instead, the standard method involves thermodynamic integration of the whole solid to a reference system, such as an Einstein crystal. This is a big change (involving the full extensive free energy), and so it is an expensive calculation to compute the free-energy along this path. In our research we are attempting to develop efficient free-energy methods for solids.
Molecular and engineering modeling of hydrogen fluorideThis project is concerned with the development of models that can predict and characterize a broad range of properties of hydrogen fluoride (HF) and its mixtures with water. The specific aim of the project is to develop and test molecular models appropriate for these systems, and to use the observed microscopic behaviors to formulate better engineering models. Hydrogen fluoride is an excellent candidate on which to base this development, and for many reasons:
Object-oriented programming in molecular simulationThis is less a research than a development project driven by an interest in applying molecular simulation as a teaching tool. More on that topic is available on the teaching page. We are investigating the feasibility of using the molecular simulation teaching tools as a basis for constructing extensible, object-oriented research-quality molecular simulation codes. Etomica is the name of the API and development environment that we are constructing. We now have a page devoted to describing and disseminating it. |