1 / 19

CEM 888: Molecular Modeling: Applications for Experimentalists

CEM 888: Molecular Modeling: Applications for Experimentalists What can theory do for the practicing chemist? What do we wish it could do? A brief progress report with examples on the performance and practical utility of theoretical tools for non-theoreticians The Evolving Roles of Theory

paul2
Download Presentation

CEM 888: Molecular Modeling: Applications for Experimentalists

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. CEM 888: Molecular Modeling: Applications for Experimentalists What can theory do for the practicing chemist? What do we wish it could do? A brief progress report with examples on the performance and practical utility of theoretical tools for non-theoreticians

  2. The Evolving Roles of Theory Pattern-matching–>Explanatory–> Exploratory–>Predictive–>Prescriptive? • Organize data around common themes • Develop physical insights into data • Interpolate and extrapolate from data to analogous systems • Predict results for proposed new expts • Identify and propose new experiments

  3. Wish List (add items as desired) Zeroth order wishes: • Theory should always get the answer right • Should be applicable to “real” molecules, materials, and situations • Should teach us something so we can better think about science (who wants to always need a supercomputer in their back pocket?)

  4. Chemistry Wishes: Structure • Internal: distances, angles, dihedrals • Rotational constants and symmetry • Conformational preferences • Solvation sphere • Lattice packing • Sensitivity to environment (solvent, pressure, fields, etc.)

  5. Chemistry Wishes: Energetics • Heats of formation • Isomer relative energies • “Effect” energies: strain, aromaticity, solvation, etc. • Heats of fusion, vaporization; heat capacities • Ionization potentials/e– affinities • UPS/XPS, E-chem, redox reagent chemistry • Bond strengths, atomization energies • MS, reaction paths

  6. Chemistry Wishes: Observables • Melting and boiling points • Density, viscosity • Refractive index, optical rotation, CD/ORD • Solubilities in various solvents • pKa (or other ionic dissociation abilities) • Dipole, polarizability • Magnetic susceptibility and e–-e– coupling

  7. Chemistry Wishes: Spectroscopies • UV-vis-IR/Raman-µwave absorption/emission • Peak positions • Transition intensities, rates (Abs, ISC, Emission, Internal conv.) • NMR • Chemical shifts,  • Spin-spin couplings, JAB • Relaxation times • NOE intensities • Conformational dynamics • EPR • Spin densities • Hyperfine coupling constants

  8. Chemistry Wishes: Reactivity and Mechanism • Activation parameters • Transition state structures and reaction paths • Conformational interconversion barriers (NMR) • Isotope effects • Solvent effects on all (structure, energetics, etc.)

  9. Chemistry Wishes: New Insights • Charge allocation to atoms--meaningful? • Why aren’t some classically valid-looking structures stable? • Are orbitals or resonance structures “real”? • How about “ring currents” • “Steric” vs. “Electronic” effects • VSEPR?

  10. Performance of the Methods: Structure • “Ordinary” compounds--simple hydrides AHn • LiH, BeH2, BH3, BH4–, CH4, NH3, NH4+, OH2, FH • NaH, MgH2, AlH3, AlH4–, SiH4, PH3, PH4+, SH2, ClH • Simple A-B bonded systems HmA-BHn • H3C-CH3, H2N-NH2, HO-OH, F-F • H3C-NH2, H3C-OH, H3C-F, etc. • Multiple bonded AB systems ABHn

  11. Beyond Minima: Reactions • Potential Energy Surfaces • Transition Structures • Reaction Paths • Transition States • Reaction Rates

  12. Generic Textbook Reaction PathUsually for unimolecular processes • How are stationary points and reaction path defined? • Are they unique and independent of coordinate system? • What is a “Reaction Coordinate” anyway, in 3n-6 dimensions?

  13. Simplest P.E. Curve: Diatomics • Diatomic dissociation is familiar • Linear structure: 3n-5 =1 mode • Reaction Coord. is uniquely defined as r(A…B)

  14. Generic Reaction PathCommon for bimolecular processes • In gas phase, two fragments always “stick” together a bit. • TS may be above or below fragment totals. • Minima may have inverse E order e.g. H3O–

  15. Rxn Paths: TS Vibrational Modes Vibrations in the TS for the (degenerate) SN2 attack of Cl– on CH3Cl. Note: all but the Reaction Coord motion are positive “ordinary” vibrations From Anwar G. Baboul and H. Bernhard Schlegel, “Improved method for calculating projected frequencies along a reaction path” J. Chem. Phys.1997, 107, 9413-9417.

  16. The Transition Structure • Stationary Point (i.e. gradients ~ 0) • One imaginary frequency (“Nimag=1”) • Locate by minimizing gradient (dE/dxi) • Structure is independent of coordinate system, nuclear masses • Is this structure the transition “state”? • How to get close enough for local optimiz’n

  17. Searching for the MEP or IRC • At the TS only, the reaction coord is well defined, as the mode with the imaginary freq. • From TS, follow steepest descent at each step. • Reaction path points not independently defined; path curves (i.e. rxn coord makeup varies) • Search scheme, step size, intermediate optimization are all important.

  18. Sample PES

  19. Anglada, J. M.; Besalú, E.; Bofill, J. M.; Crehuet, R.J. Comput. Chem.2001, 22, 387-406. Figure 9. The Müller-Brown potential surface. Dashed line, energy contours. Solid line is the reduced potential surface defined as gx0, gy=0. The black circles are the stationary points, minimum, M1 and M2, transition state, TS1. The empty circles are the starting point, P, and the turning point, TP. The black dots are the different points evaluated by the algorithm; see text for more details.

More Related