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 #369

July 2007

News and Notes

Santosa named next IMA director

We are delighted to announce that Fadil Santosa has been appointed the next director of the Institute for Mathematics and its Applications (IMA), beginning in summer 2008. Please see the press release for more information.

Solicitation of Proposals for IMA Participating Institutions Conferences: All faculty members of the Participating Institutions of the Institute for Mathematics and its Applications are invited to submit proposals for the 2008-2009 IMA Participating Institution conferences. Deadline: November 15, 2007. Please see Solicitation of Proposals for more information.

IMA Events

PI Summer Graduate Program

Applicable Algebraic Geometry

July 23 - August 10, 2007

Organizers: Laura Felicia Matusevich (Texas A & M University), Frank Sottile (Texas A & M University), Thorsten Theobald (Johann Wolfgang Goethe-Universität Frankfurt)

Classical and Quantum Approaches in Molecular Modeling

July 23 - August 3, 2007

Organizers: Eric Cances (Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS)), Giovanni Ciccotti (Università di Roma "La Sapienza"), Benedict Leimkuhler (University of Edinburgh), Nicola Marzari (Massachusetts Institute of Technology), Yousef Saad (University of Minnesota Twin Cities), Gustavo E. Scuseria (Rice University), Robert D. Skeel (Purdue University), Mark E. Tuckerman (New York University)
Schedule

Wednesday, July 4

All DayIndependence Day. The IMA is closed.

Thursday, July 5

11:15a-12:15pTim Hardy, Wayne State College
Real algebraic geometry tutorial: Distance matrices
Lind Hall 409 RAG

Monday, July 23

