Web: http://www.ima.umn.edu | Email: ima-staff@ima.umn.edu | Telephone: (612) 624-6066 | Fax: (612) 626-7370
Additional newsletters available at http://www.ima.umn.edu/newsletters

IMA Newsletter #391

May 2009

2008-2009 Program

Mathematics and Chemistry

See http://www.ima.umn.edu/2008-2009 for a full description of the 2008-2009 program on Mathematics and Chemistry.

News and Notes

IMA Events

IMA Workshop

MOLCAS

May 4-8, 2009

Organizers: Roland Lindh (Lund University)

IMA Workshop

Symposium on Spatio-Temporal Reaction-Diffusion Phenomena

May 12, 2009

Organizers: James W. Evans (Iowa State University), Arnd Scheel (University of Minnesota Twin Cities)

IMA Tutorial

Methods of Molecular Simulation

May 15-16, 2009

Organizers: Benedict Leimkuhler (University of Edinburgh), Berend Smit (University of California)

IMA Annual Program Year Workshop

Molecular Simulations: Algorithms, Analysis, and Applications

May 18-22, 2009

Organizers: Weinan E (Princeton University), Ron Elber (University of Texas), Daan Frenkel (University of Cambridge), Tamar Schlick (New York University)
Schedule

Friday, May 1

10:45am-11:15amCoffee breakLind Hall 400
1:25pm-2:25pmComputational challenges in cancer therapeuticsPaul Rejto (Pfizer Global R&D)Vincent Hall 570 IPS

Monday, May 4

All DayArrival day
From 2:00pm - 6:00pm Registration and hardware/software verification
(Roland Lindh and Valera Veryazov)
Lind Hall 305
SW5.4-8.09
10:45am-11:15amCoffee breakLind Hall 400

Tuesday, May 5

8:30am-8:45amWelcome to the IMAFadil Santosa (University of Minnesota)Lind Hall 305 SW5.4-8.09
8:45am-9:00amPresentation of the workshopRoland Lindh (Lund University)Lind Hall 305 SW5.4-8.09
9:00am-10:00amLecture 1: Overview of MOLCAS Valera Veryazov (Lund University)Lind Hall 305 SW5.4-8.09
10:00am-10:30amCoffee breakLind Hall 400 SW5.4-8.09
10:30am-11:15amLecture 2: Basis sets, CD/RIRoland Lindh (Lund University)Lind Hall 305 SW5.4-8.09
11:15am-12:00pmLecture 3: Use of symmetry
HF, DFT
Per-Åke Malmqvist (Lund University)
Valera Veryazov (Lund University)
Lind Hall 305 SW5.4-8.09
12:00pm-1:30pmLunch SW5.4-8.09
1:30pm-3:00pmLab 1: Code installation, testing, tailoring, tools...Lind Hall 305 SW5.4-8.09
3:00pm-3:30pmCoffee breakLind Hall 400 SW5.4-8.09
3:30pm-4:30pmLecture 4: Multiconfigurational methods: CASSCF/CASPT2Roland Lindh (Lund University)Lind Hall 305 SW5.4-8.09
4:30pm-5:30pmLab 2: Point energy calculations HF, DFTLind Hall 305 SW5.4-8.09
5:30pm-7:00pmLab 3: Point energy calculations CASSCF, CASPT2Lind Hall 305 SW5.4-8.09

Wednesday, May 6

8:30am-9:00amCoffee Lind Hall 400 SW5.4-8.09
9:00am-10:00amLecture 5: Geometry optimizations with and without constraintsRoland Lindh (Lund University)Lind Hall 305 SW5.4-8.09
10:00am-10:30amCoffee breakLind Hall 400 SW5.4-8.09
10:30am-12:30pmLab 4: Structure optimizations (minima, Hessian) Lind Hall 305 SW5.4-8.09
12:30pm-2:00pmLunch SW5.4-8.09
2:00pm-3:00pmLecture 6: Relativistic calculationsPer-Åke Malmqvist (Lund University)Lind Hall 305 SW5.4-8.09
3:00pm-3:30pmCoffee breakLind Hall 400 SW5.4-8.09
3:30pm-5:00pmLab 5: Structure optimizations (transition states, reactions, geometrical constraints)Lind Hall 305 SW5.4-8.09
5:00pm-5:30pmLecture 7: MOLCAS for advanced usersValera Veryazov (Lund University)Lind Hall 305 SW5.4-8.09
5:30pm-7:00pmLab 6: MOLCAS as a development platform. Implementing codes in MOLCASLind Hall 305 SW5.4-8.09

Thursday, May 7

8:30am-9:00amCoffeeLind Hall 400 SW5.4-8.09
9:00am-10:00amLecture 8: Excited statesLuis Serrano-Andrés (University of Valencia)Lind Hall 305 SW5.4-8.09
10:00am-10:30amCoffee breakLind Hall 400 SW5.4-8.09
10:30am-11:30amLab 7: Excited states (SA-CASSCF, CASPT2) Lind Hall 305 SW5.4-8.09
11:30am-12:30pmLab 8: Excited states and symmetry Lind Hall 305 SW5.4-8.09
12:30pm-2:00pmLunch SW5.4-8.09
1:00pm-2:00pmSpecial seminar: Protein folding rates from a dynamical systems perspectiveDaniel Dix (University of South Carolina)Lind Hall 409
2:00pm-3:00pmLab 8: cont.Lind Hall 305 SW5.4-8.09
3:00pm-3:30pmCoffee breakLind Hall 400 SW5.4-8.09
3:30pm-5:00pmLab 9: Excited states (optimizations, crossings, minimum energy paths)Lind Hall 305 SW5.4-8.09
5:00pm-6:30pmLab 10: Exited states CASSI, SOC Lind Hall 305 SW5.4-8.09
7:00pm-8:30pmWorkshop Dinner at PagodaPagoda Restaurant
1417 4th St. SE
Minneapolis, MN
612-378-4710
SW5.4-8.09

Friday, May 8

8:30am-9:00amCoffee Lind Hall 400 SW5.4-8.09
9:00am-9:30amLecture 9: QM/MM in Molcas Roland Lindh (Lund University)Lind Hall 305 SW5.4-8.09
9:30am-10:00amLecture 10: Solvent models Luis Serrano-Andrés (University of Valencia)Lind Hall 305 SW5.4-8.09
10:00am-10:30amCoffee breakLind Hall 400 SW5.4-8.09
10:30am-11:30amLab 11: Solvent models (ground and excited states) Lind Hall 305 SW5.4-8.09
11:30am-12:30pmLecture 11: Handling problems in MOLCAS (convergence, orbital selection, etc.)Luis Serrano-Andrés (University of Valencia)Lind Hall 305 SW5.4-8.09
12:30pm-2:00pmLunch SW5.4-8.09
2:00pm-2:30pmConclusions Roland Lindh (Lund University)Lind Hall 305 SW5.4-8.09
2:30pm-2:40pmGroup photoLind Hall 305 SW5.4-8.09
2:40pm-3:10pmCoffee and handing out of diplomasLind Hall 400 SW5.4-8.09
3:10pm-4:10pmOpen Lab (voluntary presence)Lind Hall 305 SW5.4-8.09

Monday, May 11

10:45am-11:15amCoffee breakLind Hall 400
3:30pm-4:30pmThe Truhlar group seminar: Density matrix treatment of dissipative dynamics in molecules and at surfaces David A. Micha (University of Florida)283 Kolthoff Hall 3:30

Tuesday, May 12

10:45am-11:15amCoffee breakLind Hall 400
11:15am-12:15pmMotor mediated deformation of microtubules in living cells and gliding assaysErkan Tüzel (University of Minnesota)Lind Hall 305 PS
12:00pm-1:00pmLunch snacks
(for symposium participants only)
Lind Hall 400 SW5.12.09
1:00pm-1:40pmCoherent structures in spatially extended systemsArnd Scheel (University of Minnesota)Lind Hall 409 SW5.12.09
1:40pm-1:50pmBreakLind Hall 400 SW5.12.09
1:50pm-2:30pmLattice reaction-diffusion equation analysis for autocatalysisDa-Jiang Liu (Ames Laboratory)Lind Hall 409 SW5.12.09
2:30pm-2:40pmBreakLind Hall 400 SW5.12.09
2:40pm-3:20pmStability of spikes: from autocatalysis to precipitationAlin Pogan (University of Minnesota)Lind Hall 409 SW5.12.09
3:20pm-3:40pmBreakLind Hall 400 SW5.12.09
3:40pm-4:20pmNonlinear oscillations in a closed chemical reactorYongfeng Li (University of Minnesota)Lind Hall 409 SW5.12.09
4:20pm-4:30pmBreakLind Hall 400 SW5.12.09
4:30pm-5:10pmFrom molecular-level reaction models to continuum reaction-diffusion equationsJames W. Evans (Iowa State University)Lind Hall 409 SW5.12.09

Wednesday, May 13

10:45am-11:15amCoffee breakLind Hall 400
11:15am-12:15pmNon-ergodicity of the Nosé-Hoover dynamicsFrédéric Legoll (École Nationale des Ponts-et-Chaussées)Lind Hall 409 SMC

Thursday, May 14

10:45am-11:15amCoffee breakLind Hall 400
11:15am-12:15pmAdaptive methods for efficient sampling. Applications in molecular dynamicsTony Lelievre (Ecole Nationale des Ponts et Chaussees)Lind Hall 409 SMC
1:00pm-2:00pmReading group presentation: A graph theoretical data structure for indexing molecular internal coordinates, and IMIMOL, a computer program for 3D molecule building and geometric analysis Daniel Dix (University of South Carolina)Lind Hall 409

Friday, May 15

8:15am-8:45amRegistration and coffeeEE/CS 3-176 T5.15-16.09
8:50am-9:00amWelcome to the IMAFadil Santosa (University of Minnesota)EE/CS 3-180 T5.15-16.09
9:00am-9:45amLecture 1: Introduction to molecular simulation.
An overview of molecular simulation and its application.
Berend Smit (University of California, Berkeley)EE/CS 3-180 T5.15-16.09
9:45am-10:30amLecture 2: Molecular dynamics: foundations I.
Dynamical issues, ergodicity, and averages at constant energy.
Benedict Leimkuhler (University of Edinburgh)EE/CS 3-180 T5.15-16.09
10:30am-11:00amBreakEE/CS 3-176 T5.15-16.09
11:00am-11:45amLecture 3: Molecular dynamics: foundations II.
Geometric integrators and applications to molecular dynamics.
Benedict Leimkuhler (University of Edinburgh)EE/CS 3-180 T5.15-16.09
11:45am-2:00pmLunch break T5.15-16.09
2:00pm-2:45pmLecture 4: Basics of statistical mechanics.
Description of statistical mechanical ensembles and implications for simulation.
Berend Smit (University of California, Berkeley)EE/CS 3-180 T5.15-16.09
2:45pm-3:30pmSpecial Topic 1: Quasicontinuum methods for coupling atomistic and continuum modelsMitchell Luskin (University of Minnesota)EE/CS 3-180 T5.15-16.09
3:30pm-4:00pmBreakEE/CS 3-176 T5.15-16.09
4:00pm-4:45pmSpecial Topic 2: Multiple time-scale methodsFrédéric Legoll (École Nationale des Ponts-et-Chaussées)EE/CS 3-180 T5.15-16.09

Saturday, May 16

8:30am-9:00amCoffeeEE/CS 3-176 T5.15-16.09
9:00am-9:45amLecture 5: Control of the ensemble.
Controlling temperature and pressure in molecular dynamics.
Benedict Leimkuhler (University of Edinburgh)EE/CS 3-180 T5.15-16.09
9:45am-10:30amLecture 6: Free energy calculation.
Methods for coarse-graining a complex energy surface, reaction coordinates
Robert D. Skeel (Purdue University)EE/CS 3-180 T5.15-16.09
10:30am-11:00amBreakEE/CS 3-176 T5.15-16.09
11:00am-11:45amLecture 7: Monte Carlo methods.
Methods for sampling based on Monte Carlo and variants.
Berend Smit (University of California, Berkeley)EE/CS 3-180 T5.15-16.09
11:45am-2:00pmLunch break T5.15-16.09
2:00pm-2:45pmLecture 8: Methods for polymer models.
Coarse grained models for mesoscale simulation, and methodology such as dissipative particle dynamics.
Berend Smit (University of California, Berkeley)EE/CS 3-180 T5.15-16.09
2:45pm-3:30pmSpecial topic 3. Computation of force fields for molecular simulationRobert Krasny (University of Michigan)EE/CS 3-180 T5.15-16.09
3:30pm-4:00pmBreakEE/CS 3-176 T5.15-16.09
4:00pm-4:45pmSpecial Topic 4. Adaptive methods for free energy computationTony Lelievre (Ecole Nationale des Ponts et Chaussees)EE/CS 3-180 T5.15-16.09

Monday, May 18

