1
|
Durlak P, Latajka Z. Investigations of the hydrogen bond in the crystals of tropolone and thiotropolone via car‐parrinello and path integral molecular dynamics. J Comput Chem 2018; 40:671-687. [DOI: 10.1002/jcc.25753] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2018] [Revised: 10/21/2018] [Accepted: 10/22/2018] [Indexed: 12/12/2022]
Affiliation(s)
- Piotr Durlak
- Faculty of ChemistryUniversity of Wrocław Wrocław 50‐383 Poland
| | | |
Collapse
|
2
|
Durlak P, Latajka Z. Car–Parrinello and Path Integral Molecular Dynamics Study of the Proton Transfer in the Intramolecular Hydrogen Bonds in the Ketohydrazone–Azoenol System. J Phys Chem B 2018; 122:7862-7873. [DOI: 10.1021/acs.jpcb.8b04883] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Piotr Durlak
- Faculty of Chemistry, University of Wrocław, 14 F. Joliot-Curie Street, 50-383 Wrocław, Poland
| | - Zdzisław Latajka
- Faculty of Chemistry, University of Wrocław, 14 F. Joliot-Curie Street, 50-383 Wrocław, Poland
| |
Collapse
|
3
|
Kong X, Brinkmann A, Terskikh V, Wasylishen RE, Bernard GM, Duan Z, Wu Q, Wu G. Proton Probability Distribution in the O···H···O Low-Barrier Hydrogen Bond: A Combined Solid-State NMR and Quantum Chemical Computational Study of Dibenzoylmethane and Curcumin. J Phys Chem B 2016; 120:11692-11704. [PMID: 27782387 DOI: 10.1021/acs.jpcb.6b08091] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report a combined solid-state (1H, 2H, 13C, 17O) NMR and plane-wave density functional theory (DFT) computational study of the O···H···O low-barrier hydrogen bonds (LBHBs) in two 1,3-diketone compounds: dibenzoylmethane (1) and curcumin (2). In the solid state, both 1 and 2 exist in the cis-keto-enol tautomeric form, each exhibiting an intramolecular LBHB with a short O···O distance (2.435 Å in 1 and 2.455 Å in 2). Whereas numerous experimental (structural and spectroscopic) and computational studies have been reported for the enol isomers of 1,3-diketones, a unified picture about the proton location within an LBHB is still lacking. This work reports for the first time the solid-state 17O NMR data for the O···H···O LBHBs in 1,3-diketones. The central conclusion of this work is that detailed information about the probability density distribution of the proton (nuclear zero-point motion) across an LBHB can be obtained from a combination of solid-state NMR and plane-wave DFT computations (both NMR parameter calculations and ab initio molecular dynamics simulations). We propose that the precise proton probability distribution across an LBHB should provide a common basis on which different and sometimes seemingly contradicting experimental results obtained from complementary techniques, such as X-ray diffraction, neutron diffraction, and solid-state NMR, can be reconciled.
Collapse
Affiliation(s)
- Xianqi Kong
- Department of Chemistry, Queen's University , 90 Bader Lane, Kingston, Ontario, Canada K7L 3N6
| | - Andreas Brinkmann
- Measurement Science and Standards, National Research Council Canada , 1200 Montreal Road, M-40, Ottawa, Ontario, Canada K1A 0R6
| | - Victor Terskikh
- Department of Chemistry, Queen's University , 90 Bader Lane, Kingston, Ontario, Canada K7L 3N6.,Department of Chemistry, University of Ottawa , 10 Marie Curie Private, Ottawa, Ontario, Canada K1N 6N5
| | | | - Guy M Bernard
- Department of Chemistry, University of Alberta , Edmonton, Alberta, Canada T6G 2G2
| | - Zhuang Duan
- Department of Chemistry, University of Alberta , Edmonton, Alberta, Canada T6G 2G2
| | - Qichao Wu
- Department of Chemistry, University of Alberta , Edmonton, Alberta, Canada T6G 2G2
| | - Gang Wu
- Department of Chemistry, Queen's University , 90 Bader Lane, Kingston, Ontario, Canada K7L 3N6
| |
Collapse
|
4
|
Dračínský M, Bouř P, Hodgkinson P. Temperature Dependence of NMR Parameters Calculated from Path Integral Molecular Dynamics Simulations. J Chem Theory Comput 2016; 12:968-73. [PMID: 26857802 DOI: 10.1021/acs.jctc.5b01131] [Citation(s) in RCA: 47] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The influence of temperature on NMR chemical shifts and quadrupolar couplings in model molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory (DFT) calculations of shielding and electric field gradient (EFG) tensors. An approach based on convoluting calculated shielding or EFG tensor components with probability distributions of selected bond distances and valence angles obtained from DFT-PIMD simulations at several temperatures is used to calculate the temperature effects. The probability distributions obtained from the quantum PIMD simulations, which includes nuclear quantum effects, are significantly broader and less temperature dependent than those obtained with conventional DFT molecular dynamics or with 1D scans through the potential energy surface. Predicted NMR observables for the model systems were in excellent agreement with experimental data.
Collapse
Affiliation(s)
- Martin Dračínský
- Institute of Organic Chemistry and Biochemistry , Flemingovo nám. 2, 16610 Prague, Czech Republic
| | - Petr Bouř
- Institute of Organic Chemistry and Biochemistry , Flemingovo nám. 2, 16610 Prague, Czech Republic
| | - Paul Hodgkinson
- Department of Chemistry, Durham University , South Road, DH1 3LE Durham, United Kingdom
| |
Collapse
|
5
|
Durlak P, Mierzwicki K, Latajka Z. Investigations of the Very Short Hydrogen Bond in the Crystal of Nitromalonamide via Car–Parrinello and Path Integral Molecular Dynamics. J Phys Chem B 2013; 117:5430-40. [DOI: 10.1021/jp312473b] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Piotr Durlak
- Faculty
of Chemistry, University of Wrocław,
14 F. Joliot-Curie Street, 50-383 Wrocław, Poland
| | - Krzysztof Mierzwicki
- Faculty
of Chemistry, University of Wrocław,
14 F. Joliot-Curie Street, 50-383 Wrocław, Poland
| | - Zdzisław Latajka
- Faculty
of Chemistry, University of Wrocław,
14 F. Joliot-Curie Street, 50-383 Wrocław, Poland
| |
Collapse
|
6
|
Li X, Iyengar SS. Quantum Wavepacket Ab Initio Molecular Dynamics for Extended Systems. J Phys Chem A 2011; 115:6269-84. [DOI: 10.1021/jp112389m] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaohu Li
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, IN 47405
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, IN 47405
| |
Collapse
|
7
|
Zhao X, Harbison GS. Deuterium quadrupolar tensors of L-histidine hydrochloride monohydrate-d7. J Phys Chem B 2007; 110:25059-65. [PMID: 17149930 DOI: 10.1021/jp0600965] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A deuterium NMR study at 14 T of a single crystal of L-histidine hydrochloride monohydrate has determined the deuteron quadrupole coupling constants CQ, asymmetry parameters eta, and electric field gradient orientation for the imidazolium and primary ammonium groups and for a water of crystallization. The imidazolium deuterons, which have very long relaxation times, have quite different coupling constants, reflecting different hydrogen bonding but nearly identical orientations, with the most distinct principal axis in both cases nearly parallel to the N-D vector. The -ND3 groups undergo 3-fold hops about the C-N bond axis and have typical quadrupole couplings; the D2O undergoes 2-fold hops, leading to a tensor with a large asymmetry parameter. With appropriate corrections for vibrational averaging, density functional cluster calculations give an excellent fit to the imidazolium tensor magnitudes and orientations.
Collapse
Affiliation(s)
- Xingang Zhao
- Department of Chemistry, University of Nebraska at Lincoln, Lincoln, Nebraska 68588-0304, USA
| | | |
Collapse
|