Free energy of hydrophobic hydration:A molecular dynamics study of noble gases in water

The potential utility and limitations of two methods to determine free energy differences from molecular dynamics simulations (MD) are studied. The computation of the free energy of hydration of the inert gases serves as a simple but illustrative example. Good results are obtained for the inert gases from a perturbation treatment, using a reference ensemble obtained from a MD simulation of a cavity in water, if these atoms are comparable in size to the cavity and the calculated free energy differences are small. This limits the applicability of the perturbation treatment of a small number of cases. Larger free energy differences can be obtained with reasonable accuracy from MD simulations with continuously changing interaction parameters. This integration method is more generally applicable, but makes an additional simulation necessary

Molecular dynamics simulations of free energy and conformational transition rates of calix[4]arene in chloroform

In a previous article we introduced a reaction coordinate based on the unstable normal mode at the saddle point of the potential energy surface. We here calculate the free-energy distribution along this coordinate for the isomerization of calix[4]arene in vacuo and in chloroform using umbrella sampling, with one umbrella covering the entire range of the reaction coordinate. An excellent first guess at this umbrella is obtained by performing a normal-mode analysis at various points along the reaction path. The isomerization rate constant of this reaction is determined using the reactive flux method and is found to be in good agreement with experimental data. The rate was found to be independent of the location of the transition state, as it should be

GronOR:Massively parallel and GPU-accelerated non-orthogonal configuration interaction for large molecular systems

GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented

ASCR/HEP Exascale Requirements Review Report

This draft report summarizes and details the findings, results, and recommendations derived from the ASCR/HEP Exascale Requirements Review meeting held in June, 2015. The main conclusions are as follows. 1) Larger, more capable computing and data facilities are needed to support HEP science goals in all three frontiers: Energy, Intensity, and Cosmic. The expected scale of the demand at the 2025 timescale is at least two orders of magnitude -- and in some cases greater -- than that available currently. 2) The growth rate of data produced by simulations is overwhelming the current ability, of both facilities and researchers, to store and analyze it. Additional resources and new techniques for data analysis are urgently needed. 3) Data rates and volumes from HEP experimental facilities are also straining the ability to store and analyze large and complex data volumes. Appropriately configured leadership-class facilities can play a transformational role in enabling scientific discovery from these datasets. 4) A close integration of HPC simulation and data analysis will aid greatly in interpreting results from HEP experiments. Such an integration will minimize data movement and facilitate interdependent workflows. 5) Long-range planning between HEP and ASCR will be required to meet HEP's research needs. To best use ASCR HPC resources the experimental HEP program needs a) an established long-term plan for access to ASCR computational and data resources, b) an ability to map workflows onto HPC resources, c) the ability for ASCR facilities to accommodate workflows run by collaborations that can have thousands of individual members, d) to transition codes to the next-generation HPC platforms that will be available at ASCR facilities, e) to build up and train a workforce capable of developing and using simulations and analysis to support HEP scientific research on next-generation systems.Comment: 77 pages, 13 Figures; draft report, subject to further revisio

The Entropy Production Fluctuation Theorem and the Nonequilibrium Work Relation for Free Energy Differences

There are only a very few known relations in statistical dynamics that are valid for systems driven arbitrarily far-from-equilibrium. One of these is the fluctuation theorem, which places conditions on the entropy production probability distribution of nonequilibrium systems. Another recently discovered far-from-equilibrium expression relates nonequilibrium measurements of the work done on a system to equilibrium free energy differences. In this paper, we derive a generalized version of the fluctuation theorem for stochastic, microscopically reversible dynamics. Invoking this generalized theorem provides a succinct proof of the nonequilibrium work relation.Comment: 6 pages, 5 figures. Typos removed, references added. Includes an expanded introductio

Effect of the integration method on the accuracy and computational efficiency of free energy calculations using thermodynamic integration

Although calculations of free energy using molecular dynamics simulations have gained significant importance in the chemical and biochemical fields, they still remain quite computationally intensive. Furthermore, when using thermodynamic integration, numerical evaluation of the integral of the Hamiltonian with respect to the coupling parameter may introduce unwanted errors in the free energy. In this paper, we compare the performance of two numerical integration techniques-the trapezoidal and Simpson's rules and propose a new method, based on the analytic integration of physically based fitting functions that are able to accurately describe the behavior of the data. We develop and test our methodology by performing detailed studies on two prototype systems, hydrated methane and hydrated methanol, and treat Lennard-Jones and electrostatic contributions separately. We conclude that the widely used trapezoidal rule may introduce systematic errors in the calculation, but these errors are reduced if Simpson's rule is employed, at least for the electrostatic component. Furthermore, by fitting thermodynamic integration data, we are able to obtain precise free energy estimates using significantly fewer data points (5 intermediate states for the electrostatic component and 11 for the Lennard-Jones term), thus significantly decreasing the associated computational cost. Our method and improved protocol were successfully validated by computing the free energy of more complex systems hydration of 2-methylbutanol and of 4-nitrophenol-thus paving the way for widespread use in solvation free energy calculations of drug molecules

Ocular disease in patients with ANCA-positive vasculitis

Anti-neutrophil cytoplasmic antibody (ANCA)-positive vasculitis—the term recently applied to Wegener's granulomatosis—is a rare multi-system inflammation characterized by necrotizing granulomas and vasculitis. We investigated the ocular manifestations of this disease in a group of patients drawn from five inflammatory eye disease clinics across the United States. Of 8,562 persons with ocular inflammation, 59 individuals were diagnosed with ANCA-positive vasculitis; 35 males and 21 females, aged 16 to 96 years, were included in this study. Ocular diagnoses were scleritis (75.0%), uveitis (17.9%), and other ocular inflammatory conditions (33.9%) including peripheral ulcerative keratitis and orbital pseudotumor. Mean duration of ocular disease was 4.6 years. Oral corticosteroids and other systemic immunosuppressive agents were used by 85.7% and 78.5% of patients, respectively. Over time, patients with ANCA-positive vasculitis experienced 2.75-fold higher mortality than other patients with inflammatory eye disease