All DayMorning Chair: Daan Frenkel (University of Cambridge)
Afternoon Chair: Giovanni Ciccotti (Università di Roma "La Sapienza")
W5.18-22.09
8:15am-9:00amRegistration and coffeeEE/CS 3-180 W5.18-22.09
9:00am-9:15amWelcome to the IMAFadil Santosa (University of Minnesota)EE/CS 3-180 W5.18-22.09
9:15am-9:55amCoupled atomistic-continuum methods for fluidsWeiqing Ren (New York University)EE/CS 3-180 W5.18-22.09
9:55am-10:15amDiscussionEE/CS 3-180 W5.18-22.09
10:15am-10:45amBreakEE/CS 3-176 W5.18-22.09
10:45am-11:25amDynamics and thermodynamics of large conformational rearrangements in confined polymersDmitrii E. Makarov (University of Texas)EE/CS 3-180 W5.18-22.09
11:25am-11:45amDiscussionEE/CS 3-180 W5.18-22.09
11:45am-2:00pmLunch W5.18-22.09
2:00pm-2:40pmCoarse-graining molecular dynamics for crystalline solids Xiantao Li (Pennsylvania State University)EE/CS 3-180 W5.18-22.09
2:40pm-3:00pmDiscussionEE/CS 3-180 W5.18-22.09
3:00pm-3:40pmAdaptive kinetic Monte Carlo for simulating dynamics in atomic systems Graeme Henkelman (University of Texas)EE/CS 3-180 W5.18-22.09
3:40pm-4:00pmDiscussionEE/CS 3-176 W5.18-22.09
4:00pm-6:30pmReception and Poster Session
Poster submissions welcome from all participants
Instructions
Lind Hall 400 W5.18-22.09
Molecular dynamics simulation of Liquid methyl chloride with a treecode Ewald summation method Henry A. Boateng (University of Michigan)
Pathwise accuracy and ergodicity of metropolized integrators for SDEsNawaf Bou-Rabee (New York University)
Molecular modelling with help from invariant theory and computer algebraBastiaan J. Braams (Emory University)
The string method as a dynamical systemMaria Kourkina Cameron (Courant Institute of Mathematical Sciences)
Investigating the importance of quantum mechanical recrossing effects in the rate coefficient of a condensed-phase reaction Rosana Collepardo-Guevara (University of Oxford)
Atomistic and heterogeneous multiscale modeling of CO-oxidation on metal(100) surfaces: From nanoscale ordering to mesoscale frontsJames W. Evans (Iowa State University)
The role of Arg517 in defining the fidelity of DNA polymerase lambdaMeredith Catherine Foley (New York University)
Interface method and Green's function based Poisson–Boltzmann equation solver and interface technique based molecular dynamicsWeihua Geng (University of Michigan)
Competing quantum effects in liquid waterScott Habershon (University of Oxford)
The interplay of AAA+ molecular machines and sliding clamps at the DNA replication forkIvaylo N. Ivanov (University of California, San Diego)
Exotic self-assembly scenarios in in two-dimensional dipolar mixtures with genetic algorithmsGerhard Kahl (Technische Universität Wien)
Hard-core colloids forming stable crystals: Clusters - columns - lamellae - compactGerhard Kahl (Technische Universität Wien)
Computational alanine scanning of the alpha-lactamase inhibitor protein and TEM-1 alpha-lactamase complexHiqmet Kamberaj (University of Minnesota)
Yuk Sham (University of Minnesota)
A Cartesian treecode for screened Coulomb interactionsRobert Krasny (University of Michigan)
Analysis of folding pathways in model peptides Krzysztof Kuczera (University of Kansas)
Membrane channels: Diffusion in confined geometries and entropic barriersJuan C. Latorre (Freie Universität Berlin)
Kinetic Monte Carlo simulation of the Yttria Stabilized Zirconia (YSZ) fuel cellKah Chun Lau (George Washington University)
Effective dynamics for reaction coordinates Frédéric Legoll (École Nationale des Ponts-et-Chaussées)
Flexible-boundary QM/MM: Partial charge transfer between the QM and MM subsystemsHai Lin (University of Colorado)
An image-based reaction field method for electrostatic interactions in molecular dynamics simulations of aqueous solutionsYuchun Lin (University of North Carolina - Charlotte)
The impact of calcium on the dynamics and energetics of the complex between calmodulin and anthrax edema factorThérèse E. Malliavin (Centre National de la Recherche Scientifique (CNRS))
Fast path integral simulations using ring polymer contractionThomas Markland (University of Oxford)
Density matrix treatment of optical response with combined instantaneous and delayed dissipation: Adsorbates on solid surfacesDavid A. Micha (University of Florida)
Density-dependent analysis of nonequilibrium paths improves free energy estimatesDavid Do Le Minh (National Institutes of Health)
Some improvements of the ART method for finding transition pathways on potential energy surfacesKimiya Minoukadeh (Ecole Nationale des Ponts et Chaussees)
Simulations of phase transitions for fluids in pores using dynamic mean field theoryPeter A. Monson (University of Massachusetts)
What is a good thermostat? Emad Noorizadeh (University of Edinburgh)
Dynamical pathways to inactivation in the KcsA potassium channel Albert C. Pan (University of Chicago)
Concurrent triple-scale simulation of liquid waterMatej Praprotnik (National Institute of Chemistry)
Molecular dynamics simulations of binding of metal catalyst nanoparticles to carbon nanostructures by using ReaxFFCarlos F. Sanz-Navarro (Autonomous University of Barcelona)
Liquid-liquid critical point in supercooled siliconSrikanth Sastry (Jawaharlal Nehru Centre for Advanced Scientific Research(JNCASR))
Optimal fuzzy aggregation of Markov chainsChristof Schütte (Freie Universität Berlin)
Eric Vanden-Eijnden ()
Maximum flux transition pathJuanfang Shen (Purdue University)
Robert D. Skeel (Purdue University)
Ruijun Zhao (Purdue University)
Counting for rigidity, flexibility and extensions via the pebble game algorithm – hinge predictions and allosteryAdnan Sljoka (York University)
Statistical behavior of the self-guided dynamics algorithmsDaniel B. Smith (University of Pittsburgh)
Adaptive importance sampling strategiesGabriel Stoltz (École Nationale des Ponts-et-Chaussées (ENPC))
Models of Chemical potential and the concept of atomistic embeddingSteven M. Valone (Los Alamos National Laboratory)
What ice can teach us about water interactions: a critical comparison of the performance of different water modelsCarlos Vega de las Heras (Universidad Complutense de Madrid)
Microscopic flow around a diffusing particle Rodolphe Vuilleumier (École Normale Supérieure)
Projected dynamics from molecular dynamics simulations Rodolphe Vuilleumier (École Normale Supérieure)

Tuesday, May 19

All DayMorning Chair: Ron Elber (University of Texas)
Afternoon Chair: Tamar Schlick (New York University)
W5.18-22.09
8:30am-9:00amCoffeeEE/CS 3-176 W5.18-22.09
9:00am-9:40amModeling organic-inorganic interfaces across different scales: Modeling concepts and ApplicationsLuigi Delle Site (Max-Planck Institut für Polymerforschung)EE/CS 3-180 W5.18-22.09
9:40am-10:00amDiscussionEE/CS 3-180 W5.18-22.09
10:00am-10:40amBenard convection by dynamical nonequilibrium molecular dynamicsGiovanni Ciccotti (Università di Roma "La Sapienza")EE/CS 3-180 W5.18-22.09
10:40am-11:00amDiscussionEE/CS 3-180 W5.18-22.09
11:00am-11:30amBreakEE/CS 3-176 W5.18-22.09
11:30am-12:10pmFree-energy methods to study conformational transitions and ligand binding to flexible targetsFrancesco Luigi Gervasio (Eidgenössische TH Zürich)EE/CS 3-180 W5.18-22.09
12:10pm-12:30pmDiscussionEE/CS 3-180 W5.18-22.09
12:30pm-2:00pmLunch W5.18-22.09
2:00pm-2:40pmEnergy landscapes of clusters, glasses, and biomoleculesDavid John Wales (University of Cambridge)EE/CS 3-180 W5.18-22.09
2:40pm-3:00pmDiscussionEE/CS 3-180 W5.18-22.09
3:00pm-3:15pmBreakEE/CS 3-176 W5.18-22.09
3:15pm-3:55pmAlgorithms for millisecond-scale simulation on Anton, a specialized molecular Dynamics machineRon Dror (D. E. Shaw Research)EE/CS 3-180 W5.18-22.09
3:55pm-4:15pmDiscussionEE/CS 3-180 W5.18-22.09

Wednesday, May 20

All Day1:00pm-3:00pm Chair: Benedict Leimkuhler (University of Edinburgh)
3:30pm-5:30pm Chair: Tamar Schlick (New York University)
W5.18-22.09
10:45am-11:15amCoffee breakLind Hall 400
11:15am-12:15pmReading group presentation: A graph theoretical data structure for indexing molecular internal coordinates, and IMIMOL, a computer program for 3D molecule building and geometric analysis Daniel Dix (University of South Carolina)Lind Hall 409
1:00pm-1:40pmFolding landscape of RNA from single molecule trajectories Devarajan Thirumalai (University of Maryland)EE/CS 3-180 W5.18-22.09
1:40pm-2:00pmDiscussionEE/CS 3-180 W5.18-22.09
2:00pm-2:40pmMarkov state models for biomolecular dynamics: Theory and applicationsChristof Schütte (Freie Universität Berlin)EE/CS 3-180 W5.18-22.09
2:40pm-3:00pmDiscussionEE/CS 3-180 W5.18-22.09
3:00pm-3:30pmBreakEE/CS 3-176 W5.18-22.09
3:30pm-4:10pmCoarse-grained molecular simulation of soft matter and interfacesMichael Philip Allen (University of Warwick)EE/CS 3-180 W5.18-22.09
4:10pm-4:30pmDiscussionEE/CS 3-180 W5.18-22.09
4:30pm-5:10pmHybrid field-based models of cell membranes and proteins for signaling applicationsRavi Radhakrishnan (University of Pennsylvania)EE/CS 3-180 W5.18-22.09
5:10pm-5:30pmDiscussionEE/CS 3-180 W5.18-22.09
5:30pm-5:40pmGroup photo W5.18-22.09

Thursday, May 21

All DayMorning Chair: Weinan E (Princeton University)
Afternoon Chair: Christof Schütte (Freie Universität Berlin)
W5.18-22.09
8:30am-9:00amCoffeeEE/CS 3-176 W5.18-22.09
9:00am-9:40amProtein path sampling: Bring your algorithm and your modelDaniel M. Zuckerman (University of Pittsburgh)EE/CS 3-180 W5.18-22.09
9:40am-10:00amDiscussionEE/CS 3-180 W5.18-22.09
10:00am-10:40amGenerating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniquesPeter Bolhuis (Universiteit van Amsterdam)EE/CS 3-180 W5.18-22.09
10:40am-11:00amDiscussionEE/CS 3-180 W5.18-22.09
11:00am-11:30amBreakEE/CS 3-176 W5.18-22.09
11:30am-12:10pmTheory and modeling of reactive eventsEric Vanden-Eijnden (New York University)EE/CS 3-180 W5.18-22.09
12:10pm-12:30pmDiscussionEE/CS 3-180 W5.18-22.09
12:30pm-2:00pmLunch W5.18-22.09
2:00pm-2:40pmComputing slow reaction rates in molecular systemsEric Felix Darve (Stanford University)EE/CS 3-180 W5.18-22.09
2:40pm-3:00pmDiscussionEE/CS 3-180 W5.18-22.09
3:00pm-3:15pmBreakEE/CS 3-176 W5.18-22.09
3:15pm-3:55pmRheological properties of polymer melt between rapidly oscillating plates: an application of multiscale modelingRyoichi Yamamoto (Kyoto University)EE/CS 3-180 W5.18-22.09
3:55pm-4:15pmDiscussionEE/CS 3-180 W5.18-22.09
6:30pm-8:30pmWorkshop dinnerPagoda Restaurant
1417 4th St. SE
Minneapolis, MN
612-378-4710
W5.18-22.09

Friday, May 22

All DayChair: Robert D. Skeel (Purdue University) W5.18-22.09
8:30am-9:00amCoffeeEE/CS 3-176 W5.18-22.09
9:00am-9:40amMilestoningRon Elber (University of Texas)EE/CS 3-180 W5.18-22.09
9:40am-10:00amDiscussionEE/CS 3-180 W5.18-22.09
10:00am-10:30amBreakEE/CS 3-176 W5.18-22.09
10:30am-11:10amFlexible control of the ensembleBenedict Leimkuhler (University of Edinburgh)EE/CS 3-180 W5.18-22.09
11:10am-11:30amDiscussionEE/CS 3-180 W5.18-22.09
11:30am-11:40amClosing remarksEE/CS 3-180 W5.18-22.09
1:00pm-2:00pmAdiabatic switching for degenerate ground statesGabriel Stoltz (École Nationale des Ponts-et-Chaussées (ENPC))Lind Hall 409 SMC

Monday, May 25

All DayMemorial Day. The IMA is closed.

Tuesday, May 26

10:45am-11:15amCoffee breakLind Hall 400
11:15am-12:15pmThe Coulomb barrier in excitonic charge separationXiaoyang Zhu (University of Minnesota)Lind Hall 409 PS

Wednesday, May 27

10:45am-11:15amCoffee breakLind Hall 400

Thursday, May 28

10:45am-11:15amCoffee breakLind Hall 400

Friday, May 29

