To reveal the biochemical mechanism of pyruvic acid under the bioelectric currents, the effect of the external electric field on the intermolecular double-proton cooperativity migration was investigated for the dimer, accompanied by a comparison of the conformation transformation. The results indicate that the influences of the external electric fields on the transition state of the proton migration, i.e., the changes in the molecular structures, atoms in molecules, surface electrostatic potentials, and potential barriers, are closely related to the directions of the fields. The spatial order parameters of the conformation transformation of pyruvic acid can be quickly converged through the umbrella sampling and parameter averaging, thereby obtaining the potential energy surface of the potential of mean force. Under the strong external electric field, the barrier heights of the double-proton cooperativity migration are close to the free energy change of the conformation transformation, indicating the competition between the proton cooperativity transfer and conformation transformation. This study is of importance for revealing furthermore the essence of the competition between the proton migration and conformation transformation of the biomolecules stimulated by the bioelectric currents.
Methods
The effects of external electric fields on the double-proton cooperativity migration reaction pathway, potential barrier, and rate constant of the pyruvic acid dimer were studied using the M06-2X/6-311++G**, M06-2X/aug-cc-pVTZ, and CCSD(T)/6-311++G(2d,p) methods by using Gaussian 09 packages. A chaotic finite temperature string model with enhanced sampling for molecular conformation transformation was constructed, and the free energy change of conformation transformation under external electric field was studied by the finite temperature string (FTS) method with the CHARMM22 force field from NAMD software packages and PLUMED program.
{"title":"Kinetics of hydrogen migration and string of conformation transformation for pyruvic acid dimer under external electric field","authors":"Rui-zhi Feng, Xiao-jiang Li, Ai-chuan Liu, Jian Kuo, Su-rong Yan, Wei Zhang","doi":"10.1007/s00894-025-06560-z","DOIUrl":"10.1007/s00894-025-06560-z","url":null,"abstract":"<div><h3>Context</h3><p>To reveal the biochemical mechanism of pyruvic acid under the bioelectric currents, the effect of the external electric field on the intermolecular double-proton cooperativity migration was investigated for the dimer, accompanied by a comparison of the conformation transformation. The results indicate that the influences of the external electric fields on the transition state of the proton migration, i.e., the changes in the molecular structures, atoms in molecules, surface electrostatic potentials, and potential barriers, are closely related to the directions of the fields. The spatial order parameters of the conformation transformation of pyruvic acid can be quickly converged through the umbrella sampling and parameter averaging, thereby obtaining the potential energy surface of the potential of mean force. Under the strong external electric field, the barrier heights of the double-proton cooperativity migration are close to the free energy change of the conformation transformation, indicating the competition between the proton cooperativity transfer and conformation transformation. This study is of importance for revealing furthermore the essence of the competition between the proton migration and conformation transformation of the biomolecules stimulated by the bioelectric currents.</p><h3>Methods</h3><p>The effects of external electric fields on the double-proton cooperativity migration reaction pathway, potential barrier, and rate constant of the pyruvic acid dimer were studied using the M06-2X/6-311++G**, M06-2X/aug-cc-pVTZ, and CCSD(T)/6-311++G(2d,p) methods by using Gaussian 09 packages. A chaotic finite temperature string model with enhanced sampling for molecular conformation transformation was constructed, and the free energy change of conformation transformation under external electric field was studied by the finite temperature string (FTS) method with the CHARMM22 force field from NAMD software packages and PLUMED program.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 12","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-11-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145449732","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-29DOI: 10.1007/s00894-025-06548-9
Lisha Lin, Jing Li
Context
The advancement of novel sensing materials significantly enhances the precision of early disease diagnosis and treatment. This study focuses on the fabrication of nano-sensing materials, specifically Janus arsenene/PtSeS (stacking-I and stacking-II), and examines their performance in strain sensing. The calculation results indicate that both stacking configurations exhibit excellent thermodynamic stability, with formation energies of − 3.27 eV and − 3.35 eV, respectively. Electronic structure analysis shows that both are indirect bandgap semiconductors, with bandgap values of 1.229 eV and 0.845 eV, respectively, and charge transfer occurs at the interface (net transfer amounts of 0.26 |e| and 0.28 |e|, respectively). Under applied tensile or compressive strain, ranging from 0 to 4%, the bandgaps of both heterojunctions exhibit a gradual decrease or increase. In the visible light spectrum, the absorption coefficient of the heterojunctions approaches approximately 1.6 × 105 cm−1. These research outcomes are expected to facilitate the integration of these materials into medical diagnostic devices, including endoscopes.
Methods
All density functional theory calculations were conducted utilizing the CASTEP software package. The electron–electron exchange interaction was addressed through the Perdew-Burke-Ernzerhof functional, which operates within the framework of the generalized gradient approximation. To mitigate the known issue of band gap underestimation associated with the PBE functional, the Heyd-Scuseria-Ernzerhof hybrid functional was employed to enhance the electronic structure and optical absorption properties of the arsenene/PtSeS heterojunction.
{"title":"Strain regulation of the optoelectronic characteristics of arsenene/PtSeS heterostructure: promoting smart diagnosis","authors":"Lisha Lin, Jing Li","doi":"10.1007/s00894-025-06548-9","DOIUrl":"10.1007/s00894-025-06548-9","url":null,"abstract":"<div><h3>Context</h3><p>The advancement of novel sensing materials significantly enhances the precision of early disease diagnosis and treatment. This study focuses on the fabrication of nano-sensing materials, specifically Janus arsenene/PtSeS (stacking-I and stacking-II), and examines their performance in strain sensing. The calculation results indicate that both stacking configurations exhibit excellent thermodynamic stability, with formation energies of − 3.27 eV and − 3.35 eV, respectively. Electronic structure analysis shows that both are indirect bandgap semiconductors, with bandgap values of 1.229 eV and 0.845 eV, respectively, and charge transfer occurs at the interface (net transfer amounts of 0.26 |e| and 0.28 |e|, respectively). Under applied tensile or compressive strain, ranging from 0 to 4%, the bandgaps of both heterojunctions exhibit a gradual decrease or increase. In the visible light spectrum, the absorption coefficient of the heterojunctions approaches approximately 1.6 × 10<sup>5</sup> cm<sup>−1</sup>. These research outcomes are expected to facilitate the integration of these materials into medical diagnostic devices, including endoscopes.</p><h3>Methods</h3><p>All density functional theory calculations were conducted utilizing the CASTEP software package. The electron–electron exchange interaction was addressed through the Perdew-Burke-Ernzerhof functional, which operates within the framework of the generalized gradient approximation. To mitigate the known issue of band gap underestimation associated with the PBE functional, the Heyd-Scuseria-Ernzerhof hybrid functional was employed to enhance the electronic structure and optical absorption properties of the arsenene/PtSeS heterojunction.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145399445","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-28DOI: 10.1007/s00894-025-06519-0
Cleuton de Souza Silva
Context
The indium compounds were investigated both with and without spin–orbit correction, which highlighted the importance of this correction for indium compounds. The methods used in this work showed good results for enthalpy of formation when compared with other experimental results and studies in the literature with CCSD(T)/CBS. However, for bond dissociation, the results of the studied methods were highly divergent and not considered very effective when compared with experimental data from the literature. Notably, the G4-T-V2, G4-T-V3, and G4-T-V4 methods demonstrated good results for enthalpy of formation. These methods are particularly noteworthy because they utilise three energy coefficient adjustments and do not rely on higher-level correction.
Methods
We used the G4MP2-XK method along with 30 other composite methods developed by Martin and Semidalas.
{"title":"Using spin–orbit DFT module (SODFT) in indium compounds: enthalpies of formation and bond dissociation","authors":"Cleuton de Souza Silva","doi":"10.1007/s00894-025-06519-0","DOIUrl":"10.1007/s00894-025-06519-0","url":null,"abstract":"<div><h3>Context</h3><p>The indium compounds were investigated both with and without spin–orbit correction, which highlighted the importance of this correction for indium compounds. The methods used in this work showed good results for enthalpy of formation when compared with other experimental results and studies in the literature with CCSD(T)/CBS. However, for bond dissociation, the results of the studied methods were highly divergent and not considered very effective when compared with experimental data from the literature. Notably, the G4-T-V2, G4-T-V3, and G4-T-V4 methods demonstrated good results for enthalpy of formation. These methods are particularly noteworthy because they utilise three energy coefficient adjustments and do not rely on higher-level correction.</p><h3>Methods</h3><p>We used the G4MP2-XK method along with 30 other composite methods developed by Martin and Semidalas.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145385624","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-28DOI: 10.1007/s00894-025-06539-w
Sumeyra Can, Mehmet Cinar, Ozlem Baris
Context
The structural and spectroscopic properties of drug groups can be determined using spectroscopic techniques and quantum chemical calculations. In this study, the structural and spectroscopic properties of nicardipine, a potent calcium channel blocker belonging to the dihydropyridine drug group, were investigated. For this purpose, FT-IR, Raman, NMR, and UV–Vis spectra of nicardipine hydrochloride were recorded. The ground state geometry optimization performed under vacuum conditions showed excellent agreement with X-ray diffraction (XRD) data, with root-mean-square deviations of bond lengths and bond angles calculated as 0.017 Å and 2.2°, respectively. Vibrational spectrum analysis revealed all characteristic modes, and theoretical B3LYP predictions largely agreed with experimental data. NMR spectra recorded in chloroform solvent confirmed the molecular structure of nicardipine. A comparison with previously reported data revealed significant differences in the proton NMR results, while the carbon NMR data showed consistency. The gauge-invariant atomic orbital (GIAO) method, commonly used for NMR spectrum prediction, yielded results within an acceptable error range despite minor inconsistencies arising from differences between experimental and computational conditions. In the UV–Vis spectrum of nicardipine HCl, an absorption peak corresponding to π-π* excitation was observed at 238 nm.
Method
The geometric optimizations and vibration spectra of the nicardipine drug molecule were performed using the Gaussian 09 program with the B3LYP functional and the 6–311 + + G(d,p) basis set. Visualizations were performed using the GaussView 5.0 interface program. NMR spectra were recorded in chloroform solvent using the GIAO method. Molecular electrostatic potential (MEP) maps and Mulliken atomic charges were analyzed to gain deeper insight into the electronic properties.
{"title":"Detailed investigation of the spectroscopic features of nicardipine hydrochloride using experimental and quantum chemical methods","authors":"Sumeyra Can, Mehmet Cinar, Ozlem Baris","doi":"10.1007/s00894-025-06539-w","DOIUrl":"10.1007/s00894-025-06539-w","url":null,"abstract":"<div><h3>Context</h3><p>The structural and spectroscopic properties of drug groups can be determined using spectroscopic techniques and quantum chemical calculations. In this study, the structural and spectroscopic properties of nicardipine, a potent calcium channel blocker belonging to the dihydropyridine drug group, were investigated. For this purpose, FT-IR, Raman, NMR, and UV–Vis spectra of nicardipine hydrochloride were recorded. The ground state geometry optimization performed under vacuum conditions showed excellent agreement with X-ray diffraction (XRD) data, with root-mean-square deviations of bond lengths and bond angles calculated as 0.017 Å and 2.2°, respectively. Vibrational spectrum analysis revealed all characteristic modes, and theoretical B3LYP predictions largely agreed with experimental data. NMR spectra recorded in chloroform solvent confirmed the molecular structure of nicardipine. A comparison with previously reported data revealed significant differences in the proton NMR results, while the carbon NMR data showed consistency. The gauge-invariant atomic orbital (GIAO) method, commonly used for NMR spectrum prediction, yielded results within an acceptable error range despite minor inconsistencies arising from differences between experimental and computational conditions. In the UV–Vis spectrum of nicardipine HCl, an absorption peak corresponding to π-π* excitation was observed at 238 nm.</p><h3>Method</h3><p>The geometric optimizations and vibration spectra of the nicardipine drug molecule were performed using the Gaussian 09 program with the B3LYP functional and the 6–311 + + G(d,p) basis set. Visualizations were performed using the GaussView 5.0 interface program. NMR spectra were recorded in chloroform solvent using the GIAO method. Molecular electrostatic potential (MEP) maps and Mulliken atomic charges were analyzed to gain deeper insight into the electronic properties.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145385632","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-27DOI: 10.1007/s00894-025-06545-y
Samantha E. Knoth, Daniel Tunega, Adelia J. A. Aquino
Context
Exo-tetrahydrodicyclopentadiene (exo-TCD) is a key component of Jet Propellant-10 (JP-10), a high-density hydrocarbon fuel extensively used in aerospace applications. The addition of aluminum particles enhances fuel performance and reactivity, making the understanding of initial decomposition pathways crucial. This study used density functional theory (DFT) calculations to investigate the initial hydrogen abstraction reactions in the decomposition of exo-TCD, with emphasis on radical formation processes. A significant aspect of this work is the role of the γ-Al2O3 surface in facilitating these reaction pathways, especially considering surface defects modeled by removing hydrogen from active hydroxyl groups. Five known active hydroxyl sites on γ-Al2O3 (A_Ia, A_Ib, A_IIA, B_IIb, and B_III) were used to construct complexes with exo-TCD. The formed complexes are primarily van der Waals interactions, with energies ranging from −11 to −20 kcal/mol and no substantial energy differences between configurations. The results indicate that hydrogen abstraction from the R4 site of exo-TCD is the most energetically favorable, owing to the molecular structure. Surface defects can boost reactivity by facilitating hydrogen abstraction, as seen in spontaneous H transfer to the active A_Ib site and low energetic barrier to the transition state of the H-abstraction of the B_IIb site. These findings improve the understanding of TCD decomposition and the catalytic role of γ-Al2O3, aiding the development of better propulsion fuels and energetic materials.
Methods
The calculations used the Perdew–Burke–Ernzerhof PBE exchange–correlation functional with split-valence polarization (SVP) and triple-zeta valence polarization (TZVP) basis sets, combined with the resolution of identity (RI) method to accelerate four-center electron repulsion integrals. The PBE results were benchmarked with the hybrid meta-GGA functional M06-2X. Dispersion correction D3 was applied throughout. All computations were performed using the Turbomole program.
{"title":"Surface interactions and radical generation in TCD decomposition: a DFT approach","authors":"Samantha E. Knoth, Daniel Tunega, Adelia J. A. Aquino","doi":"10.1007/s00894-025-06545-y","DOIUrl":"10.1007/s00894-025-06545-y","url":null,"abstract":"<div><h3>Context</h3><p>Exo-tetrahydrodicyclopentadiene (exo-TCD) is a key component of Jet Propellant-10 (JP-10), a high-density hydrocarbon fuel extensively used in aerospace applications. The addition of aluminum particles enhances fuel performance and reactivity, making the understanding of initial decomposition pathways crucial. This study used density functional theory (DFT) calculations to investigate the initial hydrogen abstraction reactions in the decomposition of exo-TCD, with emphasis on radical formation processes. A significant aspect of this work is the role of the γ-Al<sub>2</sub>O<sub>3</sub> surface in facilitating these reaction pathways, especially considering surface defects modeled by removing hydrogen from active hydroxyl groups. Five known active hydroxyl sites on γ-Al<sub>2</sub>O<sub>3</sub> (A_Ia, A_Ib, A_IIA, B_IIb, and B_III) were used to construct complexes with exo-TCD. The formed complexes are primarily van der Waals interactions, with energies ranging from −11 to −20 kcal/mol and no substantial energy differences between configurations. The results indicate that hydrogen abstraction from the R4 site of exo-TCD is the most energetically favorable, owing to the molecular structure. Surface defects can boost reactivity by facilitating hydrogen abstraction, as seen in spontaneous H transfer to the active A_Ib site and low energetic barrier to the transition state of the H-abstraction of the B_IIb site. These findings improve the understanding of TCD decomposition and the catalytic role of γ-Al<sub>2</sub>O<sub>3</sub>, aiding the development of better propulsion fuels and energetic materials.</p><h3>Methods</h3><p>The calculations used the Perdew–Burke–Ernzerhof PBE exchange–correlation functional with split-valence polarization (SVP) and triple-zeta valence polarization (TZVP) basis sets, combined with the resolution of identity (RI) method to accelerate four-center electron repulsion integrals. The PBE results were benchmarked with the hybrid meta-GGA functional M06-2X. Dispersion correction D3 was applied throughout. All computations were performed using the Turbomole program.\u0000</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12559037/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145375764","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-27DOI: 10.1007/s00894-025-06533-2
Isaac O. M. Magalhães, Nailton Martins Rodrigues, Daniel F. Scalabrini Machado, José A. S. Laranjeira, Julio R. Sambrano, João B. L. Martins
Context
H(_{varvec{2}})S and SO(_{varvec{2}}) are harmful gases found in nature and from industrial production. This study investigated their adsorption within IRMOF-1 cavities by means of periodic density functional theory. QTAIM-based IGMH isosurface and adsorption energy analyses provided insight into the SO(_{varvec{2}})(^{varvec{...}})H(_{varvec{2}})S interactions, where SO(_{varvec{2}}) is directed to the inorganic moiety. The mixed gases show significant interactions between SO(_{varvec{2}}) and H(_{varvec{2}})S. Regardless of the density functional approximation chosen, SO(_{varvec{2}}) exhibited a stronger adsorption preference over H(_{varvec{2}})S. B3LYP BSSE and dispersion-corrected interaction energies were consistently stronger than HSE06, whereas the B3LYP band gap presents larger values than HSE06. Grand Canonical Monte Carlo (GCMC) simulations further confirmed that SO(_{varvec{2}}) preferentially occupies the IRMOF-1 pores, with higher pressures enhancing adsorption.
Methods
Periodic optimizations were performed using the density functional approximations (DFAs) with B3LYP and HSE06 hybrid functionals, incorporating D3 dispersion and BSSE counterpoise corrections. Two basis sets were employed to analyze the B3LYP and HSE06 electronic and adsorption results: the pob-TZVP basis set and the modified m-6-311G(d) basis set. GCMC simulations were conducted to evaluate adsorption isotherms, while QTAIM and IGMH were used to account for the interaction behavior.
{"title":"Adsorption of H(_{2})S and SO(_{2}) in IRMOF-1: a computational study using DFT and GCMC simulations","authors":"Isaac O. M. Magalhães, Nailton Martins Rodrigues, Daniel F. Scalabrini Machado, José A. S. Laranjeira, Julio R. Sambrano, João B. L. Martins","doi":"10.1007/s00894-025-06533-2","DOIUrl":"10.1007/s00894-025-06533-2","url":null,"abstract":"<div><h3>Context</h3><p>H<span>(_{varvec{2}})</span>S and SO<span>(_{varvec{2}})</span> are harmful gases found in nature and from industrial production. This study investigated their adsorption within IRMOF-1 cavities by means of periodic density functional theory. QTAIM-based IGMH isosurface and adsorption energy analyses provided insight into the SO<span>(_{varvec{2}})</span> <span>(^{varvec{...}})</span>H<span>(_{varvec{2}})</span>S interactions, where SO<span>(_{varvec{2}})</span> is directed to the inorganic moiety. The mixed gases show significant interactions between SO<span>(_{varvec{2}})</span> and H<span>(_{varvec{2}})</span>S. Regardless of the density functional approximation chosen, SO<span>(_{varvec{2}})</span> exhibited a stronger adsorption preference over H<span>(_{varvec{2}})</span>S. B3LYP BSSE and dispersion-corrected interaction energies were consistently stronger than HSE06, whereas the B3LYP band gap presents larger values than HSE06. Grand Canonical Monte Carlo (GCMC) simulations further confirmed that SO<span>(_{varvec{2}})</span> preferentially occupies the IRMOF-1 pores, with higher pressures enhancing adsorption.</p><h3>Methods</h3><p>Periodic optimizations were performed using the density functional approximations (DFAs) with B3LYP and HSE06 hybrid functionals, incorporating D3 dispersion and BSSE counterpoise corrections. Two basis sets were employed to analyze the B3LYP and HSE06 electronic and adsorption results: the pob-TZVP basis set and the modified m-6-311G(d) basis set. GCMC simulations were conducted to evaluate adsorption isotherms, while QTAIM and IGMH were used to account for the interaction behavior.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145375704","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-27DOI: 10.1007/s00894-025-06546-x
Anibal Sierraalta, Rafael Añez, David S. Coll
Context
This work investigates the initial stages of the reaction in the conversion of chloromethane to light olefins (CMTO) on NaZSM-5 zeolite. Density functional theory (DFT) calculations indicate that the overall reaction is endothermic in bulk but exothermic on the surface. The formation of the :CH2 species is highly endothermic and has a high activation barrier. The migration of the CH3 group between oxygen atoms within the zeolite structure is energetically less demanding than the formation of the :CH2 species. Finally, the formation of ZOCH2CH3 from ZOCH3 + CH3Cl occurs through a two-step mechanism that includes a migration of the CH3 group with energy barriers of +95.2 and +148.2 kJ mol−1. The results obtained with ONIOM2 are qualitatively comparable to those obtained with VASP in terms of energy and geometry.
Method
DFT calculations were carried out using Gaussian 09 and VASP code. For cluster models, ONIOM2(DF:PM3) was used with the ωB97x-D, PBE, and M062X density functionals (DF) along with the relativistic Stevens effective core potential and its corresponding basis set. QST2, IRC, and frequency calculations are performed to characterize the transition states. For periodic calculations, structure optimization was carried out using the PBE functional with the D3 dispersion of Grimme, a plane-wave basis set, and the PAW representation of the atomic cores. The electronic structures of the systems were analyzed using PDOS.
{"title":"New considerations on the first steps for the conversion of chloromethane to olefins from a DFT perspective","authors":"Anibal Sierraalta, Rafael Añez, David S. Coll","doi":"10.1007/s00894-025-06546-x","DOIUrl":"10.1007/s00894-025-06546-x","url":null,"abstract":"<div><h3>Context</h3><p>This work investigates the initial stages of the reaction in the conversion of chloromethane to light olefins (CMTO) on NaZSM-5 zeolite. Density functional theory (DFT) calculations indicate that the overall reaction is endothermic in bulk but exothermic on the surface. The formation of the :CH<sub>2</sub> species is highly endothermic and has a high activation barrier. The migration of the CH<sub>3</sub> group between oxygen atoms within the zeolite structure is energetically less demanding than the formation of the :CH<sub>2</sub> species. Finally, the formation of ZOCH2CH3 from ZOCH3 + CH3Cl occurs through a two-step mechanism that includes a migration of the CH<sub>3</sub> group with energy barriers of +95.2 and +148.2 kJ mol<sup>−1</sup>. The results obtained with ONIOM2 are qualitatively comparable to those obtained with VASP in terms of energy and geometry.</p><h3>Method</h3><p>DFT calculations were carried out using Gaussian 09 and VASP code. For cluster models, ONIOM2(DF:PM3) was used with the ωB97x-D, PBE, and M062X density functionals (DF) along with the relativistic Stevens effective core potential and its corresponding basis set. QST2, IRC, and frequency calculations are performed to characterize the transition states. For periodic calculations, structure optimization was carried out using the PBE functional with the D3 dispersion of Grimme, a plane-wave basis set, and the PAW representation of the atomic cores. The electronic structures of the systems were analyzed using PDOS.\u0000</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145375746","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-27DOI: 10.1007/s00894-025-06547-w
Cassia R. Rambo, Samuel F. C. Paula, André L. M. Porto, Isac G. Rosset
Context
Hydroxylated derivatives of progesterone are biologically relevant compounds obtained through regio- and stereoselective biotransformations, but determining their stereochemical configuration is challenging due to subtle NMR chemical shift differences. In this study, monohydroxylated (prog-15-OH) and dihydroxylated (prog-7,15-di-OH) progesterone derivatives produced by Penicillium oxalicum CBMAI 1996 were analyzed. Quantum NMR calculations combined with statistical validation allowed the unambiguous identification of the 15β-OH and 7β,15β-di-OH isomers, with DP4+ probabilities of 100% and 99.8%, respectively. These findings confirm the stereochemical outcome of enzymatic hydroxylation and demonstrate the reliability of integrating computational NMR with experimental spectroscopy for resolving stereochemical ambiguities in complex steroidal frameworks.
Methods
Conformational searches were performed using the MMFF94 force field, followed by DFT geometry optimizations at the B3LYP/6-311G(2d,p) level with IEFPCM (chloroform). NMR shielding tensors were computed by the GIAO method at the mPW1PW91/6-311G(d,p) level with PCM (chloroform), and chemical shifts were referenced to TMS. Experimental and calculated 1H/13C shifts were compared using correlation coefficients (R2), mean absolute deviation (MAD), root mean square deviation (RMSD), and DP4+ probability analysis. All calculations were performed with Gaussian 16.
{"title":"Stereochemical determination of hydroxylated derivatives of progesterone using DFT NMR calculations and DP4+ probability","authors":"Cassia R. Rambo, Samuel F. C. Paula, André L. M. Porto, Isac G. Rosset","doi":"10.1007/s00894-025-06547-w","DOIUrl":"10.1007/s00894-025-06547-w","url":null,"abstract":"<div><h3>Context</h3><p>Hydroxylated derivatives of progesterone are biologically relevant compounds obtained through regio- and stereoselective biotransformations, but determining their stereochemical configuration is challenging due to subtle NMR chemical shift differences. In this study, monohydroxylated (prog-15-OH) and dihydroxylated (prog-7,15-di-OH) progesterone derivatives produced by <i>Penicillium oxalicum</i> CBMAI 1996 were analyzed. Quantum NMR calculations combined with statistical validation allowed the unambiguous identification of the 15β-OH and 7β,15β-di-OH isomers, with DP4+ probabilities of 100% and 99.8%, respectively. These findings confirm the stereochemical outcome of enzymatic hydroxylation and demonstrate the reliability of integrating computational NMR with experimental spectroscopy for resolving stereochemical ambiguities in complex steroidal frameworks.</p><h3>Methods</h3><p>Conformational searches were performed using the MMFF94 force field, followed by DFT geometry optimizations at the B3LYP/6-311G(2d,p) level with IEFPCM (chloroform). NMR shielding tensors were computed by the GIAO method at the mPW1PW91/6-311G(d,p) level with PCM (chloroform), and chemical shifts were referenced to TMS. Experimental and calculated <sup>1</sup>H/<sup>13</sup>C shifts were compared using correlation coefficients (<i>R</i><sup>2</sup>), mean absolute deviation (MAD), root mean square deviation (RMSD), and DP4+ probability analysis. All calculations were performed with Gaussian 16.\u0000</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145375720","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-27DOI: 10.1007/s00894-025-06554-x
Qiong Wu, Zhihao Zheng, Siyang Cheng, Juan Xu, Zusheng Hang, Weihua Zhu
The adsorption and detection of explosive and hazardous compounds are highly important for preventing accidents. In this study, two new synthesized cyclo[n]carbons (n = 10, C10. n = 14, C14) were employed as new sensors to adsorb six explosive (N3H, NH2NO2) or hazardous compounds (CO, HCHO, HCN, NH3). Based on the DFT calculations, first of all, it is found that all six adsorbates are adsorbed inside the conjugated ring of of C10 and C14 with the adsorption energies ranged from 15.7 to 38.9 kJ/mol, showing the typical weak-to-moderate physical adsorption characteristic. Then, C14 indicated the stronger adsorption toward six small molecules than C10 in general. And, both C10 and C14 showed the best selectivity toward the adsorption of NH2NO2 than others. Plus, the adsorption mechanisms to the same adsorbate are different sometimes for two cyclo[n]carbons, due to the difference in the shape, size, and type of chemical bonds of them. Finally, the short recovery time ranged from 8.7 fs to 0.024 ms, and significant changes in the spectrum also supports the high potential for using C10 and C14 as sensors for adsorbing and detecting these six explosive or hazardous compounds.
All of the calculations on the structures were carried out by using the Gaussian 09 software at the M06-2X/6-311G(d,p), M06-2X-D3/6-311G(d,p), and M06-2X-D3/6-311++G(d,p) level. In addition, further calculations on the properties and interactions were performed by using the Multiwfn software.
{"title":"Adsorption of explosive and hazardous compounds by cyclo[10]carbon and cyclo[14]carbon: A DFT study","authors":"Qiong Wu, Zhihao Zheng, Siyang Cheng, Juan Xu, Zusheng Hang, Weihua Zhu","doi":"10.1007/s00894-025-06554-x","DOIUrl":"10.1007/s00894-025-06554-x","url":null,"abstract":"<p>The adsorption and detection of explosive and hazardous compounds are highly important for preventing accidents. In this study, two new synthesized cyclo[n]carbons (<i>n</i> = 10, C<sub>10</sub>. <i>n</i> = 14, C<sub>14</sub>) were employed as new sensors to adsorb six explosive (N<sub>3</sub>H, NH<sub>2</sub>NO<sub>2</sub>) or hazardous compounds (CO, HCHO, HCN, NH<sub>3</sub>). Based on the DFT calculations, first of all, it is found that all six adsorbates are adsorbed inside the conjugated ring of of C<sub>10</sub> and C<sub>14</sub> with the adsorption energies ranged from 15.7 to 38.9 kJ/mol, showing the typical weak-to-moderate physical adsorption characteristic. Then, C<sub>14</sub> indicated the stronger adsorption toward six small molecules than C<sub>10</sub> in general. And, both C<sub>10</sub> and C<sub>14</sub> showed the best selectivity toward the adsorption of NH<sub>2</sub>NO<sub>2</sub> than others. Plus, the adsorption mechanisms to the same adsorbate are different sometimes for two cyclo[n]carbons, due to the difference in the shape, size, and type of chemical bonds of them. Finally, the short recovery time ranged from 8.7 fs to 0.024 ms, and significant changes in the spectrum also supports the high potential for using C<sub>10</sub> and C<sub>14</sub> as sensors for adsorbing and detecting these six explosive or hazardous compounds.</p><p>All of the calculations on the structures were carried out by using the Gaussian 09 software at the M06-2X/6-311G(d,p), M06-2X-D3/6-311G(d,p), and M06-2X-D3/6-311++G(d,p) level. In addition, further calculations on the properties and interactions were performed by using the Multiwfn software.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145375762","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-10-23DOI: 10.1007/s00894-025-06531-4
Manoj Vishwakarma, Ranjeet Vishwakarma, Kritika Garg, Ayon Bhattacharjee, Debanjan Bhattacharjee
Context
This article focuses on the temperature-dependent Raman spectroscopy analysis of a 7.O4O.7 dimeric liquid crystal, a single molecule. This dimeric compound consists of two symmetric monomers linked by a spacer. This compound exhibits two liquid crystalline phases: SmF and SmA. The article includes the shifting of peak position, variations in integral intensity, and changes in full width at half maximum (FWHM) with temperature fluctuations. The variations in Raman spectral characteristics with temperature provide evidence of charge displacement during phase transitions of the liquid crystal. The orientational and vibrational freedoms of molecules related to Raman bands in relation to temperature variation are also examined. Density functional theoretical approaches were employed to optimise the compound and to compare the Raman spectra at room temperature with experimental Raman data. The experimental and theoretical Raman results align exactly at room temperature, indicating that the basis set and functional employed in this investigation were valid for the 7.O4O.7 dimeric liquid crystal molecule analysed. The assignment of vibrational modes, observed in both experimental and DFT-generated data, was based on potential energy distribution (PED) utilising vibrational energy distribution analysis (VEDA) calculations.
Methods
All theoretical calculations conducted for the study utilised the Gaussian 09 software package. Gauss View 06 and VEDA 4 were employed for the analysis and visualisation of the compounds. The compound had been tuned for optimal performance using the B3LYP functional with the 6-31G (d,p) basis set. Experimental Raman spectroscopy was conducted using the Horiba Lab RAM HR Evolution Raman Spectrometer to obtain temperature-dependent Raman spectra.
{"title":"Exploration of a 7.O4O.7 dimeric liquid crystal using temperature-dependent Raman spectroscopy and comparison with density functional theory","authors":"Manoj Vishwakarma, Ranjeet Vishwakarma, Kritika Garg, Ayon Bhattacharjee, Debanjan Bhattacharjee","doi":"10.1007/s00894-025-06531-4","DOIUrl":"10.1007/s00894-025-06531-4","url":null,"abstract":"<div><h3>Context</h3><p>This article focuses on the temperature-dependent Raman spectroscopy analysis of a 7.O4O.7 dimeric liquid crystal, a single molecule. This dimeric compound consists of two symmetric monomers linked by a spacer. This compound exhibits two liquid crystalline phases: SmF and SmA. The article includes the shifting of peak position, variations in integral intensity, and changes in full width at half maximum (FWHM) with temperature fluctuations. The variations in Raman spectral characteristics with temperature provide evidence of charge displacement during phase transitions of the liquid crystal. The orientational and vibrational freedoms of molecules related to Raman bands in relation to temperature variation are also examined. Density functional theoretical approaches were employed to optimise the compound and to compare the Raman spectra at room temperature with experimental Raman data. The experimental and theoretical Raman results align exactly at room temperature, indicating that the basis set and functional employed in this investigation were valid for the 7.O4O.7 dimeric liquid crystal molecule analysed. The assignment of vibrational modes, observed in both experimental and DFT-generated data, was based on potential energy distribution (PED) utilising vibrational energy distribution analysis (VEDA) calculations.</p><h3>Methods</h3><p>All theoretical calculations conducted for the study utilised the Gaussian 09 software package. Gauss View 06 and VEDA 4 were employed for the analysis and visualisation of the compounds. The compound had been tuned for optimal performance using the B3LYP functional with the 6-31G (d,p) basis set. Experimental Raman spectroscopy was conducted using the Horiba Lab RAM HR Evolution Raman Spectrometer to obtain temperature-dependent Raman spectra.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145342748","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}