FUNCTIONALISATION, HYDRATION & DYNAMICS
|
|
- Horace Henderson
- 5 years ago
- Views:
Transcription
1 FUNCTIONALISATION, HYDRATION & DYNAMICS MATERIALS STUDIO TUTORIAL WAMMBAT - DAY 4
2 Our main goal: MD simulation of glucose in hydrated & functionalised silica nanopore It is not a rocket science and you can correct any outcome right away as it OCCURS. ~Axel Kohlmeyer
3 GAME PLAN 1. Introduction to Dynamics: Glucose dynamics in a water box. MSD Calculations: a) Preparation of structures & geometry optimization, Hydration in Amorphous Cell Package, NVT ensmble, b) Hydration (lower density), NPT ensemble, NVT ensemble, NVE ensemble and MSD calculation.
4 GAME PLAN 2. Glucose adsorption on functionalised silica surface: Structures preparation (silica surface, functionalisation), Hydration, NVT ensemble and MSD calculation.
5 GAME PLAN 3. Glucose dynamics in functionalised silica pore due to MSD calculation: Preparation of structures, - pore multiplication, - functionalisation, - glucose insertion, Hydration, Dynamics.
6 Before we start A bit of theory...
7 MSD
8 MSD <r 2 > T 2 T T 0.5
9 MSD caging regime diffusive regime balistic regime
10 MSD
11 MSD MSD Time...
12 MSD MSD Time...
13 MSD MSD Time...
14 MSD MSD Time...
15 MSD MSD Time NO DATA
16 MSD MSD Time NO DATA NO DATA
17 MSD MSD Time
18 MSD
19 MSD
20 CVFF The Consistent-Valence ForceField is a generalized valence forcefield. Parameters are provided for amino acids, water, and a variety of other functional groups.
21 CVFF pot = D b [1 e α b b 0 b ሿ 2 + H θ (θ θ 0 ሻ 2 + H φ [1 s cos nφ θ φ ൧ + H X X 2 X energy of deformation of bond lengths, bond angles, torsion angles and out-of-plane interactions + F bb b b 0 b b 0 + F θθ θ θ 0 θ θ 0 b b θ θ + F bθ b b 0 θ θ 0 + F φθθ cos θ θ 0 θ θ 0 + b θ φ X X F XX XX vibrational frequencies and the dynamic properties of molecules + σ ε [ r r 12-2 r r 6 ] + q iq j εr ij non-bond interactions
22 CVFF pot = D b [1 e α b b 0 b ሿ 2 + H θ (θ θ 0 ሻ 2 + H φ [1 s cos nφ θ φ ൧ + H X X 2 X + F bb b b 0 b b 0 + F θθ θ θ 0 θ θ 0 b b θ θ + F bθ b b 0 θ θ 0 + F φθθ cos θ θ 0 θ θ 0 + b θ φ X X F XX XX + σ ε [ r r 12-2 r r 6 ] + q iq j εr ij
23 COMPASS Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies is the first ab initio forcefield that enables accurate and simultaneous prediction of gasphase properties (structural, conformational, vibrational, etc.) and condensed-phase properties (equation of state, cohesive energies, etc.) for a broad range of molecules and polymers. It is also the first high quality forcefield to consolidate parameters of organic and inorganic materials.
24 COMPASS COMPASS COMPASS II Atom Types Forcefield Terms
25 MD software There are lots of MD programs, may not as graphically sophisticated as Materials Studio, but professional likewise. MD scientist's work is always divided into several steps and often executed with different software, for example: - building crystal (with or without the pore), - hydration, - optimalisation, - dynamics, - results analysis and presentation (visualisation). Each program propose different opportunities and in some way is hardly replaceable in it's special field.
26 MD software Materials Studio allows user to perform all the steps, but sometimes other specific parameters or functions are necessary and usage of more adjustable software is vital. List of paramters that might be set and possibilities of changes is often more available than in commercial software, but user is demanded to know much more about how program works. Major part of MD software is OpenSource and required changes are possible straight from the source. fragment of LAMMPS source code
27 Geometry optimization Each system must be optimized before your proper simulation. Every algorithm leads potential energy function to its minimum. You must be aweare that potential energy is rather multidimensional hypersurface than 3D function (seen on the pic.). Initial energy Global or local minimum
28 Geometry optimization Example of convergence criteria: 1. Convergence of the function (energy) f x k+1 f(x k ሻ < ε f 2. Convergence of the variable (position) x k+1 x k < ε x 3. Covergence of the gradient (force) Initial energy Global or local minimum f x k+1 f(x k ሻ < ε g
29 Geometry optimization Examples of most popular energy minimization algorithms: - conjugate gradients, - steepest descent. And also some damping dynamics algorithms such as: - quickmin, - fire. Initial energy Global or local minimum
30 Potential energy Geometry optimization 60º 109.5º 145º 180º 215º 215.5º 300º Angle
31 PART 1 - Introduction The first task is to create the whole system: molecule in water. Usually it is a very time-consuming step during the whole research process. This example is very simple, but imagine some more complicated biological systems as presented on the picture. Each step may be done in few different approaches. One of them is our goal: hydration.
32 Periodic Boundary Conditions
33 Goal: glucose in water no PBC Glucose PBC addition Water addition
34 Goal: glucose in water no PBC No problems with non-bonded interactions Huge number of atoms Expensive PBC Atom may interact with itself Small number of atoms Cheaper
35 Example: Hydration by packmol Packmol is an opensource project developed for quick biosystem building. It does not have it's own graphical interface and should be run in system shell. Our group logs into BEM cluster and performs calculations there in command line by Secure Shell (SSH). During WAMMBAT we have already used BEM cluster by Materials Studio gateway due to the hospitality of WCSS team.
36 Example: Hydration by packmol How it works? 1. User defines box sizes and number of molecules (each kind). 2. Packmol guesses positions of all the molecules as they were balls of selected radius (tolerance). 3. Packmol tries to calculate better coordinates to fit all the particles into selected box and stops while tolerance is reached. Final structure is set of geometrically fitted molecules, so it's geometry should be optimized with proper forcefield.
37 After packmol hydration Geometry Optimization includes forces and optimize geometry NPT ensemble density (volume) equilibration NVT temperaturę equilibration
38 Example: Amorphous Cell - Packing geometry optimization + next step geometry optimization + next step
39 Example: Amorphous Cell - Packing after a while..
40 Pros and cons of both methods parameter Packmol Amorphous Cell - Packing Handling Difficulty low-medium very easy Graphical Interface text-mode (unix systems) module built in MS Computation Efficiency very efficient even in really huge systems only small systems Geometry Optimization no automatic Density Control user defines number of molecules --> system needs equilibration of density automatic
41 Queueing An imprortant issue is PBS - the queue system. Any job before running on BEM (or any other cluster) has to be queued. User must define computing resources (number of parallel processors, RAM memory, time of reservation) and some additional parametrs, then is given specified priority by the queue system. Materials Studio gateway does it automatically.
42 Queueing User Gathering results Queue system Hardware
43 Queueing
44 Queueing User Job Priority Queue Hardware - Memory - Nodes - Max time - length - available slots - resources (nodes, memory) - user priority
45 Queueing 12 nodes 1/12 in use Requested resources 12 nodes Find available slots 11/12 3/12 8/12 1/12 12/12 12/12 12/12 12/12 12/12 12/12 12/12 12/12 12/12 12/12 12/12 7/12 12/12 12/12 12/12 12/12
46 WCSS
47 Goal: glucose in water Glucose PBC addition Water addition
48 Clean up the scene Tousands of constantly opening windows during work in Materials Studio may be very annoying, so please remember to uncheck Automatically view output checkbox in more menu in Job Control tab before you start any job.
49 Clean up the scene The most important issue is not to be lost in data mess. First of all we'll manage all the files created through the last days. - Open Materials Studio and your wammbat project. - Create a new folder in the main location: Day 4 - Dynamics and Create new folders Task 1, Task 2, Task 3 inside it.
50 Clean up the scene Find and paste proper files from last days: Task 1 Task 2 Task 3 glucose (Day 1) glucose+.xsd (Day 3) glucose+.xsd (Day 3) Si-2NH3.xsd (Day 3) functionalised silica pore (Day 2)
51 Clean up the scene For both crystal structures change representations to Constraint and check if crystal part is fixed. functionalised silica pore (Day 2) Si-2NH3.xsd (Day 3)
52 Structures preparation Open a new 3D Atomistic document in Task1 folder. Name it water cvff. Draw Oxygen atom using Sketch atom and update hydrogens.
53 Structures preparation Before water molecule will be ready to use it should be optimized in forcefield that will be used afterwards. Set Geometry Optimization task and More... options as presented. Open Forcite Calculation.
54 Structures preparation Materials Studio modules (Forcite, Amorphous Cell and others) "remember" the last used properties. You can open different document and be sure the settings are set exactly the same as in the previous document. If you see: "set as presented" on the tutorial page, it means for 90% that we used this module before and it does not require any changes.
55 Structures preparation In Energy tab set cvff Forcefield, Forcefield assigned charges, Medium quality, Atom based Electrostatic and van der Waals. In Job Control tab set one core.
56 Downloading your results To download results after a job is completed rightclick on the job, select Actions --> Download Results
57 Structures preparation Without closing Forcite window open glucose molecule and RUN optimization again.
58 Structures preparation You can rename and move optimized Water and glucose documents to the Task1 folder and delete other files.
59 Structures preparation To build a box: Open Build --> Crystal --> Build Crystal... Set Lengths in Lattice Parameters tab to 20 Angstroms and click Build. Concentration of glucose in this box: 0.2 M. C = n V = mol dm 3 = M
60 Hydration in Amorphous Cell To pack water inside the box, open: Amorphous Cell --> Calculation module. Select Packing task, Density to 1,0 g/cc, Output to 4 frames and set water opt as a Molecule. Open Options menu.
61 Hydration in Amorphous Cell Set options as presented. In previous menu you set 4 frames for output. Sometimes system after packing may have extreme energy and explode. Usually temperature, pressure and density is reached gradually to prevent such situations. In this tutorial we try to save time and prepare 4 structures in case of explosion.
62 Hydration in Amorphous Cell In Energy tab set Ewald in Electrostatic summation During this task MS cannot divide your box into a few smaller boxes, so you cannot use more cores.
63 NVT ensemble System is ready to carry on the dynamics calculations. Open Forcite Calculation. Choose Dynamics and open More... menu.
64 NVT ensemble In Dynamics tab set: NVT ensemble, Random velocities, 298 K, 1,0 fs time step and 130 ps total simulation time. To obtain a detailed plot, set an output frame every 200 steps. In Thermostat tab choose Nose thermostat.
65 NVT ensemble Go back to main menu, and in Energy tab set PPPM electrostatic summation and Ewald van der Waals. Remember to set maximum available core number in Job Control properties. Then click Run.
66 Ewald summation method Ewald summation is an algorithm wich effciently calculates the interactions between periodic mirrors of atoms or molecules. Van der Waals interactions deacys like 1 r 6 so there is no porblem usually periodic mirrors are over the cutoff. In case of Coulomb electrostatic interactions k q iq i_mirror r i i_mirror potential deacys slow. Ewald proposed 1 + which states for short range, rapidly varying function r r and long range flat function. The f r = erfc r = 2 π x exp t 2 dt. = f r r 1 f r
67 Ewald summation method Reciprocal space Real space r r Each charge is neutralised with neutralising charge Another contribution counteracts neutralisation Real space r
68 Particle-particle particle-mesh In the particle-particle particle-mesh method (PPPM, P 3 M) some simplifications are made. Beyond some distance R c particles (atoms) are discetized charge density gives a contribution to neighbouring grid nodes. For close particels (atoms) interactions are calcualted by common Ewald summation. Mesh grid
69 Ewald vs. PPPM (P3M) Ewald summation Summation over every atom Uses Fourier transform Particle Particle Particle mesh Summation over a discretized grid of charge denisty Uses fast Fourier transform Faster for smaller systems Significantly faster for large systems (N > 10000ሻ # Atoms Ewald (c. time) P3M (c. time)
70 But in real work...
71 Dynamics Calculation stages to define temperature and density NPT NVT NVE equilibration of density (volume changes) equilibration of temperature final dynamics
72 Thermostating & Barostating Simulation system is dynamical. In order to preserve the system specific temperature or pressure we employ in the calculation the concept of thermostat and barostat. Conceptionaly these are separate systems with defined pressure an temperature connected to our simulation balls. Mathematicaly we simply add two Newtons equations, additionaly equations defining temperature and pressure.
73 Berendsen Thermostat The velocities are scaled at each step, such that the rate of change of temperature is proportional to the difference in temperature. Scaling factor:
74 Berendsen Barostat. In barostat the mechanism is similar: Scaling factor:
75 NPT V nt/p During simulation ensemble uses barostat and thermostat to equilibrate temperature and pressure. Varying, operating both: P and T results in changes of volume of the simulation box. In other words, our goal is to reach density of water about 1.0 g/cc, that is a property of our system in 298 K and 1 bar.
76 NPT Box length running average
77 NPT Density
78 NVT High kinetic energy equilibrated temperature Running average Random Velocities High potential Energy Temperature Strong thermostat work
79 NVE
80 New system In Task 1 folder open glucose cvff.xsd document. Prepare Packing exactly the same as last time except density, set it to 0.8 g/cc. Click Run.
81 NPT ensemble System is ready to carry on dynamics. Open Forcite Calculation. Choose Dynamics and open More... menu.
82 NPT ensemble In Dynamics tab choose NPT ensemble, 1,013e-4 GPa pressure, 1,0 fs time step, 8 ps total time and output every 20 steps. Thermostat tab leave without changes. Set Berendsen barostat.
83 NPT ensemble Go back to previous menu and set PPPM Electrostatic and Ewald vdv. Remember to set maximum core number in Job Control tab. Run dynamics.
84 NPT ensemble Open plots and check if simulation time was sufficient to equilibrate density in your box.
85 NVT ensemble Open Forcite Calculation. Choose Dynamics and open More... menu. Check Restart box. Forcite will take final velocities from NPT dynamics.
86 NVT ensemble In Dynamics tab set NVT ensemble, 8,0 ps total simulation time and output every 20 frames. Thermostat leave without changes. Move back to the main Forcite menu.
87 NVT ensemble Rest of the tabs leave without changes and Run dynamics.
88 NVT ensemble Open plots and check if simulation time was sufficient to equilibrate temperature inside your box.
89 NVE ensemble Open Forcite Calculation. Choose Dynamics and open More... menu.
90 NVE ensemble As you see, in NPT we set Dynamics, Thermostat and Barostat properties. In NVT there was no Barostat, and now in NVE ensemble there is neither Barostat nor Thermostat. Set NVE ensemble and 130ps in total simulation time. Set output every 140 steps.
91 NVE ensemble Rest of the tabs leave without changes and Run dynamics.
92 Calculate MSD Materials Studio store all data about dynamics, so different calculations might be performed any time it is needed. Now, we will perform MSD analyse for both tasks. Make sure if you have opened proper xtd document before you start analysing.
93 Calculate MSD Materials Studio store all data about dynamics, so different calculations might be performed any time it is needed. First of all select glucose by shift+double click. Then Edit --> Sets
94 Calculate MSD Make sure you have opened proper file and open Forcite --> Analysis and choose mean square displacement (MSD) Set origin to the begining of simulation and using arrows set length frames to about the 2/5 of simulation time. Now Analyze.
95 Calculate MSD Diffusive regime
96 Calculate MSD Calculations resulted with two new documents - plot and data table. Open Table and select B column. Open Statistics --> --> Model Building --> Multiple Linear Regression
97 Calculate MSD Select Time and click OK. Here is your result.
98 Calculate MSD D[ m2 A 2 your result [ s ሿ = ps ሿ
99 Calculate MSD Do it again with second simulation results!
100 You didn't receive ? :(
101 Confined Diffusion
102 Diffusion parameter - concentration - thermostating - forcefield - system size - randomness in three dimensions - temperature - accuracy -...
103 Diffusion Parameter Really, really short simulation Proper order of magnitude! Result similar to reality!
104 It's high time we started our second task!
105 PART 2 - Glucose adsorption We will repeat Monte Carlo project using Dynamics instead. We will also hydrate our system. Due to the periodic boundary conditions and very small scale we can simplify the pore to the silica layer with water above.
106 Goal: Box with glucose and water Box Box + glucose Water addition
107 Prepare silica layer There are some structures in Task 2 folder already: - glucose+.xsd, - Si-2NH3.xsd.
108 Structures preparation Open a new 3D Atomistic document in Task2 folder. Draw Oxygen atom using Sketch atom and update hydrogens.
109 Structures preparation Open Forcite Calculation. Set Geometry Optimization task.
110 Structures preparation In Energy tab set COMPASS II Forcefield.
111 Structures preparation Paste glucose into the box and move it to a proper destination using toolbar.
112 Solvent surface 1A ball radius 2A radius vdw surface vdw radius 2A solvation surface
113 Solvent Surface exclude too-close (high energy) locations exclude crystal inside
114 Prepare hydration surface To create hydration surface open: Tools --> Atom Volumes & Surfaces and choose Task: Solvent surfaces. Then Create surface.
115 Prepare hydration surface The last step is to inverse hydration surface to set high values inside the surface, that results afterwards in packing the surface's outher part. Select isosurface bu double click, right click and choose Display Style --> Isosurface tab Select High values inside checkbox.
116 Pack water inside the box Amorphous Cell Packing module packs molecules one after another and optimalise whole system geometry. To pack water inside the box, open: Amorphous Cell --> Calculation module. Select Packing task, Density as 1,0 g/cc, and water Molecule. Open Task --> more menu and check Pack in isosurface enclosed volume. Open Options menu.
117 Pack water inside the box Packing parameters are probably set by default, but make sure they are set correctly. Open Options menu and select: parameter Torsions Temerature value calculate automatically 298 K Loading steps 1000 Check energies Optimize geometry select select
118 Pack water inside the box Energy and Job Control tabs are already set. You should only change forcefield to COMPASS II. Now you can Run packing.
119 Pack water inside the box Select isosurface and delete it. It is no longer neccesary. Now, you are ready to run dynamics calculations.
120 NVT ensemble Open Forcite Calculation. Choose Dynamics and open More... menu.
121 NVT ensemble To observe diffusion the simulation time has to be of the order of nanoseconds. Now set the Total simulation time as 5,0 picoseconds to examine whole process. Afterwards, calculations will continue for the whole night. To equilibrate temperature time about 120ps is far enough. Afterwards, dynamics will be continued for 1ns in NVE ensamble and results of the analysis will be conducted tommorow.
122 NVT ensemble Go back to previous menu, check Energy tab and Run dynamics.
123 Dynamics As you noticed before, each job has it's own folder. Check what files they contain. Open the temperature plot. It presents that the thermostat works in a proper way. The temperature optimization will be continued. To start it please use prepared script that will run NVT and NVE simulations one after another.
124 Continue dynamics We equilibrated temperature and system is ready for NVE ensemble, but the rest of dynamics we will continue in the end of today laboratories.
125 PART 3 - Glucose dynamics The last task is to functionalise and hydrate complete silica pore. The dynamics, that we will start today will continue the whole night and will be analysed tomorrow.
126 Goal: Pore with glucose and water silica preparation functionalization water addition
127 Clean up the scene again Check if you recognise all the subfolders names in your catalog tree. Change names if necessary. The way you order your project should be suitable for easy navigation for you, but we recommend to follow us:
128 Put glucose inside the pore To put glucose inside the pore select deprotonated glucose, that was optimized in COMPASSII forcefield and paste it into your pore. Select glucose by shift+click and move using moving toolbar. -NH3 -COO -NH3 -COO Be careful to put glucose in between ("z" direction) NH3 groups and COO- groups and in the center of the pore tunnel. : -NH3 (the same)
129 Prepare hydration surface It is time we did hydration surface exactly the same as in Task2: Open functionalised silica document from Task 3 folder. Open Tools --> Atom Volumes & Surfaces and choose Solvent surfaces task. Then Create surface.
130 Prepare hydration surface The last step is to inverse hydration surface to set high values inside the surface, that results afterwards in packing the surface's outher part. Select isosurface bu double click, right click and choose Display Style - -> Isosurface tab Select High values inside checkbox.
131 Pack water inside the box To pack water inside the box, open: Amorphous Cell --> Calculation module. Select Packing task, Density as 1,0 g/cc, and water Molecule. Open Task --> more menu and check Pack in isosurface enclosed volume.
132 Pack water inside the box Open Task --> more menu and check Pack in isosurface enclosed volume. In Energy tab set COMPASS II Forcefield, Charges Forcefield assigned, Quality Medium and Ewald Electrostatic summation method. Now you can Run packing. Remember to delete isosurface while packing is finished.
133 Dynamics Our structure is to big to be calculated as fast as before. The dynamics will continue for the whole night. In Forcite Calculation --> Dynamics set NVE ensemble, 1 fs timestep, 500 ps total time, and frame output every 500 frames. In Job Control tab select 2 cores. Now, open Task 2 ready structure. In Dynamics options set 1000ps total time and Run again.
134 CAUTION! Our calculations will be continued so please DO NOT log out or shut down your computer till tomorrow!
135 To be continued...
Silica surface - Materials Studio tutorial. CREATING SiO 2 SURFACE
Silica surface - Materials Studio tutorial CREATING SiO 2 SURFACE Our goal surface of SiO2 6.948 Ǻ Import structure The XRD experiment gives us such parameters as: lattice parameters, symmetry group and
More informationMolecular modeling for functional materials design in chemical engineering and more
Molecular modeling for functional materials design in chemical engineering and more Dr inż. Łukasz Radosiński Group of Bioprocess and Biomedical Engineering Wroclaw University of Science and Technology
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
More informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
More informationAPBS electrostatics in VMD - Software. APBS! >!Examples! >!Visualization! >! Contents
Software Search this site Home Announcements An update on mailing lists APBS 1.2.0 released APBS 1.2.1 released APBS 1.3 released New APBS 1.3 Windows Installer PDB2PQR 1.7.1 released PDB2PQR 1.8 released
More informationChem 253. Tutorial for Materials Studio
Chem 253 Tutorial for Materials Studio This tutorial is designed to introduce Materials Studio 7.0, which is a program used for modeling and simulating materials for predicting and rationalizing structure
More informationWinmostar tutorial LAMMPS Melting point V X-Ability Co,. Ltd. 2017/8/17
Winmostar tutorial LAMMPS Melting point V7.025 X-Ability Co,. Ltd. question@winmostar.com Contents Configure I. Build solid phase II. Equilibration of solid phase III. Equilibration of liquid phase IV.
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationExercises for Windows
Exercises for Windows CAChe User Interface for Windows Select tool Application window Document window (workspace) Style bar Tool palette Select entire molecule Select Similar Group Select Atom tool Rotate
More informationFROM ATOM TO POLYMER. creating and manipulating models of chemical compounds In BIOVIA Materials studio
FROM ATOM TO POLYMER creating and manipulating models of chemical compounds In BIOVIA Materials studio What can be built in Materials Studio? Single structures What can be built in Materials Studio? Homopolymers
More informationJASS Modeling and visualization of molecular dynamic processes
JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St Petersburg State University, Physics faculty, Department of Computational Physics Supervisor PhD Stepanova Margarita
More informationWinmostar tutorial LAMMPS Polymer modeling V X-Ability Co,. Ltd. 2017/7/6
Winmostar tutorial LAMMPS Polymer modeling V7.021 X-Ability Co,. Ltd. question@winmostar.com 2017/7/6 Contents Configure I. Register a monomer II. Define a polymer III. Build a simulation cell IV. Execute
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationAthena Visual Software, Inc. 1
Athena Visual Studio Visual Kinetics Tutorial VisualKinetics is an integrated tool within the Athena Visual Studio software environment, which allows scientists and engineers to simulate the dynamic behavior
More informationDiffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study
Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron
More informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationCreating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner
Table of Contents Creating a Pharmacophore Query from a Reference Molecule & Scaffold Hopping in CSD-CrossMiner Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Features
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More information4 Materials and methods
14 Materials and methods 4 Materials and methods This chapter gives an overview on the materials and methods used in this work: In Chapter 4.1, the raw materials are presented as well as the crystallization
More informationCPMD Tutorial Atosim/RFCT 2009/10
These exercices were inspired by the CPMD Tutorial of Axel Kohlmeyer http://www.theochem.ruhruni-bochum.de/ axel.kohlmeyer/cpmd-tutor/index.html and by other tutorials. Here is a summary of what we will
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationWinmostar tutorial Gromacs Interfacial Tension V X-Ability Co., Ltd. 2018/01/15
Winmostar tutorial Gromacs Interfacial Tension V8.007 X-Ability Co., Ltd. question@winmostar.com Summary In this tutorial, we will calculate density distribution and interfacial tension between liquid-liquid
More informationIntroduction to Spark
1 As you become familiar or continue to explore the Cresset technology and software applications, we encourage you to look through the user manual. This is accessible from the Help menu. However, don t
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationAppendix A Prototypes Models
Appendix A Prototypes Models This appendix describes the model of the prototypes used in Chap. 3. These mathematical models can also be found in the Student Handout by Quanser. A.1 The QUANSER SRV-02 Setup
More informationMultiphysics Modeling
11 Multiphysics Modeling This chapter covers the use of FEMLAB for multiphysics modeling and coupled-field analyses. It first describes the various ways of building multiphysics models. Then a step-by-step
More informationON SITE SYSTEMS Chemical Safety Assistant
ON SITE SYSTEMS Chemical Safety Assistant CS ASSISTANT WEB USERS MANUAL On Site Systems 23 N. Gore Ave. Suite 200 St. Louis, MO 63119 Phone 314-963-9934 Fax 314-963-9281 Table of Contents INTRODUCTION
More informationMaterials that you may find helpful when working through this exercise
Detailed steps illustrating how to use VASP on the Suns in Fitz 177 For use in lab: 11/10/2009 (Original file by Dr. Rachel Getman, 11/18/2007. Editted for use by Dorrell McCalman 11/09/2009.) Note on
More informationApplications of Molecular Dynamics
June 4, 0 Molecular Modeling and Simulation Applications of Molecular Dynamics Agricultural Bioinformatics Research Unit, Graduate School of Agricultural and Life Sciences, The University of Tokyo Tohru
More information2D Clausius Clapeyron equation
2D Clausius Clapeyron equation 1/14 Aim: Verify the Clausius Clapeyron equation by simulations of a 2D model of matter Model: 8-4 type potential ( Lennard-Jones in 2D) u(r) = 1 r 8 1 r 4 Hard attractive
More information2D Clausius Clapeyron equation
2D Clausius Clapeyron equation 1/14 Aim: Verify the Clausius Clapeyron equation by simulations of a 2D model of matter Model: 8-4 type potential ( Lennard-Jones in 2D) u(r) = 1 r 8 1 r 4 Hard attractive
More informationMERGING (MERGE / MOSAIC) GEOSPATIAL DATA
This help guide describes how to merge two or more feature classes (vector) or rasters into one single feature class or raster dataset. The Merge Tool The Merge Tool combines input features from input
More informationQuantumWise. QuantumWise is now part of Synopsys
Table of Contents Table of Contents NiSi2 Si interface Create the NiSi2/Si device The screening region Increase the length of the central region Set-up the calculation for the undoped device Dope the device
More informationUniversal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information
Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information B. Shadrack Jabes, Dusan Bratko, and Alenka Luzar Department of Chemistry,
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationPhysics E-1ax, Fall 2014 Experiment 3. Experiment 3: Force. 2. Find your center of mass by balancing yourself on two force plates.
Learning Goals Experiment 3: Force After you finish this lab, you will be able to: 1. Use Logger Pro to analyze video and calculate position, velocity, and acceleration. 2. Find your center of mass by
More informationST-Links. SpatialKit. Version 3.0.x. For ArcMap. ArcMap Extension for Directly Connecting to Spatial Databases. ST-Links Corporation.
ST-Links SpatialKit For ArcMap Version 3.0.x ArcMap Extension for Directly Connecting to Spatial Databases ST-Links Corporation www.st-links.com 2012 Contents Introduction... 3 Installation... 3 Database
More informationThe Geodatabase Working with Spatial Analyst. Calculating Elevation and Slope Values for Forested Roads, Streams, and Stands.
GIS LAB 7 The Geodatabase Working with Spatial Analyst. Calculating Elevation and Slope Values for Forested Roads, Streams, and Stands. This lab will ask you to work with the Spatial Analyst extension.
More information2008 Biowerkzeug Ltd.
2008 Biowerkzeug Ltd. 1 Contents Summary...3 1 Simulation...4 1.1 Setup...4 1.2 Output...4 2 Settings...5 3 Analysis...9 3.1 Setup...9 3.2 Input options...9 3.3 Descriptions...10 Please note that we cannot
More informationTutorial 12 Excess Pore Pressure (B-bar method) Undrained loading (B-bar method) Initial pore pressure Excess pore pressure
Tutorial 12 Excess Pore Pressure (B-bar method) Undrained loading (B-bar method) Initial pore pressure Excess pore pressure Introduction This tutorial will demonstrate the Excess Pore Pressure (Undrained
More informationCalculating Bond Enthalpies of the Hydrides
Proposed Exercise for the General Chemistry Section of the Teaching with Cache Workbook: Calculating Bond Enthalpies of the Hydrides Contributed by James Foresman, Rachel Fogle, and Jeremy Beck, York College
More informationIntroduc)on to IQmol: Part I.!!! Shirin Faraji, Ilya Kaliman, and Anna Krylov
Introduc)on to IQmol: Part I!!! Shirin Faraji, Ilya Kaliman, and Anna Krylov! 1 Resources! Written by Dr. Andrew Gilbert Keep yourself up to date with IQmol website: http://iqmol.org! IQmol Youtube channel:
More informationModel Building An Introduction to Atomistic Simulation
Materials and Modelling MPhil 2006-07 COURSE MP3: MONTE CARLO AND MOLECULAR DYNAMICS COMPUTING CLASS 1 Model Building An Introduction to Atomistic Simulation Wednesday 22 nd November 2006 14.00 16.00 1
More informationPerforming a Pharmacophore Search using CSD-CrossMiner
Table of Contents Introduction... 2 CSD-CrossMiner Terminology... 2 Overview of CSD-CrossMiner... 3 Searching with a Pharmacophore... 4 Performing a Pharmacophore Search using CSD-CrossMiner Version 2.0
More informationTutorial. Getting started. Sample to Insight. March 31, 2016
Getting started March 31, 2016 Sample to Insight CLC bio, a QIAGEN Company Silkeborgvej 2 Prismet 8000 Aarhus C Denmark Telephone: +45 70 22 32 44 www.clcbio.com support-clcbio@qiagen.com Getting started
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationWinmostar tutorial LAMMPS Polymer Annealing V X-Ability Co., Ltd. 2018/01/15
Winmostar tutorial LAMMPS Polymer Annealing V8.007 X-Ability Co., Ltd. question@winmostar.com 2018/01/15 Summary In this tutorial we will calculate glass transition temperature from the cooling process
More informationIntroduction to Hartree-Fock calculations in Spartan
EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,
More informationCHARTING THE HEAVENS USING A VIRTUAL PLANETARIUM
Name Partner(s) Section Date CHARTING THE HEAVENS USING A VIRTUAL PLANETARIUM You have had the opportunity to look at two different tools to display the night sky, the celestial sphere and the star chart.
More informationGeneral Chemistry Lab Molecular Modeling
PURPOSE The objectives of this experiment are PROCEDURE General Chemistry Lab Molecular Modeling To learn how to use molecular modeling software, a commonly used tool in chemical research and industry.
More informationv Prerequisite Tutorials GSSHA WMS Basics Watershed Delineation using DEMs and 2D Grid Generation Time minutes
v. 10.1 WMS 10.1 Tutorial GSSHA WMS Basics Creating Feature Objects and Mapping Attributes to the 2D Grid Populate hydrologic parameters in a GSSHA model using land use and soil data Objectives This tutorial
More informationLearning to Use Scigress Wagner, Eugene P. (revised May 15, 2018)
Learning to Use Scigress Wagner, Eugene P. (revised May 15, 2018) Abstract Students are introduced to basic features of Scigress by building molecules and performing calculations on them using semi-empirical
More informationTable of Contents. Table of Contents Exploring graphene. Build a graphene sheet Build a CNT Transmission spectrum of a GNR Twisted nanoribbon
Table of Contents Table of Contents Exploring graphene Build a graphene sheet Build a CNT Transmission spectrum of a GNR Twisted nanoribbon Transmission spectrum Möbius nanoribbon Buckling a graphene sheet
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationMolecular Modeling and Conformational Analysis with PC Spartan
Molecular Modeling and Conformational Analysis with PC Spartan Introduction Molecular modeling can be done in a variety of ways, from using simple hand-held models to doing sophisticated calculations on
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
More informationProject 4/5 - Molecular dynamics part II: advanced study of Lennard-Jones fluids, deadline December 1st (noon)
Format for delivery of report and programs The format of the project is that of a printed file or hand-written report. The programs should also be included with the report. Write only your candidate number
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationComputer Simulation of Shock Waves in Condensed Matter. Matthew R. Farrow 2 November 2007
Computer Simulation of Shock Waves in Condensed Matter Matthew R. Farrow 2 November 2007 Outline of talk Shock wave theory Results Conclusion Computer simulation of shock waves Shock Wave Theory Shock
More informationAssignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans
Assignment A02: Geometry Definition: File Formats, Redundant Coordinates, PES Scans In Assignments A00 and A01, you familiarized yourself with GaussView and G09W, you learned the basics about input (GJF)
More informationhydrated Nafion-117 for fuel cell application
A molecular dynamics simulation study of oxygen within hydrated Nafion-117 for fuel cell application Jeffrey P. Fuller, Giuseppe F. Brunello, Seung Soon Jang School of Materials Science and Engineering
More informationTutorial I: IQ MOL and Basic DFT and MP2 Calculations 1 / 30
Tutorial I: IQ MOL and Basic DFT and MP2 Calculations Q-Chem User Workshop, Denver March 21, 2015 1 / 30 2 / 30 Introduction to IQMOL DFT and MP2 Calculations 3 / 30 IQMOL and Q-CHEM IQMOL is an open-source
More informationOPSIAL Manual. v Xiaofeng Tan. All Rights Reserved
OPSIAL Manual v1.0 2016 Xiaofeng Tan. All Rights Reserved 1. Introduction... 3 1.1 Spectral Calculator & Fitter (SCF)... 3 1.2 Automated Analyzer (AA)... 3 2. Working Principles and Workflows of OPSIAL...
More informationWORKING WITH DMTI DIGITAL ELEVATION MODELS (DEM)
WORKING WITH DMTI DIGITAL ELEVATION MODELS (DEM) Contents (Ctrl-Click to jump to a specific page) Manipulating the DEM Step 1: Finding the DEM Tiles You Need... 2 Step 2: Importing the DEM Tiles into ArcMap...
More informationComputational Chemistry Lab Module: Conformational Analysis of Alkanes
Introduction Computational Chemistry Lab Module: Conformational Analysis of Alkanes In this experiment, we will use CAChe software package to model the conformations of butane, 2-methylbutane, and substituted
More informationCHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued
CHE3935 Lecture 4 Quantum Mechanical Simulation Methods Continued 1 OUTLINE Review Introduction to CPMD MD and ensembles The functionals of density functional theory Return to ab initio methods Binding
More informationThis tutorial is intended to familiarize you with the Geomatica Toolbar and describe the basics of viewing data using Geomatica Focus.
PCI GEOMATICS GEOMATICA QUICKSTART 1. Introduction This tutorial is intended to familiarize you with the Geomatica Toolbar and describe the basics of viewing data using Geomatica Focus. All data used in
More informationJaguar DFT Optimizations and Transition State Searches
Jaguar DFT Optimizations and Transition State Searches Density Functional Theory (DFT) is a quantum mechanical (QM) method that gives results superior to Hartree Fock (HF) in less computational time. A
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationTechnical Note Calculations of Orbital Overlap Range Function EDR( r ; d) and Overlap Distance D(r )using Multiwfn
Technical Note Calculations of Orbital Overlap Range Function EDR( r ; d) and Overlap Distance D(r )using Multiwfn Abstract The orbital overlap range function EDR( r; d) (J. Chem. Phys. 2014, 141, 144104)
More informationRelated topics Velocity, acceleration, force, gravitational acceleration, kinetic energy, and potential energy
Related topics Velocity, acceleration, force, gravitational acceleration, kinetic energy, and potential energy Principle A mass, which is connected to a cart via a silk thread, drops to the floor. The
More informationUsing Web-Based Computations in Organic Chemistry
10/30/2017 1 Using Web-Based Computations in Organic Chemistry John Keller UAF Department of Chemistry & Biochemistry The UAF WebMO site Practical aspects of computational chemistry theory and nomenclature
More informationCS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C.
CS 273 Prof. Serafim Batzoglou Prof. Jean-Claude Latombe Spring 2006 Lecture 12 : Energy maintenance (1) Lecturer: Prof. J.C. Latombe Scribe: Neda Nategh How do you update the energy function during the
More informationWorking with Digital Elevation Models in ArcGIS 8.3
Working with Digital Elevation Models in ArcGIS 8.3 The homework that you need to turn in is found at the end of this document. This lab continues your introduction to using the Spatial Analyst Extension
More informationSimple Harmonic Motion
Introduction Simple Harmonic Motion The simple harmonic oscillator (a mass oscillating on a spring) is the most important system in physics. There are several reasons behind this remarkable claim: Any
More informationAssignment 1: Molecular Mechanics (PART 2 25 points)
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard September 2, 2015 Assignment 1: Molecular Mechanics (PART 2 25 points) In this assignment, you will perform some additional molecular mechanics calculations
More informationOECD QSAR Toolbox v.4.1. Tutorial illustrating new options for grouping with metabolism
OECD QSAR Toolbox v.4.1 Tutorial illustrating new options for grouping with metabolism Outlook Background Objectives Specific Aims The exercise Workflow 2 Background Grouping with metabolism is a procedure
More informationM E R C E R W I N WA L K T H R O U G H
H E A L T H W E A L T H C A R E E R WA L K T H R O U G H C L I E N T S O L U T I O N S T E A M T A B L E O F C O N T E N T 1. Login to the Tool 2 2. Published reports... 7 3. Select Results Criteria...
More informationExample Resistive Heating
Example Resistive Heating SOLVED WITH COMSOL MULTIPHYSICS 3.5a COPYRIGHT 2008. All right reserved. No part of this documentation may be photocopied or reproduced in any form without prior written consent
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationExperiment 14 It s Snow Big Deal
Experiment 14 It s Snow Big Deal OUTCOMES After completing this experiment, the student should be able to: use computer-based data acquisition techniques to measure temperatures. draw appropriate conclusions
More informationAssignment 1: Molecular Mechanics (PART 1 25 points)
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard August 19, 2015 Assignment 1: Molecular Mechanics (PART 1 25 points) In this assignment, you will perform some molecular mechanics calculations using the
More informationSimulation of molecular systems by molecular dynamics
Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction
More informationPreparing a PDB File
Figure 1: Schematic view of the ligand-binding domain from the vitamin D receptor (PDB file 1IE9). The crystallographic waters are shown as small spheres and the bound ligand is shown as a CPK model. HO
More informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More informationCalculating NMR Chemical Shifts for beta-ionone O
Calculating NMR Chemical Shifts for beta-ionone O Molecular orbital calculations can be used to get good estimates for chemical shifts. In this exercise we will calculate the chemical shifts for beta-ionone.
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationThe micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction
ISBN 978-1-84626-081-0 Proceedings of the 2010 International Conference on Application of Mathematics and Physics Volume 1: Advances on Space Weather, Meteorology and Applied Physics Nanjing, P. R. China,
More informationSimple Harmonic Motion
Physics Topics Simple Harmonic Motion If necessary, review the following topics and relevant textbook sections from Serway / Jewett Physics for Scientists and Engineers, 9th Ed. Hooke s Law (Serway, Sec.
More informationVirtual Cell Membrane Potential Tutorial IV
Virtual Cell Membrane Potential Tutorial IV Creating the BioModel Creating the Application!" Application I -Studying voltage changes in a compartmental model!" Application II - Studying voltage, sodium,
More informationStandards-Based Quantification in DTSA-II Part II
Standards-Based Quantification in DTSA-II Part II Nicholas W.M. Ritchie National Institute of Standards and Technology, Gaithersburg, MD 20899-8371 nicholas.ritchie@nist.gov Introduction This article is
More information41. Sim Reactions Example
HSC Chemistry 7.0 41-1(6) 41. Sim Reactions Example Figure 1: Sim Reactions Example, Run mode view after calculations. General This example contains instruction how to create a simple model. The example
More informationTutorial 8 Raster Data Analysis
Objectives Tutorial 8 Raster Data Analysis This tutorial is designed to introduce you to a basic set of raster-based analyses including: 1. Displaying Digital Elevation Model (DEM) 2. Slope calculations
More informationTutorial 23 Back Analysis of Material Properties
Tutorial 23 Back Analysis of Material Properties slope with known failure surface sensitivity analysis probabilistic analysis back analysis of material strength Introduction Model This tutorial will demonstrate
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationVCell Tutorial. Building a Rule-Based Model
VCell Tutorial Building a Rule-Based Model We will demonstrate how to create a rule-based model of EGFR receptor interaction with two adapter proteins Grb2 and Shc. A Receptor-monomer reversibly binds
More information