10:45am-11:15amCoffee breakLind Hall 400
Abstracts
Michael Philip Allen (University of Warwick) Coarse-grained molecular simulation of soft matter and interfaces
Abstract: In this talk I shall describe two kinds of soft matter system in which the interfacial behaviour is of interest. Firstly, the adsorption of polymers and peptides on solid surfaces is of growing interest because of the possible applications in biotechnology. We have recently improved a Monte Carlo simulation method for studying lattice polymers near to surfaces, and I shall present some recent results. Secondly, I shall discuss the simulation of confined liquid crystals, with special emphasis on hard platelets. These idealized systems may be related to experiments on colloidal suspensions.
Henry A. Boateng (University of Michigan) Molecular dynamics simulation of Liquid methyl chloride with a treecode Ewald summation method
Abstract: We will present results on the application of the treecode method to the computation of the real space part of the Ewald sum in molecular dynamics simulation. We show that for reasonable parameters, the treecode speeds up the Ewald sum for large systems and reproduces the structural and dynamical properties of the liquid within error bars.
Peter Bolhuis (Universiteit van Amsterdam) Generating, sampling, and analyzing the dynamical pathways of folding proteins with rare event techniques
Abstract: All-atom explicit solvent molecular dynamics simulations of (un)folding and other conformational changes in proteins remain a great challenge because of the long time scales involved. These long times are due to high free energy barriers between metastable states. Rare event simulation methodology, such as Replica Exchange (REMD) and Metadynamics can overcome these free energy barriers, and give thermodynamic equilibrium properties. However, they do not preserve the dynamical pathways needed to gain insight into kinetics and mechanism of conformational change. Transition Path Sampling (TPS) is a technique designed for harvesting an ensemble of unbiased dynamical transition paths, and allows computation of the rates and extraction of the reaction coordinates. In this presentation I will discuss the various rare event methods their applicability and limitations. I will illustrate these methods on several examples: the folding of the 20-residue Trp-cage, the folding of a 37-residue WW-domain and the unfolding step in the photocycle of the 125-residue bacterial sensor Photoactive Yellow protein (PYP). For each of these systems we have studied the mechanism, reaction coordinates, folding rates, and/or barriers heights. This work shows that by employing a toolbox of computational rare event techniques we can overcome the timescale gap inherent in Molecular Dynamics of rare events and gain more insight in complex transitions such as protein folding and unfolding.
Nawaf Bou-Rabee (New York University) Pathwise accuracy and ergodicity of metropolized integrators for SDEs
Abstract: Metropolized integrators for ergodic stochastic differential equations (SDE) are exhibited which (i) are ergodic with respect to the (known) equilibrium distribution of the SDE and (ii) approximate pathwise the solutions of the SDE on finite time intervals. Both these properties are demonstrated including precise strong error estimates. It is also shown that the Metropolized integrator retains these properties even in situations where the drift in the SDE is nonglobally Lipschitz, and vanilla explicit integrators for SDEs typically become unstable and fail to be ergodic.
Bastiaan J. Braams (Emory University) Molecular modelling with help from invariant theory and computer algebra
Abstract: Invariant theory concerns an algebra K[V]G: the polynomials on a vector space V over a field K that are invariant under the action of a group G that is represented on V. In the case of a finite group acting on a finite-dimensional vector space over a field of characteristic zero, and in other cases not of present concern, this algebra is finitely generated in the Cohen-Macaulay form of a finite-dimensional free module over a polynomial ring. Generators of that ring (primary invariants) and a basis for the module (secondary invariants) may be found by computer algebra. Molecular modelling encompasses many things: molecular dynamics for the study of reaction and isomerization processes, diffusion Monte Carlo for ground state properties, path-integral Monte Carlo for thermal properties, and direct diagonalization of a configuration-interaction matrix for calculation of a ro-vibrational spectrum. A key tool in molecular modelling for any specific system is an analytical, fitted potential energy surface: the potential energy of the system as a function of the spatial configuration of nuclei, as is obtained from first principles by solution of the electronic Schroedinger equation. In the work to be described here molecular modelling meets invariant theory and computer algebra in the construction of potential energy surfaces. A polynomial model in suitable variables is employed and the model is explicitly invariant under the full molecular permutation symmetry group. The MAGMA computer algebra system is employed to obtain the generators. Using such a representation we are able to construct quite routinely highly accurate potential energy surfaces for systems of N=7 or more atoms in full (3N-6) dimensionality, and with use of a many-body expansion also for larger systems such as water clusters.
Maria Kourkina Cameron (Courant Institute of Mathematical Sciences) The string method as a dynamical system
Abstract: This project is devoted to the theoretical analysis of the zero-temperature string method to identify minimum energy paths (MEP's) on a given energy landscape. The MEP's are curves connecting local minima on the energy landscape which are everywhere tangent to the gradient of the potential except possibly at critical points. In practice, MEP's are mountain pass curves that play a special role e.g. in the context of rare reactive events that occur when one considers a steepest descent dynamics on the potential perturbed by a small random noise. The string method aims at identifying MEP's by moving each point of the curve by steepest descent on the energy landscape. Here we address the question of whether such a curve evolution necessarily converges to a MEP. Surprisingly, the answer is no, for an interesting reason: MEP's may not be isolated, in the sense that there may be families of them that can be continuously deformed into one another. This degeneracy is related to the presence of critical points of Morse index 2 or higher along the MEP. In this paper, we elucidate this issue and characterize completely the limit set of a curve evolving by the string method. We establish rigorously that the limit set of such a curve is again a curve when the MEP's are isolated. We also show under the same hypothesis that the string evolution converges to an MEP. However, we also identify and classify situations where the limit set is not a curve and may contain higher-dimensional parts. We present a collection of examples where the limit set of a path contains a 2D region, a 2D surface, or a 3D body. In some of our examples the evolving path wanders around without converging to its limit set. In other examples it fills a region, converging to its limit set which is not an MEP.
Giovanni Ciccotti (Università di Roma "La Sapienza") Benard convection by dynamical nonequilibrium molecular dynamics
Abstract: We present the need (especially critical for multiscale simulations) of extending standard equilibrium Molecular Dynamics to nonequilibrium situations. We show how this is possible in full rigor by use of the nonlinear Kubo-Onsager relation, connecting dynamical nonequilibrium averages to equilibrium (more genarally, initial distribution) averages over suitable time-evolved variables. To show the power of the method we apply it to study the onset of the convective circulation in liquids establishing as stationary asymptotic state the convective rolls know also as Benard cells.
Rosana Collepardo-Guevara (University of Oxford) Investigating the importance of quantum mechanical recrossing effects in the rate coefficient of a condensed-phase reaction
Abstract: We have used the ring polymer molecular dynamics method to study the Azzouz-Borgis model for proton transfer between phenol (AH) and trimethylamine (B) in liquid methyl chloride. When the A-H distance is used as the reaction coordinate, the ring polymer trajectories are found to exhibit multiple recrossings of the transition state dividing surface and to give a rate coefficient that is smaller than the quantum transition state theory value by an order of magnitude. This is to be expected on kinematic grounds for a heavy-light-heavy reaction when the light atom transfer coordinate is used as the reaction coordinate, and it clearly precludes the use of transition state theory with this reaction coordinate. As has been shown previously for this problem, a solvent polarization coordinate defined in terms of the expectation value of the proton transfer distance in the ground adiabatic quantum state provides a better reaction coordinate with less recrossing. These results are discussed in light of the wide body of earlier theoretical work on the Azzouz-Borgis model and the considerable range of previously reported values for its proton and deuteron transfer rate coefficients.
Eric Felix Darve (Stanford University) Computing slow reaction rates in molecular systems
Abstract: Joint work with Amirali Kia and Jose Solomon. The Mori-Zwanzig formalism is an effective tool to derive differential equations describing the evolution of a small number of resolved variables. We will present its application to the derivation of generalized Langevin equations and generalized non-Markovian Fokker-Planck equations. We show how long time scales rates and meta-stable basins can be extracted from these equations. Numerical algorithms are proposed to discretize these equations. An important aspect is the numerical solution of the orthogonal dynamics equation which is a partial differential equation in a high dimensional space. We propose efficient numerical methods to solve this orthogonal dynamics equation. In addition, we present a projection formalism of the Mori-Zwanzig type that is applicable to discrete maps. Numerical applications are presented from the field of protein modeling.
Luigi Delle Site (Max-Planck Institut für Polymerforschung) Modeling organic-inorganic interfaces across different scales: Modeling concepts and Applications
Abstract: The study of the interaction between (large) molecules and inorganic (metal) surfaces is a typical multiscale problem. In fact while the specific chemistry emerges at the quantum scale with the interaction of the electrons of the surface with the valence electrons of the molecule, the molecular conformations are instead the expression of a large scale statistical behavior. The delicate interplay between these two aspects, intimately linked during the adsorption process, gives rise to the interface properties. I will show the application of a sequential modeling procedure which links the quantum to the classical atomistic scale in a consistent way. Within this framework the ideal approach would be to identify and use only those degrees of freedom strictly required by the problem. To achieve this goal I will illustrate the first step of the development of a spatially adaptive simulation scheme (AdResS) which allows to switch from the atomistic to a coarse-grained description on demand.
Daniel Dix (University of South Carolina) Reading group presentation: A graph theoretical data structure for indexing molecular internal coordinates, and IMIMOL, a computer program for 3D molecule building and geometric analysis
Abstract: First a mathematical structure related to partially ordered sets and iterated line graphs is defined and shown to be the appropriate index set for the usual valence internal coordinates of molecules, such as proteins. Second a computer program, IMIMOL, based on this formalism will be demonstrated. The program allows the building of three dimensional (static) molecular models by gluing together pieces. It also allows the determination of molecular shapes which satisfy particular geometric criteria. Using it, for example, an `alpha helix' can be found by solving a system of equations.
Daniel Dix (University of South Carolina) Special seminar: Protein folding rates from a dynamical systems perspective
Abstract: The protein folding process is modeled by a classical mechanical system. Particles for the nuclei of protein atoms as well as solvent atoms must be included, and are subject to a molecular mechanics force field, as well as containing forces. Individual trajectories do not exhibit protein folding in the large time limit. A probability density function on phase space, representing an initial unfolded ensemble, will evolve into an intensely oscillatory function, and the time dependent expectation of an observable quantity with respect to this evolved density can exhibit irreversible relaxation behavior. The problem of determining the rate of this relaxation is shown to be the `rate of decay of correlations ' problem in the ergodic theory of smooth flows. This problem has only been solved (by others) in certain uniformly mixing systems, and the methods do not easily generalize to cover the protein folding case. Thus we abstract theoretical features of the protein folding case and investigate the rate of decay of correlations problem in a highly simplified model system retaining those features. The model has a new weaker mechanism causing correlation decay, which can be observed even on periodic orbits of the flow.
Daniel Dix (University of South Carolina) Reading group presentation: A graph theoretical data structure for indexing molecular internal coordinates, and IMIMOL, a computer program for 3D molecule building and geometric analysis
Abstract:
Ron Dror (D. E. Shaw Research) Algorithms for millisecond-scale simulation on Anton, a specialized molecular Dynamics machine
Abstract: Anton, a recently completed special-purpose machine for molecular dynamics simulations, achieves a speedup of several orders of magnitude over state-of-the-art general-purpose systems by mapping novel algorithms to massively parallel application-specific hardware. This talk will describe algorithms and numerical techniques used by Anton, including some developed especially for Anton and others adapted to exploit its architecture. These algorithms underlie Anton's ability to avoid communication bottlenecks even at high levels of parallelism, to effectively map computation onto fast, hardwired arithmetic pipelines, and to produce accurate results efficiently and deterministically.
Ron Elber (University of Texas) Milestoning
Abstract: Milestoning is a theory for stochastic coarse graining in space and time and an algorithm to compute kinetics and reaction dynamics of complex molecular systems. I will discuss the theory and algorithm properties, implementation, and an application to the folding of a helix and to the recovery stroke of the protein myosin. I will review the difficulties of simulating functions of biological molecules on relevant time scales and suggest that Milestoning addresses some of these problems. Folding pathways of a helical peptide with an energy landscape that is expected to be a funnel is illustrated. The myosin example shows coupling between local molecular events and global transitions which extends to the millisecond time scale.
James W. Evans (Iowa State University) Atomistic and heterogeneous multiscale modeling of CO-oxidation on metal(100) surfaces: From nanoscale ordering to mesoscale fronts
Abstract: A long-standing challenge has been to develop realistic atomistic-level models of key surface reactions on single-crystal metal catalyst surfaces. We have developed realistic multi-site lattice gas models for these systems incorporating DFT energetics, and thus accurately describing the complex nanoscale ordering of reactants [1]. Model behavior is determined from KMC simulation. A second challenge is to incorporate these atomistic models into a multiscale treatment of spatio-temporal behavior on a much larger scale (of microns). This is achieved using our Heterogeneous Coupled Lattice Gas (HCLG) simulation strategy [2], together with a precise description of collective diffusion in these mixed reactant adlayers. The result is a "first-principles" description of reaction-diffusion front propagation [3]. [1] e.g., Liu, Evans, J. Chem. Phys. 124 (2006) 154705. [2] Tammaro, Sabella, Evans, J. Chem. Phys. 103 (1995) 10277. [3] Liu, Evans, Surf. Sci. (2009) doi:10.1016/j.susc.2008.10.058
Meredith Catherine Foley (New York University) The role of Arg517 in defining the fidelity of DNA polymerase lambda
Abstract: The maintenance of the genetic information in a cell is essential for its survival. DNA polymerases play a central role in this process since they replicate and repair DNA. These enzymes must be versatile and specific to accommodate the four combinations of the Watson-Crick base pairs while preventing errors such as incorrect nucleotide insertion. DNA polymerase (pol) lambda is a mammalian enzyme involved in the base excision repair (BER) and non-homologous end-joining (NHEJ) DNA repair pathways. It is from the X-family of polymerases and is shaped like a hand with finger, palm, and thumb subdomains. It also has 8-kDa and BRCT domains that aid its participation in the repair pathways. Pol lambda has an unusual fidelity profile since it produces many deletion errors and relatively few base substitution errors. Substantial experimental and computational work on a related enzyme, pol beta, has elucidated the conformational changes induced upon its binding the correct and incorrect nucleotide. These motions involve a large-scale rearrangement of the thumb and subtle motions of active-site protein residues that serve as "gate-keepers" in regulating pol beta's transition to an active state. Binding of the incorrect nucleotide hampers these transitions. Here we describe dynamics simulations of pol lambda/DNA complexes to uncover the important motions and rearrangements prior to chemistry that account for pol lambda's specific error profile. As in pol beta, a series of active-site protein residues are hypothesized to modulate the assembly of the active site. However, DNA motion is also involved. Arg517 emerges as a key residue for discriminating against incorrect nucleotide incorporation. This residue also plays an important role in the stabilization of the active position of the DNA even when the DNA is misaligned and poised for a deletion error.
Weihua Geng (University of Michigan) Interface method and Green's function based Poisson–Boltzmann equation solver and interface technique based molecular dynamics
Abstract: Implicit solvent models, which treat the solvent as a macroscopic continuum while admitting a microscopic atomic description for the biomolecule, are efficient multiscale approaches to complex, large scale biological systems. We summarize our recent advances in mathematical methods for the Poisson-Boltzmann (PB) Model. These include rigorous mathematical treatments of interface jump conditions, surface geometric singularities, charge singularities, and associated electrostatic force computation for molecular dynamics.
Francesco Luigi Gervasio (Eidgenössische TH Zürich) Free-energy methods to study conformational transitions and ligand binding to flexible targets
Abstract: Large scale conformational rearrangements represent both a challenge and an opportunity for rational drug design. Exploring the conformational space of a target protein with sufficient detail is computationally very demanding and most current docking programs are very limited in this respect. If it were possible, however, it could open the avenue to the design of more selective drug candidates. Here we show how atomistic molecular dynamics together with methods developed to accelerate rare events can be used to study both kinase conformational transitions and ligand binding. Using a new sampling method which is able to find the low free energy channel between an initial and final state[1] we determined the atomistic dynamics of the open-to-closed movement of the cyclin dependent kinase 5 (CDK5). [4] The inactivation movement has a two-step mechanism. First alpha-C-helix rotates about 45 allowing the interaction between Glu51 and Arg149, then T-loop has to refold to assume the closed conformation. We found out the key role of Arg149 in the closure movement as the link residue between the N-lobe and the C-lobe, allowing a concerted movement of the two lobes. We do not find evidence for the DFG out flip. A complementary method, Metadynamics,[2,3] was used to study the undocking path of a congeneric series of ligands to CDK2. Also in this case the large scale dynamics of the kinase plays a fundamental role. 1. D. Branduardi, F. L. Gervasio and M. Parrinello. Rare events in molecular systems: From A to B in Free Energy Space J. Chem. Phys., 126, 054103, 2007. 2. F. L. Gervasio*, A. Laio and M. Parrinello. Flexible Docking in Solution Using Metadynamics J. Am. Chem. Soc., 127(8): 2600-2607, 2005. 3. A. Laio and F. L. Gervasio* Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science Rep. Prog. Phys. , 71, 126601 (22p), 2008. 4. A. Berteotti, A. Cavalli, D. Branduardi, F. L. Gervasio*, M. Recanatini and M. Parrinello, Protein Conformational Transitions: The Closure Mechanism of a Kinase Explored by Atomistic Simulations nism of CDK5. J. Am. Chem. Soc.,131, 244-250, 2009.
Scott Habershon (University of Oxford) Competing quantum effects in liquid water
Abstract: Numerous simulation studies comparing the classical and quantum molecular dynamics of simple empirical liquid water models have suggested that quantum mechanical zero point energy and tunnelling increase the rates of translational and orientational dynamics by a factor of about 1.4. In this work, we suggest that the magnitude of this effect has been overestimated, principally as a result of the use of either rigid water models or models in which intramolecular flexibility is described by simple harmonic functions. We have developed a new simple point charge model for liquid water, q-TIP4P/F, in which O-H stretches are described by Morse-type functions. This model correctly describes the liquid structure, translational diffusion and dipole absorption frequencies in quantum (path integral-based) simulations, and also gives a good description of the temperature-dependence of the liquid density. By comparing classical and quantum simulations, we find that quantum effects increase the rates of translation diffusion and orientational relaxation by a factor of 1.15, an effect which is smaller than observed in all previous simulations, regardless of quantum simulation method or water model. We attribute the small quantum effect in our water model to two competing phenomena. First, zero point energy weakens the intermolecular hydrogen-bonding network resulting in a less viscous liquid and faster diffusion. Second, intramolecular zero point energy changes the average water monomer geometry, resulting in a larger molecular dipole moment and stronger intermolecular interactions. By comparison to other simulations of liquid water models, we suggest that the small quantum effect in our model is associated with its ability to reproduce a single broad O-H stretching band in the infra-red dipole absorption spectrum.
Graeme Henkelman (University of Texas) Adaptive kinetic Monte Carlo for simulating dynamics in atomic systems
Abstract: I will present a computational method for simulating the dynamics of atomic systems on time scales much longer than can be accessed with classical dynamics. Possible reaction mechanisms available to the system are found by exploring the potential energy surface from minima to find nearby saddle points. Reaction rates are then calculated using harmonic transition state theory, and the system is propagated in time according to the kinetic Monte Carlo algorithm. The method can be run in parallel, and a distributed computing system has allowed us to simulate dynamics on metal surfaces and in grain boundaries over experimental time scales using an embedded atom potential. The algorithm is efficient enough to model the evolution of systems with ab-initio forces as well, for which I will show a few examples, including metal cluster formation on oxides and catalytic reactions on metal surfaces.
Ivaylo N. Ivanov (University of California, San Diego) The interplay of AAA+ molecular machines and sliding clamps at the DNA replication fork
Abstract: Replication of chromosomal DNA is accomplished by a dynamic protein assembly termed the replisome. Within the replisome, the function of replicative DNA polymerases is to rapidly and faithfully duplicate the cell’s genetic material prior to cell division. During replication proliferating cell nuclear antigen (PCNA, sliding clamp) serves as an accessory protein whose role is to topologically tether DNA polymerase to DNA. PCNA also acts as a scaffold in the recruitment of proteins involved in cell-cycle control and DNA repair. Mounting recent experimental evidence has pointed to the involvement of sliding clamps in almost every aspect of DNA metabolism. Yet many mechanistic details of how these fascinating proteins are loaded onto DNA and function within the machinery of the replisome have remained unknown. To perform its vital functions the clamp has to be opened and resealed at DNA primer-template junctions by the action of a clamp loader ATPase – replication factor C (RFC). Recent computational work on the RFC/PCNA complex will be presented. It established that upon clamp opening the complex undergoes a large conformational rearrangement, leading to the formation of an extended interface between the surface of the sliding clamp and RFC. An interface complementary to ring-open PCNA transforms the free energy landscape underlying the closed- to open state transition, effectively trapping PCNA in an open conformation. Thus, careful comparison of free energy profiles for clamp opening in the presence and absence of RFC has allowed us to substantiate the role of the clamp loader in the initial stage of the clamp-loading cycle.
Gerhard Kahl (Technische Universität Wien) Exotic self-assembly scenarios in in two-dimensional dipolar mixtures with genetic algorithms
Abstract: Joint work with Julia Fornleitner, Federica LoVerso, and Christos N. Likos. Genetic algorithms (GAs) are powerful tools, applied in many fields of physics as they offer an efficient strategy in high-dimensional optimisation problems. We have developed GA-techniques that allow us to study equilibrium crystal structures in an unbiased and parameter-free way. Recently we have refined our approach and extended it to binary mixtures of dipolar particles in a two-dimensional geometry, realisable in experiments via, e.g. magnetic colloids on a pendant water droplet and exposed to an external field or, equivalently, via polystyrene particles floating at an oil-water-interface. Several non-trivial and exotic structures are discovered and their dependence on concentration and size ratio of the particle species is discussed.
Gerhard Kahl (Technische Universität Wien) Hard-core colloids forming stable crystals: Clusters - columns - lamellae - compact
Abstract: Joint work with Gernot J. Pauschenwein. The success of experiments to design colloidal particles that interact like classical hard spheres reintroduced lively work on systems interacting via a hard core and various attached potential tails. We investigated the equilibrium crystal structures at zero temperature of hard--core systems interacting additionally either via soft shoulder (polymer--grafted colloids) or Yukawa interactions (equally charged colloids in salty solvent). Our investigations are based on optimisation strategies that use ideas of genetic algorithms. The problem of excluding a priori structures with overlapping cores from search space was solved by introducing a particular parametrisation of all crystal lattices. This parametrisation drastically improved the performance of the genetic algorithm. The application of this approach to the soft shoulder potential revealed a large number of equilibrium structures in dependence of ambient pressure, ranging from clusters over columns and lamellae to finally compact structures. The close--packed structure with the lowest energy per particle can, in dependence of the shoulder width, be one with a different stacking of hexagonal layers than both fcc (ABC) or hcp (AB), as the genetic algorithm revealed and was proven theoretically. For hard core Yukawa systems the genetic algorithm gives, as expected, fcc or bcc depending on pressure and interaction range over the vast majority of parameter space. However our investigation gives evidence that the phase transformation from bcc to fcc at high pressure is happening through a continuous Bain transformation, giving rise to a small regime of centered tetragonal symmetry.
Hiqmet Kamberaj (University of Minnesota), Yuk Sham (University of Minnesota) Computational alanine scanning of the alpha-lactamase inhibitor protein and TEM-1 alpha-lactamase complex
Abstract: Protein association plays an essential role in many biological processes. Understanding the general principles governing the protein-protein interactions at the molecular level is of great importance. In this study, we present a computational based alanine scanning method that is used to elucidate the interactions between alpha-lactamase and the alpha-lactamase inhibitor protein (BLIP). The advantage of such approach is that it provides a direct quantification of the specific per residue interactions between the two proteins which can be directly validated with experimental studies. Such analysis provides insight behind the structural and energetic basis behind the molecular recognition process. Here, we show the association between BLIP and TEM-1 is governed by clusters of "hot spots" residues, which co-operatively modulate the formation of the BLIP-TEM-1 complex. We also show that mutations in each cluster do not influence the binding energy contribution of the residues comprising other "hot spots" clusters, in agreement with experimental results. The method presented here have broad range application to study effect of mutation in drug binding as well as de-novel protein design.
Robert Krasny (University of Michigan) A Cartesian treecode for screened Coulomb interactions
Abstract: A treecode algorithm is presented for evaluating the electrostatic potential in charged particle systems undergoing screened Coulomb interactions in 3D. The method uses a far-field Taylor expansion in Cartesian coordinates to compute particle-cluster interactions. The Taylor coefficients are evaluated using new recurrence relations which permit efficient computation of high order approximations. Two types of clusters are considered, uniform cubes and adapted rectangular boxes. The treecode error, CPU time, and memory usage are reported and compared with direct summation for randomly distributed particles inside a cube, on the surface of a sphere, and on an 8-sphere configuration. For a given order of Taylor approximation, the treecode CPU time scales as O(NlogN) and the memory usage scales as O(N), where N is the number of particles in the system. Results show that the treecode is well suited for non-homogeneous particle distributions as in the sphere and 8-sphere test cases. This is joint work with Peijun Li (Purdue) and Hans Johnston (U. Mass - Amherst), and has appeared in J. Comput. Phys., vol. 228 (2009) 3858-3868.
Robert Krasny (University of Michigan) Special topic 3. Computation of force fields for molecular simulation
Abstract: No Abstract
Krzysztof Kuczera (University of Kansas) Analysis of folding pathways in model peptides
Abstract: Molecular dynamics and replica-exchange simulations are performed for a model peptide, NAcetyl-ALA5-Amide, using several popular force fields. We obtain helix populations, melting curves, folding and nucleation times. The distributions of conformer populations are used to measure folding cooperativity. Finally, a statistical analysis of the sample of helix-coil transition paths is performed, uncovering interesting microscopic properties of the folding process.
Juan C. Latorre (Freie Universität Berlin) Membrane channels: Diffusion in confined geometries and entropic barriers
Abstract: In this work we present a mathematical approach to the study of diffusing particles confined to narrow geometries, such as water and ion channels in cell membranes. In this model, a strong potential (such as an electric field) drives the particle to remain close to the center of the channel, so that the dynamics in the direction parallel to the channel may be described by an averaged equation with an effective Fixman potential driving the motion of the particle. The Fixman potential depends on the particular geometry of the channel, accounting for entropic barriers (bottlenecks) whose effective heights are proportional to temperature. We compare transport properties such as flux rate and effective diffusivity when the typical barrier is either entropic or energetic.
Kah Chun Lau (George Washington University) Kinetic Monte Carlo simulation of the Yttria Stabilized Zirconia (YSZ) fuel cell
Abstract: A Kinetic Monte Carlo (KMC) model is developed to simulate non-symmetrically the cathode side of a Yttria Stabilized Zirconia (YSZ) fuel cell, in order to translate experimental, and ultimately theoretical rates into an atomistic model of the fuel cell. The KMC model consists of a set of several electrochemical reaction rates, adopted from experiments and first-principles calculations. The KMC simulations are used to model these simultaneously occurring events, and to determine potential limitations in cathode/YSZ performance. The focus of this work is ionic current density (J), studied as a function of various physical parameters: oxygen partial pressure, external applied bias voltage, operating temperature, dopant concentration, relative permittivity of YSZ, and geometrical features of the YSZ electrolyte. This simple model can be used as a baseline to translate elementary electrochemical reaction rates into atomistic simulations of a working solid oxide fuel cell cathodes, pertinent to the complete set of experimental operating conditions.
Frédéric Legoll (École Nationale des Ponts-et-Chaussées) Non-ergodicity of the Nosé-Hoover dynamics
Abstract: The Nosé-Hoover dynamics is a deterministic method that is commonly used to sample the canonical Gibbs measure. This dynamics extends the physical Hamiltonian dynamics by the addition of a thermostat variable, that is coupled nonlinearly with the physical variables. The validity of the method depends on the dynamics being ergodic. It has been numerically observed for a long time that such a thermostat, applied to some model problems (including the one-dimensional harmonic oscillator), is actually not ergodic. In this work, we first show that, for some multidimensional systems, the averaged dynamics, obtained in the limit of infinite thermostat mass, has many invariants, thus giving theoretical support for either non-ergodicity or slow ergodization. Next, in the case of one-dimensional Hamiltonian systems, we go further and prove non-ergodicity of the thermostat for large (but finite) thermostat masses. Numerical experiments will illustrate the theoretical results.
Frédéric Legoll (École Nationale des Ponts-et-Chaussées) Special Topic 2: Multiple time-scale methods
Abstract: No Abstract
Frédéric Legoll (École Nationale des Ponts-et-Chaussées) Effective dynamics for reaction coordinates
Abstract: Computing canonical averages is a standard task in molecular dynamics. The main difficulty comes from the existence of metastable states. In order to facilitate sampling, and to obtain a better understanding of the system, low-dimensional reaction coordinates are often introduced. The statistics of these reaction coordinates are completely described by the associated free energy. However, the link between free energy barriers and dynamical informations, such as transition rates from one well to another, is less clear. In this work, we design an effective, low-dimensional, dynamics on the reaction coordinate, which is a good approximation of the reference evolution, computed from the dynamics of the complete system. The accuracy of the effective dynamics is supported by error estimators. We also show, on a simple system, that our dynamics correctly reproduces the residence times in the wells.
Benedict Leimkuhler (University of Edinburgh) Flexible control of the ensemble
Abstract: For molecular dynamics to be useful for recovering macroscopically-relevant information (e.g. thermodynamic averages), it must be regulated by auxiliary control laws. Typically these would represent thermal or pressure controls. In this talk I will discuss a general family of stochastic-dynamic methods, including gentle thermostats which mildly perturb averaged dynamics and some very aggressive thermostats intended to enhance sampling of the canonical ensemble.
Benedict Leimkuhler (University of Edinburgh) Lecture 2: Molecular dynamics: foundations I.
Dynamical issues, ergodicity, and averages at constant energy.
Abstract: No Abstract
Benedict Leimkuhler (University of Edinburgh) Lecture 3: Molecular dynamics: foundations II.
Geometric integrators and applications to molecular dynamics.
Abstract: No Abstract
Benedict Leimkuhler (University of Edinburgh) Lecture 5: Control of the ensemble.
Controlling temperature and pressure in molecular dynamics.
Abstract: No Abstract
Tony Lelievre (Ecole Nationale des Ponts et Chaussees) Adaptive methods for efficient sampling. Applications in molecular dynamics
Abstract: One aim of molecular dynamics simulations is to sample Boltzmann-Gibbs measures associated to some potentials in high dimensional spaces, to compute macroscopic quantities (such as chemical reaction constants, or diffusions constants) by statistical means in the canonical (NVT) ensemble. Numerical methods typically rely on ergodic limits for processes solution to well-chosen stochastic differential equations (SDEs). The main difficulty comes from existence of metastable states in which the stochastic processes remain for long time: this may slow down dramatically the convergence of the ergodic limit. We present a class of adaptive importance sampling methods which enable fast exploration of the configurational space, by modifying the potential seen by the particles (the associated SDE becomes non-homogeneous and nonlinear). These methods accelerate the longtime convergence while they allow to obtain, in the longtime limit, the quantities of practical interest (the marginal law associated to the slow variables in the system). We propose a proof of convergence for some of these methods, based on entropy techniques. References: - T. Lelievre, M. Rousset and G. Stoltz, Computation of free energy profiles with parallel adaptive dynamics, Journal of Chemical Physics 126, 134111 (2007) - T. Lelievre, M. Rousset and G. Stoltz, Long-time convergence of an Adaptive Biasing Force method, Nonlinearity, 21, 1155-1181 (2008) - T. Lelievre, A general two-scale criteria for logarithmic Sobolev inequalities, to appear in Journal of Functional Analysis
Tony Lelievre (Ecole Nationale des Ponts et Chaussees) Special Topic 4. Adaptive methods for free energy computation
Abstract: No Abstract
Xiantao Li (Pennsylvania State University) Coarse-graining molecular dynamics for crystalline solids
Abstract: I will discuss a coarse-grained model of molecular dynamics in crystalline solids. Such system usually has an underlying lattice structure, and the empirical potentials are of quite simple forms. Thus it provides a nice setting in which coarse-graining (CG) methods can be considered. The goal of the CG method is to reduce the atomic degrees of freedom, and arrive at an effective model in which only a small number of atoms are explicitly involved.
Hai Lin (University of Colorado) Flexible-boundary QM/MM: Partial charge transfer between the QM and MM subsystems
Abstract: A prominent limitation of the commonly-used QM/MM methods is that no partial transfer is allowed between the QM and MM subsystems. A challenging example is the Eigen cation (Fig. 1), where one models the inner moiety H3O+ by QM and the surrounding three H2O molecules by MM. Both mutual polarization and partial charge transfer are expected between the QM and MM subsystems. This work aims to go beyond the above limit by incorporating partial charge transfer between the QM and MM subsystems, which (we hope) could lead to more accurate description for the QM/MM electrostatic interactions.
Yuchun Lin (University of North Carolina - Charlotte) An image-based reaction field method for electrostatic interactions in molecular dynamics simulations of aqueous solutions
Abstract: Joint work with Andrij Baumketnera, Wei Caib, Shaozhong Dengb, Donald Jacobsa, and Zhenli Xub. A new solvation model is proposed for simulations of biomolecules in aqueous solutions that combines the strengths of explicit and implicit solvent representations. Solute molecules are placed in spherical cavities filled with explicit water, thus providing microscopic detail where it is needed. Solvent outside of the cavities is replaced with a dielectric continuum whose effect on the solute is modeled through the reaction field corrections. With this explicit/implicit model, the electrostatic potential represents a solute molecule in an infinite bath of solvent, thus avoiding unphysical interactions between periodic images of the solute commonly used in the lattice-sum explicit solvent simulations. For improved computational efficiency, our model employs an accurate and efficient multiple-image method to compute reaction fields together with the fast multipole-expansion technique for the direct Coulomb interactions. To minimize the surface effects, periodic boundary conditions are employed for non-electrostatic interactions. The proposed model is applied to study liquid water. The effect of main geometric parameters of the model, which include the size of the cavity, the number of charge images used to compute reaction field and the thickness of the buffer layer, is investigated in comparison with the particle-mesh Ewald simulations as a reference. An optimal set of parameters is obtained that allows for a faithful representation of many structural, dielectric and dynamic properties of the simulated water, while maintaining manageable a computational cost. With predicable increasing accuracy of the multiple-image charge representation of the reaction field, it is concluded that the proposed model achieves convergence with only one image charge in the case of pure water. Future applications to pKa calculations, conformational sampling of solvated biomolecules and electrolyte solutions are briefly discussed.

a: Department of Physics and Optical Science, University of North Carolina at Charlotte, Charlotte, NC 28223, United States b: Department of Mathematics and Statistics, University of North Carolina at Charlotte, Charlotte, NC 28223, United States

Roland Lindh (Lund University) Presentation of the workshop
Abstract: No Abstract
Roland Lindh (Lund University) Lecture 2: Basis sets, CD/RI
Abstract: No Abstract
Roland Lindh (Lund University) Lecture 4: Multiconfigurational methods: CASSCF/CASPT2
Abstract: No Abstract
Roland Lindh (Lund University) Lecture 5: Geometry optimizations with and without constraints
Abstract: No Abstract
Roland Lindh (Lund University) Lecture 9: QM/MM in Molcas
Abstract: No Abstract
Roland Lindh (Lund University) Conclusions
Abstract: No Abstract
Mitchell Luskin (University of Minnesota) Special Topic 1: Quasicontinuum methods for coupling atomistic and continuum models
Abstract: No Abstract
Dmitrii E. Makarov (University of Texas) Dynamics and thermodynamics of large conformational rearrangements in confined polymers
Abstract: In this talk I will discuss some of the problems arising in connection with the dynamics of polymers under confinement (such as proteins traversing transmembrane pores). Confinement often partitions a polymer's configuration space into distinct basins of attraction separated by large free energy barriers. Because of the predominantly entropic nature of such barriers, transitions between different basins often cannot be characterized in terms of dominant pathways. I will particularly focus on the problem of a polymer reversal inside a long pore and compare exact rate calculations with predictions from polymer scaling theory and with various approximate theories.
Thérèse E. Malliavin (Centre National de la Recherche Scientifique (CNRS)) The impact of calcium on the dynamics and energetics of the complex between calmodulin and anthrax edema factor
Abstract: Joint work with E. Laine (1), J.D. Yoneda (2), and A. Blondel (1). Among the toxins secreted by Bacillus anthracis the edema factor EF, an adenylate cyclase, provokes severe cellular dysfunction by accumulating cAMP from ATP. EF is activated by calmodulin (CaM), involved in many calcium signaling pathways. The stability of the EF-CaM complex depends on the level of calcium bound to CaM while the architecture of the complex loaded with 2, 3 or 4 Ca2+ ions remains practically unchanged. That is why modeling the electrostatic effect of Calcium through EF- CaM structure is challenging. Here, we aim at describing the calcium-induced changes in EF-CaM dynamics and energetics through a consensual view of its residue network organization. The analysis of molecular dynamics (MD) simulations of EF-CaM with 0,2 and 4 Ca2+ ions helped characterize CaM conformational plasticity and led to a model of the EF- CaM interaction, in which CaM acts as a spring that maintains EF in an open active conformation (Laine et al., 2008). The complex was then analyzed (Laine et al.,2009) in order to extract simplified features describing its energetics and dynamics. The dynamics were analyzed using the LFA (Local Feature Analysis) (Zhang and Wrigers, 2006) and the generalized correlations (Lange and Grubmuller, 2006). These two approaches permitted to determine communication paths along the network of EF residues, from the helical domain to the domain CB. Beside, a new definition of the complex domains arose from the LFA analysis. This new definition, along with an approach derived from the energetic dependency maps (Hamacher et al, 2006), permitted to draw a simplified model of the energetic influences inside the complex. The variations of dynamics and energetics according to the level of Calcium complexation provide phenomenological reasons for the stability of the complex EF-CaM, by showing a disruption of the communication path or an imbalance of the energetic influences. The computation of various dynamical covariances and energetic dependency maps from the MD trajectories further raised the concept of residue network connectedness. This connectedness quality provides a frame for unifying the dynamics and energetics of the complex and a criterion for assessing its stability. . Hamacher, K., J. Trylska, and J. McCammon. 2006. Dependency map of proteins in the small ribosomal subunit. PLoS Comput. Biol. 2:e10. Laine E, Blondel A, Malliavin TE (2009) Dynamics and energetics: a consensus analysis of the impact of calcium on EF-CaM protein complex. Biophys J 96:1249-1263. Laine E, Yoneda JD, Blondel A, Malliavin TE (2008) The conformational plasticity of calmodulin upon calcium complexation gives a model of its interaction with the oedema factor of Bacillus anthracis. Proteins 71:1813-1829. Lange OF, Grubmuller, H (2006) Generalized correlation for biomolecular dynamics. Proteins 62:1053-1061. Zhang Z, Wriggers, W (2006) Local feature analysis: a statistical theory for reproducible essential dynamics of large macromolecules. Proteins 64:391-403. (1) Unite de Bioinformatique Structurale, Institut Pasteur, 28, rue du Dr. Roux, F-75015 Paris, France (2) Programa de Posgraduac¸ao em Quimica Organica, Universidade Federal Fluminense, Outeiro de S. Joao Batista s/n, 24020150, Niteroi, RJ, Brazil.
Per-Åke Malmqvist (Lund University), Valera Veryazov (Lund University) Lecture 3: Use of symmetry
HF, DFT
Abstract: No Abstract
Per-Åke Malmqvist (Lund University) Lecture 6: Relativistic calculations
Abstract: No Abstract
Thomas Markland (University of Oxford) Fast path integral simulations using ring polymer contraction
Abstract: A quantum simulation of an imaginary time path integral typically requires around n times more computational effort than the corresponding classical simulation, where n is the number of ring polymer beads (or imaginary time slices) used in the calculation. However, this estimate neglects the fact that the potential energies of many systems can be decomposed into a sum of rapidly-varying short-range and slowly-varying long-range contributions. For such systems, the computational effort of the path integral simulation can be reduced considerably by evaluating the long-range forces on a contracted ring polymer with fewer beads than are needed to evaluate the short-range forces. Using this approach it is shown that, for a typical application to a flexible model of liquid water, near classical simulation times can be achieved while still obtaining the exact path integral result.
David A. Micha (University of Florida) Density matrix treatment of optical response with combined instantaneous and delayed dissipation: Adsorbates on solid surfaces
Abstract: A density operator treatment has been developed for the dissipative quantum dynamics of molecular systems in condensed matter. It incorporates instantaneous and delayed rates of dissipation originating in electronic and vibrational motions of the molecular environment, respectively. The theory is applied here to the optical response of the adsorbates CO/Cu(001) and Ag/Si(111).
David Do Le Minh (National Institutes of Health) Density-dependent analysis of nonequilibrium paths improves free energy estimates
Abstract: When a system is driven out of equilibrium by a time-dependent protocol that modifies the Hamiltonian, it follows a nonequilibrium path. Samples of these paths can be used in nonequilibrium work theorems to estimate equilibrium quantities, such as free energy differences. Here, we consider analyzing paths generated with one protocol using another one. It is posited that analysis protocols which minimize the lag, the difference between the nonequilibrium and the instantaneous equilibrium densities, will reduce the dissipation of reprocessed trajectories and lead to better free energy estimates. Indeed, when minimal lag analysis protocols based on exactly soluble propagators or relative entropies are applied to several test cases, substantial gains in the accuracy and precision of estimated free energy differences are observed.
Kimiya Minoukadeh (Ecole Nationale des Ponts et Chaussees) Some improvements of the ART method for finding transition pathways on potential energy surfaces
Abstract: The Activation-Relaxation Technique nouveau (ARTn) is an eigenvector following method for systematic search of saddle points and transition pathways on a given potential energy surface. We propose a variation of this method aiming at improving the efficiency of the local convergence close to the saddle point. We prove the convergence and robustness of this new algorithm. The efficiency of the method is tested in the case of point defects in body centered cubic iron.
Peter A. Monson (University of Massachusetts) Simulations of phase transitions for fluids in pores using dynamic mean field theory
Abstract: We consider the mean field kinetic equations describing the relaxation dynamics of a lattice model of a fluid confined in a porous material. The dynamical theory embodied in these equations can be viewed as a mean field approximation to a Kawasaki dynamics Monte Carlo simulation of the system, as a theory of diffusion, or as a dynamical density functional theory. The solutions of the kinetic equations for long times coincide with the solutions of the static mean field equations for the inhomogeneous lattice gas. The approach is applied to a lattice gas model of a fluid confined in a finite length slit pore open at both ends and in contact with the bulk fluid at a temperature where capillary condensation and hysteresis occur. The states emerging dynamically during irreversible changes in the chemical potential are compared with those obtained from the static mean field equations for states associated with a quasistatic progression up and down the adsorption/desorption isotherm. In the capillary transition region the dynamics involves the appearance of undulates (adsorption) and liquid bridges (adsorption and desorption) which are unstable in the static mean field theory in the grand ensemble for the open pore but which are stable in the static mean field theory in the canonical ensemble for an infinite pore.
Emad Noorizadeh (University of Edinburgh) What is a good thermostat?
Abstract: Molecular dynamics perturbs Hamiltonian dynamics in order to compute correct ensemble averages. On the other hand computation of dynamical averages such as autocorrelation functions requires the molecular dynamics trajectory to be close to a microcanonical dynamics. This suggest that, ideally we would like to have a small growth of the perturbation and a fast rate of convergence to the Boltzmann-Gibbs measure. Here we show how to achieve such a gentle thermostat by combining a deterministic method (e.g. Nose-Hoover) with Langevin dynamics. We use the concept of hypoellipticity of the forward Kolmogorov equation, and Hörmander's condition, to investigate the existence of a unique invariant measure, which implies ergodicity. We compare different methods in numerical experiment. This is a joint work with Ben Leimkuhler (School of Mathematics, University of Edinburgh) and Florian Theil (Mathematics Institute, University of Warwick). Reference:
[1] B. Leimkuler, E. Noorizadeh and F. Theil, J. Stat. Phys. 135, 261-277 (2009)
Albert C. Pan (University of Chicago) Dynamical pathways to inactivation in the KcsA potassium channel
Abstract: Ion channels gate the passage of ions through cell membranes in response to external stimuli. In the case of the archetypal potassium ion channel, KcsA, which opens and closes upon a change in pH, recent experimental results have demonstrated the existence of two physical gates: an intracellular gate and a gate at the selectivity filter. Lowering the pH opens the intracellular gate allowing ions to pass. After the intracellular gate opens, however, the channel can still inactivate by constricting the selectivity filter and impeding the flow of ions even though the bottom gate remains open. We use a combination of advanced simulation techniques such as umbrella sampling, the string method with swarms of trajectories, free energy perturbation theory and Markov State models to probe the detailed mechanism of this physiologically important process.
Matej Praprotnik (National Institute of Chemistry) Concurrent triple-scale simulation of liquid water
Abstract: We present a triple-scale simulation of liquid water by concurrently coupling the atomic, coarse-grained, and continuum descriptions of the liquid. The triple-scale scheme, which is shown to correctly describe the hydrodynamics, successfully sorts out the problem of large molecule insertion in the hybrid particle-continuum simulations of molecular liquids. The presented approach allows for efficient grand-canonical molecular dynamics simulations of open systems involving an explicit solvent, e.g., liquid water.
Ravi Radhakrishnan (University of Pennsylvania) Hybrid field-based models of cell membranes and proteins for signaling applications
Abstract: We discuss theoretical and computational methodologies for quantitatively describing how cell-membrane topologies are actively mediated and manipulated by intracellular protein assemblies. Such scenarios are ubiquitous in intracellular trafficking mechanisms, i.e., active transport mechanisms characterized by vesicle nucleation and budding of the cell membrane orchestrated by protein-interaction networks. We will describe the development and application of the kinetic Monte Carlo time-dependent Ginzburg Landau (KMC-TDGL) algorithm for unified dynamics of curvature-inducing proteins and membranes. We will also describe the surface-evolution method which predicts minimum energy conformations of highly curved axis-symmetric membrane structures and an extension of the TDGL method which relaxes the assumption of axis-symmetry. Together, these methods enable us to quantify thermal effects in membrane mediated events. We will briefly focus on an application to clathrin-dependent endocytosis.
Paul Rejto (Pfizer Global R&D) Computational challenges in cancer therapeutics
Abstract: Recent advances have dramatically advanced our understanding of cancer at the molecular level. In turn, new therapeutic agents that target specific molecular defects in cancer have been developed, though cancer remains a significant health threat. Following an introduction to the molecular biology of cancer, a statistical approach to distinguish driver mutations from passengers based on non-random clustering will be discussed. Next, approaches to pharmaceutical intervention will be reviewed and an integrated approach to link targeted therapeutics with specific patient populations will be shown. Throughout, open questions will be presented with a focus on problems of potential mathematical interest.
Weiqing Ren (New York University) Coupled atomistic-continuum methods for fluids
Abstract: I will present a seamless multiscale method for the study of fluid systems with unknown constitutive relations and/or unknown boundary conditions. The multiscale method aims at capturing the macroscale behavior of the fluid system; this is done by combining the conservation laws for the mass and momentum (the macro model) and the molecular dynamics (the micro model). The macro model provides the necessary constraint for the micro model and the micro model supplies the missing data (the stress or the boundary conditions) needed in the macro model. In the multiscale method, the macro and micro models evolve simultaneously using different time steps, and they exchange data at every step. The micro model uses its own appropriate (micro) time step. The macro model uses a macro time step but runs at a slower pace than required by accuracy and stability considerations in order for the micro model to relax. One of the advantages of this method is that it does not require reinitialization of the micro model at each macro time step or each macro iteration step. In this talk, I will discuss the algorithm of the multiscale method, the error analysis, and its application to complex fluids.
Fadil Santosa (University of Minnesota) Welcome to the IMA
Abstract: No Abstract
Carlos F. Sanz-Navarro (Autonomous University of Barcelona) Molecular dynamics simulations of binding of metal catalyst nanoparticles to carbon nanostructures by using ReaxFF
Abstract: Carbon nanofibers have attractive properties as support materials for catalyst nanoparticles. Nonetheless, there is still a need for more insights into the actual atomic structure of a metal particle interacting with carbon nanostructures. Detailed information on the local structure of adhered nanoclusters is still inaccessible to direct experimental imaging. In contrast, molecular dynamics (MD) simulations can provide atomistic structural information of a nanoscale system. In addition, the Reax force field (ReaxFF) aims to achieve a good compromise between accuracy and computational efficiency. We present novel results of MD simulations of platinum and nickel clusters adsorbed on several different carbon supports by using the ReaxFF. This work seeks to give insights into recent experiments that have found a high correlation between carbon nanofiber-induced strain and catalytic activity of Ni clusters. We show how some metal atoms are detached from the surface of the nanoclusters over time scales unreachable by electronic-structure methods. The adatom migration is compensated by a rearrangement of the atomic structure of the cluster. We also show differences in the bond length distribution of the cluster when adsorbed to carbon nanocones with different radii as well as differences found between Pt and Ni clusters. In addition, we brief on preliminary results on diffusion of oxygen, hydrogen and carbon atoms over the adsorbed cluster.
Srikanth Sastry (Jawaharlal Nehru Centre for Advanced Scientific Research(JNCASR)) Liquid-liquid critical point in supercooled silicon
Abstract: Liquid silicon has previously been shown to exhibit a liquid-liquid transition in the supercooled state at zero pressure, using computer simulations employing the Stillinger-Weber (SW) potential. The associated liquid-liquid (LL) critical point lies at negative pressures. Computer simulation evidence of such a negative pressure critical point is presented, along with characterization of structural and dynamical properties of the liquid in the vicinity of the critical point. Preliminary results are presented on crystal nucleation, which near the liquid-liquid critical point will be governed by the interplay of critical density fluctuations and nearly tetrahedral, crystal-like, local order in the liquid that grows with decreasing density.
Christof Schütte (Freie Universität Berlin) Markov state models for biomolecular dynamics: Theory and applications
Abstract: The talk will center around "conformation dynamics" and "rare events."
Christof Schütte (Freie Universität Berlin), Eric Vanden-Eijnden Optimal fuzzy aggregation of Markov chains
Abstract: Joint work with Marco Sarich (Freie Universität Berlin). We map a discrete-time Markov Chain onto another Markov Chain describing transitions between some clusters only. Every such aggregated Markov chain and affiliation function can be lifted again onto the full state space to define the so-called lifted transition matrix. The optimal aggregated Markov chain and affiliation function can then be determined by minimizing some appropriately defined distance between the lifted transition matrix and the transition matrix of the original chain.
Luis Serrano-Andrés (University of Valencia) Lecture 8: Excited states
Abstract: No Abstract
Luis Serrano-Andrés (University of Valencia) Lecture 10: Solvent models
Abstract: No Abstract
Luis Serrano-Andrés (University of Valencia) Lecture 11: Handling problems in MOLCAS (convergence, orbital selection, etc.)
Abstract: No Abstract
Juanfang Shen (Purdue University), Robert D. Skeel (Purdue University), Ruijun Zhao (Purdue University) Maximum flux transition path
Abstract: Given two metastable states A and B of a biomolecular system, the problem is to compute a representative transition path by which the mechanism of reaction can be analyzed. It is often necessary to find such a path in collective variable space for large biomolecular system. The maximum flux transition path (MFTP) is defined as a path in collective variable space that crosses each isocommittor at a point which (locally) has the highest crossing rate of distinct reactive trajectories. (Here the committor is defined to be the probability that a trajectory at that point will reach B before A in collective variable space.) Such an algorithm and computational performance will be discussed.
Robert D. Skeel (Purdue University) Lecture 6: Free energy calculation.
Methods for coarse-graining a complex energy surface, reaction coordinates
Abstract: No Abstract
Adnan Sljoka (York University) Counting for rigidity, flexibility and extensions via the pebble game algorithm – hinge predictions and allostery
Abstract: In rigidity theory, specifically combinatorial rigidity, one can simply count vertices and edges (constraints) in a graph and its subgraphs to determine the rigidity and flexibility of a corresponding framework. The 6|V|–6 counting condition for 3-dimensional body-hinge structures (modulo molecular conjecture), and a fast `pebble game' algorithm which tracks the underlying count in the multigraph, have led to the development of the program FIRST, for rapid analysis of the rigidity and flexibility of proteins and other large macromolecular structures such as viral capsids. An important recent application of FIRST is the introduction of the program FRODA which is a Monte Carlo based geometric simulation technique, which uses the predictions of rigid and flexible regions from FIRST to simulate the motions of proteins. FRODA can run up to a million times faster than the Molecular Dynamics Simulations. We will give a brief description and illustration of the 6|V| – 6 pebble game algorithm. We further extend the pebble game algorithm to quantify the relative degrees of freedom of a specified region (core) in the multigraph and identify the regions that are relevant as constraints with respect to the core and those regions that are irrelevant and that can be ignored both for focus and for speed of simulations. With these new extensions (algorithms) we have developed a novel approach in detecting hinge motions between protein domains. We can accurately locate hinge points using only a single protein structure. Other application of this work is to detect regions required for allostery; and the rigidity impact of mutation studies. We present examples of these new extensions along with biological applications.
Berend Smit (University of California, Berkeley) Lecture 1: Introduction to molecular simulation.
An overview of molecular simulation and its application.
Abstract: No Abstract
Berend Smit (University of California, Berkeley) Lecture 4: Basics of statistical mechanics.
Description of statistical mechanical ensembles and implications for simulation.
Abstract: No Abstract
Berend Smit (University of California, Berkeley) Lecture 7: Monte Carlo methods.
Methods for sampling based on Monte Carlo and variants.
Abstract: No Abstract
Berend Smit (University of California, Berkeley) Lecture 8: Methods for polymer models.
Coarse grained models for mesoscale simulation, and methodology such as dissipative particle dynamics.
Abstract: No Abstract
Daniel B. Smith (University of Pittsburgh) Statistical behavior of the self-guided dynamics algorithms
Abstract: The self-guided dynamics algorithms are a pair of techniques to increase the speed of sampling from the possible conformations for a macromolecule. Wu and Wang (1998), and later Wu and Brooks (2003), proposed a numerical scheme that introduced a biasing term based on a running average of past states of the system into the equations of motion to increase the rate of barrier-crossings. In the present work, an analytical form for the guiding force is derived from the numerical methods. That analytical form is used to study the sampling properties of the two algorithms. It is shown that, up to a separability assumption, the systems do not sample from the Boltzmann distribution, and the magnitudes of the error are estimated. Finally, a new thermostat for the self-guided algorithms that conserves the Boltzmann distribution is proposed.
Gabriel Stoltz (École Nationale des Ponts-et-Chaussées (ENPC)) Adiabatic switching for degenerate ground states
Abstract: The Gell-Mann and Low switching allows to transform eigenstates of an unperturbed Hamiltonian H0 into eigenstates of the modified Hamiltonian H0 + V. This switching can be performed when the initial eigenstate is not degenerate, under some gap conditions with the remainder of the spectrum. We show here how to extend this approach to the case when the ground state of the unperturbed Hamiltonian is degenerate. More precisely, we prove that the switching procedure can still be performed when the initial states are eigenstates of the finite rank self-adjoint operator cP0 V cP0, where cP0 is the projection onto the degenerate eigenspace of H0.
Gabriel Stoltz (École Nationale des Ponts-et-Chaussées (ENPC)) Adaptive importance sampling strategies
Abstract: I study from a mathematical viewpoint a nonlinear stochastic dynamics (called the Adaptive Biasing Force method), which allows to adaptively bias a metastable dynamics by computing the free energy associated with a chosen slow direction. Some applications in computational statistical physics as well as Bayesian statistics are presented.
Devarajan Thirumalai (University of Maryland) Folding landscape of RNA from single molecule trajectories
Abstract: The sequence dependent folding landscapes of nucleic acid hairpins reflect much of the complexity of biomolecular folding. Recently, mechanical folding trajectories, generated using single molecule force clamp experiments by attaching semiflexible polymers to the ends of hairpins have been used to infer their folding landscapes. Using simulations and theory, we study the effect of the dynamics of the attached handles on the RNA free energy profile F(zm), where zm is the molecular extension of the hairpin. Accurate measurements of F(zm) requires stiff polymers with small L/lp, where L is the contour length of the handle, and lp is the persistence length. Paradoxically, reliable estimates of hopping rates can only be made using flexible handles. We show that F(zm) at an external tension fm, the force (f) at which the folded and unfolded states are equally populated, in conjunction with Kramers' theory, can provide accurate estimates of the force-dependent hopping rates in the absence of handles at arbitrary values of f. Our theoretical framework shows that zm is a good reaction coordinate for nucleic acid hairpins under tension.
Erkan Tüzel (University of Minnesota) Motor mediated deformation of microtubules in living cells and gliding assays
Abstract: Microtubules are often viewed as mechanically rigid compressive struts that help maintain cell shape and aid in the transport of cellular cargo by serving as tracks for the molecular motors. However, individual microtubules in living cells are often highly bent, and even though the different force mechanisms that contribute to microtubule deformations are known, it is not clear how all these mechanisms act together in a given cell type. Recent experiments on LLC-PK1 epithelial cells strongly suggest that F-actin dynamics and polymerization of microtubules play a minor role, and that the dominant mechanism is anterograde transport via molecular motors. Surprisingly, quantitative analysis of these deformations using curvature distributions exhibit striking similarities to in vitro gliding assays. Motivated by these experiments, both in vivo and in vitro, we have modeled the deformation of microtubules under the influence of molecular motor forces by coarse-grained simulations. In the simulations, microtubules are modeled as semi-flexible polymers with rigid bond constraints embedded in a solvent. Molecular motors exert forces on the microtubules, and walk along microtubule tracks according to their known force-velocity relations, bind and unbind stochastically. Simulation results support our experimental findings and further elucidate on the interplay between molecular motors and passive cross- linkers. Our results suggest that molecular motors are not necessarily just cargo carriers, but can play a dynamic role in the deformation and positioning the microtubule array.
Steven M. Valone (Los Alamos National Laboratory) Models of Chemical potential and the concept of atomistic embedding
Abstract: The concept of embedding some fragment in a larger system or reservoir appears everywhere in our scientific deliberations. Examples cover such diverse situations as an atom in a crystal, a functional group in a molecule, a molecular wire between a source and sink, or a protein in aqueous solution. Virtually every atomistic potential energy surface is built on one embedding strategy or another. This situation presents a difficulty for systems where charges are being transferred among fragments. Governance of transfer of charge is embodied in the concept of the chemical potential. When an embedded fragment and a reservoir can be readily identified, transfer of electrons between them can be described by an open ensemble. However, when interactions are strong as happens in crystals and molecules, the ensemble picture is not strictly applicable. Nevertheless, we still think in terms of open systems. Customary atomistic models use a linear model of chemical potential that is highly problematic. Here one strategy is presented to accommodate the open system point of view, while overcoming several limitations. That strategy focuses on decomposition into fragments of the many-body electronic hamiltonian itself, rather than decomposing the ground-state energy. The strategy produces a new, nonlinear model of chemical potential that has numerous ramifications for designing atomic potentials that can account for charge transfer.
Eric Vanden-Eijnden (New York University) Theory and modeling of reactive events
Abstract: In the first part of the talk, I will explain why we may need to go beyond the standard framework of transition state theory (TST) to describe activated processes and reactive events, and I will present another framework, termed transition path theory (TPT), that permits to do that. Unlike TST, which gives mainly an expression for the rate of the reactive event, TPT describes more fully the statistical properties of the reactive trajectories (i.e. those trajectory by which the reactive event occurs), in particular in terms of their probability density function and their probability current. In the second part of the talk, I will describe how TPT can be use to design and/or improve numerical methods for computing the pathways and rate of reactive events. I will focus in particular on the string method and milestoning.
Carlos Vega de las Heras (Universidad Complutense de Madrid) What ice can teach us about water interactions: a critical comparison of the performance of different water models
Abstract: Coauthors: J. L. F. Abascal, M. M. Conde and J. L. Aragones. The performance of several popular water models (TIP3P, TIP4P, TIP5P and TIP4P/2005) is analysed. For that purpose the predictions for ten different properties of water are investigated, namely: 1. vapour-liquid equilibria (VLE) and critical temperature; 2. surface tension; 3. densities of the different solid structures of water (ices); 4. phase diagram; 5. melting point properties; 6. maximum in density at room pressure and thermal coefficients α and κT; 7. structure of liquid water and ice; 8. equation of state at high pressures; 9. diffusion coefficient; 10. dielectric constant. For each property, the performance of each model is analysed in detail with a critical discussion of the possible reason of the success or failure of the model. A final judgement on the quality of these models is provided. TIP4P/2005 provides the best description of almost all properties of the list, with the only exception of the dielectric constant. In the second position, TIP5P and TIP4P yield an overall similar performance, and the last place with the poorest description of the water properties is provided by TIP3P. The ideas leading to the proposal and design of the TIP4P/2005 are also discussed in detail. TIP4P/2005 is probably close to the best description of water that can be achieved with a non polarizable model described by a single Lennard-Jones (LJ) site and three charges.
Valera Veryazov (Lund University) Lecture 1: Overview of MOLCAS
Abstract: No Abstract
Valera Veryazov (Lund University) Lecture 7: MOLCAS for advanced users
Abstract: No Abstract
Rodolphe Vuilleumier (École Normale Supérieure) Microscopic flow around a diffusing particle
Abstract: The fundamental model for diffusion is the Stokes-Einstein approach. It relates the diffusion coefficient to the friction coefficient which is in turn extracted from the hydrodynamic friction felt by a moving sphere in a liquid environment. The hydrodynamic flow around the sphere, commonly called Stokes flow, depends on the boundary conditions on the sphere. Two boundary conditions are usually considered: stick boundary conditions (no velocity at contact), as in the orginal approach, or slip boundary conditions (finite tangential velocity). As these two boundary conditions lead to different diffusion coefficient as a function of particle size, the usual way to investigate these from numerical simulations has been to compute the diffusion coefficient as a function of particle size. This leads for diffusion of Lennard-Jones particles to the suggestion of slip boundary conditions. Here, we want to present a direct calculation of the microscopic velocity field around a diffusing particle from numerical simulations. This allow for comparison between the atomic model and the hydrodynamics approach. It is first demonstrated that the hydrodynamics flow is well recovered after only a few atomic radius from the tagged particle. However, two effects of the velocity fluctuations of the diffusing particle are evidenced. The boundary conditions are shown to include a finite normal velocity at contact, which would seem at first in contradiction with the non-penetrability of the particles but is an effect of fluctuations. Then, the flux of momentum in the flow is not determined solely by viscosity but also from the diffusion of the tagged particle. This gives new insight in the diffusion process in liquids and in particular on the role of fluctuations at the atomic level. We will also try to show how this method could be also used in order to provide a more mechanistic description of the diffusion processes, in terms similar to reaction paths and transition states.
Rodolphe Vuilleumier (École Normale Supérieure) Projected dynamics from molecular dynamics simulations
Abstract: It is often desirable to reduce the description of an atomic system to the dynamics of only a few relevant, usually, slow variables. Projection methods offer a rigourous framework to perform such reduction of variables. The main outcome of the method is the definition of a memory or friction kernel for the reduced dynamics. The memory kernel is equal to the time correlation of a so-called random force acting on the relevant variables. The evolution of the random force is not governed by the natural microscopic dynamics but by a projected dynamics. We have tried to illustrate this projected dynamics through the design of a method to directly compute the projected correlation functions from Molecular Dynamics simulations. Not only did this allow us to compute the memory kernel for the self-diffusion of Lennard-Jones particles but il also allowed for the decomposition of the memory kernel into short range and long range interactions contributions. Finally, long time behavior of the projected dynamics is investigated.
David John Wales (University of Cambridge) Energy landscapes of clusters, glasses, and biomolecules
Abstract: Coarse-graining the potential energy surface into the basins of attraction of local minima, provides a computational framework for investigating structure, dynamics and thermodynamics in molecular science. Steps between local minima form the basis for global optimisation via basin-hopping and for calculating thermodynamic properties using the superposition approach and basin-sampling. To treat global dynamics we must include transition states of the potential energy surface, which link local minima via steepest-descent paths. We may then apply the discrete path sampling method, which provides access to rate constants for rare events. In large systems the paths between minima with unrelated structures may involve hundreds of stationary points of the potential energy surface. New algorithms have been developed for both geometry optimisation and making connections between distant local minima, which allow us to treat such systems. Applications will be presented for a wide variety of atomic and molecular clusters, glass-formers and biomolecules. Results for folding, misfolding and aggregation of peptides and proteins illustrate how experimental time and length scales can now be addressed for such systems.
Ryoichi Yamamoto (Kyoto University) Rheological properties of polymer melt between rapidly oscillating plates: an application of multiscale modeling
Abstract:

The behavior of supercooled polymer melt composed of short chains with ten beads between rapidly oscillating plates is simulated by using a hybrid simulation of molecular dynamics and computational fluid dynamics. The flow profiles of polymer melt near an oscillating plate are quite different from those of Newtonian fluid. The viscous boundary layer of the melt is much thinner than that of the Newtonian fluid due to the shear thinning of the melt. Three different rheological regimes, i.e., the viscous fluid, viscoelastic liquid, and viscoelastic solid regimes, form over the oscillating plate according to the local Deborah numbers. The melt behaves as a viscous fluid when ftR < 1, and the crossover between the liquid-like and solid-like regime takes place around ftR = 1 (where f is the angular frequency of the plate and tR and ta are the Rouse and the alpha relaxation time, respectively).

Xiaoyang Zhu (University of Minnesota) The Coulomb barrier in excitonic charge separation
Abstract: When a molecular material is electronically excited by a photon, the Coulomb attraction between the excited electron and the hole gives rise to an atomic-H-like quasi-particle called an exciton. The bound electron-hole pair also forms across a material interface, such as the donor/acceptor (D/A) interface in an organic heterojunction solar cell; the results are charge-transfer (CT) excitons. In a conventional p-n junction cell, the exciton binding energy is very small and there is a built-in potential to ensure charge separation. In contrast, there is no a priori a built-in potential in a solar cell based on organic molecules and polymers. Charge separation requires an energetic driving force provided by the differences in the electronic levels at the D/A interface. From typical dielectric constants of organic semiconductors and sizes of conjugated molecules, one can estimate that the Coulomb energy of a CT exciton across a D/A interface is more than one order of magnitude greater than k_B T at room temperature. We can also estimate this by solving the Schrodinger equation in a dielectric continuum model. The results are a series of atomic-H like states with the binding energy of the 1s state on the order of a few hundred meV. How can the e-h pair escape this Coulomb trap? To answer this question, we carry out two experiments: femtosecond time-resolved two-photon photoemission (TR-2PPE) spectroscopy to probe the physical nature of the CT excitons and vibrational spectroscopy to probe charge transfer and electric field at buried interfaces. In the first experiment, we use a crystalline pentacene thin film as a model system and excite an electron above the surface. When there is only the excess electron, the electron is bound to the surface by the positive polarization cloud in the pentacene molecules; this attractive potential leads to a bound electronic state called an image potential state. When a positive hole is present on a pentacene molecule, the electron is attracted to the surface by both the positive hole and the polarization cloud; the result of this combined attractive potential is approximately the CT exciton referenced to the image band. We observe a series of CT excitons with binding energies <=0.5 eV below the image band minimum. Given the large binding energy of the lowest CT exciton state (1s), we conclude that hot CT exciton states must be involved in charge separation. The second experiment probes electric field at the D/A interface. This field may come from direct charge transfer, from unintentional doping, or from charge redistribution or polarization. Experimentally, the presence of charge transfer and interfacial field in a bulk heterojunction can be probed by vibrational spectroscopy and the Stark effect. We apply this approach to a model bulk heterojunctions systems of electron donating polymers and a soluble form of C60. We find that the electric field for at the D/A interface is as large as 4x10^6 V/cm in the wrong direction for charge separation. This removal of this Coulomb barrier occurs upon photoinduced charge separation in the presence of sufficient phase separation and local crystallinity. These findings suggest key design principles for successful organic bulk heterojunction solar cells.
Daniel M. Zuckerman (University of Pittsburgh) Protein path sampling: Bring your algorithm and your model
Abstract: Conformational changes in proteins are intrinsic to almost every function, but sampling such changes remains one of the most difficult bio-computational problems. The challenge lies on two levels. First, algorithms which can efficiently focus computer time on the correct path ensemble are necessary. Second, even if perfect path-sampling algorithms were available, intrinsic "fast" timescales of protein motions are too costly for atomistic models; coarse-grained models are therefore required in most systems of interest. The group's work on both these issues will be described. Algorithmically, the "weighted ensemble" method is generalized and set on a firm statistical footing. Toward the goal of determining the most chemically realistic models which can be fully sampled, we describe a novel class of semi-atomistic models, based on pre-calculated residue libraries.
Visitors in Residence
Carlo Adamo École Nationale Supérieure de Chimie de Paris 4/3/2009 - 6/30/2009
Michael Albright St. Cloud State University 5/15/2009 - 5/16/2009
Michael Philip Allen University of Warwick 5/16/2009 - 5/22/2009
Donald G. Aronson University of Minnesota 9/1/2002 - 8/31/2009
Manuel Athenes Commissariat à l'Energie Atomique (CEA)-Centre d'Études Nucléaires de Grenoble (CEN) 5/16/2009 - 5/22/2009
Katie Beutel University of Minnesota 5/15/2009 - 5/16/2009
Brian Bies Washington University 5/26/2009 - 8/7/2009
Henry A. Boateng University of Michigan 5/14/2009 - 5/23/2009
Peter Bolhuis Universiteit van Amsterdam 5/16/2009 - 5/21/2009
Stephen Bond University of Illinois at Urbana-Champaign 5/14/2009 - 5/23/2009
Nawaf Bou-Rabee New York University 5/17/2009 - 5/22/2009
Bastiaan J. Braams Emory University 5/14/2009 - 5/23/2009
Peter Brune University of Chicago 9/8/2008 - 6/9/2009
Sun Young Bu University of North Carolina 2/1/2009 - 5/30/2009
Maria-Carme T. Calderer University of Minnesota 9/1/2008 - 6/30/2009
Hannah Callender University of Minnesota 9/1/2007 - 8/31/2009
Maria Kourkina Cameron Courant Institute of Mathematical Sciences 5/17/2009 - 5/22/2009
Alessandro Cembran University of Minnesota 5/15/2009 - 5/16/2009
Xianjin Chen University of Minnesota 9/1/2008 - 8/31/2010
Giovanni Ciccotti Università di Roma "La Sapienza" 5/16/2009 - 5/22/2009
Ilaria Ciofini École Nationale Supérieure de Chimie de Paris 4/3/2009 - 6/30/2009
Rosana Collepardo-Guevara University of Oxford 5/17/2009 - 5/22/2009
Grazia Cottone Università di Palermo 5/17/2009 - 5/23/2009
Eric Felix Darve Stanford University 5/17/2009 - 5/21/2009
Andrew Davies St. Cloud State University 5/15/2009 - 5/16/2009
Luigi Delle Site Max-Planck Institut für Polymerforschung 5/6/2009 - 5/30/2009
Ajitha Devarajan Iowa State University 5/4/2009 - 5/8/2009
Daniel Dix University of South Carolina 1/1/2009 - 6/30/2009
Ron Dror D. E. Shaw Research 5/17/2009 - 5/22/2009
Olivier Dubois University of Minnesota 9/3/2007 - 8/31/2009
Allison Lee Dzubak University of Minnesota 5/3/2009 - 5/8/2009
Allison Lee Dzubak University of Minnesota 5/15/2009 - 5/16/2009
Weinan E Princeton University 5/20/2009 - 5/22/2009
Robert S. Eisenberg Rush University Medical Center 5/16/2009 - 5/22/2009
Ron Elber University of Texas 5/17/2009 - 5/22/2009
James W. Evans Iowa State University 3/15/2009 - 5/22/2009
Luke Fiedler University of Minnesota 5/4/2009 - 5/8/2009
Meredith Catherine Foley New York University 5/17/2009 - 5/22/2009
Christopher Fraser University of Chicago 8/27/2008 - 6/30/2009
Daan Frenkel University of Cambridge 5/16/2009 - 5/22/2009
Alexander Gaenko Iowa State University 5/3/2009 - 5/8/2009
Weihua Geng University of Michigan 5/17/2009 - 5/22/2009
Francesco Luigi Gervasio Eidgenössische TH Zürich 5/16/2009 - 5/22/2009
Scott Habershon University of Oxford 5/16/2009 - 5/22/2009
Amir Haji-Akbari University of Michigan 5/14/2009 - 5/23/2009
Jaebeom Han University of Minnesota 5/15/2009 - 5/16/2009
Carsten Hartmann Freie Universität Berlin 5/16/2009 - 5/22/2009
Neale Haugen St. Cloud State University 5/15/2009 - 5/16/2009
Graeme Henkelman University of Texas 5/17/2009 - 5/22/2009
Mark S. Herman University of Minnesota 9/1/2008 - 8/31/2010
Peter Hinow University of Minnesota 9/1/2007 - 8/21/2009
Jingfang Huang University of North Carolina 12/30/2008 - 5/31/2009
Yunkyong Hyon University of Minnesota 9/1/2008 - 8/31/2010
Olexandr Isayev Case Western Reserve University 5/17/2009 - 5/23/2009
Leah Yakeerah Isseroff Princeton University 5/3/2009 - 5/8/2009
Ivaylo N. Ivanov University of California, San Diego 5/17/2009 - 5/22/2009
Mark Iwen University of Minnesota 9/1/2008 - 8/31/2010
Alexander Izzo Bowling Green State University 9/1/2008 - 6/30/2009
Srividhya Jeyaraman University of Minnesota 9/1/2008 - 8/31/2010
Lijian Jiang University of Minnesota 9/1/2008 - 8/31/2010
Segun Jung New York University 5/17/2009 - 5/22/2009
Gerhard Kahl Technische Universität Wien 5/16/2009 - 5/22/2009
Hiqmet Kamberaj University of Minnesota 5/18/2009 - 5/22/2009
Dalal Kanan Princeton University 5/3/2009 - 5/8/2009
Markus Keel University of Minnesota 7/21/2008 - 6/30/2009
Svetlana Kotochigova Temple University 5/3/2009 - 5/8/2009
Robert Krasny University of Michigan 5/13/2009 - 5/23/2009
David Krisiloff Princeton University 5/3/2009 - 5/8/2009
Krzysztof Kuczera University of Kansas 5/17/2009 - 5/22/2009
Laura Largo Escudero University of Minnesota 5/4/2009 - 5/8/2009
Juan C. Latorre Freie Universität Berlin 5/17/2009 - 5/22/2009
Kah Chun Lau George Washington University 5/17/2009 - 5/22/2009
Claude Le Bris CERMICS 9/11/2008 - 5/30/2009
Chiun-Chang Lee National Taiwan University 8/26/2008 - 8/24/2010
Young-Ju Lee Rutgers University 5/14/2009 - 5/23/2009
Frédéric Legoll École Nationale des Ponts-et-Chaussées 5/5/2009 - 5/23/2009
Benedict Leimkuhler University of Edinburgh 5/3/2009 - 5/23/2009
Tony Lelievre Ecole Nationale des Ponts et Chaussees 5/12/2009 - 5/22/2009
Hannah Ruth Leverentz University of Minnesota 5/4/2009 - 5/8/2009
Kun Li Pennsylvania State University 5/14/2009 - 5/23/2009
Xiantao Li Pennsylvania State University 5/17/2009 - 5/22/2009
Yongfeng Li University of Minnesota 9/1/2008 - 8/31/2010
Peilin Liao Princeton University 5/3/2009 - 5/8/2009
Russ Lidberg St. Cloud State University 5/15/2009 - 5/16/2009
Hai Lin University of Colorado 3/1/2009 - 5/30/2009
Pinsker Lin University of Minnesota 5/15/2009 - 5/16/2009
Tai-Chia Lin National Taiwan University 8/23/2008 - 7/31/2009
Yuchun Lin University of North Carolina - Charlotte 5/17/2009 - 5/22/2009
Roland Lindh Lund University 5/2/2009 - 5/9/2009
Chun Liu University of Minnesota 9/1/2008 - 8/31/2010
Da-Jiang Liu Ames Laboratory 5/11/2009 - 5/13/2009
Gang Lu California State University 5/17/2009 - 5/23/2009
Mitchell Luskin University of Minnesota 9/1/2008 - 6/30/2009
Dmitrii E. Makarov University of Texas 5/17/2009 - 5/22/2009
Thérèse E. Malliavin Centre National de la Recherche Scientifique (CNRS) 5/16/2009 - 5/21/2009
Per-Åke Malmqvist Lund University 5/2/2009 - 5/8/2009
Manjeera Mantina University of Minnesota 5/4/2009 - 5/8/2009
Paul Maragakis D. E. Shaw Research 5/17/2009 - 5/22/2009
Luca Maragliano University of Chicago 5/16/2009 - 5/23/2009
Thomas Markland University of Oxford 5/17/2009 - 5/22/2009
Vasileios Maroulas University of Minnesota 9/1/2008 - 8/31/2010
Kevin W. Mclaughlin University of Wisconsin - River Falls 1/6/2009 - 6/30/2009
David A. Micha University of Florida 5/1/2009 - 5/31/2009
Thomas F Miller III California Institute of Technology 5/17/2009 - 5/22/2009
David Do Le Minh National Institutes of Health 5/17/2009 - 5/22/2009
Kimiya Minoukadeh Ecole Nationale des Ponts et Chaussees 5/14/2009 - 5/23/2009
Peter A. Monson University of Massachusetts 5/17/2009 - 5/22/2009
Abdul Rehaman Moughal Shahi Université de Genève 2/2/2009 - 6/30/2009
Kwangho Nam Harvard University 5/17/2009 - 5/22/2009
Emad Noorizadeh University of Edinburgh 5/15/2009 - 5/22/2009
Christoph Ortner University of Oxford 4/1/2009 - 5/3/2009
Victor B. Oyeyemi Princeton University 5/3/2009 - 5/8/2009
Albert C. Pan University of Chicago 5/17/2009 - 5/22/2009
Tom Pandsack University of Minnesota 5/15/2009 - 5/16/2009
Apirak Payaka University of Minnesota 5/15/2009 - 5/16/2009
Huo-Lei Peng Ohio State University 5/4/2009 - 5/8/2009
Malgorzata Peszynska Oregon State University 5/17/2009 - 5/20/2009
Alin Pogan University of Minnesota 5/12/2009 - 5/12/2009
Matej Praprotnik National Institute of Chemistry 5/15/2009 - 5/22/2009
Ravi Radhakrishnan University of Pennsylvania 5/17/2009 - 5/22/2009
Neeraj Rai University of Minnesota 5/4/2009 - 5/8/2009
Neeraj Rai University of Minnesota 5/18/2009 - 5/22/2009
Paul Rejto Pfizer Global R&D 4/30/2009 - 5/2/2009
Weiqing Ren New York University 5/17/2009 - 5/22/2009
Karsten Reuter Max-Planck-Institut 5/16/2009 - 5/20/2009
Andrew Ritzmann Princeton University 5/3/2009 - 5/8/2009
Adrian E. Roitberg University of Florida 5/17/2009 - 5/22/2009
Mikhail Ryazantsev Bowling Green State University 5/3/2009 - 5/8/2009
Daniel Sadowsky University of Minnesota 5/4/2009 - 5/8/2009
Anna Sankari Lund University 5/2/2009 - 5/9/2009
Fadil Santosa University of Minnesota 7/1/2008 - 6/30/2010
Carlos F. Sanz-Navarro Autonomous University of Barcelona 5/17/2009 - 5/23/2009
Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research(JNCASR) 5/17/2009 - 5/22/2009
Arnd Scheel University of Minnesota 9/1/2008 - 6/30/2009
Tamar Schlick New York University 5/18/2009 - 5/22/2009
Christof Schütte Freie Universität Berlin 5/18/2009 - 5/23/2009
L. Ridgway Scott University of Chicago 9/1/2008 - 6/30/2009
Tsvetanka Sendova University of Minnesota 9/1/2008 - 8/31/2010
Luis Serrano-Andrés University of Valencia 5/2/2009 - 5/12/2009
Sami Shahin Southwest Minnesota State University 5/18/2009 - 5/21/2009
Yuk Sham University of Minnesota 9/1/2008 - 6/30/2009
Juanfang Shen Purdue University 5/17/2009 - 5/22/2009
Robert D. Skeel Purdue University 3/29/2009 - 6/27/2009
Adnan Sljoka York University 5/15/2009 - 5/22/2009
Berend Smit University of California, Berkeley 5/14/2009 - 5/16/2009
Daniel B. Smith University of Pittsburgh 5/17/2009 - 5/24/2009
Andrew M. Stein University of Minnesota 9/1/2007 - 5/1/2009
Gabriel Stoltz École Nationale des Ponts-et-Chaussées (ENPC) 5/16/2009 - 5/22/2009
Huan Sun Pennsylvania State University 2/5/2009 - 5/31/2009
Panida Surawatanawong Texas A & M University 5/3/2009 - 5/9/2009
Devarajan Thirumalai University of Maryland 5/19/2009 - 5/22/2009
Donald G. Truhlar University of Minnesota 9/1/2008 - 6/30/2009
Erkan Tüzel University of Minnesota 9/1/2007 - 8/7/2009
Steven M. Valone Los Alamos National Laboratory 5/7/2009 - 6/30/2009
Eric Vanden-Eijnden New York University 5/21/2009 - 5/22/2009
Carlos Vega de las Heras Universidad Complutense de Madrid 5/17/2009 - 5/22/2009
Maddalena Venturoli New York University 5/17/2009 - 5/22/2009
Valera Veryazov Lund University 5/2/2009 - 5/9/2009
Bess Vlaisavljevich University of Minnesota 5/4/2009 - 5/8/2009
Rodolphe Vuilleumier École Normale Supérieure 5/17/2009 - 5/24/2009
David John Wales University of Cambridge 5/16/2009 - 5/21/2009
Bo Wang University of Minnesota 5/4/2009 - 5/8/2009
Zhian Wang University of Minnesota 9/1/2007 - 8/31/2009
Stephen Wiggins University of Bristol 1/10/2009 - 6/30/2009
Wei Xiong University of Minnesota 9/1/2008 - 8/31/2010
Xiang Xu Pennsylvania State University 1/26/2009 - 6/1/2009
Ryoichi Yamamoto Kyoto University 5/17/2009 - 5/22/2009
Ke Yang University of Minnesota 5/4/2009 - 5/8/2009
Ke Yang University of Minnesota 5/18/2009 - 5/22/2009
Xiao Yi University of Minnesota 5/15/2009 - 5/16/2009
Xiao Yi University of Minnesota 5/18/2009 - 5/22/2009
Darrin M. York University of Minnesota 5/18/2009 - 5/22/2009
Chensong Zhang Pennsylvania State University 5/15/2009 - 5/23/2009
Chunfeng Zhao University of Calgary 5/17/2009 - 5/23/2009
Ruijun Zhao Purdue University 5/17/2009 - 5/22/2009
Guishan Zheng Harvard University 5/17/2009 - 5/23/2009
Jingjing Zheng University of Minnesota 5/4/2009 - 5/8/2009
Weigang Zhong University of Minnesota 9/1/2008 - 8/31/2010
Yan Zhou University of Minnesota 5/15/2009 - 5/16/2009
Xiaoyang Zhu University of Minnesota 5/26/2009 - 5/26/2009
Daniel M. Zuckerman University of Pittsburgh 5/19/2009 - 5/21/2009
Legend: Postdoc or Industrial Postdoc Long-term Visitor

IMA Affiliates:
Arizona State University, Boeing, Corning Incorporated, ExxonMobil, Ford, General Motors, Georgia Institute of Technology, Honeywell, IBM, Indiana University, Iowa State University, Kent State University, Korea Advanced Institute of Science and Technology (KAIST), Lawrence Livermore National Laboratory, Lockheed Martin, Los Alamos National Laboratory, Medtronic, Michigan State University, Michigan Technological University, Microsoft Research, Mississippi State University, Motorola, Northern Illinois University, Ohio State University, Pennsylvania State University, Purdue University, Rice University, Rutgers University, Sandia National Laboratories, Schlumberger Cambridge Research Laboratories, Schlumberger-Doll, Seoul National University, Siemens, Telcordia, Texas A & M University, University of Central Florida, University of Chicago, University of Cincinnati, University of Delaware, University of Houston, University of Illinois at Urbana-Champaign, University of Iowa, University of Kentucky, University of Maryland, University of Michigan, University of Minnesota, University of Notre Dame, University of Pittsburgh, University of Tennessee, University of Wisconsin, University of Wyoming, US Air Force Research Laboratory, Wayne State University, Worcester Polytechnic Institute