Q-Chem 4.0 User's Manual

Total Page:16

File Type:pdf, Size:1020Kb

Q-Chem 4.0 User's Manual Q-Chem User's Manual Version 4.0.1 September, 2012 Version 4.0.1 September, 2012 Q-Chem User's Guide This version was edited by: Dr. Andrew Gilbert, Dr. Emil Proynov, and Prof. Anna Krylov Version 4.0 was edited by: Dr. Emil Proynov, Dr. Jing Kong, and Prof. John Herbert with contributions from those listed in the New Features Section 1.6.2 Version 3.2 was edited by: Dr. Yihan Shao with contributions from: Dr. Nick Besley (Partial Hessian) Dr. David Casanova (SF-XCIS) Dr. Jeng-Da Chai (Variations of !B97 functional) Dr. Deborah Crittenden (Wigner intracule) Dr. Evgeny Epifanovsky (Coupled-cluster parallelization) Prof. Steve Gwaltney (Onsager) Prof. John Herbert (LRC-DFT) Prof. Cherri Hsu (Electron transfer analysis) Dr. Rustam Khaliullin (ALMO, EDA, CTA) Dr. Ester Livshits (BNL functional) Dr. Alek Marenich (SM8) Prof. Young Min Rhee (SOS-CIS(D), SOS-CIS(D0)) Prof. David Sherrill (DFT-D) Dr. Vitalii Vanovschi (Distributed multipole analysis) Prof. Troy van Voorhis (Constrained DFT, Onsager, RCA) Dr. Lee Woodcock (QM/MM hessian) 3 Version 3.1 was edited by: Dr. Andrew Gilbert with contributions from: Dr. Greg Beran (Coupled-cluster active space methods) Prof. Dan Chipman and Dr. Shawn T. Brown (SS(V)PE solvation model) Dr. Laszlo Fusti-Molnar (Fourier Transform Coulomb Method) Prof. Martin Head-Gordon (Auxiliary bases, SOS MP2, perfect and imperfect pairing) Prof. John Herbert (Ab initio dynamics, Born-Oppenheimer dynamics) Dr. Jing Kong (Fast XC calculations) Prof. Anna Krylov (EOM methods) Dr. Joerg Kussman and Prof. Dr. Christian Ochsenfeld (Linear scaling NMR and optical properties) Dr. Ching Yeh Lin (Anharmonic Corrections) Dr. Rohini Lochan (SOS and MOS-MP2) Prof. Vitaly Rassolov (Geminal Models) Dr. Ryan Steele (Dual basis methods) Dr. Yihan Shao (Integral algorithm improvements, QM/MM and improved TS finder) Versions 3.1 and 3.2 are revisions and expansions based on version 2.1, which was written by: Dr. Jeremy Dombroski Prof. Martin Head-Gordon Dr. Andrew Gilbert Published by: Customer Support: Q-Chem, Inc. Telephone: (412) 687-0695 5001 Baum Blvd Facsimile: (412) 687-0698 Suite 690 email: [email protected] Pittsburgh, PA 15213 website: http://www.q-chem.com Q-Chem is a trademark of Q-Chem, Inc. All rights reserved. The information in this document applies to version 4.0.1 of Q-Chem. This document version generated on October 1, 2012. © Copyright 2012 Q-Chem, Inc. This document is protected under the U.S. Copyright Act of 1976 and state trade secret laws. Unauthorized disclosure, reproduction, distribution, or use is prohibited and may violate federal and state laws. Contents 1 Introduction 16 1.1 About This Manual . 16 1.2 Chapter Summaries . 17 1.3 Contact Information . 17 1.3.1 Customer Support . 17 1.4 Q-Chem, Inc. 18 1.5 Company Mission . 18 1.6 Q-Chem Features . 18 1.6.1 New Features in Q-Chem 4.0.1 . 18 1.6.2 New Features in Q-Chem 4.0 . 19 1.6.3 New Features in Q-Chem 3.2 . 22 1.6.4 New Features in Q-Chem 3.1 . 23 1.6.5 New Features in Q-Chem 3.0 . 24 1.6.6 Summary of Features Prior to Q-Chem 3.0 . 26 1.7 Current Development and Future Releases . 28 1.8 Citing Q-Chem ..................................... 28 References and Further Reading . 28 2 Installation 30 2.1 Q-Chem Installation Requirements . 30 2.1.1 Execution Environment . 30 2.1.2 Hardware Platforms and Operating Systems . 30 2.1.3 Memory and Hard Disk . 31 2.2 Installing Q-Chem ................................... 32 2.3 Q-Chem Auxiliary files ($QCAUX ) ......................... 32 2.4 Q-Chem Runtime Environment Variables . 32 2.5 User Account Adjustments . 33 2.6 Further Customization . 33 2.6.1 .qchemrc and Preferences File Format . 34 2.6.2 Recommendations . 34 2.7 Running Q-Chem ................................... 35 2.7.1 Running Q-Chem in parallel . 36 2.8 IQmol Installation Requirements . 37 2.9 Testing and Exploring Q-Chem ............................ 39 3 Q-Chem Inputs 40 3.1 IQmol .......................................... 40 3.2 General Form . 40 3.3 Molecular Coordinate Input ($molecule)....................... 41 CONTENTS 6 3.3.1 Reading Molecular Coordinates From a Previous Calculation . 43 3.3.2 Reading Molecular Coordinates from Another File . 44 3.4 Cartesian Coordinates . 44 3.4.1 Examples . 44 3.5 Z -matrix Coordinates . 45 3.5.1 Dummy Atoms . 47 3.6 Job Specification: The $rem Array Concept . 47 3.7 $rem Array Format in Q-Chem Input . 48 3.8 Minimum $rem Array Requirements . 49 3.9 User-Defined Basis Sets ($basis)............................ 49 3.10 Comments ($comment)................................. 49 3.11 User-Defined Pseudopotentials ($ecp)......................... 49 3.12 User-defined Parameters for DFT Dispersion Correction ($empirical dispersion) . 50 3.13 Addition of External Charges ($external charges).................. 50 3.14 Intracules ($intracule) ................................. 50 3.15 Isotopic Substitutions ($isotopes)........................... 50 3.16 Applying a Multipole Field ($multipole field) .................... 50 3.17 Natural Bond Orbital Package ($nbo) ........................ 51 3.18 User-Defined Occupied Guess Orbitals ($occupied and $swap occupied virtual) . 51 3.19 Geometry Optimization with General Constraints ($opt) . 51 3.20 Polarizable Continuum Solvation Models ($pcm) .................. 51 3.21 Effective Fragment Potential calculations ($efp fragmentsand $efp params) . 52 3.22 SS(V)PE Solvation Modeling ($svp and $svpirf ) .................. 52 3.23 Orbitals, Densities and ESPs on a Mesh ($plots) .................. 52 3.24 User-Defined van der Waals Radii ($van der waals) . 52 3.25 User-Defined Exchange-Correlation Density Functionals ($xc functional) . 52 3.26 Multiple Jobs in a Single File: Q-Chem Batch Job Files . 53 3.27 Q-Chem Output File . 55 3.28 Q-Chem Scratch Files . 55 4 Self-Consistent Field Ground State Methods 56 4.1 Introduction . 56 4.1.1 Overview of Chapter . 56 4.1.2 Theoretical Background . 57 4.2 Hartree{Fock Calculations . 60 4.2.1 The Hartree-Fock Equations . 60 4.2.2 Wavefunction Stability Analysis . 61 4.2.3 Basic Hartree-Fock Job Control . 62 4.2.4 Additional Hartree-Fock Job Control Options . 65 4.2.5 Examples . 68 4.2.6 Symmetry . 69 4.3 Density Functional Theory . 70 4.3.1 Introduction . 70 4.3.2 Kohn-Sham Density Functional Theory . 71 4.3.3 Exchange-Correlation Functionals . 72 4.3.4 Long-Range-Corrected DFT . 77 4.3.5 Nonlocal Correlation Functionals . 87 4.3.6 DFT-D Methods . 89 4.3.7 XDM DFT Model of Dispersion . 91 CONTENTS 7 4.3.8 DFT-D3 Methods . 95 4.3.9 Double-Hybrid Density Functional Theory . 98 4.3.10 Asymptotically Corrected Exchange-Correlation Potentials . 103 4.3.11 DFT Numerical Quadrature . 104 4.3.12 Angular Grids . 105 4.3.13 Standard Quadrature Grids . 105 4.3.14 Consistency Check and Cutoffs for Numerical Integration . 106 4.3.15 Basic DFT Job Control . 107 4.3.16 Example . 112 4.3.17 User-Defined Density Functionals . 112 4.4 Large Molecules and Linear Scaling Methods . 115 4.4.1 Introduction . 115 4.4.2 Continuous Fast Multipole Method (CFMM) . 116 4.4.3 Linear Scaling Exchange (LinK) Matrix Evaluation . 118 4.4.4 Incremental and Variable Thresh Fock Matrix Building . 119 4.4.5 Incremental DFT . 120 4.4.6 Fourier Transform Coulomb Method . 122 4.4.7 Multiresolution Exchange-Correlation (mrXC) Method . 124 4.4.8 Examples . 126 4.5 SCF Initial Guess . 127 4.5.1 Introduction . 127 4.5.2 Simple Initial Guesses . 127 4.5.3 Reading MOs from Disk . 129 4.5.4 Modifying the Occupied Molecular Orbitals . 129 4.5.5 Basis Set Projection . 131 4.5.6 Examples . 133 4.6 Converging SCF Calculations . 135 4.6.1 Introduction . 135 4.6.2 Basic Convergence Control Options . 135 4.6.3 Direct Inversion in the Iterative Subspace (DIIS) . 137 4.6.4 Geometric Direct Minimization (GDM) . 139 4.6.5 Direct Minimization (DM) . 140 4.6.6 Maximum Overlap Method (MOM) . 140 4.6.7 Relaxed Constraint Algorithm (RCA) . 141 4.6.8 Examples . 143 4.7 Dual-Basis Self-Consistent Field Calculations . ..
Recommended publications
  • MODELLER 10.1 Manual
    MODELLER A Program for Protein Structure Modeling Release 10.1, r12156 Andrej Saliˇ with help from Ben Webb, M.S. Madhusudhan, Min-Yi Shen, Guangqiang Dong, Marc A. Martı-Renom, Narayanan Eswar, Frank Alber, Maya Topf, Baldomero Oliva, Andr´as Fiser, Roberto S´anchez, Bozidar Yerkovich, Azat Badretdinov, Francisco Melo, John P. Overington, and Eric Feyfant email: modeller-care AT salilab.org URL https://salilab.org/modeller/ 2021/03/12 ii Contents Copyright notice xxi Acknowledgments xxv 1 Introduction 1 1.1 What is Modeller?............................................. 1 1.2 Modeller bibliography....................................... .... 2 1.3 Obtainingandinstallingtheprogram. .................... 3 1.4 Bugreports...................................... ............ 4 1.5 Method for comparative protein structure modeling by Modeller ................... 5 1.6 Using Modeller forcomparativemodeling. ... 8 1.6.1 Preparinginputfiles . ............. 8 1.6.2 Running Modeller ......................................... 9 2 Automated comparative modeling with AutoModel 11 2.1 Simpleusage ..................................... ............ 11 2.2 Moreadvancedusage............................... .............. 12 2.2.1 Including water molecules, HETATM residues, and hydrogenatoms .............. 12 2.2.2 Changing the default optimization and refinement protocol ................... 14 2.2.3 Getting a very fast and approximate model . ................. 14 2.2.4 Building a model from multiple templates . .................. 15 2.2.5 Buildinganallhydrogenmodel
    [Show full text]
  • Correlation.Pdf
    INTRODUCTION TO THE ELECTRONIC CORRELATION PROBLEM PAUL E.S. WORMER Institute of Theoretical Chemistry, University of Nijmegen, Toernooiveld, 6525 ED Nijmegen, The Netherlands Contents 1 Introduction 2 2 Rayleigh-SchrÄodinger perturbation theory 4 3 M¿ller-Plesset perturbation theory 7 4 Diagrammatic perturbation theory 9 5 Unlinked clusters 14 6 Convergence of MP perturbation theory 17 7 Second quantization 18 8 Coupled cluster Ansatz 21 9 Coupled cluster equations 26 9.1 Exact CC equations . 26 9.2 The CCD equations . 30 9.3 CC theory versus MP theory . 32 9.4 CCSD(T) . 33 A Hartree-Fock, Slater-Condon, Brillouin 34 B Exponential structure of the wavefunction 37 C Bibliography 40 1 1 Introduction These are notes for a six hour lecture series on the electronic correlation problem, given by the author at a Dutch national winterschool in 1999. The main purpose of this course was to give some theoretical background on the M¿ller-Plesset and coupled cluster methods. Both computational methods are available in many quantum chemical \black box" programs. The audi- ence consisted of graduate students, mostly with an undergraduate chemistry education and doing research in theoretical chemistry. A basic knowledge of quantum mechanics and quantum chemistry is pre- supposed. In particular a knowledge of Slater determinants, Slater-Condon rules and Hartree-Fock theory is a prerequisite of understanding the following notes. In Appendix A this theory is reviewed briefly. Because of time limitations hardly any proofs are given, the theory is sketchily outlined. No attempt is made to integrate out the spin, the theory is formulated in terms of spin-orbitals only.
    [Show full text]
  • Supporting Information
    Electronic Supplementary Material (ESI) for RSC Advances. This journal is © The Royal Society of Chemistry 2020 Supporting Information How to Select Ionic Liquids as Extracting Agent Systematically? Special Case Study for Extractive Denitrification Process Shurong Gaoa,b,c,*, Jiaxin Jina,b, Masroor Abroc, Ruozhen Songc, Miao Hed, Xiaochun Chenc,* a State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Beijing, 102206, China b Research Center of Engineering Thermophysics, North China Electric Power University, Beijing, 102206, China c Beijing Key Laboratory of Membrane Science and Technology & College of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, PR China d Office of Laboratory Safety Administration, Beijing University of Technology, Beijing 100124, China * Corresponding author, Tel./Fax: +86-10-6443-3570, E-mail: [email protected], [email protected] 1 COSMO-RS Computation COSMOtherm allows for simple and efficient processing of large numbers of compounds, i.e., a database of molecular COSMO files; e.g. the COSMObase database. COSMObase is a database of molecular COSMO files available from COSMOlogic GmbH & Co KG. Currently COSMObase consists of over 2000 compounds including a large number of industrial solvents plus a wide variety of common organic compounds. All compounds in COSMObase are indexed by their Chemical Abstracts / Registry Number (CAS/RN), by a trivial name and additionally by their sum formula and molecular weight, allowing a simple identification of the compounds. We obtained the anions and cations of different ILs and the molecular structure of typical N-compounds directly from the COSMObase database in this manuscript.
    [Show full text]
  • GROMACS: Fast, Flexible, and Free
    GROMACS: Fast, Flexible, and Free DAVID VAN DER SPOEL,1 ERIK LINDAHL,2 BERK HESS,3 GERRIT GROENHOF,4 ALAN E. MARK,4 HERMAN J. C. BERENDSEN4 1Department of Cell and Molecular Biology, Uppsala University, Husargatan 3, Box 596, S-75124 Uppsala, Sweden 2Stockholm Bioinformatics Center, SCFAB, Stockholm University, SE-10691 Stockholm, Sweden 3Max-Planck Institut fu¨r Polymerforschung, Ackermannweg 10, D-55128 Mainz, Germany 4Groningen Biomolecular Sciences and Biotechnology Institute, University of Groningen, Nijenborgh 4, NL-9747 AG Groningen, The Netherlands Received 12 February 2005; Accepted 18 March 2005 DOI 10.1002/jcc.20291 Published online in Wiley InterScience (www.interscience.wiley.com). Abstract: This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simula- tions. The package includes about 100 utility and analysis programs.
    [Show full text]
  • Effective Hamiltonians Derived from Equation-Of-Motion Coupled-Cluster
    Effective Hamiltonians derived from equation-of-motion coupled-cluster wave-functions: Theory and application to the Hubbard and Heisenberg Hamiltonians Pavel Pokhilkoa and Anna I. Krylova a Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482 Effective Hamiltonians, which are commonly used for fitting experimental observ- ables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not triv- ial. In this contribution, we apply Bloch's formalism to equation-of-motion coupled- cluster (EOM-CC) wave functions to rigorously derive effective Hamiltonians in the Bloch's and des Cloizeaux's forms. We report the key equations and illustrate the theory by examples of systems with electronic states of covalent and ionic characters. We show that the Hubbard and Heisenberg Hamiltonians are extracted directly from the so-obtained effective Hamiltonians. By making quantitative connections between many-body states and simple models, the approach also facilitates the analysis of the correlated wave functions. Artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed. I. INTRODUCTION Coarse graining is commonly used in computational chemistry and physics. It is ex- ploited in a number of classic models serving as a foundation of modern solid-state physics: tight binding[1, 2], Drude{Sommerfeld's model[3{5], Hubbard's [6] and Heisenberg's[7{9] Hamiltonians. These models explain macroscopic properties of materials through effective interactions whose strengths are treated as model parameters. The values of these pa- rameters are determined either from more sophisticated theoretical models or by fitting to experimental observables.
    [Show full text]
  • Chem Compute Quickstart
    Chem Compute Quickstart Chem Compute is maintained by Mark Perri at Sonoma State University and hosted on Jetstream at Indiana University. The Chem Compute URL is https://chemcompute.org/. We will use Chem Compute as the frontend for running electronic structure calculations with The General Atomic and Molecular Electronic Structure System, GAMESS (http://www.msg.ameslab.gov/gamess/). Chem Compute also provides access to other computational chemistry resources including PSI4 and the molecular dynamics packages TINKER and NAMD, though we will not be using those resource at this time. Follow this link, https://chemcompute.org/gamess/submit, to directly access the Chem Compute GAMESS guided submission interface. If you are a returning Chem Computer user, please log in now. If your University is part of the InCommon Federation you can log in without registering by clicking Login then "Log In with Google or your University" – select your University from the dropdown list. Otherwise if this is your first time working with Chem Compute, please register as a Chem Compute user by clicking the “Register” link in the top-right corner of the page. This gives you access to all of the computational resources available at ChemCompute.org and will allow you to maintain copies of your calculations in your user “Dashboard” that you can refer to later. Registering also helps track usage and obtain the resources needed to continue providing its service. When logged in with the GAMESS-Submit tabs selected, an instruction section appears on the left side of the page with instructions for several different kinds of calculations.
    [Show full text]
  • Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development
    processes Review Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development Outi M. H. Salo-Ahen 1,2,* , Ida Alanko 1,2, Rajendra Bhadane 1,2 , Alexandre M. J. J. Bonvin 3,* , Rodrigo Vargas Honorato 3, Shakhawath Hossain 4 , André H. Juffer 5 , Aleksei Kabedev 4, Maija Lahtela-Kakkonen 6, Anders Støttrup Larsen 7, Eveline Lescrinier 8 , Parthiban Marimuthu 1,2 , Muhammad Usman Mirza 8 , Ghulam Mustafa 9, Ariane Nunes-Alves 10,11,* , Tatu Pantsar 6,12, Atefeh Saadabadi 1,2 , Kalaimathy Singaravelu 13 and Michiel Vanmeert 8 1 Pharmaceutical Sciences Laboratory (Pharmacy), Åbo Akademi University, Tykistökatu 6 A, Biocity, FI-20520 Turku, Finland; ida.alanko@abo.fi (I.A.); rajendra.bhadane@abo.fi (R.B.); parthiban.marimuthu@abo.fi (P.M.); atefeh.saadabadi@abo.fi (A.S.) 2 Structural Bioinformatics Laboratory (Biochemistry), Åbo Akademi University, Tykistökatu 6 A, Biocity, FI-20520 Turku, Finland 3 Faculty of Science-Chemistry, Bijvoet Center for Biomolecular Research, Utrecht University, 3584 CH Utrecht, The Netherlands; [email protected] 4 Swedish Drug Delivery Forum (SDDF), Department of Pharmacy, Uppsala Biomedical Center, Uppsala University, 751 23 Uppsala, Sweden; [email protected] (S.H.); [email protected] (A.K.) 5 Biocenter Oulu & Faculty of Biochemistry and Molecular Medicine, University of Oulu, Aapistie 7 A, FI-90014 Oulu, Finland; andre.juffer@oulu.fi 6 School of Pharmacy, University of Eastern Finland, FI-70210 Kuopio, Finland; maija.lahtela-kakkonen@uef.fi (M.L.-K.); tatu.pantsar@uef.fi
    [Show full text]
  • Bornologically Isomorphic Representations of Tensor Distributions
    Bornologically isomorphic representations of distributions on manifolds E. Nigsch Thursday 15th November, 2018 Abstract Distributional tensor fields can be regarded as multilinear mappings with distributional values or as (classical) tensor fields with distribu- tional coefficients. We show that the corresponding isomorphisms hold also in the bornological setting. 1 Introduction ′ ′ ′r s ′ Let D (M) := Γc(M, Vol(M)) and Ds (M) := Γc(M, Tr(M) ⊗ Vol(M)) be the strong duals of the space of compactly supported sections of the volume s bundle Vol(M) and of its tensor product with the tensor bundle Tr(M) over a manifold; these are the spaces of scalar and tensor distributions on M as defined in [?, ?]. A property of the space of tensor distributions which is fundamental in distributional geometry is given by the C∞(M)-module isomorphisms ′r ∼ s ′ ∼ r ′ Ds (M) = LC∞(M)(Tr (M), D (M)) = Ts (M) ⊗C∞(M) D (M) (1) (cf. [?, Theorem 3.1.12 and Corollary 3.1.15]) where C∞(M) is the space of smooth functions on M. In[?] a space of Colombeau-type nonlinear generalized tensor fields was constructed. This involved handling smooth functions (in the sense of convenient calculus as developed in [?]) in par- arXiv:1105.1642v1 [math.FA] 9 May 2011 ∞ r ′ ticular on the C (M)-module tensor products Ts (M) ⊗C∞(M) D (M) and Γ(E) ⊗C∞(M) Γ(F ), where Γ(E) denotes the space of smooth sections of a vector bundle E over M. In[?], however, only minor attention was paid to questions of topology on these tensor products.
    [Show full text]
  • Solution: the First Element of the Dual Basis Is the Linear Function Α
    Homework assignment 7 pp. 105 3 Exercise 2. Let B = f®1; ®2; ®3g be the basis for C defined by ®1 = (1; 0; ¡1) ®2 = (1; 1; 1) ®3 = (2; 2; 0): Find the dual basis of B. Solution: ¤ The first element of the dual basis is the linear function ®1 such ¤ ¤ ¤ that ®1(®1) = 1; ®1(®2) = 0 and ®1(®3) = 0. To describe such a function more explicitly we need to find its values on the standard basis vectors e1, e2 and e3. To do this express e1; e2; e3 through ®1; ®2; ®3 (refer to the solution of Exercise 1 pp. 54-55 from Homework 6). For each i = 1; 2; 3 you will find the numbers ai; bi; ci such that e1 = ai®1 + bi®2 + ci®3 (i.e. the coordinates of ei relative to the ¤ ¤ ¤ basis ®1; ®2; ®3). Then by linearity of ®1 we get that ®1(ei) = ai. Then ®2(ei) = bi, ¤ and ®3(ei) = ci. This is the answer. It can also be reformulated as follows. If P is the transition matrix from the standard basis e1; e2; e3 to ®1; ®2; ®3, i.e. ¡1 t (®1; ®2; ®3) = (e1; e2; e3)P , then (P ) is the transition matrix from the dual basis ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¡1 t e1; e2; e3 to the dual basis ®1; ®2; a3, i.e. (®1; ®2; a3) = (e1; e2; e3)(P ) . Note that this problem is basically the change of coordinates problem: e.g. ¤ 3 the value of ®1 on the vector v 2 C is the first coordinate of v relative to the basis ®1; ®2; ®3.
    [Show full text]
  • Vectors and Dual Vectors
    Vectors By now you have a pretty good experience with \vectors". Usually, a vector is defined as a quantity that has a direction and a magnitude, such as a position vector, velocity vector, acceleration vector, etc. However, the notion of a vector has a considerably wider realm of applicability than these examples might suggest. The set of all real numbers forms a vector space, as does the set of all complex numbers. The set of functions on a set (e.g., functions of one variable, f(x)) form a vector space. Solutions of linear homogeneous equations form a vector space. We begin by giving the abstract rules for forming a space of vectors, also known as a vector space. A vector space V is a set equipped with an operation of \addition" and an additive identity. The elements of the set are called vectors, which we shall denote as ~u, ~v, ~w, etc. For now, you can think of them as position vectors in order to keep yourself sane. Addition, is an operation in which two vectors, say ~u and ~v, can be combined to make another vector, say, ~w. We denote this operation by the symbol \+": ~u + ~v = ~w: (1) Do not be fooled by this simple notation. The \addition" of vectors may be quite a different operation than ordinary arithmetic addition. For example, if we view position vectors in the x-y plane as \arrows" drawn from the origin, the addition of vectors is defined by the parallelogram rule. Clearly this rule is quite different than ordinary \addition".
    [Show full text]
  • Parameterizing a Novel Residue
    University of Illinois at Urbana-Champaign Luthey-Schulten Group, Department of Chemistry Theoretical and Computational Biophysics Group Computational Biophysics Workshop Parameterizing a Novel Residue Rommie Amaro Brijeet Dhaliwal Zaida Luthey-Schulten Current Editors: Christopher Mayne Po-Chao Wen February 2012 CONTENTS 2 Contents 1 Biological Background and Chemical Mechanism 4 2 HisH System Setup 7 3 Testing out your new residue 9 4 The CHARMM Force Field 12 5 Developing Topology and Parameter Files 13 5.1 An Introduction to a CHARMM Topology File . 13 5.2 An Introduction to a CHARMM Parameter File . 16 5.3 Assigning Initial Values for Unknown Parameters . 18 5.4 A Closer Look at Dihedral Parameters . 18 6 Parameter generation using SPARTAN (Optional) 20 7 Minimization with new parameters 32 CONTENTS 3 Introduction Molecular dynamics (MD) simulations are a powerful scientific tool used to study a wide variety of systems in atomic detail. From a standard protein simulation, to the use of steered molecular dynamics (SMD), to modelling DNA-protein interactions, there are many useful applications. With the advent of massively parallel simulation programs such as NAMD2, the limits of computational anal- ysis are being pushed even further. Inevitably there comes a time in any molecular modelling scientist’s career when the need to simulate an entirely new molecule or ligand arises. The tech- nique of determining new force field parameters to describe these novel system components therefore becomes an invaluable skill. Determining the correct sys- tem parameters to use in conjunction with the chosen force field is only one important aspect of the process.
    [Show full text]
  • Dmol Guide to Select a Dmol3 Task 1
    DMOL3 GUIDE MATERIALS STUDIO 8.0 Copyright Notice ©2014 Dassault Systèmes. All rights reserved. 3DEXPERIENCE, the Compass icon and the 3DS logo, CATIA, SOLIDWORKS, ENOVIA, DELMIA, SIMULIA, GEOVIA, EXALEAD, 3D VIA, BIOVIA and NETVIBES are commercial trademarks or registered trademarks of Dassault Systèmes or its subsidiaries in the U.S. and/or other countries. All other trademarks are owned by their respective owners. Use of any Dassault Systèmes or its subsidiaries trademarks is subject to their express written approval. Acknowledgments and References To print photographs or files of computational results (figures and/or data) obtained using BIOVIA software, acknowledge the source in an appropriate format. For example: "Computational results obtained using software programs from Dassault Systèmes Biovia Corp.. The ab initio calculations were performed with the DMol3 program, and graphical displays generated with Materials Studio." BIOVIA may grant permission to republish or reprint its copyrighted materials. Requests should be submitted to BIOVIA Support, either through electronic mail to [email protected], or in writing to: BIOVIA Support 5005 Wateridge Vista Drive, San Diego, CA 92121 USA Contents DMol3 1 Setting up a molecular dynamics calculation20 Introduction 1 Choosing an ensemble 21 Further Information 1 Defining the time step 21 Tasks in DMol3 2 Defining the thermostat control 21 Energy 3 Constraints during dynamics 21 Setting up the calculation 3 Setting up a transition state calculation 22 Dynamics 4 Which method to use?
    [Show full text]