1
|
Rusakov YY, Rusakova IL. Getaway from the Geometry Factor Error in the Molecular Property Calculations: Efficient pecG- n ( n = 1, 2) Basis Sets for the Geometry Optimization of Molecules Containing Light p Elements. J Chem Theory Comput 2024. [PMID: 39075034 DOI: 10.1021/acs.jctc.4c00772] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
The basis of all molecular property quantum chemical calculations is the correct equilibrium geometry. In this paper, new efficient pecG-n (n = 1, 2) basis sets for the geometry optimization of molecules containing hydrogen and p elements of 2-3 periods are proposed. These basis sets were optimized via the property-energy consistent (PEC) algorithm directed to the minimization of the molecular energy gradient relative to the bond lengths. New basis sets are compact and give equilibrium geometries of very high quality, which is comparable to that provided by considerably larger energy-optimized basis sets. The equilibrium geometries obtained with the pecG-n (n = 1, 2) basis sets and the other basis sets of diverse quality were tested in the CCSD calculations of different second-order molecular properties, including NMR shielding constants, static polarizabilities, and static magnetizabilities. As a result, new basis sets have demonstrated far superior performance as compared to the other energy-optimized basis sets of the same or close sizes commonly used at the geometry optimization stage.
Collapse
Affiliation(s)
- Yuriy Yu Rusakov
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| | - Irina L Rusakova
- Siberian Branch of the Russian Academy of Sciences, A. E. Favorsky Irkutsk Institute of Chemistry, Favorsky St. 1, 664033 Irkutsk, Russian Federation
| |
Collapse
|
2
|
Tasdemir HU. Effects of intramolecular hydrogen bonding on nuclear magnetic resonance, electron paramagnetic resonance and molecular docking studies: Mexiletine molecule. J Mol Model 2024; 30:41. [PMID: 38228865 DOI: 10.1007/s00894-024-05838-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Accepted: 01/09/2024] [Indexed: 01/18/2024]
Abstract
CONTEXT In this study, the molecular structure of the mexiletine molecule was investigated. Since the Mexiletine molecule is a drug active ingredient, its molecular structure and spectroscopic properties are important. The effects of intramolecular hydrogen bonding on Nuclear Magnetic Resonance Parameters (NMR), Electron Paramagnetic Resonance (EPR) parameters and molecular docking studies were examined in the mexiletine molecule. The effects of intramolecular hydrogen bonding on EPR parameters and molecular docking studies are the most important steps for this study. METHOD Conformational space scanning required for molecular structure calculations was carried out with the Molecular Mechanic Force Field method. DFT method with 6-311 + + G(d,p) basis set level was used to obtain the most stable structure among the conformations. NMR parameters (1H and 13C chemical shift values) were also performed using the same basis set as the DFT method. The radicals created to calculate the Electron Paramagnetic Resonance parameters were modeled using the DFT/B3LYP/6-311 + + G(d,p) method basis set level. Molecular Docking studies were carried out with the Autodock vina program.
Collapse
Affiliation(s)
- Halil Ugur Tasdemir
- Ahmet Kelesoglu Education Faculty, Physics Department, Necmettin Erbakan University, Konya, Turkey.
| |
Collapse
|
3
|
Soliman UA. Computational Studies, NMR and Ir Spectral Analysis, Normal Coordinate Analysis, and Thermodynamic Properties of 2-Fluoro-4-Pyridineboronic Acid. J STRUCT CHEM+ 2020. [DOI: 10.1134/s0022476620030075] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
4
|
Mohamed TA, Soliman UA, Zoghaib WM. Raman, infrared and NMR spectra, vibrational assignments and quantum mechanical calculations of centrosymmetric 3,6-Dicholoro-1,2,4,5-tetrazine. J Mol Struct 2019. [DOI: 10.1016/j.molstruc.2018.10.047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
5
|
Mohamed TA, Soliman UA, Shaaban IA, Zoghaib WM, Wilson LD. Raman, infrared and NMR spectral analysis, normal coordinate analysis and theoretical calculations of 5-(methylthio)-1,3,4-thiadiazole-2(3H)-thione and its thiol tautomer. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2015; 150:339-349. [PMID: 26056985 DOI: 10.1016/j.saa.2015.05.039] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/21/2015] [Revised: 04/28/2015] [Accepted: 05/04/2015] [Indexed: 06/04/2023]
Abstract
Raman (3400-100 cm(-1)) and infrared (4000-200 cm(-1)) spectra of 5-(methylthio)-1,3,4-thiadiazole-2(3H)-thione (C3H4N2S3; MTT) were measured in the solid state, and the (1)H/(13)C NMR spectra were obtained in DMSO-d6. Initially, twelve structures were proposed as a result of thiol-thione tautomerism and the internal rotation about the C-S bonds. The energies and vibrational frequencies of the optimized structures were calculated using the 6-31G(d) basis set with the methods of MP2 and DFT/B3LYP with Gaussian 98 quantum calculations. Additionally, (1)H/(13)C NMR chemical shifts were predicted for the thiol (structure 5) and thione (structure 9) tautomers by means of B3LYP/6-311+G(d,p) calculations utilizing the GIAO approximation and the PCM solvation model. After complete relaxation of twelve candidate isomers, the thione tautomer (structure 9) was favored owing to its low energy and its predicted real spectral frequencies. These results agree with the recorded infrared and Raman results, in addition to the observed/calculated (1)H and (13)C NMR spectra. Aided by normal coordinate analysis and potential energy distributions (PEDs), complete vibrational assignments have been proposed for all observed fundamentals for the thione tautomer. With the aid of MP2/6-31G(d) potential surface scans, CH3, CH3S, and SH barriers to internal rotations were estimated with the optimized structural parameters from the MP2 method with the 6-31G(d) basis set. The results are discussed herein and compared with similar model compounds whenever appropriate.
Collapse
Affiliation(s)
- Tarek A Mohamed
- Department of Chemistry, Al-Azhar University (Men's Campus), Nasr City 11884, Cairo, Egypt.
| | - Usama A Soliman
- Department of Chemistry, Al-Azhar University (Men's Campus), Nasr City 11884, Cairo, Egypt
| | - Ibrahim A Shaaban
- Department of Chemistry, Al-Azhar University (Men's Campus), Nasr City 11884, Cairo, Egypt
| | - Wajdi M Zoghaib
- Department of Chemistry, Sultan Qaboos University, P.O. Box 36, Al Khod, Muscat, Oman
| | - Lee D Wilson
- Department of Chemistry, 110 Science Place (Rm. 156), University of Saskatchewan, Saskatoon, SK S7N 5C9, Canada
| |
Collapse
|
6
|
Prabavathi N, Nilufer A, Krishnakumar V. Vibrational spectroscopic (FT-IR and FT-Raman) studies, natural bond orbital analysis and molecular electrostatic potential surface of Isoxanthopterin. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2013; 114:101-113. [PMID: 23751224 DOI: 10.1016/j.saa.2013.05.013] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/05/2013] [Revised: 05/03/2013] [Accepted: 05/06/2013] [Indexed: 06/02/2023]
Abstract
The FTIR and FT-Raman spectra of Isoxanthopterin have been recorded in the region 4000-450 and 4000-100 cm(-1), respectively. The optimized geometry, frequency and intensity of the vibrational bands of Isoxanthopterin were obtained by the density functional theory (DFT) using 6-311++G(d,p) basis set. The harmonic vibrational frequencies were scaled and compared with experimental values. The observed and the calculated frequencies are found to be in good agreement. The (1)H and (13)C nuclear magnetic resonance chemical shifts of the molecule were also calculated using the gauge independent atomic orbital (GIAO) method. The UV-visible spectrum was also recorded and compared with the theoretical values. The calculated HOMO and LUMO energies show that charge transfer occurs within molecule. The first order hyperpolarizability (β0), related properties (β, α0 and Δα) and the Mulliken charges of the molecule were also computed using DFT calculations. Stability of the molecule arising from hyperconjugative interactions, charge delocalization have been analyzed using natural bond orbital (NBO) analysis. The results show that charge in electron density (ED) in the σ* and π* antibonding orbitals and second order delocalization energies (E2) confirms the occurrence of intramolecular charge transfer (ICT) within the molecule. Information about the charge density distribution of the molecule and its chemical reactivity has been obtained by mapping molecular electrostatic potential surface. In addition, the non-linear optical properties were discussed from the dipole moment values and excitation wavelength in the UV-visible region.
Collapse
Affiliation(s)
- N Prabavathi
- Department of Physics, Sri Sarada College for Women (Autonomous), Salem, India.
| | | | | |
Collapse
|
7
|
Prabavathi N, Nilufer A, Krishnakumar V. Molecular structure, vibrational, UV, NMR, hyperpolarizability, NBO and HOMO-LUMO analysis of Pteridine2,4-dione. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2012; 99:292-302. [PMID: 23085999 DOI: 10.1016/j.saa.2012.09.003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2012] [Revised: 08/22/2012] [Accepted: 09/04/2012] [Indexed: 06/01/2023]
Abstract
The FTIR and FT-Raman spectra of Pteridine2,4-dione has been recorded in the region 4000-450 and 4000-100 cm(-1), respectively. The tautomeric stability, optimized geometry, frequency and intensity of the vibrational bands of Pteridine2,4-dione were obtained by the density functional theory (DFT) using 6-311++G(d,p) basis set. The harmonic vibrational frequencies were calculated and the scaled values have been compared with experimental FTIR and FT-Raman spectra. The observed and the calculated frequencies are found to be in good agreement. The (1)H and(13)C NMR spectra chemical shifts of the molecule were also calculated using the gauge independent atomic orbital (GIAO) method. The theoretical UV-Vis spectrum of the compound using CIS method and the electronic properties, such as HOMO and LUMO energies, were performed by time-dependent DFT (TD-DFT) approach. The calculated HOMO and LUMO energies show that charge transfer occurs within molecule. The first order hyperpolarizability (β(0)) of these novel molecular system and related properties (β, α(0) and Δα) of Pteridine2,4-dione are calculated using DFT/6-311++G (d,p) method on the finite-field approach. The Mulliken charges, the values of electric dipole moment (μ) of the molecule were computed using DFT calculations. The change in electron density (ED) in the σ(*) antibonding orbitals and stabilization energies E(2) have been calculated by natural bond (NBO) analysis to give clear evidence of stabilization originating in the hyper conjugation of hydrogen-bonded interactions.
Collapse
Affiliation(s)
- N Prabavathi
- Department of Physics, Sri Sarada College for Women (Autonomous), Salem 636 016, India.
| | | | | |
Collapse
|
8
|
Mohamed TA, Hassan AM, Soliman UA, Zoghaib WM, Husband J, Abdelall MM. Infrared, Raman and NMR spectra, conformational stability, normal coordinate analysis and B3LYP calculations of 5-amino-4-cyano-3-(methylthio)-1H-pyrazole-1-carbothioamide. J Mol Struct 2011. [DOI: 10.1016/j.molstruc.2010.11.008] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
9
|
Mohamed TA, Shabaan IA, Zoghaib WM, Husband J, Farag RS, Alajhaz AENM. Tautomerism, normal coordinate analysis, vibrational assignments, calculated IR, Raman and NMR spectra of adenine. J Mol Struct 2009. [DOI: 10.1016/j.molstruc.2009.09.040] [Citation(s) in RCA: 45] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
10
|
Kupka T. H2O, H2, HF, F2 and F2O nuclear magnetic shielding constants and indirect nuclear spin-spin coupling constants (SSCCs) in the BHandH/pcJ-n and BHandH/XZP Kohn-Sham limits. MAGNETIC RESONANCE IN CHEMISTRY : MRC 2009; 47:959-970. [PMID: 19681101 DOI: 10.1002/mrc.2500] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
Good performance of segmented contracted basis sets XZP, where X = D, T, Q and 5, for obtaining H(2)O, H(2), HF, F(2) and F(2)O nuclear isotropic shielding constants in the BHandH Kohn-Sham basis set limit was shown. The results of two- and three-parameter complete basis set limit extrapolation schemes were compared with experimental results, earlier literature data and benchmark ab initio results. Similar convergence patterns of shieldings obtained from calculations using general purpose XZP basis sets and from polarization-consistent basis sets pcS-n and pcJ-n, where n = 0, 1, 2, 3 and 4, designed to accurately predict magnetic properties were observed. On the contrary, the SSCCs were more sensitive to the XZP basis set size and generally less accurate than those estimated using pcJ-n basis set family. The BHandH density functional markedly outperforms B3LYP method in predicting heavy atom shieldings and SSCCs values in the studied systems.
Collapse
Affiliation(s)
- Teobald Kupka
- Department of Chemistry, University of Opole, 48, Oleska Street, 45-052 Opole, Poland.
| |
Collapse
|
11
|
Cakmak I. GIAO calculations of chemical shifts in enantiometrically pure 1-trifluoromethyl tetrahydroisoquinoline alkaloids. ACTA ACUST UNITED AC 2005. [DOI: 10.1016/j.theochem.2004.12.004] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
12
|
Kupka T, Ruscic B, Botto RE. Hartree-Fock and density functional complete basis-set (CBS) predicted nuclear shielding anisotropy and shielding tensor components. SOLID STATE NUCLEAR MAGNETIC RESONANCE 2003; 23:145-167. [PMID: 12763561 DOI: 10.1016/s0926-2040(02)00020-6] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/24/2023]
Abstract
The nuclear shielding anisotropy and shielding tensor components calculated using the hybrid density functional B3PW91 are reported for a model set of compounds comprised of N2, NH3, CH4, C2H4, HCN and CH3CN. An estimation of density functional theory (DFT) and Hartree-Fock complete basis-set limit (CBS) parameters from a 2 (3) point exact fit vs. least-squares fit was obtained with the cc-pVxZ and aug-cc-pVxZ basis sets (x=D, T, Q, 5, 6). Both Hartree-Fock- and DFT-predicted CBS shielding anisotropies and shielding tensor components of the model molecules were in reasonable agreement with available experimental data. The utility of using a limited CBS approach for calculating accurate anisotropic shielding parameters of larger molecules as complementary methods to solid-state NMR is proposed.
Collapse
Affiliation(s)
- Teobald Kupka
- Chemistry Division, Argonne National Laboratory, 9700 S. Cass Ave., Argonne, IL 60439, USA
| | | | | |
Collapse
|
13
|
Kupka T, Wrzalik R, Pasterna G, Pasterny K. Theoretical DFT and experimental Raman and NMR studies on thiophene, 3-methylthiophene and selenophene. J Mol Struct 2002. [DOI: 10.1016/s0022-2860(02)00177-1] [Citation(s) in RCA: 58] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
14
|
Pasterny K, Wrzalik R, Kupka T, Pasterna G. Theoretical and experimental vibrational studies on liquid thiophene and its acetonitrile solution. J Mol Struct 2002. [DOI: 10.1016/s0022-2860(02)00264-8] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
15
|
Blicharska B, Kupka T. Theoretical DFT and experimental NMR studies on uracil and 5-fluorouracil. J Mol Struct 2002. [DOI: 10.1016/s0022-2860(02)00171-0] [Citation(s) in RCA: 56] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
16
|
Kupka T, Gerothanassis I, Demetropoulos I. Density functional study of a model amide. Prediction of formamide geometry, dipole moment, IR harmonic vibration ν CO and GIAO NMR shieldings. ACTA ACUST UNITED AC 2000. [DOI: 10.1016/s0166-1280(00)00441-3] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
17
|
|