All DayThemes for Birds of a Feather Sessions will be specified on Day 1. SP7.23-8.3.07
8:30a-9:15aRegistration and coffeeEE/CS 3-176 SP7.23-8.3.07
9:15a-9:30aWelcome and introductionDouglas N. Arnold (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07
9:30a-10:30aIntroduction to molecular dynamicsRobert D. Skeel (Purdue University)EE/CS 3-180 SP7.23-8.3.07
10:30a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-12:00pStatistical mechanics and molecular dynamicsMark E. Tuckerman (New York University)EE/CS 3-180 SP7.23-8.3.07
12:00p-2:00pLunch SP7.23-8.3.07
2:00p-3:00pDynamical equations and numerical integratorsBenedict Leimkuhler (University of Edinburgh)EE/CS 3-180 SP7.23-8.3.07
3:00p-3:30pCoffee SP7.23-8.3.07
3:30p-5:30pHands-on computer session on MD simulation (parallel session)

VMD and NAMD software packages/tutorials

Eric Barth (Kalamazoo College)
Stephen Bond (University of Illinois at Urbana-Champaign)
Robert D. Skeel (Purdue University)
Chris R. Sweet (University of Notre Dame)
EE/CS 3-180 hands-on computer session

EE/CS 3-230 parallel session
SP7.23-8.3.07
3:30p-4:00pIntegrators for highly oscillatory Hamiltonian systems: an homogenization approachFrederic Legoll (École Nationale des Ponts-et-Chaussées)EE/CS 3-230 SP7.23-8.3.07
4:00p-4:30pAnalysis and computational studies of the ergodicity of the Nose-Hoover thermostatMitchell Luskin (University of Minnesota Twin Cities)EE/CS 3-230 SP7.23-8.3.07

Tuesday, July 24

9:00a-9:30aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:30a-10:30aFree energy calculations and the potential of mean forceMark E. Tuckerman (New York University)EE/CS 3-180 SP7.23-8.3.07
10:30a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-12:00pMolecular samplingRobert D. Skeel (Purdue University)EE/CS 3-180 SP7.23-8.3.07
12:00p-2:00pLunch SP7.23-8.3.07
2:00p-3:00pConstraints and coarse grainingGiovanni Ciccotti (Università di Roma "La Sapienza")EE/CS 3-180 SP7.23-8.3.07
3:00p-3:30pCoffeeEE/CS 3-176 SP7.23-8.3.07
3:30p-5:00pHands-on computer session (parallel session)

VMD and NAMD software packages/tutorials

Eric Barth (Kalamazoo College)
Stephen Bond (University of Illinois at Urbana-Champaign)
Robert D. Skeel (Purdue University)
Chris R. Sweet (University of Notre Dame)
EE/CS 3-180 hands-on computer session

EE/CS 3-230 parallel session
SP7.23-8.3.07
3:30p-4:00pTowards a mathematical justification of kinetic theoryFlorian Theil (University of Warwick)EE/CS 3-230 SP7.23-8.3.07
4:00p-4:30pAnalysis of a lattice model for phase transitions and derivation of kinetic relationsJohannes Zimmer (University of Bath)EE/CS 3-230 SP7.23-8.3.07
5:00p-6:30pReception and Poster SessionLind Hall 400 SP7.23-8.3.07
Molecular modelling the structure and dynamics of alginate oligosaccharidesHoda Abdel-Aal Bettley (University of Manchester)
Method for determination of Hubbard model phase diagram from optical lattice experiments by two parameter scalingVivaldo L. Campo (University of Minnesota Twin Cities)
Absolute entropy and free energy of free and bound states of a mobile loop of alpha-amylase using the hypothetical scanning methodSrinath Cheluvaraja (University of Pittsburgh)
A comparison of maximum likelihood and weighted residual approximations to the potential of mean forceEric C. Cyr (University of Illinois at Urbana-Champaign)
Objective structures and their applicationsKaushik Dayal (University of Minnesota Twin Cities)
Modelling of local defects in crystalsAmélie Deleurence (École Nationale des Ponts-et-Chaussées (ENPC))
A systematic method to explore possible silicon tip structures used in AFMSeyed-Alireza Ghasemi (Universität Basel)
Multilevel summation method for Coulomb interactionsDavid J. Hardy (University of Illinois at Urbana-Champaign)
Conformational reinvestigation of two cyclic pentapeptides: to a generic approach in drug developmentPieter Hendrickx (University of Ghent (UG))
An ab initio molecular dynamics simulation of solid CL-20: mechanism and kinetics of thermal decompositionOlexandr Isayev (Jackson State University)
Numerical method for solving stochastic differential equations with non-Gaussian noiseChangho Kim (Korea Advanced Institute of Science and Technology (KAIST))
Particle-Scaling function (P3S) algorithm for electrostatic problems in free boundary conditionsAlexey Neelov (Universität Basel)
Temperature-regulated microcanonical dynamicsEmad Noorizadeh (University of Edinburgh)
A Bell-Evans-Polanyi principle for molecular dynamics trajectories and its implications for global optimizationShantanu Roy (Universität Basel)
Uniqueness of the density-to-potential mapping in excited-state density-functional theoryPrasanjit Samal (University of Minnesota Twin Cities)
Precision problems in density functional development for better molecular modelingMichael Teter (Cornell University)
Mesoscopic model for the fluctuating hydrodynamics of binary and ternary mixtures Erkan Tüzel (North Dakota State University)
New numerical algorithms and software for minimizing biomolecular potential energy functionsDexuan Xie (University of Wisconsin)

Wednesday, July 25

9:00a-9:30aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:30a-10:30aTransition pathways of rare eventsEric Vanden-Eijnden (New York University)EE/CS 3-180 SP7.23-8.3.07
10:30a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-12:00pRate constants Giovanni Ciccotti (Università di Roma "La Sapienza")
Mark E. Tuckerman (New York University)
EE/CS 3-180 SP7.23-8.3.07
12:15p-2:00pLunch SP7.23-8.3.07
2:00p-3:00pAdaptive methods for free energy computation and coarse-graining strategies Eric Felix Darve (Stanford University)EE/CS 3-180 SP7.23-8.3.07
3:00p-3:30pCoffeeEE/CS 3-176 SP7.23-8.3.07
3:30p-4:25p"Birds of a feather" sessions:

Data Abstractions
(Robert D. Skeel, coordinator)
EECS 3-180

Coarse Graining
(Claude Le Bris, coordinator)
EECS 3-111

EECS 3-180
EECS 3-111
SP7.23-8.3.07
4:35p-5:30p"Birds of a feather" sessions:

Thermostatting
(Stephen Bond and Frederic Legoll, coordinators)
EECS 3-180

Collective Variables
(Giovanni Ciccotti and Eric Vanden-Eijnden), coordinators
EECS 3-111

EECS 3-180
EECS 3-111
SP7.23-8.3.07

Thursday, July 26

8:45a-9:15aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:15a-9:55aMilestoning, extending time scale of molecular simulationsRon Elber (University of Texas)EE/CS 3-180 SP7.23-8.3.07
9:55a-10:35aNanomechanics of biomoleculesChristof Schuette (Freie Universität Berlin)EE/CS 3-180 SP7.23-8.3.07
10:35a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-11:20aA reduced stochastic model for shock and detonation wavesGabriel Stoltz (École Nationale des Ponts-et-Chaussées (ENPC))EE/CS 3-180 SP7.23-8.3.07
11:40a-12:00pEfficient implementation of adaptively biased molecular dynamicsCeleste Sagui (North Carolina State University)EE/CS 3-180 SP7.23-8.3.07
12:00p-12:15pGroup Photo SP7.23-8.3.07
12:15p-2:00pLunch SP7.23-8.3.07
2:00p-2:40pHybrid Monte Carlo methodsJesus A. Izaguirre (University of Notre Dame)EE/CS 3-180 SP7.23-8.3.07
2:40p-3:00pNormal mode partitioning of Langevin dynamics for the simulation of biomoleculesChris R. Sweet (University of Notre Dame)EE/CS 3-180 SP7.23-8.3.07
3:00p-3:30pCoffeeEE/CS 3-176 SP7.23-8.3.07
3:30p-3:50pEfficient modeling of solvent environments Michael Feig (Michigan State University)EE/CS 3-180 SP7.23-8.3.07
3:50p-4:10pCluster optimization in protein dockingJulie C. Mitchell (University of Wisconsin)EE/CS 3-180 SP7.23-8.3.07
4:10p-4:30pUmbrella sampling for non-equilibrium processesAaron Dinner (University of Chicago)EE/CS 3-180 SP7.23-8.3.07
4:30p-5:00pSecond ChancesEE/CS 3-180 SP7.23-8.3.07
6:30p-8:30pGroup dinner at KikugawaKikugawa at Riverplace
43 Main St. SE. Minneapolis, MN (612) 378-3006
SP7.23-8.3.07

Friday, July 27

8:45a-9:15aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:15a-9:55aCoarse-graining the free energy of atomistic systems: a mathematical approachClaude Le Bris (École Nationale des Ponts-et-Chaussées (ENPC))EE/CS 3-180 SP7.23-8.3.07
9:55a-10:15aEstimating accuracy in classical molecular simulationStephen Bond (University of Illinois at Urbana-Champaign)EE/CS 3-180 SP7.23-8.3.07
10:15a-10:35aAdaptive methods for molecular dynamics Tony Lelievre (Ecole Nationale des Ponts et Chaussees)EE/CS 3-180 SP7.23-8.3.07
10:35a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-11:20aBoundary conditions for non-equilibrium molecular dynamics simulationXiantao Li (Pennsylvania State University)EE/CS 3-180 SP7.23-8.3.07
11:20a-11:40aQuantum reaction ratesRaymond Kapral (University of Toronto)EE/CS 3-180 SP7.23-8.3.07
11:40a-2:00pLunch SP7.23-8.3.07
2:00p-2:40pAccelerated molecular dynamics methodsArthur F. Voter (Los Alamos National Laboratory) SP7.23-8.3.07
2:40p-3:00pCoffeeEE/CS 3-176 SP7.23-8.3.07
3:00p-3:20pMD modelling of primary damage production in displacement cascadesRoman Voskoboynikov (University of Cambridge)EE/CS 3-180 SP7.23-8.3.07
3:20p-3:40pThermal boundary conditions for molecular dynamics simulationsSimon Gill (University of Leicester)EE/CS 3-180 SP7.23-8.3.07
3:40p-4:00pProtein materials balance strength, energy dissipation and robustness by selecting nanopatterned, hierarchical featuresMarkus J. Buehler (Massachusetts Institute of Technology)EE/CS 3-180 SP7.23-8.3.07
4:00p-4:30pSecond ChancesEE/CS 3-180 SP7.23-8.3.07

Saturday, July 28

8:30a-9:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:00a-9:40aQuantum Drude models : A complete description of polarization and dispersionGlenn Martyna (IBM Corporation)EE/CS 3-180 SP7.23-8.3.07
9:40a-10:00aEffective normal modes from finite temperature molecular dynamics simulationsRodolphe Vuilleumier (Université de Paris VI (Pierre et Marie Curie))EE/CS 3-180 SP7.23-8.3.07
10:00a-10:20aDirect calculation of interfacial free energy using molecular simulation Brian Laird (University of Kansas)EE/CS 3-180 SP7.23-8.3.07
10:20a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-11:20aAspects of nonautonomous molecular dynamicsMichel Cuendet (Swiss Institute of Bioinformatics)EE/CS 3-180 SP7.23-8.3.07
11:20a-12:00pSecond chances and closing remarksEE/CS 3-180 SP7.23-8.3.07

Sunday, July 29

All DayNo activity scheduled. SP7.23-8.3.07

Monday, July 30

9:00a-9:30aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:30a-10:45aIntroduction to quantum chemistry and ab initio calculations of interatomic potentials Eric Cances (Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS))EE/CS 3-180 SP7.23-8.3.07
10:45a-11:15aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:15a-12:30pNumerical algorithms for density functional theoryYousef Saad (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07
12:30p-2:30aLunch SP7.23-8.3.07
2:30p-3:30pDensity functional theory for periodic systemsEric Cances (Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS))EE/CS 3-180 SP7.23-8.3.07
3:30p-4:00pCoffeeEE/CS 3-176 SP7.23-8.3.07
4:00p-5:00pNumerical algorithms for density functional theory (continued)Yousef Saad (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07

Tuesday, July 31

9:00a-9:30aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:30a-10:45aDensity-functional practice Nicola Marzari (Massachusetts Institute of Technology)EE/CS 3-180 SP7.23-8.3.07
10:45a-11:15aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:15a-12:30pElectronic correlations and Hubbard approachesMatteo Cococcioni (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07
12:30p-2:30aLunch SP7.23-8.3.07
2:30p-3:30pTutorial on static calculationsMatteo Cococcioni (University of Minnesota Twin Cities)
Ismaila Dabo (Massachusetts Institute of Technology)
EE/CS 3-180 SP7.23-8.3.07
3:30p-4:00pCoffeeEE/CS 3-176 SP7.23-8.3.07
4:00p-5:30pTutorial on Car-Parrinello molecular dynamicsNicola Marzari (Massachusetts Institute of Technology)
Arash A. Mostofi (University of Cambridge)
EE/CS 3-180 SP7.23-8.3.07

Wednesday, August 1

8:30a-9:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:00a-10:00aReal space pseudopotentials applied to nanoscale systemsJames R. Chelikowsky (University of Texas)EE/CS 3-180 SP7.23-8.3.07
10:00a-10:15aCoffeeEE/CS 3-176 SP7.23-8.3.07
10:15a-11:15aDensity-functional theory and its generalizations: legendre transform, constrained search, open problemsPaul W. Ayers (McMaster University)EE/CS 3-180 SP7.23-8.3.07
11:15a-11:30aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:30a-12:00pTBAMatt Challacombe (Los Alamos National Laboratory)EE/CS 3-180 SP7.23-8.3.07
12:00p-12:30pNovel materials for quantum computingNicholas M. Harrison (Imperial College London)EE/CS 3-180 SP7.23-8.3.07
12:30p-2:30aLunch SP7.23-8.3.07
2:30p-3:00pTheoretical description of electrons in singloe molecule magnetsErnest R. Davidson (University of Washington)EE/CS 3-180 SP7.23-8.3.07
3:00p-4:00pA consistent, linear-response approach to LDA+U Matteo Cococcioni (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07
4:00p-4:30pReal-space finite difference method for O(N) first-principles molecular dynamics with plane waves accuracyJean-Luc Fattebert (Lawrence Livermore National Laboratory)EE/CS 3-180 SP7.23-8.3.07
4:30p-4:45pGroup Photo SP7.23-8.3.07
4:45p-6:15pPoster SessionLind Hall 400 SP7.23-8.3.07
Molecular modelling the structure and dynamics of alginate oligosaccharidesHoda Abdel-Aal Bettley (University of Manchester)
ethod for determination of Hubbard model phase diagram from optical lattice experiments by two parameter scalingVivaldo L. Campo (University of Minnesota Twin Cities)
Real-space corrections for electrostatic interactions in periodic boundary conditionsIsmaila Dabo (Massachusetts Institute of Technology)
Objective structures and their applicationsKaushik Dayal (University of Minnesota Twin Cities)
Modelling of local defects in crystalsAmélie Deleurence (École Nationale des Ponts-et-Chaussées (ENPC))
A systematic method to explore possible silicon tip structures used in AFMSeyed-Alireza Ghasemi (Universität Basel)
Conformational reinvestigation of two cyclic pentapeptides: to a generic approach in drug developmentPieter Hendrickx (University of Ghent (UG))
An ab initio molecular dynamics simulation of solid CL-20: mechanism and kinetics of thermal decompositionOlexandr Isayev (Jackson State University)
Numerical method for solving stochastic differential equations with non-Gaussian noiseChangho Kim (Korea Advanced Institute of Science and Technology (KAIST))
Particle-Scaling function (P3S) algorithm for electrostatic problems in free boundary conditionsAlexey Neelov (Universität Basel)
A Bell-Evans-Polanyi principle for molecular dynamics trajectories and its implications for global optimizationShantanu Roy (Universität Basel)
Uniqueness of the density-to-potential mapping in excited-state density-functional theoryPrasanjit Samal (University of Minnesota Twin Cities)
Local exchange potentials: A mathematical viewpointGabriel Stoltz (École Nationale des Ponts-et-Chaussées (ENPC))
Precision problems in density functional development for better molecular modelingMichael Teter (Cornell University)
Mesoscopic model for the fluctuating hydrodynamics of binary and ternary mixturesErkan Tüzel (North Dakota State University)
New numerical algorithms and software for minimizing biomolecular potential energy functionsDexuan Xie (University of Wisconsin)

Thursday, August 2

8:30a-9:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:00a-10:00aWavelets for electronic structure calculations and electrostatic problemsStefan Goedecker (Universität Basel)EE/CS 3-180 SP7.23-8.3.07
10:00a-10:15aCoffeeEE/CS 3-176 SP7.23-8.3.07
10:15a-11:15aExchange and correlation in electronic systems: the hole story Axel D. Becke (Dalhousie University)EE/CS 3-180 SP7.23-8.3.07
11:15a-11:30aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:30a-12:00pEfficient Kohn-Sham density functional calculations using the Gaussian and plane waves approachJürg Hutter (Universität Zürich)EE/CS 3-180 SP7.23-8.3.07
12:00p-12:30pTBADavid Langreth (Rutgers University)EE/CS 3-180 SP7.23-8.3.07
12:30p-2:30pLunch SP7.23-8.3.07
2:30p-3:00pLinear-scaling density-functional calculations with plane-wavesArash A. Mostofi (University of Cambridge)EE/CS 3-180 SP7.23-8.3.07
3:00p-3:30pA Linear-scaling AO-based MP2 method for large molecules by rigorous integral estimatesChristian Ochsenfeld (Eberhard-Karls-Universität Tübingen)EE/CS 3-180 SP7.23-8.3.07
3:30p-4:00pCoffeeEE/CS 3-176 SP7.23-8.3.07
4:00p-5:30pSecond Chances session on fast algorithms EE/CS 3-180 SP7.23-8.3.07
6:30p-8:30pGroup dinner at Caspian BistroCaspian Bistro 2418 University Ave SE Minneapolis, MN 55414
(612) 623-1113
SP7.23-8.3.07

Friday, August 3

8:30a-9:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
9:00a-9:30aDealing with spatial regions Andreas Savin (Université de Paris VI (Pierre et Marie Curie))EE/CS 3-180 SP7.23-8.3.07
9:30a-10:00aKohn-Sham methods for implicit density functionalsViktor N. Staroverov (University of Western Ontario)EE/CS 3-180 SP7.23-8.3.07
10:00a-10:30aQM/MM studies on enzymesWalter Thiel ( Max-Planck-Institut für Kohlenforschung)EE/CS 3-180 SP7.23-8.3.07
10:30a-11:00aCoffeeEE/CS 3-176 SP7.23-8.3.07
11:00a-11:30aTBADonald Truhlar (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07
11:30a-12:00pTBAY. Alexander Wang (University of British Columbia)EE/CS 3-180 SP7.23-8.3.07
12:00p-12:30pTBARenata Wentzcovitch (University of Minnesota Twin Cities)EE/CS 3-180 SP7.23-8.3.07
12:30p-2:30aLunch SP7.23-8.3.07
2:30p-3:00pOrbital-free embedding potential: properties, approximations, and the use in computer simulations to couple quantum chemical and classical levels of descriptionTomasz A. Wesolowski (Université de Genève)EE/CS 3-176 SP7.23-8.3.07
3:00p-4:30pSecond Chances session on DFT EE/CS 3-180 SP7.23-8.3.07

Wednesday, August 8

All DayWorkshop Outline: Posing of problems by the 6 industry mentors. Half-hour introductory talks in the morning followed by a welcoming lunch. In the afternoon, the teams work with the mentors. The goal at the end of the day is to get the students to start working on the projects.EE/CS 3-180 MM8.8-17.07
9:00a-9:30aCoffee and RegistrationEE/CS 3-176 MM8.8-17.07
9:30a-9:40aWelcome and IntroductionDouglas N. Arnold (University of Minnesota Twin Cities)
Richard J. Braun (University of Delaware)
Fernando Reitich (University of Minnesota Twin Cities)
Arnd Scheel (University of Minnesota Twin Cities)
EE/CS 3-180 MM8.8-17.07
9:40a-10:00aTeam 1: Supersonic designNatalia Alexandrov (NASA Langley Research Center)EE/CS 3-180 MM8.8-17.07
10:00a-10:20aTeam 2: 802.11 WLAN MAC layer modelingRadu V. Balan (Siemens Corporate Research, Inc.)EE/CS 3-180 MM8.8-17.07
10:20a-10:40aTeam 3: Associating earth-orbiting objects detected by astronomical telescopesGary B. Green (The Aerospace Corporation)EE/CS 3-180 MM8.8-17.07
10:40a-11:00aBreakEE/CS 3-176 MM8.8-17.07
11:00a-11:20aTeam 4: High dimensional, nonlinear, non-convex optimization problems in the area of aircraft and vehicle designJohn R. Hoffman (Lockheed Martin Missiles and Space Company, Inc.)EE/CS 3-180 MM8.8-17.07
11:20a-11:40aTeam 5: Size and shape comparisons from noisy, unlabeled, incomplete configurations of landmarks in three-dimensional spaceMark A. Stuff (General Dynamics Advanced Information Systems)EE/CS 3-180 MM8.8-17.07
11:40a-12:00pTeam 6: Wavelength assignment and conversion in optical networkingLisa Zhang (Lucent Technologies Bell Laboratories)EE/CS 3-180 MM8.8-17.07
12:00p-1:30pLunchLind Hall 400 MM8.8-17.07
1:30p-4:30pafternoon - start work on projectsBreak-out Rooms MM8.8-17.07

Thursday, August 9

All DayStudents work on the projects. Mentors guide their groups through the modeling process, leading discussion sessions, suggesting references, and assigning work.Break-out Rooms MM8.8-17.07

Friday, August 10

All DayStudents work on the projects. Mentors guide their groups through the modeling process, leading discussion sessions, suggesting references, and assigning work.Break-out Rooms MM8.8-17.07
Abstracts
Tim Hardy, Wayne State College
Real algebraic geometry tutorial: Distance matrices
Abstract: "distance matrix" is a matrix containing the distances, taken pairwise, of a set of points. I shall trace the history of distance matrices during the last 150 years. I shall develop distance matrices from an inner product. I shall present some applications and at least one generalization.
Hoda Abdel-Aal Bettley (University of Manchester) Molecular modelling the structure and dynamics of alginate oligosaccharides
Abstract: Joint work with Richard A. Bryce. Alginate copolymers are a key component of the extracellular polymeric substances (EPS) matrix of bacterial microorganisms such as P. aeruginosa, and occur in high abundance in nature. Alginate heteropolysaccharides comprise alternating blocks of alpha-(1->4)-linked L-guluronate and beta-(1->4)-linked D-mannuronate. To understand the microscopic behaviour and interactions of these flexible acidic sugars within the EPS matrix, a suitable molecular-level model is required. We derive a molecular mechanical force field for the two uronic acids, with validation against available experimental data. We then construct a detailed study of the structure and dynamics of alginate chains. We explore alginate models of increasing complexity, from disaccharides to single- and double-stranded oligomer helices, employing the techniques of molecular dynamics and replica-exchange molecular dynamics. The condensed phase behaviour of these systems is discussed, including the role of counterions and the implications for interaction with other constituents of the EPS matrix.
Hoda Abdel-Aal Bettley (University of Manchester) Molecular modelling the structure and dynamics of alginate oligosaccharides
Abstract: Same abstract as the 7/24 poster session.
Natalia Alexandrov (NASA Langley Research Center) Team 1: Supersonic design
Abstract: Designing affordable, efficient, quiet supersonic passenger aircraft has been under investigation for many years. Obstacles to designing such aircraft are also many, both in fundamental physics and in computational science and engineering. The problem of design is multidisciplinary in its nature and the goals of the constituent disciplines that govern the behavior of an aircraft are often at odds. In particular, aircraft that yields low sonic boom may not be attractive aerodynamically, while aerodynamically optimized aircraft may produce unacceptable sonic boom. One of the essential difficulties in using direct optimization methods to design for low boom and low drag is in modeling the design problem. For instance, it is not clear what objective functions to use.

Some Early Boom Shaping Developments (Ferri, 1969)

This project will use simple aerodynamic and sonic boom models to examine modeling of the design problem itself. We will attempt to establish a meaningful direct functional dependence between the shape of the aircraft and aerodynamic and noise quantities of interest by studying the sensitivity of these quantities to changes in shape. We will experiment with several direct multiobjective optimization problem formulations. References:

  1. Seebass, R., Argrow, B.; "Sonic Boom Minimization Revisited", AIAA Paper 98-2956
  2. Shepherd, K.P., Sullivan, B.M.; "A Loudness Calculation Procedure Applied to Shaped Sonic Booms", NASA Technical Paper 3134, 1991
  3. Carlson, H.W., Maglieri, D.J.; "Review of Sonic Boom Generation Theory and Prediction Methods", J. Acoust. Soc. Amer., 51, pp. 675-685 (1972)
  4. Alonso, J.J., Kroo, I.M., Jameson, A. "Advanced Algorithms for Design and Optimization of Quiet Supersonic Platform", 40th AIAA Aerospace Sciences Meeting and Exhibit, AIAA Paper 2002-0144, Reno, NV, January 2002
  5. Raymer, D.P.; "Aircraft Design: a Conceptual Approach", Third Edition, AIAA, 1999

Prerequisites:

Required: Scientific computing skills (Matlab or Fortran 90/95 or C), 1 semester in nonlinear optimization

Desired: Some background in statistical modeling, numerical analysis, multiobjective optimization

Keywords: multidisciplinary optimization, supersonic design, low boom, aerodynamic optimization

Paul W. Ayers (McMaster University) Density-functional theory and its generalizations: legendre transform, constrained search, open problems
Abstract: The quantum many-electron problem is easy in principle (solve the N-electron Schrödinger equation) and hard in practice (because the cost of numerical methods typically grows exponentially with the number of variables). However, there are simplifying features. First, the dimensionality can be reduced because electronic Hamiltonians contain only 1-body and 2-body terms. (This leads to reduced density-matrix methods.) Second, the dimensionality can be reduced because electrons are identical particles: if you know everything about one electron, then you know everything about all of the electrons. (This leads to electron-propagator theory and density-functional theory.) There is a “catch.” Reducing the number of dimensions leads to other problems associated with approximating the energy functional and/or associated with restricting the domain of the variational procedure. Two powerful techniques for resolving these difficulties are the Legendre transform and constrained-search formulations of density functional theory. This talk will discuss these formulations, and show how they can be extended to define "generalized" density-functional theories. I'll conclude with some of my favorite open problems in density-functional theory.
Radu V. Balan (Siemens Corporate Research, Inc.) Team 2: 802.11 WLAN MAC layer modeling
Abstract: 802.11 Wireless Local Area Networks (WLANs) have become as ubiquitous as Internet access for personal computers. The basic unit of a WLAN is composed of one Access Point (AP), and several mobile stations (STAs), all forming a Basic Service Set (BSS). A typical WLAN setup is depicted in Figure 1.

Figure 1: A typical Basic Service Set (BSS), with one AP, and several mobile stations.

The IEEE Standard governing WLANs describes two modes of operation: Distributed Coordination Function (DCF), and Point Coordination Function (PCF). By and large, chipset manufacturers implement only the DCF mode, and compatibility testing is done for this mode exclusively. The DCF is a contention-based mechanism where each wireless device (AP, or STA) competes for air time. More specifically, the 802.11 standard is implemented as follows:

  1. At regular intervals (typically hundreds of ms) the AP broadcasts a beacon signal, which resets all devices internal clocks;
  2. Assume a transmission opportunity ended at time t0. Depending on the status of the internal Backoff counter (BCK) of the device, the following actions can take place:
    • If BCKr=0, and there is no activity on air for DIFS (Distributed Inter Frame Spacing = 50us in 802.11b) time, then station starts transmitting its data packet;
    • If BCK=0 and during the DIFS period there is activity on air then device generates a random BCK between 0 and CW-1 (initially CW=CWmin = 16, in 802.11b); Then the following rules apply:
      • For each slot time (Ts) of medium inactivity, the BCK decrements;
      • The countdown is stopped whenever medium is busy, and the the countdown is resumed only after a supplemental AIFS (Arbitration Inter Frame Spacing, =DIFS in 802.11b) wait;
      • When BCK reaches 0, the device transmits its packet data;
      • If receiver (AP, or STA) receives successfully the packet, then it sends back on the air an Acknowledgement (ACK) frame, after a SIFT (Short Inter Frame Spacing = 10us) period after transmitter finishes its transmission;
      • If transmitter receives the ACK correctly, then it assumes data was received correctly, and transmission ends; On the other hand, if ACK is not broadcast, or the transmitter does not receive correctly the ACK, then it assumes the transmission was not successful, and the following rules apply:
        1. If current number of retransmissions has not reached a max threshold, then increment the Number of Retransmissions counter
        2. If CW<= CWmax, then CW doubles;
        3. A new random BCK is generated between 0 and CW-1;
        4. Transmission process is restarted from step b. above
    • When transmission ends, a post-backoff mechanism is implemented, by which a random BCK between 0 and CWmin-1 is generated, and a virtual countdown process is started obeying b.i and b.ii above.
These (somewhat simplified) rules govern the behavior of 802.11 devices. A big challenge in WLAN research is in modeling such a system. The purpose of this research group is to advance the current state-of-the-art model to allow for effective network control algorithm design.

Description of the problem

Basically there are two distinct regimes, completely opposite from one another:

  1. Deterministic Regime: when no collision happen, and the initial MAC instance time are sufficiently far apart, then transmission happens in a deterministic mode. Such a case may happen when only voice data (such as VoIP stations are connected to the AP), or periodic transmitting stations are present. The deterministic regime analysis gives an upper bound on system performance;
  2. Stochastic Regime: once collisions happen, or medium is detected busy during a packet arrival, the random generation of a Backoff counter happens, and the contention-based mechanism kicks in.
The deterministic regime is used to compute maximal performance of a WLAN. In such a case, performance may be superior even to the PCF mode, where AP acts as a transmission controller. However, in highly loaded networks, collisions are quite frequent, and the stochastic regime is more likely. Several works proposed stochastic models for this regime. Each work concerned one feature or another of network behavior. Bianchi [1] was the first to propose the use of Markov Chain in modeling the saturation regime of a WLAN. Since his paper, several others considered saturation, and non-saturation modeling of WLANs, increasing the model complexity, and taking into account more phenomena observed in experimental setups. In particular [2] represents a relatively good stochastic model for several regimes of WLAN. A somewhat refined diagram is presented in Figure 2. However the current state-of-the-art model is not sufficient for several reasons:
  1. It does not take into account the deterministic regime, nor does the performance converge to that upper bound;
  2. The Markovianity assumption is not always justifiable; is it possible to introduce a deterministic-stochastic hybrid model?
  3. Subsequent improvements of the standard are not yet captured by the current model; in particular the 802.11n draft introduces new MAC mechanisms.

The goal of this research group is to address one or more of the issues above. Ideally, students should have:

- familiarity with basic stochastic modeling concepts (such as Markov chains)
- familiarity with use of network simulation software (such as ns2);
- familiarity with time-series data analysis software (such as perl, Matlab);

Figure 2: A Markov Chain Model for a WLAN device.

Bibliography

[1] G.Bianchi, Performance Analysis of the IEEE 802.11 Distributed Coordination Function, IEEE Journal on Selected Areas of Communications, 18 (3), 2000, 535-547.

[2] P.E.Engelstad and O.N.Osterbo, Non-Saturation and Saturation Analysis of IEEE 802.11e EDCA with Starvation Prediction, MSWiM.05: Proceedings of the 8th ACM international symposium on Modeling, analysis and simulation of wireless and mobile systems, Montreal, Canada, 2005.

Axel D. Becke (Dalhousie University) Exchange and correlation in electronic systems: the hole story
Abstract: Exchange and correlation effects in electronic systems are rigorously related to a two-electron function called the exchange-correlation "hole". Modelling of the hole in real space is a powerful route to development and refinement of exchange and correlation functionals in DFT. We have developed real-space models of all correlation types of importance in chemical physics (dynamical, nondynamical, and dispersion) and these will be reviewed.
Stephen Bond (University of Illinois at Urbana-Champaign) Estimating accuracy in classical molecular simulation
Abstract: In classical molecular dynamics, the accuracy of numerical methods should be measured with respect to statistical averages, due to the lack of a meaningful exact trajectory. In this talk I will survey some results from backward error analysis and show how (under certain assumptions) these results can be applied to compute estimates of the error in averages from molecular dynamics simulations. Results from several test problems will be explored including examples from constant temperature molecular dynamics using a Nosé thermostat.
Markus J. Buehler (Massachusetts Institute of Technology) Protein materials balance strength, energy dissipation and robustness by selecting nanopatterned, hierarchical features
Abstract: Deformation and fracture are fundamental phenomena with major implications on the stability and reliability of machines, buildings and biological systems. All deformation processes begin with erratic motion of individual atoms around flaws or defects that quickly evolve into formation of macroscopic fractures as chemical bonds rupture rapidly, eventually compromising the integrity of the entire structure. However, most existing theories of fracture treat matter as a continuum, neglecting the existence of atoms or nanoscopic features. Clearly, such a description is questionable. Here we discuss an atomistic approach to describe such processes using ultra large-scale molecular dynamics (MD) simulation implemented supercomputers. MD provides unparalleled insight into the complex atomic-scale deformation processes, linking nano to macro, without relying on empirical input, since all atomic interaction parameters can be derived from fundamental quantum chemical theories. We demonstrate how MD can be used within a multi-scale simulation framework to predict the elastic and fracture properties of hierarchical protein materials, marvelous examples of structural designs that balance a multitude of tasks, representing some of the most sustainable material solutions that integrate structure and function across the scales. Breaking the material into its building blocks enables us to perform systematic studies of how microscopic design features influence the mechanical behavior at larger scales. We review studies of collagen – Nature’s super-glue, spider silk – a natural fiber that can reach the strength of a steel cable, as well as intermediate filaments – an important class of structural proteins responsible for the mechanical integrity of cells, which, if flawed, can cause serious diseases such as the rapid aging disease progeria. The common ground of these examples is the significance of the material properties at large deformation, its alteration under stress, presence of defects or the effect of variation of environmental conditions. Our studies elucidate intriguing material concepts that enable to balance strength, energy dissipation and robustness by selecting nanopatterned, hierarchical features.
Vivaldo L. Campo (University of Minnesota Twin Cities) Method for determination of Hubbard model phase diagram from optical lattice experiments by two parameter scaling
Abstract: We propose[1] an experimental scheme to obtain the phase diagram of the Hubbard model using cold atoms in optical lattices.[2] The scheme is based on measuring the total energy for a sequence of trapping potentials with different profile and is independent of di- mensionality. Its essential ingredient is a two-parameter scaling pro- cedure that combines a variant of the familiar finite-size scaling with a nontrivial additional ’finite-curvature scaling’ necessary to reach the homogeneous limit. We illustrate the viability of the scheme in the one-dimensional case, using simulations based on the Bethe Ansatz local-density approximation as a substitute for experimental data, and show that the filling corresponding to the Mott transition can be determined with better than 3% accuracy.
Vivaldo L. Campo (University of Minnesota Twin Cities) ethod for determination of Hubbard model phase diagram from optical lattice experiments by two parameter scaling
Abstract: Same abstract as the 7/24 poster session.
Eric Cances (Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS)) Introduction to quantum chemistry and ab initio calculations of interatomic potentials
Abstract: Quantum Chemistry aims at understanding the properties of matter through the modelling of its behaviour at a subatomic scale, where matter is described as an assembly of nuclei and electrons. At this scale, the equation that rules the interactions between these constitutive elements is the Schrdinger equation. It can be considered (except in few special cases notably those involving relativistic phenomena or nuclear reactions) as a universal model for at least three reasons. First it contains all the physical information of the system under consideration so that any of the properties of this system can be deduced in theory from the Schrdinger equation associated to it. Second, the Schrdinger equation does not involve any empirical parameter, except some fundamental constants of Physics (the Planck constant, the mass and charge of the electron, ...); it can thus be written for any kind of molecular system provided its chemical composition, in terms of natures of nuclei and number of electrons, is known. Third, this model enjoys remarkable predictive capabilities, as confirmed by comparisons with a large amount of experimental data of various types. Unfortunately, the Schrdinger equation cannot be directly simulated, except for very small chemical systems. It indeed reads as a time-dependent 3(M+N)-dimensional partial differential equation, where M is the number of nuclei and N the number of the electrons in the system under consideration. On the basis of asymptotic and semiclassical limit arguments, it is however often possible to approximate the Schrdinger dynamics by the so-called Born-Oppenheimer dynamics, in which nuclei behave as classical point-like particles. The internuclei (or interatomic) potential can be computed ab initio, by solving the time-independent electronic Schrdinger equation. The latter equation is a 3N-dimensional partial differential equation (it is in fact a spectral problem), for which several approximation methods are available. The main of them are the wavefunction methods and the Density Functional Theory (DFT). They lead in particular to the Hartree-Fock model and to the Kohn-Sham model, respectively. In this introductory lecture, I will first present the electronic Schrdinger equation, and show how to derive from this equation the Hartree-Fock and Kohn-Sham models. Although obtained from totally different approaches, these models have similar mathematical structures. They read as constrained optimization problems, whose Euler-Lagrange equations are nonlinear eigenvalue problems. In the second part of the lecture, I will introduce the Linear Combination of Atomic Orbitals (LCAO) discretization method, and briefly discuss two important numerical issues: ensuring self-consistent field (SCF) convergence, and reducing computational scaling.
Eric Cances (Centre d'Enseignement et de Recherche en Mathématiques et Calcul Scientifique (CERMICS)) Density functional theory for periodic systems
Abstract: By means of thermodynamic limit arguments, the Kohn-Sham model can be extended to the case when the nuclear configuration is periodic. The resulting periodic Kohn-Sham model is very useful to understand the electronic structure of perfect crystals, but also to simulate the condensed phase by means of artificial periodic boundary conditions. In the latter case, the periodic cell may contain a very large number of nuclei and electrons. In the first part of the lecture, I will present the Bloch theory of periodic Schrdinger operators, and recall how the band structure of their spectrum can be related to the insulating or conducting character of the underlying material. I will also present the concept of Wannier functions, which is of growing importance in various applications (computation of dielectric and magnetic properties, local correlation methods, electronic structure of crystals with defects, ...). In the second part of the lecture, I will focus on the discretization of the periodic Kohn-Sham model in plane-wave (PW) basis sets.
James R. Chelikowsky (University of Texas) Real space pseudopotentials applied to nanoscale systems
Abstract: One of the most challenging issues in materials physics is to predict the properties of matter at the nanoscale. In this size regime, new structural and electronic properties exist that resemble neither the atomic, nor solid state. These altered properties can have profound technological implications. Theoretical methods to address such issues face formidable challenges. Nanoscale systems may contain thousands of electrons and atoms, and often possess little symmetry. I will illustrate some recent advances in this area based on new computational methods and apply these techniques to systems ranging from clusters of a few dozen atoms to quantum dots containing thousands of atoms. Recent publications: Y. Zhou, Y. Saad, M.L. Tiago, and J.R. Chelikowsky: "Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration,'' Phys. Rev. E 74, 066704 (2006). M.L. Tiago, Y. Zhou, M.M.G. Alemany, Y. Saad, J.R. Chelikowsky: "The Evolution of Magnetism in Iron from the Atom to the Bulk,'' Phys. Rev. Lett. 97, 147201 (2006). M. Lopez del Puerto, M.L. Tiago, and J.R. Chelikowsky: "Excitonic effects and optical properties of passivated CdSe clusters,'' Phys. Rev. Lett. 97, 096401 (2006).
Srinath Cheluvaraja (University of Pittsburgh) Absolute entropy and free energy of free and bound states of a mobile loop of alpha-amylase using the hypothetical scanning method
Abstract: A new simulation method is studied which is capable of calculating the absolute entropy and free energy of a microstate. The method is based on calculating the probability of a member in the statistical ensemble by growing the configuration in a step-by-step manner and expressing it as a product of conditional probabilities, each of which is calculated by performing simulations at every step. This is a potentially exact method with the errors arising only from inadequate sampling. The method is applied to a small protein loop in implicit solvent and the results obtained agree quite well with those obtained from other (approximate) methods. Unlike other methods, all long range interactions are taken into account in the reconstruction procedure.
Giovanni Ciccotti (Università di Roma "La Sapienza") Constraints and coarse graining
Abstract: Holonomic constraints as coarse graining of (large) molecules. The dynamical evolution of a system subjected to holonomic constraints (SHAKE and beyond): Lagrange & Hamilton eq.s of first type; the standard (unstable) solution; Shake , i.e. an evolution satisfying exactly the constraints at every time step. The statistical measure and its associated evolution equation (Liouville equation) for mechanical systems with holonomic constraints: nonhamiltonian systems and their statistical behavior; the eq.s of motions of a mechanical systems subjected to constraints as a nonhamiltonian system; generalized Liouville eq. and statistical equilibrium probability.
Giovanni Ciccotti (Università di Roma "La Sapienza"), Mark E. Tuckerman (New York University) Rate constants
Abstract: The linear response theory formulation of a rate constan as a suitable time correlation function; its representation via the transition state theory (TST) term and the transmission coefficient. Calculation of the TST term via the blue-moon approach. i.e. impose a reaction coordinate as a constraint, compute its associated free energy and, by exponentiation, compute the associated probability. Use the unbiased constrained ensemble to generate a proper ensemble of reactive trajectories and compute the transmission coeffcient. The product of the TST term times the transmission coefficient provides the rate constant.
Matteo Cococcioni (University of Minnesota Twin Cities) Electronic correlations and Hubbard approaches
Abstract: Electronic correlation play a very important role in the physical properties of many materials characterized by very localized electronic states (as transition-metals or rare-earths compounds). Unfortunately, standard approximations to Density Functional Theory (like the Local Density Approximation or the Generalized Gradient Approximation) are not accurate enough to give a proper description of correlation effects. While multi-configurational Quantum-Chemistry methods are too expensive for most systems but very small molecules, several corrections to DFT have been introduced to give a better description of correlated states without giving up the relatively low computational cost of DFT calculations. After a general introduction on electronic correlation, one of the simplest corrective approaches, named LDA+U, will be described along with its implementation in a plane-wave pseudopotential DFT code.
Matteo Cococcioni (University of Minnesota Twin Cities), Ismaila Dabo (Massachusetts Institute of Technology) Tutorial on static calculations
Abstract: A hands-on tutorial on DFT total energy and force calculations will cover optimization of input parameters through convergence tests, electronic relaxation and structural optimization. The simulations will be performed using the plane-wave pseudopotential implementation of DFT contained in the Quantum-ESPRESSO package.
Matteo Cococcioni (University of Minnesota Twin Cities) A consistent, linear-response approach to LDA+U
Abstract: Hubbard U-corrected DFT functionals have been very successful in describing several strongly-correlated systems for which "standard" approximations to DFT fail. Unfortunately no explicit expression exists for the effective electronic interaction parameter (the Hubbard U) contained in the corrective ("+U") functional and semiempirical estimates have been often used to determine its value. In this talk, after a general introduction to the LDA+U method, I will present our linear response approach to the evaluation of the Hubbard U [1]. Within this approach the on-site electronic coupling is computed from the response of the considered system to a shift in the potential acting on its correlated atomic states. Specifically, it is evaluated as the difference between the inverse of the bare and fully interacting response matrices. The U we obtain thus corresponds to the effective (atomically averaged) interaction between electrons that are located on the same site. In this way the strength of the "+U" correction is consistently evaluated from the same DFT scheme we aim to correct; the LDA+U is transformed in a completely ab-initio method with no need for any empirical evaluation of the effective coupling. The results are also largely independent on the choice of the localized orbitals: the same occupation matrix that enters the expression of the "+U" correction is consistently used to compute the effective interaction parameter. With this approach we successfully studied the structural, electronic, chemical and electrochemical properties of several transition metals compounds. Examples of applications will include minerals in the Earth's interior [1], cathode materials for next-generation lithium-ion batteries [2] and catalysis reactions on molecules [3]. [1] M. Cococcioni and S. de Gironcoli, PRB (2005). [2] F. Zhou, M. Cococcioni, A. C. Marianetti, D. Morgan and G. Ceder, PRB (2004). [3] H. J. Kulik, M. Cococcioni, D. Scherlis and N. Marzari, PRL (2007).
Michel Cuendet (Swiss Institute of Bioinformatics) Aspects of nonautonomous molecular dynamics
Abstract: Performing constant temperature (and pressure) molecular dynamics simulation requires adequate thermostating (and barostating) algorithms, which couple additional degrees of freedom to the physical phase space variables. Thermostating is particularly important when external work is performed on the system in order to induce a nonequilibrium transition between two states. It appeared recently that nonequilibrium statistical mechanical relations such as the Jarzynski identity and the Crooks theorem naturally emerge as a property of nonautonomous thermostated dynamical systems. This however imposes an additional condition on the form of the thermostating dynamics, which is not required in the equilibrium case and is not fulfilled by all state-of-the-art thermostats.
Eric C. Cyr (University of Illinois at Urbana-Champaign) A comparison of maximum likelihood and weighted residual approximations to the potential of mean force
Abstract: The potential of mean force (PMF) describes the change in free energy along a reaction coordinate and determines the strength and likelihood of association in molecular systems. Current methods, like Thermodynamic Integration (TI) and WHAM, use piecewise linear and piecewise constant approximations of the PMF. We propose two methods that allow the use of arbitrary basis functions, one based on maximum likelihood estimation and another on the method of weighted residuals. Numerical comparisons with TI and WHAM are performed.
Ismaila Dabo (Massachusetts Institute of Technology) Real-space corrections for electrostatic interactions in periodic boundary conditions
Abstract: Joint work with Boris Kozinsky (Department of Physics, MIT), Nicholas E. Singh-Miller, and Nicola Marzari (Department of Materials Science and Engineering, MIT). We address periodic-image errors arising from the use of periodic boundary conditions to describe systems that do not exhibit full three- dimensional periodicity. We show that the difference between the periodic potential, straightforwardly obtained from a Fourier transform, and the exact potential can be characterized analytically. In light of this observation, we present an efficient real-space method to correct periodic-image errors, demonstrating that exponential convergence of the energy with respect to cell size can be achieved in practical periodic boundary-condition calculations. Comparing the method with existing schemes, we find that it is particularly advantageous for studying charged systems and systems exhibiting partial periodicity.
Eric Felix Darve (Stanford University) Adaptive methods for free energy computation and coarse-graining strategies
Abstract: Free energy computation is one of the main goals of molecular dynamics simulations. Numerically, these calculations are made difficult by the presence of free energy barriers separating meta-stable basins. Several methods have been developed to address this issue such as the adaptive biasing force (ABF), in which an external force is progressively adapted. Upon convergence, the reaction coordinate along which the free energy is computed exhibits a diffusive behavior in a flat free energy profile, leading to a rapid convergence of the computed averages and a rapid reduction of statistical errors. Such methods can serve as the starting point for computing coarse grained models of large proteins in which the number of degree of freedom is greatly reduced. We will show the role the free energy plays in formulating these models. In this case as well, ABF can be applied to yield faster convergence.
Ernest R. Davidson (University of Washington) Theoretical description of electrons in singloe molecule magnets
Abstract: Single molecule magnets are usually based on transition metals with partially filled d shells. When several metal centers are involved this leads to molecules with many single occupied orbitals coupled into an intermediate spin state. Conventional methods of quantum chemistry are not able to deal with this situation, so the Heisenberg model hamiltonian is often used with parameters estimated from DFT calculations. Practical as well as logical problems with this approach will be discussed. Some of the difficulties with treating exchange in DFT for open shell systems will be presented.
Kaushik Dayal (University of Minnesota Twin Cities) Objective structures and their applications
Abstract: Objective structures [1] are a recent classification of atomic and molecular structures that includes numerous systems of current interest. Examples are nanotubes and biological virus components. While these structures are not crystalline, there are many strong analogies to crystals. These analogies can be exploited to generalize computational techniques developed in the context of crystals to these new systems. Further, they may provide strategies for self-assembly and structure determination. Some examples of current work exploiting these ideas will be presented. These include finding the normal modes of these structures (in analogy with phonons), Objective Molecular Dynamics in analogy with conventional MD, and DFT and other first-principles techniques in the Objective setting. Recent work on finding all possible objective structures will be described. [1] R. D. James, Objective Structures, J. Mech. Phys. Solids, 54, 2354-2390 (2006).
Kaushik Dayal (University of Minnesota Twin Cities) Objective structures and their applications
Abstract: Same abstract as the 7/24 poster session.
Amélie Deleurence (École Nationale des Ponts-et-Chaussées (ENPC)) Modelling of local defects in crystals
Abstract: Work in collaboration with Eric Cancès (CERMICS-ENPC, France) and Mathieu Lewin (Department of Mathematics, Université de Cergy, France) We present mathematical results obtained for a new mean-field model dedicated to the description of interacting electrons in crystals with local defects. We work with a reduced Hartree-Fock model, obtained from the usual Hartree-Fock model by neglecting the exchange term. Then we deduce a new variational model for computing the perturbation in a basis of precomputed maximally localized Wannier functions of the reference perfect crystal. Finally we show some numerical results in which we have applied the variational approximation and used a one-dimensional model with Yukawa interaction potential.
Amélie Deleurence (École Nationale des Ponts-et-Chaussées (ENPC)) Modelling of local defects in crystals
Abstract: Same abstract as the 7/24 poster session.
Aaron Dinner (University of Chicago) Umbrella sampling for non-equilibrium processes
Abstract: Many systems of significant fundamental and applied interest are irreversible. These include, but are not limited to, living systems, chemical reactors, systems with driven flows of matter and energy, and photoactivated systems. For theoretical studies of such non-equilibrium processes, the steady-state distribution is of central importance because it enables calculation of static averages of observables for comparison to experimental measurements. For systems at equilibrium, low probability states can be explored efficiently in simulations with umbrella sampling methods, in which biasing potentials that are functions of one or more order parameters are used to enhance sampling of selected regions of phase space. What complicates extending umbrella sampling to simulations of non-equilibrium processes is that, by definition, they do not obey detailed balance (microscopic reversibility). As such, one must account for the fact that the steady-state probability of observing particular values of the order parameters can be determined by a balance of flows in phase space through different possible transitions. In this talk, I will describe the first algorithm for enforcing equal sampling of different regions of phase space in an ergodic system arbitrarily far from equilibrium, which enables its steady-state probability distribution to be determined with high accuracy. The efficiency of the algorithm will be demonstrated by applying it to a model of a genetic toggle switch which evolves irreversibly according to a continuous time Monte Carlo procedure.
Ron Elber (University of Texas) Milestoning, extending time scale of molecular simulations
Abstract: A new computational algorithm to extend time scales of atomically detailed simulations is illustrated. The algorithm (Milestoning) is based on partitioning the dynamics to a sequence of trajectories between “milestones” and constructing a non-Markovian model for the motion along a reaction coordinate. Besides "toy" models, two molecular systems are discussed: (i) The kinetics of a conformational transition in a blocked alanine is computed and shown to be accurate, more efficient than straightforward Molecular Dynamics by a factor of about 9, and non-exponential. (ii) The allosteric conformational transition in Scapharca hemoglobin is calculated. The results for the rate (about 10+/-9 microseconds) are in accord with experiment and are obtained about 1,000 times faster than straightforward Molecular Dynamics. No assumption of an activated process or states separated by significant barrier is made. A general scaling argument predicts linear speed-up with the number of milestones for diffusive processes, and exponential speed-up for transitions over barriers. The algorithm is also trivial to parallelize. As a side result Milestoning also produces the free energy profile along the reaction coordinate, and is able to describe non-equilibrium motions along one (or a few) degrees of freedom.
Jean-Luc Fattebert (Lawrence Livermore National Laboratory) Real-space finite difference method for O(N) first-principles molecular dynamics with plane waves accuracy
Abstract: Representing the electronic structure in Density Functional Theory (DFT) by a set of localized wave functions discretized on a real-space mesh essentially leads to a linear scaling of the computational cost with the size of the physical system. This can be achieved by formulating the DFT energy functional in terms of general non-orthogonal orbitals which are then optimized under localization constraints (spatial confinement). Multigrid preconditioning and a block version of Anderson's extrapolation scheme are used to accelerate convergence towards the ground state. For localization regions --- constraints --- large enough, one can reduce truncation error to a value smaller than discretization error and achieve the level of accuracy of a Plane Waves calculation. Accuracy is improved by allowing for flexible localization regions that can adapt to the system. This also reduces problems with local minima and enables energy conserving Born-Oppenheimer molecular dynamics simulations. Our implementation of this approach scales on hundreds of processors and becomes competitive with Plane Waves codes around 500 atoms. References:
[1] J.-L. Fattebert and F. Gygi, Phys. Rev. B 73, 115124 (2006)
[2] J.-L. Fattebert and F. Gygi, Comput. Phys. Comm. 162, 24 (2004) This work was performed under the auspices of the U.S. Department of Energy by University of California Lawrence Livermore National Laboratory under contract No. W-7405-Eng-48.
Michael Feig (Michigan State University) Efficient modeling of solvent environments
Abstract: The accurate representation of the solvent environment is an essential but also very costly component of biomolecular simulations. As a computationally less expensive alternative, implicit solvent models based on continuum electrostatic theory have become attractive alternatives to the standard explicit solvent models. Examples of biomolecular simulations with implicit aqueous solvent and implicit biological membranes are presented along with remaining mathematical challenges to fully take advantage of the much reduced model complexity.
Seyed-Alireza Ghasemi (Universität Basel) A systematic method to explore possible silicon tip structures used in AFM
Abstract: Same abstract as the 7/24 poster session.
Seyed-Alireza Ghasemi (Universität Basel) A systematic method to explore possible silicon tip structures used in AFM
Abstract: We present a systematic way to construct silicon tips by using the minima hopping method(MHM) with limited temperature. In Molecular Dynamics part of the MHM we limit the temperature to some value used in experiments. We used Tight-Binding scheme to evaluate energy and forces of silicon and hydrogen atoms in MHM. Structures with lowest energy obtained with this method are used in AFM simulation and force between tip and surface is calculated in terms of distance and some of obtained results are presented in this poster.
Simon Gill (University of Leicester) Thermal boundary conditions for molecular dynamics simulations
Abstract: A method for applying transient thermal boundary conditions to a molecular dynamics simulation is presented. The Kapitza effect at the boundaries is avoided by a diffuse feedback-controlled buffer zone. It is shown that this technique can be extended to provide a simple interface for coupling between concurrent atomistic and continuum simulations for ballistic heat transport.
Stefan Goedecker (Universität Basel) Wavelets for electronic structure calculations and electrostatic problems
Abstract: Wavelets are a systematic localized basis set that is well suited for representing Kohn-Sham orbitals. I will explain the algorithms that we are using in our new ABINIT wavelet program and show perfomance results. In the second part I will show how scaling functions can be used to solve electrostatic problems both for continuous and discrete charge distributions under various boundary conditions.
Gary B. Green (The Aerospace Corporation) Team 3: Associating earth-orbiting objects detected by astronomical telescopes
Abstract: Project description: Astronomical telescopes detect the passage of an earth-orbiting object as a streak in an image. Over a period of months, it is possible that many objects will pass through the field of view, some appearing more than once. There are estimates of 100,000 objects in orbit that might be detected by high resolution telescopes. A large field of view telescope may see 100 streaks a night. Most of these objects are space debris that pose a hazard to operational satellites. There is keen interest within the space community to discover and track all these objects.

If the telescope sensor is properly instrumented, it is possible to obtain time-tagged pairs of angles that relate the space object position to the sensor. With enough angle pairs, it is possible to estimate the position and velocity (the state) of the object, along with estimates of the uncertainties of these parameters. The workshop problem is to develop techniques to identify all the streaks made by each object. Streaks created by an object must somehow be associated with one another and disassociated from those made by other objects. One solution approach treats the state data as vectors in R6 and uses statistical clustering techniques for the association. A variation on this approach addresses physical properties of the orbits, sorting according to those least likely to change with small state variations.

Regardless of the approach, there are several interesting aspects to the problem. Automatic streak detection is required, with transform techniques of interest. Orbit mechanics are essential to effective state estimation as well as clustering techniques. In addition, traditional clustering techniques are computationally taxing. A related problem is identification of asteroids that might pose a hazard to planet earth. References:

Vallado, David A., Fundamentals of Astrodynamics and Applications, Edition 2, Microsoft Press, 2004; Milani, Andrea, "Three Short Lectures on Identifications and Orbit Determination," http://copernico.dm.unipi.it/~milani/preprints/preprint.html, 2006; Kaufman, L. and Rousseeuw, P., Finding Groups in Data - An Introduction to Cluster Analysis. Wiley Interscience 2005

Prerequisites:

Required: computing proficiency demonstrated by knowledge of at least one compiler, one semester differential equations, one semester statistics

Desired: one semester numerical analysis, familiarity with orbit mechanics and estimation theory.

Keywords: orbit mechanics, astronomical telescopes, statistical clustering

The Pan-starrs telescope on Mount Haleakela in Hawaii will be used, among other tasks, to search for asteroids. However, using its 1.4 billion pixel sensor, it will also detect earth-orbiting objects.

David J. Hardy (University of Illinois at Urbana-Champaign) Multilevel summation method for Coulomb interactions
Abstract: The multilevel summation method (MSM) computes an approximation to the pairwise interaction potential and its gradient with an amount of work that scales linearly as the size of the system. The potential is smoothly split into a sum of partial potentials of increasing range and decreasing variability with the longest-range parts interpolated from grids of increasing coarseness. Multilevel summation is especially appropriate for dynamical simulations, because it can produce continuous forces for both nonperiodic and periodic boundary conditions. A small correction to the grid-based force approximation can be used to conserve both energy and momentum.
Nicholas M. Harrison (Imperial College London) Novel materials for quantum computing
Abstract: The controlled transport of spin polarised electrons on a 1 nanometre length scale is a realistic prospect and could be the basis for new multifunctional devices with a component density an order of magnitude higher than current VLSI technology. The fundamental materials chemistry challenge is to produce a nano-structured semiconductor that is ferromagnetic at room temperature. Ideally the electronic and magnetic properties need to be robust but tunable through control of composition and structure. The results of recent high quality theoretical calculations on a number of pure carbon materials will be presented. A novel mechanism for long range magnetic coupling in extended pi-bonded systems will be discussed and documented with explicit calculations on graphene ribbons and defective graphene sheets. A putative ordered defect phase which gives rise to a semiconducting ground state that is ferromagnetic at room temperature will be presented. It will also be shown that the band gap and magnetic coupling may be controlled by varying the defect density.
Pieter Hendrickx (University of Ghent (UG)) Conformational reinvestigation of two cyclic pentapeptides: to a generic approach in drug development
Abstract: Same abstract as the 7/24 poster session.
Pieter Hendrickx (University of Ghent (UG)) Conformational reinvestigation of two cyclic pentapeptides: to a generic approach in drug development
Abstract: Joint work with I. Van den Eynde2, J.C. Martins1 and D. Tourwé2 As part of a bigger research effort which consists of the development of a generic approach for drug development, two cyclic peptides based on the active sequence of somatostatin were characterized using nuclear magnetic resonance spectroscopy (NMR). In order to obtain the 3D conformation of these molecules, conformational relevant parameters (eg. nOe, scalar couplings and amide temperature coefficients) were recorded. A simulated annealing protocol using distance restraints was used to obtain a set of conformations which were validated by analysis of the hydrogen bridges present and by recalculating the measured scalar couplings. Future efforts will consist of generating mathematical descriptions (templates) of this kind of small peptide chains as well as creating templates of non-peptidic small molecules. This must allow for an intelligent choice for drug templates upon positive screening results of the small peptide molecules. 1NMR and Structure Analysis Unit, Department of Organic Chemistry, Ghent University, Krijgslaan 281 S4 B-9000 Gent, Belgium 2Unit for Organic Studies, Department of Organic Chemistry, Vrije Universiteit Brussel, Pleinlaan 2 9G616 B-1050 Elsene, Belgium
John R. Hoffman (Lockheed Martin Missiles and Space Company, Inc.) Team 4: High dimensional, nonlinear, non-convex optimization problems in the area of aircraft and vehicle design
Abstract: Presently, when a physics motivated vehicle designer explores vehicle designs for a new concept, he is often faced with an enormous range of choices and constraints. For an example, an aircraft designer has Aircraft shape, fuel type, and engine as his main free variables. While his main constraints are dictated by the laws of physics (weight, size, power, lift, and stall). Additionally, he has his objective which is typically some combination/subset of acceleration, maneuverability, range, endurance, payload capacity (size, weight and power), max and min speeds, manufacturing cost, maintainability, reliability, development cost, takeoff length, landing length, noise footprint and other items.

I am interested in examining the following problem: Given a set of performance objectives, how does one determine the space of designs available to the designer and find the optimal designs? How does the designer best visualize this space of options? Because he doesn't want just "the" answer, he wants to understand many aspects of the answer. While I'm interested in the general vehicle design problem, we will focus on aircraft design using a baseline tool that is to be determined as a concrete example with which we can test our ideas.

Jürg Hutter (Universität Zürich) Efficient Kohn-Sham density functional calculations using the Gaussian and plane waves approach
Abstract: The Gaussian and plane waves (GPW) approach combines the description of the Kohn-Sham orbitals as a linear combination of Gaussian functions with a representation of the electron density in plane waves. The unique properties of Gaussian functions allow for a fast and accurate calculation of the density in the plane wave basis. The plane wave representation of the density leads to an easy solution of Poisson's equation and thereby a representation of the electrostatic potential. Matrix elements of this potential can be calculated using the same methods. The auxiliary representation of the density is further used in the calculation of the exchange-correlation energy and potential. The resulting approach scales O(N log N) in the number of electrons and has many additional interesting features, namely, a small prefactor, early onset of linear scaling, and a nominal quadratic scaling in the basis set size for fixed system size. The GPW method is combined with a direct optimization of the subspace of occupied Kohn-Sham orbitals using an orbital transformation (OT) method. A variation of this method has recently been implemented that only requires matrix multiplications. The method combines a small prefactor with efficient implementation on parallel computers, thereby shifting the break even point with linear scaling algorithms to much larger systems. A strategy to combine the OT method with sparse linear algebra will be outlined.
Olexandr Isayev (Jackson State University) An ab initio molecular dynamics simulation of solid CL-20: mechanism and kinetics of thermal decomposition
Abstract: Joint work with Leonid Gorb,(1 and 2) Mo. Qasim,(2) and Jerzy Leszczynski. (1) CL-20 (Octahydro-1,3,4,7,8,10-hexanitro-5,2,6-(iminomethenimino)-1H-imidazo[4,5-b]-pyrazin) is one of the most important high energetic nitramines which are used as explosives and propellants. The decomposition of CL-20 is very complicated and involves hundreds of elementary reactions. Accurate knowledge of these reactions and predictions of their kinetic parameters are important for modeling these complex processes in combustion and explosion. However, due to the energetic nature of these materials and the rapid rates of the intermediate reactions, it is difficult to monitor these individual reactions experimentally. Recent advances in first principles modeling have led to enormous progress toward understanding complex condensed phase chemical phenomena. Theoretical methods, especially accurate ab initio molecular dynamics method, provide a viable alternative to study the dynamics of these reactions. Electronic properties and the dynamics of the initial thermal decomposition step of gas-, and e-crystal phases of CL-20 are investigated using the Car-Parrinello molecular dynamics (CPMD) approach. The difference in the reaction pathways between gas- and crystal phases has been accounted for the strong intermolecular interaction into crystal lattice. The thermal rate constants have been also estimated. 1) Computational Center of Molecular Structure and Interactions, Jackson State University, Jackson, MS 39217
2) US Army ERDC, Vicksburg, MS, 39180.
Olexandr Isayev (Jackson State University) An ab initio molecular dynamics simulation of solid CL-20: mechanism and kinetics of thermal decomposition
Abstract: Same abstract as the 7/24 poster session.
Jesus A. Izaguirre (University of Notre Dame) Hybrid Monte Carlo methods
Abstract: Hybrid Monte Carlo, introduced in quantum chromodynamics, is a rigorous sampling method that eliminates the bias due to time step discretization of molecular dynamics (MD). I will discuss the method, give some example of its application, and show some recent improvements to its efficiency, particularly by discussing the Shadow Hybrid Monte Carlo method.
Raymond Kapral (University of Toronto) Quantum reaction rates
Abstract: The calculation of quantum reaction rates in condensed phase systems from the perspective of quantum-classical Liouville dynamics will be described. The rate problem will be cast in a form that involves quantum equilibrium sampling combined with quantum-classical dynamics for the evolution. Application to proton transfer in the condensed phase and polar molecule nanoclusters will be described. The results of simulations of nonadiabatic reaction rates using a surface-hopping scheme based on quantum-classical Liouville dynamics will be presented. The dynamic and structural factors that influence the rate and reaction mechanism will be discussed.
Changho Kim (Korea Advanced Institute of Science and Technology (KAIST)) Numerical method for solving stochastic differential equations with non-Gaussian noise
Abstract: We propose numerical integration schemes to solve stochastic differential equations driven by two important non-Gaussian noises – dichotomous Markov noise and Poissonian white shot noise. Our formula, which is based on an integral equation, which is equivalent to the stochastic differential equation, utilizes a discrete time approximation with fixed integration time step. We show that our integration formula approaches the Euler formula if these noises approach the Gaussian white noise. We further propose a simplified weak scheme that significantly reduces the computation time, while still satisfying the moment properties up to the required order. Our approach is readily applicable to dynamical systems driven by arbitrary types of noise, provided there exists a way to describe the random increment of the noise accumulated during each integration time step. The accuracy and efficiency of the proposed algorithms are numerically examined.
Changho Kim (Korea Advanced Institute of Science and Technology (KAIST)) Numerical method for solving stochastic differential equations with non-Gaussian noise
Abstract: Same abstract as the 7/24 poster session.
Brian Laird (University of Kansas) Direct calculation of interfacial free energy using molecular simulation
Abstract: The interfacial free energy is the work required to reversibly form a unit area of interface between two materials and is a principal quantity governing crystal growth kinetics and morphology, nucleation and wetting phenomena. In this talk I will review our recent work on the direct calculation of the interfacial free energy for a number of model systems including crystal/fluid, wall-fluid and wall-crystal systems. For the crystal-melt interface, our method of cleaving walls, based on thermodynamic integration [Phys. Rev. Lett, 85, 4751 (2000)], is of sufficient precision to resolve the often small anisotropy in the interfacial free energy, which is crucial in determining the kinetics and morpholgy of dendritic growth. Our calculations are motivated by the enormous difficulty in obtaining precise interfacial free energies by experiment. We report values of the crystal-melt interfacial free energy for the hard-sphere and Lennard-Jones systems, as well as recent results on the series of inverse-power potentials. In addition, recent work extending the method to wall-fluid and wall-crystal systems is discussed and applied to a recent controversy regarding surface-induced prefreezing in hard-spheres.
Claude Le Bris (École Nationale des Ponts-et-Chaussées (ENPC)) Coarse-graining the free energy of atomistic systems: a mathematical approach
Abstract: We present a possible approach for the computation of free energies and ensemble averages in the context of one-dimensional coarse-grained models in materials science. The approach is based upon a thermodynamic limit process, and makes use of ergodic type theorems and large deviation theory. In addition to providing a possible efficient computational strategy for ensemble averages, the approach in particular allows for assessing the validity and the accuracy of some basic assumptions of the finite-temperature quasicontinuum method. This is joint work with X. Blanc (Univ. Paris 6), F. Legoll (ENPC Paris), C. Patz (WIAS Berlin).
Frederic Legoll (École Nationale des Ponts-et-Chaussées) Integrators for highly oscillatory Hamiltonian systems: an homogenization approach
Abstract: We introduce a systematic way to construct symplectic schemes for the numerical integration of a large class of highly oscillatory Hamiltonian systems. The bottom line of our construction is to consider the Hamilton-Jacobi form of the Newton equations of motion, and to perform a two-scale expansion of the solution, for small times and high frequencies. Several options for the derivation are presented. The various integrators obtained are tested and compared to several existing algorithms. The numerical results demonstrate their efficiency. This is joint work with C. Le Bris (ENPC Paris).
Benedict Leimkuhler (University of Edinburgh) Dynamical equations and numerical integrators
Abstract: This lecture will introduce the basic types of dynamical models arising in molecular dynamics, and the formulation of structure-preserving numerical methods for integrating the equations of motion. Topics to be covered include: multiple time-scale modelling, the use of constraints, the computation of averages, the role of symplecticness and time-reversal symmetry, the properties of appropriate numerical methods, splitting methods, and SHAKE and RATTLE discretizations. The talk will conclude with a discussion of the use of backward error analysis to correct the numerical bias introduced in dynamical thermostatting methods.
Tony Lelievre (Ecole Nationale des Ponts et Chaussees) Adaptive methods for molecular dynamics