Section 2 - Input Description

Total Page:16

File Type:pdf, Size:1020Kb

Section 2 - Input Description 2-1 ( 7 Aug 02) Section 2 - Input Description This section of the manual describes the input to GAMESS. The section is written in a reference, rather than tutorial fashion. However, there are frequent reminders that more information can be found on a particular input group, or type of calculation, in the 'Further Information' section of this manual. There are also a number of examples shown in the 'Input Examples' section. The order of this section is chosen to approximate the order in which most people prepare their input ($CONTRL, $BASIS/$DATA, $GUESS, and so on). The next page contains a list of all possible input groups, in the order in which they can be found in this section. 2-2 Table of Contents $CONTRL 2-6 $MOLGRF 2-80 $SYSTEM 2-14 $STONE 2-82 $BASIS 2-16 $RAMAN 2-84 $DATA/$DATAS/$DATAL 2-19 $ALPDR 2-84 $ZMAT 2-26 $MOROKM 2-85 $LIBE 2-28 $FFCALC 2-87 $SCF 2-29 $TDHF 2-88 $SCFMI 2-33 $EFRAG 2-90 $DFT 2-34 $FRAGNAME 2-92 $MP2 2-37 $FRGRPL 2-96 $CCINP 2-39 $PRTEFP 2-97 $GUESS 2-40 $DAMP 2-98 $VEC 2-42 $DAMPGS 2-99 $MOFRZ 2-42 $PCM 2-100 $STATPT 2-43 $PCMCAV 2-104 $TRUDGE 2-46 $NEWCAV 2-105 $TRURST 2-49 $IEFPCM 2-105 $FORCE 2-49 $PCMITR 2-106 $CPHF 2-52 $DISBS 2-108 $HESS 2-52 $DISREP 2-109 $GRAD 2-53 $COSGMS 2-111 $DIPDR 2-53 $SCRF 2-111 $VIB 2-53 $ECP 2-112 $MASS 2-53 $RELWFN 2-114 $IRC 2-54 $EFIELD 2-116 $VSCF 2-57 $INTGRL 2-117 $VIBSCF 2-58 $FMM 2-118 $GRADEX 2-59 $TRANS 2-118 $DRC 2-61 $CIINP 2-120 $GLOBOP 2-64 $DET/$CIDE/$GEN/$CIGEN 2-121 $SURF 2-67 $GCILST 2-124 $LOCAL 2-69 $SODET 2-125 $TWOEI 2-72 $DRT/$CIDRT 2-126 $TRUNCN 2-73 $MCSCF 2-129 $ELMOM 2-75 $MCQDPT 2-132 $ELPOT 2-75 $CISORT 2-137 $ELDENS 2-77 $GUGEM 2-137 $ELFLDG 2-77 $GUGDIA 2-137 $POINTS 2-78 $GUGDM 2-138 $GRID 2-78 $GUGDM2 2-139 $PDC 2-78 $LAGRAN 2-140 2-3 $TRFDM2 2-140 $TRANST 2-141 2-4 name function module:routine* Molecule, basis, wavefunction specification: $CONTRL chemical control data INPUTA: START $SYSTEM computer related control data INPUTA: START $BASIS basis set INPUTB: BASISS $DATA molecule, basis set INPUTB: MOLE $ZMAT coded z-matrix ZMATRX: ZMATIN $LIBE linear bend data ZMATRX: LIBE $SCF HF-SCF wavefunction control SCFLIB: SCFIN $SCFMI SCF-MI input control data SCFMI :MIINP $DFT density functional input DFT :DFTINP $MP2 2nd order Moller-Plesset MP2: MP2INP $CCTYP coupled cluster input CCSDT: CCINP $GUESS initial orbital selection GUESS: GUESMO $VEC orbitals (formatted) GUESS: READMO $MOFRZ freezes MOs during SCF runs EFPCOV:MFRZIN Potential energy surface options: $STATPT geometry search control STATPT: SETSIG $TRUDGE nongradient optimization TRUDGE: TRUINP $TRURST restart data for TRUDGE TRUDGE: TRUDGX $FORCE hessian, normal coordinates HESS: HESSX $CPHF coupled-Hartree-Fock options CPHF: CPINP $HESS force constant matrix (formatted) HESS: FCMIN $GRAD gradient vector (formatted) HESS: EGIN $DIPDR dipole deriv. matrix (formatted) HESS: DDMIN $VIB HESSIAN restart data (formatted) HESS: HSSNUM $MASS isotope selection VIBANL:RAMS $IRC intirisic reaction path RXNCRD:IRCX $VSCF vibrational SCF and MP2 VSCF :VSCFIN $VIBSCF VSCF restart data (formatted) VSCF :VGRID $DRC dynamic reaction path DRC: DRCDRV $GLOBOP Monte Carlo global optiization GLOBOP:GLOPDR $GRADEX gradient extremal path GRADEX:GRXSET $SURF potential surface scan SURF :SRFINP Interpretation, properties: $LOCAL orbital localization control LOCAL: LMOINP $TWOEI J,K integrals (formatted) LOCCD: TWEIIN $TRUNCN localized orbital truncations EFPCOV:TRNCIN $ELMOM electrostatic moments PRPLIB: INPELM $ELPOT electrostatic potential PRPLIB: INPELP $ELDENS electron density PRPLIB: INPELD $ELFLDG electric field/gradient PRPLIB: INPELF $POINTS property calculation points PRPLIB: INPPGS $GRID property calculation mesh PRPLIB: INPPGS $PDC MEP fitting mesh PRPLIB: INPPDC $MOLGRF orbital plots PARLEY:PLTMEM $STONE distributed multipole analysis PRPPOP: STNRD $RAMAN Raman intensity RAMAN :RAMANX $ALPDR alpha polar. der. (Formatted) RAMAN :ADMIN $MOROKM Morokuma energy decomposition MOROKM:MOROIN $FFCALC finite electric field FFIELD: FFLDX $TDHF time dependent HF NLO properties TDHF: TDHFX 2-5 Solvation models: $EFRAG effective fragment potentials EFINP:EFINP $FRAGNAME specific named fragment pot. EFINP:RDSTFR $FRGRPL inter-fragment repulsion EFINP:RDDFRL $PRTEFP simplified EFP generation EFINP :PREFIN $DAMP EFP multipole screening fit CHGPEN:CGPINP $DAMPGS initial guess screening params CHGPEN:CGPINP $PCM polarizable continuum model PCM :PCMINP $PCMCAV PCM cavity generation PCM :MAKCAV $NEWCAV PCM escaped charge cavity PCM :DISREP $IEFPCM PCM integral equation form. data PCM :IEFDAT $PCMITR PCM iterative IEF input PCMIEF:ITIEFIN $DISBS PCM dispersion basis set PCMDIS:ENLBS $DISREP PCM dispersion/repulsion PCMVCH:MORETS $COSGMS conductor-like screening model COSMO :COSMIN $SCRF self consistent reaction field SCRF:ZRFINP Integral and integral modification options: $ECP effective core potentials ECPLIB:ECPPAR $RELWFN relativistic correction INPUTB:RWFINP $EFIELD external electric field PRPLIB:INPEF $INTGRL format for 2e- integrals INPUTA:START $FMM fast multipole method QMFM :QFMMIN $TRANS integral transformation TRFIN :TRANS MCSCF and CI wavefunctions, and their properties: $CIINP control over CI calculation GAMESS:WFNCI $DET determinant full CI for MCSCF ALDECI:DETINP $CIDET determinant full CI ALDECI:DETINP $GEN determinant general CI for MCSCF ALGNCI:GCIINP $CIGEN determinant general CI ALGNCI:GCIINP $GCILST general determinant list ALGNCI:GCIGEN $SODET second order determinant CI FSODCI:SOCINP $DRT distinct row table for MCSCF GUGDRT:ORDORB $CIDRT distinct row table for CI GUGDRT:ORDORB $MCSCF parameters for MCSCF MCSCF :MCSCF $MCQDPT multireference pert. Theory MCQDPT:MQREAD $CISORT integral sorting GUGSRT:GUGSRT $GUGEM Hamiltonian matrix formation GUGEM :GUGAEM $GUGDIA Hamiltonian eigenvalues/vectors GUGDGA:GUGADG $GUGDM 1e- density matrix GUGDM :GUGADM $GUGDM2 2e- density matrix GUGDM2:GUG2DM $LAGRAN CI lagrangian matrix LAGRAN:CILGRN $TRFDM2 2e- density backtransformation TRFDM2:TRF2DM $TRANST transition moments, spin-orbit TRNSTN:TRNSTX * this column is more useful to programmers than to users. 2-6 $CONTRL $CONTRL group (optional) This is a free format group specifying global switches. SCFTYP together with CITYP, MPLEVL, or CCTYP here, or DFTTYP in $DFT specifies the wavefunction. You may choose from: = RHF Restricted Hartree Fock calculation (default) = UHF Unrestricted Hartree Fock calculation = ROHF Restricted open shell Hartree-Fock. (high spin, see GVB for low spin) = GVB Generalized valence bond wavefunction or OCBSE type ROHF. (needs $SCF input) = MCSCF Multiconfigurational SCF wavefunction (this requires $DET or $DRT input) = NONE indicates a single point computation, rereading a converged SCF function. This option requires that you select CITYP=GUGA, ALDET, FSOCI, or GENCI with only RUNTYP=ENERGY or TRANSITN, and with GUESS=MOREAD. Obviously, at most one of CITYP, MPLEVL, CCTYP, or DFTTYP may be chosen in any given run. CITYP = chooses CI computation after the SCF, for any SCFTYP except UHF. = NONE skips the CI. (default) = GUGA runs the Unitary Group CI package, which requires $CIDRT input. Gradients are available only for RHF, so for other SCFTYPs, you may choose only RUNTYP=ENERGY, TRUDGE, SURFACE, FFIELD, or TRANSITN. = ALDET runs the Ames Laboratory determinant full CI package, requiring $CIDET input. Use with RUNTYP=ENERGY only. = FSOCI runs a full second order CI using determinants, with RUNTYP=ENERGY only. The input is $CIDET and $SODET. = GENCI runs a determinant CI program that permits arbitrary specification of the determinants, requiring $CIGEN input. Use with RUNTYP=ENERGY only. MPLEVL = chooses Moller-Plesset perturbation theory level, after the SCF. See $MP2 and $MCQDPT input groups. = 0 skip the MP computation (default) = 2 perform second order energy correction. MP2 (aka MBPT(2)) is only implemented for RHF, UHF, ROHF, and MCSCF wavefunctions. Gradients are available for RHF and UHF, so for the others you may pick from RUNTYP=ENERGY, TRUDGE, SURFACE, or FFIELD only. CCTYP chooses a Coupled-Cluster computation after SCF, available for SCFTYP=RHF. See also the $CCINP group. = NONE skips CC computation (default). 2-7 $CONTRL = LCCD perform a coupled-cluster calculation using the linearized coupled-cluster method with double excitations. = CCD perform a CC calculation using the coupled-cluster method with double excitations. = CCSD perform a CC calculation using the coupled-cluster method with single and double excitations. = CCSD(T) in addition to the CCSD run, the non-iterative triples corrections are calculated to give the standard CCSD[T] and CCSD(T) energies. = R-CC in addition to standard CCSD, CCSD[T], and CCSD(T) calculations, renormalized R-CCSD[T] and R-CCSD(T) calculations are performed. The cost of the renormalized calculations equals standard CCSD(T). = CR-CC in addition to CCSD, CCSD[T], CCSD(T), R-CCSD[T], and R-CCSD(T) calculations, the completely renormalized CR- CCSD[T] and CR-CCSD(T) enegies are computed. The cost of CR-CCSD[T] and CR-CCSD(T) calculations, in their noniterative triples corrections portion, is twice the standard [T] and (T) corrections. The most reasonable choices are CCSD, CCSD(T), or CR-CC. Analytic gradients are not available, so use CCTYP only for RUNTYP=ENERGY, TRUDGE, SURFACE, or maybe FFIELD. Any publication describing the results of CC calculations obtained using GAMESS should give reference to P. Piecuch, S.A. Kucharski, K. Kowalski, and M. Musial, Comput.Phys. Commun., in press, 2002 For more information about the R-CCSD(T) and CR-CCSD(T) methods, see Section 4, 'Further Information'. RUNTYP specifies the type of computation, for example at a single geometry point: = ENERGY Molecular energy. (default) = GRADIENT Molecular energy plus gradient. = HESSIAN Molecular energy plus gradient plus second derivatives, including harmonic vibrational analysis. See the $FORCE input group. multiple geometry options: = OPTIMIZE Optimize the molecular geometry using analytic energy gradients. See $STATPT. = TRUDGE Non-gradient total energy minimization.
Recommended publications
  • Quantum Mechanical Calculation of Molybdenum and Tungsten Influence on the Crm-Oxide Catalyst Acidity Oyegoke Toyese1 Fadimatu N
    Hittite Journal of Science and Engineering, 2020, 7 (4) 297–311 ISSN NUMBER: 2148–4171 DOI: 10.17350/HJSE19030000199 Quantum Mechanical Calculation of Molybdenum and Tungsten Influence on the CrM-oxide Catalyst Acidity Oyegoke Toyese1 Fadimatu N. Dabai1 Adamu Uzairur2 and Baba El-Yakubu Jibril1 1Ahmadu Bello University, Department of Chemical Engineering, Zaria, Nigeria 2Ahmadu Bello University, Department of Chemistry, Zaria, Nigeria Article History: ABSTRACT Received: 2020/06/18 Accepted: 2020/11/02 Online: 2020/12/31 emi-empirical calculations were employed to understand the effects of introducing pro- Smoters such as molybdenum (Mo) and tungsten (W) on chromium (III) oxide catalyst Correspondence to: Oyegoke Toyese, for the dehydrogenation of propane into propylene. For this purpose, we investigated CrM- Ahmadu Bello University, Chemical Engineering, Zaria, NIGERIA oxide (M = Cr, Mo, and W) catalysts. In this study, the Lewis acidity of the catalyst was ex- E-Mail: [email protected] amined using Lewis acidity parameters (Ac), including ammonia and pyridine adsorption Phone: +234 703 047 91 06 energy. The results obtained from this study of overall acidity across all sites of the catalysts studied reveal Mo-modified catalyst as the one with the least acidity while the W-modified catalyst was found to have shown the highest acidity signifies that the introduction of Mo would reduce acidity while W accelerates it. The finding, therefore, confirms tungsten (W) to be more influential and would be more promising when compared to molybdenum (Mo) due to the better avenue that is offered by W for the promotion of electron exchange and its higher acidity(s).
    [Show full text]
  • 1Jwsn As a Solution State Predictor for Tungsten-Tin Solid State Bond
    Organometallics 2005, 24, 3897-3906 3897 1 JWSn as a Solution State Predictor for Tungsten-Tin Solid State Bond Length in Sterically Crowded Tungstenocene Stannyl Complexes T. Andrew Mobley,* Rochelle Gandour, Eric P. Gillis, Kwame Nti-Addae, Rahul Palchaudhuri, Presha Rajbhandari, Neil Tomson, Angel Vargas, and Qi Zheng Noyce Science Center, Grinnell College, 1116 8th Street, Grinnell, Iowa 50112 Received December 14, 2004 A correlation between the tungsten tin solid state bond distance (rWSn) and the solution 1 state one-bond scalar coupling of tungsten to tin ( JWSn) is reported. This correlation was observed in hydrido- and chlorotungstenocene diorganostannyl chlorides (Cp2WXSnClR2,X ) H(1), X ) Cl (2)) with sterically demanding organic substituents (R ) tBu (1a, 2a), Ph (1b, 2b)). The different steric demands of the organic substituents (tBu and Ph) force the complexes 1a and 1b (or 2a and 2b) to adopt at 190 K in the solid state different, single rotameric conformations. For 1a and 2a, the tBu substituents occupy the equatorial wedge of the tungstenocene, whereas for 1b and 2b the Ph substituents are rotated to place the chlorine in the equatorial wedge. Variable-temperature multinuclear NMR and density functional theory calculations at the B3LYP/LANL2DZ level indicate the possibility of low- lying conformational isomers for complexes 1b. Introduction solid and solution state structures of stannyl complexes with both small (Me) and more sterically hindered (t- Mid-transition metal stannyl complexes represent an Bu, CH(SiMe3)2) organic substituents. In the process of 1-9 area of increasing interest due to their potential as exploring the reactive chemistry of transition metal 8,9 unique reagents and catalysts.
    [Show full text]
  • Ability of the PM3 Quantum-Mechanical Method to Model Intermolecular Hydrogen Bonding Between Neutral Molecules
    Ability of the PM3 Quantum-Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules Marcus W. Jurema and George C. Shields* Department of Chemistry, Lake Forest College, Lake Forest, Illinois 60045 Received 31 March 1992; accepted 23 July 1992 The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest.
    [Show full text]
  • And Cu(II) Complexes of Acetoacetic Acid Hydrazide
    Asian Journal of Chemistry; Vol. 25, No. 13 (2013), 7371-7376 http://dx.doi.org/10.14233/ajchem.2013.14669 Synthesis, Characterization and Quantum Chemical Studies of Some Co(II) and Cu(II) Complexes of Acetoacetic Acid Hydrazide * F.A.O. ADEKUNLE , B. SEMIRE and O.A. ODUNOLA Department of Pure and Applied Chemistry, Ladoke Akintola University of Technology, Ogbomoso, Nigeria *Corresponding author: Tel: +234 8035821847; E-mail: [email protected] (Received: 11 October 2012; Accepted: 28 June 2013) AJC-13710 New complexes of cobalt(II) and copper(II) acetoacetic acid hydrazides have been synthesized and characterized by elemental analysis, infrared and electronic reflectance spectra and room temperature magnetic susceptibility measurements. The infrared spectra of the complexes revealed the coordination to the metal ion occurs at the carbonyl oxygen and the amino nitrogen of the hydrazide moiety. This was supported by theoretical calculations. The conjoint of the electronic spectra and magnetic susceptibility measurements suggests plausible octahedral geometry for the cobalt(II) complexes while the copper(II) complexes adopt a square planar geometry. The stereochemistry of the modeled structures calculated by semi-empirical (PM3) and density functional theory (DFT) methods are also consistent with those experimentally deduced. Key Words: Synthesis, Infrared, Hydrazides, DFT, Electronic spectra. INTRODUCTION Preparation of acetoacetic acid hydrazide: Ethylaceto acetate (51 mL, 0.04 mol) was added dropwisely to hydrazine Transition metal complexes of hydrazides have been hydrate (19.4 mL 0.04 mol) in quick fit conical flask fitted intensely investigated by coordination chemists because of with a reflux condenser. The orange coloured suspension their interesting structural properties and their wide ranging obtained was refluxed for 15 min and on addition of 170 mL 1-5 applications .
    [Show full text]
  • An Hourglass Circuit Motif Transforms a Motor Program Via Subcellularly Localized Muscle Calcium Signaling and Contraction
    RESEARCH ARTICLE An hourglass circuit motif transforms a motor program via subcellularly localized muscle calcium signaling and contraction Steven R Sando1, Nikhil Bhatla1,2,3, Eugene LQ Lee1,3, H Robert Horvitz1* 1Howard Hughes Medical Institute, Department of Biology, McGovern Institute for Brain Research, Massachusetts Institute of Technology, Cambridge, United States; 2Miller Institute, Helen Wills Neuroscience Institute, Department of Molecular and Cellular Biology, University of California, Berkeley, Berkeley, United States; 3Department of Brain and Cognitive Sciences, Massachusetts Institute of Technology, Cambridge, United States Abstract Neural control of muscle function is fundamental to animal behavior. Many muscles can generate multiple distinct behaviors. Nonetheless, individual muscle cells are generally regarded as the smallest units of motor control. We report that muscle cells can alter behavior by contracting subcellularly. We previously discovered that noxious tastes reverse the net flow of particles through the C. elegans pharynx, a neuromuscular pump, resulting in spitting. We now show that spitting results from the subcellular contraction of the anterior region of the pm3 muscle cell. Subcellularly localized calcium increases accompany this contraction. Spitting is controlled by an ‘hourglass’ circuit motif: parallel neural pathways converge onto a single motor neuron that differentially controls multiple muscles and the critical subcellular muscle compartment. We conclude that subcellular muscle units enable modulatory motor control and propose that subcellular muscle contraction is a fundamental mechanism by which neurons can reshape behavior. *For correspondence: [email protected] Competing interests: The Introduction authors declare that no How animal nervous systems differentially control muscle contractions to generate the variety of flex- competing interests exist. ible, context-appropriate behaviors necessary for survival and reproduction is a fundamental prob- Funding: See page 29 lem in neuroscience.
    [Show full text]
  • Winmostar™ User Manual Release 10.7.0
    Winmostar™ User Manual Release 10.7.0 X-Ability Co., Ltd. Sep 30, 2021 Contents 1 Introduction 2 2 Installation Guide 26 3 Main Window 30 4 Basic Operation Flow 33 5 Structure Building 36 6 Main Menu and Subwindows 43 7 Remote job 183 8 Add-On 191 9 Integration with other software 200 10 Other topics 202 11 Known problems 208 12 Frequently asked questions · Troubleshooting 212 Bibliography 242 i Winmostar™ User Manual, Release 10.7.0 This manual describes the operation method of each function of Winmostar (TM). The latest version of this document is available from Official site. If you are using Winmostar (TM) for the first time, please refer to Quick Manual. If there is an uncertain point or it does not move as expected, please confirm Frequently asked questions · Troubleshooting which is updated from time to time. For specific operational procedures for each purpose, such as chemical reaction analysis and calculation of specific physical properties, see various tutorials. Contents 1 CHAPTER 1 Introduction Winmostar (TM) provides a graphical user interface that can efficiently manipulate quantum chemical calculations, first principles calculations, and molecular dynamics calculations. From the creation of the initial structure, from the calculation execution to the result analysis, you can carry out the one operation required for the simulation on Winmostar (TM). For molecular modeling function it has been confirmed to operate up to 100,000 atoms. The function of MD calculation has been confirmed in a larger system. 1.1 About quotation When announcing data created using Winmostar (TM) in academic presentations, articles, etc., please describe the Winmostar (TM) main body as follows, for example.
    [Show full text]
  • Information to Users
    INFORMATION TO USERS This manuscript has been rqjroduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter Ace, while others may be firom any type o f computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, b%inning at the upper left-hand comer and continuing from left to right in equal sections with small overlaps. Each original is also photographed in one exposure and is included in reduced form at the back of the book. Photographs included in the original manuscript have been reproduced xerographically in this copy. Ifigher quality 6” x 9” black and white photographic prints are available for any photographs or illustrations appearing in this copy for an additional charge. Contact UMI directly to order. UMI A Bell & Howell Infoimation Company 300 North Zed) Road, Ann Arbor MI 48106-1346 USA 313/761-4700 800/5214)600 i INVESTIGATIONS OF THE STRUCTURE, ENERGETICS, AND SPECTRA OF WATER CLUSTERS DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Michael David Tissandier, B.S.
    [Show full text]
  • QM/MM Study Tutorial Using Gaussview, Gaussian, and TAO Package
    QM/MM Study Tutorial using GaussView, Gaussian, and TAO package Peng Tao and H. Bernhard Schlegel Department of Chemistry, Wayne State University, 5101 Cass Ave, Detroit, Michigan 48202 [email protected] 1 Introduction This tutorial is the supporting information for the article “A Toolkit to Assist ONIOM Calculations” (TAO) by Peng Tao and H. Bernhard Schlegel. This tutorial has been developed to demonstrate the general procedure for a quantum mechanics / molecular mechanics (QM/MM) study of a biochemical system using Gaussian, GaussView and the TAO package. The example used in this tutorial is the inhibition mechanism of matrix metalloproteinase 2 (MMP2) by a selective inhibitor (4- phenoxyphenylsulfonyl)methylthiirane (SB-3CT). This QM/MM study was described in the article “Matrix Metalloproteinase 2 Inhibition: Combined Quantum Mechanics and Molecular Mechanics Studies of the Inhibition Mechanism of (4-Phenoxyphenylsulfonyl)methylthiirane and Its Oxirane Analogue.” Biochemistry 2009, 48, 9839-9847. This tutorial is designed for users who are familiar with general use of Gaussian, GaussView and Unix/Linux, and who are planning to conduct QM/MM studies of biological systems using the ONIOM method available in the Gaussian package. For more details, please refer to the user manuals and information about general usage of Gaussian, GaussView and Unix/Linux system. In this tutorial, we use the MMP2·SB-3CT complex to demonstrate ONIOM input job preparation, job monitoring, production calculations, and analysis carried out in a typical QM/MM study. In particular, the structure of the reactant complex is optimized in this tutorial. The initial structure for the complex was generated from docking and molecular dynamics (MD) studies of SB-3CT with the crystal structure of MMP2.
    [Show full text]
  • Semiempirical Methods
    Chapter 7 Semiempirical Methods 7.1 Introduction Semiempirical methods modify Hartree-Fock (HF) calculations by introducing functions with empirical parameters. These parameters are adjusted with experimen- tal conclusions to improve the quality of computation. The real cost of computation is due to the two-electron integrals in the Hamiltonian that has been simplified in this method. Semiempirical methods are based on three approximation schemes. 1. The elimination of the core electrons from the calculation. Inner electrons do not contribute towards chemical activity, which makes it pos- sible to remove the core electron functions from the Hamiltonian calculation. Normally, the entire core (the nucleus and core electrons) of atoms is replaced by a parameterized function. This has the effect of drastically reducing the com- plexity of the calculation without a major impact on the accuracy. 2. The use of the minimum number of basis sets. In this approximation, while introducing the functions of valence electrons, only the minimum required number of basis sets will be used. This technique also reduces the complexity of computation to a large extent. 3. The reduction of the number of two-electron integrals. This approximation is introduced on the basis of experimentation rather than chemical grounds. The majority of the work in ab initio calculations is in the evaluation of the two electron integrals (Coulomb and exchange). All modern semiempirical methods are based on the modified neglect of differential over- lap (MNDO) approach. In this method, parameters are assigned for different atomic types and are fitted to reproduce properties such as heats of formation, geometrical variables, dipole moments, and first ionization energies.
    [Show full text]
  • Precision Cutting Tools Powdered Metal Catalog
    MULTI FLUTE MULTI FLUTE MULTI FLUTE FINISHERS COARSE & FINE PITCH ROUGHER ROUGH / FINISH Powder Metal (PM30) EndMills 4 flute / 350 Helix / For High Temp. Alloys L2 D2 D1 L1 4 Flute Finisher EndMill - Center Cut AlTiN 4 35o Available in Exxtral Plus Coating - T4 PM30 Exxtral Plus Flutes Right Hand Spiral & Cut Coating Weldon Flats Square & Ballnose End MULTI FLUTE MULTI FINISHERS RIGHT HAND - 4 Flute D1 D2 L1 L2 AlTiN Cutting Dia Shank Dia LOC OAL Uncoated Exxtral Plus Coated Ball End Uncoated 1/2 1/2 1-1/4 3-1/4 PM3-54008 T4-PM3-54008 PM3-54108 1/2 1/2 2 4 PM3-54208 T4-PM3-54208 PM3-54308 1/2 1/2 3 5 PM3-54408 T4-PM3-54408 PM3-54508 1/2 1/2 4 6 PM3-54608 T4-PM3-54608 PM3-54708 5/8 5/8 1-5/8 3-3/4 PM3-54010 T4-PM3-54010 PM3-54110 5/8 5/8 2-1/2 4-5/8 PM3-54210 T4-PM3-54210 PM3-54310 5/8 5/8 3 5-1/8 PM3-54410 T4-PM3-54410 PM3-54510 5/8 5/8 4 6-1/8 PM3-54610 T4-PM3-54610 PM3-54710 3/4 3/4 1-5/8 3-7/8 PM3-54012 T4-PM3-54012 PM3-54112 3/4 3/4 2 4-1/4 PM3-54212 T4-PM3-54212 PM3-54312 3/4 3/4 2-1/4 4-1/2 PM3-54212S T4-PM3-54212S PM3-54312S 3/4 3/4 3 5-1/4 PM3-54412 T4-PM3-54412 PM3-54512 3/4 3/4 4 6-1/4 PM3-54612 T4-PM3-54612 PM3-54712 3/4 3/4 5 7-1/4 PM3-54612S T4-PM3-54612S PM3-54712S 1 1 2 4-1/2 PM3-54016 T4-PM3-54016 PM3-54116 1 1 3 5-1/2 PM3-54216 T4-PM3-54216 PM3-54316 1 1 4 6-1/2 PM3-54416 T4-PM3-54416 PM3-54516 1 1 5 7-1/2 PM3-54416S T4-PM3-54416S PM3-54516S 1-1/4 1-1/4 2 4-1/2 PM3-54020 T4-PM3-54020 PM3-54120 1-1/4 1-1/4 3 5-1/2 PM3-54220 T4-PM3-54220 PM3-54320 1-1/4 1-1/4 4 6-1/2 PM3-54420 T4-PM3-54420 PM3-54520 1-1/4 1-1/4 5 7-1/2 PM3-54420S T4-PM3-54420S PM3-54520S 1-1/4 1-1/4 6 8-1/2 PM3-54620 T4-PM3-54620 PM3-54720 1-1/2 1-1/4 2 4-1/2 PM3-54024 T4-PM3-54024 PM3-54124 1-1/2 1-1/4 3 5-1/2 PM3-54224 T4-PM3-54224 PM3-54324 1-1/2 1-1/4 4 6-1/2 PM3-54424 T4-PM3-54424 PM3-54524 1-1/2 1-1/4 6 8-1/2 PM3-54624 T4-PM3-54624 PM3-54724 2 PRECISION CUTTING TOOLS 1-(800) 334-9373 1-(562) 921-7898 Fax 1-(562) 926-0156 6,8 flute / 350 Helix / For High Temp.
    [Show full text]
  • Interaction of Low-Energy Ions and Atoms of Light Elements with a Fluorinated Carbon Molecular Lattice
    1508 J. Phys. Chem. A 2007, 111, 1508-1514 Interaction of Low-Energy Ions and Atoms of Light Elements with a Fluorinated Carbon Molecular Lattice Pavel V. Avramov*,†,‡ and Boris I. Yakobson§ Takasaki-branch, AdVanced Science Research Center, Japan Atomic Energy Agency, Takasaki, 370-1292, Japan, L.V. Kirensky Institute of Physics SB RAS, 660036 Krasnoyarsk, Russian Federation, and Department of Mechanical Engineering and Material Science, and Department of Chemistry, Rice UniVersity, Houston, Texas 77005 ReceiVed: September 22, 2006; In Final Form: December 20, 2006 The mechanism of interaction of low-energy atoms and ions of light elements (H, H+, He, Li, the kinetic energy of the particles 2-40 eV) with C6H6,C6F12,C60, and C60F48 molecules was studied by ab initio MD simulations and quantum-chemical calculations. It was shown that starting from 6 Å from the carbon skeleton for the “C6H6 + proton” and “C60 + proton” systems, the electronic charge transfer from the aromatic molecule to H+ occurs with a probability close to 1. The process transforms the H+ to a hydrogen atom and the neutral C6H6 and C60 molecules to cation radicals. The mechanism of interaction of low-energy protons with C6F12 and C60F48 molecules has a substantially different character and can be considered qualitatively as the interaction between a neutral molecule and a point charge. The Coulomb perturbation of the system arising from the interaction of the uncompensated proton charge with the Mulliken charges of fluorine atoms results in an inversion of the energies of the electronic states localized on the proton and on the C6F12 and C60F48 molecules and makes the electronic charge transfer energetically unfavorable.
    [Show full text]
  • CHAPTER I - Introduction * * * ************************************
    (23 Sep 99) ************************************ * * * CHAPTER I - Introduction * * * ************************************ Input Philosophy Input to PC GAMESS may be in upper or lower case. There are three types of input groups in PC GAMESS: 1. A pseudo-namelist, free format, keyword driven group. Almost all input groups fall into this first category. 2. A free format group which does not use keywords. The only examples of this category are $DATA, $ECP, $POINTS, and $STONE. 3. Formatted data. This data is never typed by the user, but rather is generated in the correct format by some earlier PC GAMESS run. All input groups begin with a $ sign in column 2, followed by a name identifying that group. The group name should be the only item appearing on the input line for any group in category 2 or 3. All input groups terminate with a $END. For any group in category 2 and 3, the $END must appear beginning in column 2, and thus is the only item on that input line. Type 1 groups may have keyword input on the same line as the group name, and the $END may appear anywhere. Because each group has a unique name, the groups may be given in any order desired. In fact, multiple occurrences of category 1 groups are permissible. * * * Most of the groups can be omitted if the program defaults are adequate. An exception is $DATA, which is always required. A typical free format $DATA group is $DATA STO-3G test case for water CNV 2 OXYGEN 8.0 STO 3 HYDROGEN 1.0 -0.758 0.0 0.545 STO 3 $END Here, position is important.
    [Show full text]