首页 > 最新文献

Journal of Molecular Modeling最新文献

英文 中文
Computational investigation of medicinal plants' active ingredient effects as a potential Zika virus treatment: molecular docking, molecular dynamics simulations, ADMET screening, and DFT. 药用植物活性成分对寨卡病毒潜在治疗作用的计算研究:分子对接、分子动力学模拟、ADMET筛选和DFT
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-04 DOI: 10.1007/s00894-025-06558-7
Farid Elbamtari, Etibaria Belghalia, Khadija Zaki, Marwa Alaqarbeh, Abdelkrim Guendouzi, M'barek Choukrad, Abdelouahid Sbai, Mohammed Bouachrine, Tahar Lakhlifi

Context: Recent outbreaks of the Zika virus (ZIKV) worldwide have underscored its growing epidemiological significance, leading to its recognition as an international health concern. The steady annual rise in ZIKV cases has transformed it into a major challenge for global public health systems. Despite ongoing efforts, the development of effective therapeutic agents against the virus remains difficult. Among the promising avenues for treatment are natural products, particularly those derived from medicinal and aromatic plants.

Methods: These substances act as reservoirs of beneficial chemical compounds that can contribute to developing effective therapies. This work used computer methods to examine 26 bioactive molecules derived from plants as potential Zika inhibitors. Baicalin, epicatechin gallate, epigallocatechin gallate, isoquercetin, and sophoroflavenone are plant-derived bioactive molecules that have demonstrated significant stability at the active site of the receptor examined (PDB code: 5TFR). They provided intense binding energies and were also stabilized at the active site of the target receptor by standard hydrogen bonds. These results were validated by molecular dynamics simulation at 500 ns. The molecules chosen to meet essential therapeutic criteria, such as those of Lipinski, have good ADMET characteristics and are not toxic. As a result, they have excellent pharmacokinetic properties and appreciable bioavailability. The findings of this research strongly suggest that these five molecules could be potential inhibitors of anti-Zika action in the future.

背景:最近在世界范围内爆发的寨卡病毒(ZIKV)凸显了其日益增长的流行病学意义,使其成为一个国际卫生问题。寨卡病毒病例每年稳步上升已使其成为全球公共卫生系统面临的重大挑战。尽管正在进行努力,但开发针对该病毒的有效治疗剂仍然很困难。有希望的治疗途径是天然产品,特别是从药用和芳香植物中提取的产品。方法:这些物质作为有益化合物的储存库,有助于开发有效的治疗方法。这项工作使用计算机方法检测了26种来自植物的生物活性分子,作为潜在的寨卡病毒抑制剂。黄芩苷、表儿茶素没食子酸酯、表没食子儿茶素没食子酸酯、异槲皮素和sophoroflavenone是植物源性的生物活性分子,在所检测的受体(PDB代码:5TFR)的活性位点上表现出显著的稳定性。它们提供了强烈的结合能,并通过标准氢键稳定在目标受体的活性位点。这些结果通过500 ns的分子动力学模拟得到了验证。选择满足基本治疗标准的分子,如Lipinski的分子,具有良好的ADMET特性并且无毒。因此,它们具有优良的药代动力学特性和可观的生物利用度。这项研究的结果强烈表明,这五种分子可能是未来抗寨卡病毒作用的潜在抑制剂。
{"title":"Computational investigation of medicinal plants' active ingredient effects as a potential Zika virus treatment: molecular docking, molecular dynamics simulations, ADMET screening, and DFT.","authors":"Farid Elbamtari, Etibaria Belghalia, Khadija Zaki, Marwa Alaqarbeh, Abdelkrim Guendouzi, M'barek Choukrad, Abdelouahid Sbai, Mohammed Bouachrine, Tahar Lakhlifi","doi":"10.1007/s00894-025-06558-7","DOIUrl":"https://doi.org/10.1007/s00894-025-06558-7","url":null,"abstract":"<p><strong>Context: </strong>Recent outbreaks of the Zika virus (ZIKV) worldwide have underscored its growing epidemiological significance, leading to its recognition as an international health concern. The steady annual rise in ZIKV cases has transformed it into a major challenge for global public health systems. Despite ongoing efforts, the development of effective therapeutic agents against the virus remains difficult. Among the promising avenues for treatment are natural products, particularly those derived from medicinal and aromatic plants.</p><p><strong>Methods: </strong>These substances act as reservoirs of beneficial chemical compounds that can contribute to developing effective therapies. This work used computer methods to examine 26 bioactive molecules derived from plants as potential Zika inhibitors. Baicalin, epicatechin gallate, epigallocatechin gallate, isoquercetin, and sophoroflavenone are plant-derived bioactive molecules that have demonstrated significant stability at the active site of the receptor examined (PDB code: 5TFR). They provided intense binding energies and were also stabilized at the active site of the target receptor by standard hydrogen bonds. These results were validated by molecular dynamics simulation at 500 ns. The molecules chosen to meet essential therapeutic criteria, such as those of Lipinski, have good ADMET characteristics and are not toxic. As a result, they have excellent pharmacokinetic properties and appreciable bioavailability. The findings of this research strongly suggest that these five molecules could be potential inhibitors of anti-Zika action in the future.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":"61"},"PeriodicalIF":2.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111786","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}
引用次数: 0
Terahertz-induced Berry curvature control of spin-selective transport in chiral DNA molecules. 手性DNA分子中自旋选择性输运的太赫兹诱导Berry曲率控制。
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-04 DOI: 10.1007/s00894-026-06637-3
Moses Udoisoh, Temitope Esther Olajide

Context: Chirality-induced spin selectivity (CISS) represents a remarkable quantum phenomenon whereby electron transmission through chiral molecules, such as DNA, becomes intrinsically spin-polarized even in the absence of magnetic fields. Despite extensive experimental verification of static CISS effects, achieving dynamic control over spin polarization remains an open challenge. Terahertz (THz) radiation offers a promising route to externally modulate molecular electronic structure on sub-picosecond timescales. In this study, we develop a theoretical model that unifies THz-driven Floquet dynamics with the spin-orbit coupling inherent to chiral DNA, thereby introducing the concept of Floquet-CISS, a light-induced regime of topologically controlled spin transport in biological helices.

Method: An effective low-energy Hamiltonian incorporating kinetic motion along the DNA helix, spin-orbit coupling, and the interaction with circularly polarized THz fields was formulated and solved using Floquet theory. The resulting quasi-energy spectra, Berry curvature, and spin polarization were numerically evaluated using plane-wave expansion and LAPACK-based diagonalization. The simulations reveal that THz fields dynamically reshape the Berry curvature, induce tunable spin-split Floquet bands, and produce helicity-dependent spin polarization exceeding 60%. These effects arise entirely from light-matter coupling without magnetic components, establishing DNA as a bio-topological spin filter capable of ultrafast, reversible spin control. The Floquet-CISS mechanism provides a theoretical blueprint for THz-programmable molecular spintronics and paves the way toward optically reconfigurable bio-quantum devices.

背景:手性诱导自旋选择性(CISS)代表了一种显著的量子现象,即电子通过手性分子(如DNA)的传输即使在没有磁场的情况下也会变得本质上自旋极化。尽管对静态CISS效应进行了大量的实验验证,但实现对自旋极化的动态控制仍然是一个开放的挑战。太赫兹(THz)辐射为在亚皮秒时间尺度上外部调制分子电子结构提供了一条很有前途的途径。在这项研究中,我们建立了一个理论模型,将太赫兹驱动的Floquet动力学与手性DNA固有的自旋轨道耦合结合起来,从而引入了Floquet- ciss的概念,这是一种光诱导的生物螺旋中拓扑控制的自旋输运机制。方法:利用Floquet理论建立了包含DNA螺旋运动、自旋-轨道耦合以及与圆极化太赫兹场相互作用的有效低能哈密顿量。利用平面波展开和基于lapack的对角化方法对得到的准能谱、Berry曲率和自旋极化进行了数值计算。仿真结果表明,太赫兹场可动态地重塑Berry曲率,诱导可调谐的自旋分裂Floquet带,并产生超过60%的螺旋相关自旋极化。这些效应完全是由光-物质耦合产生的,没有磁性成分,这使DNA成为一种生物拓扑自旋过滤器,能够超快、可逆地控制自旋。Floquet-CISS机制为太赫兹可编程分子自旋电子学提供了理论蓝图,并为光学可重构的生物量子器件铺平了道路。
{"title":"Terahertz-induced Berry curvature control of spin-selective transport in chiral DNA molecules.","authors":"Moses Udoisoh, Temitope Esther Olajide","doi":"10.1007/s00894-026-06637-3","DOIUrl":"https://doi.org/10.1007/s00894-026-06637-3","url":null,"abstract":"<p><strong>Context: </strong>Chirality-induced spin selectivity (CISS) represents a remarkable quantum phenomenon whereby electron transmission through chiral molecules, such as DNA, becomes intrinsically spin-polarized even in the absence of magnetic fields. Despite extensive experimental verification of static CISS effects, achieving dynamic control over spin polarization remains an open challenge. Terahertz (THz) radiation offers a promising route to externally modulate molecular electronic structure on sub-picosecond timescales. In this study, we develop a theoretical model that unifies THz-driven Floquet dynamics with the spin-orbit coupling inherent to chiral DNA, thereby introducing the concept of Floquet-CISS, a light-induced regime of topologically controlled spin transport in biological helices.</p><p><strong>Method: </strong>An effective low-energy Hamiltonian incorporating kinetic motion along the DNA helix, spin-orbit coupling, and the interaction with circularly polarized THz fields was formulated and solved using Floquet theory. The resulting quasi-energy spectra, Berry curvature, and spin polarization were numerically evaluated using plane-wave expansion and LAPACK-based diagonalization. The simulations reveal that THz fields dynamically reshape the Berry curvature, induce tunable spin-split Floquet bands, and produce helicity-dependent spin polarization exceeding 60%. These effects arise entirely from light-matter coupling without magnetic components, establishing DNA as a bio-topological spin filter capable of ultrafast, reversible spin control. The Floquet-CISS mechanism provides a theoretical blueprint for THz-programmable molecular spintronics and paves the way toward optically reconfigurable bio-quantum devices.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":"60"},"PeriodicalIF":2.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111807","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}
引用次数: 0
Effect of nanoalloying on dynamic thermophysical response of polycrystalline copper-tantalum 纳米合金化对多晶铜钽动态热物理响应的影响
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-026-06630-w
Mahesh Kumar Gupta, Santosh Kumar Rai, Vinay Panwar, R. P. Mahapatra, Abhishek Tevatia

Context

In the present work, the investigation of polycrystalline nanomaterials has been extended to a specific nanoalloy of copper and tantalum having a 9:1 atomic concentration. The study aims to analyze the influence of temperature and average grain size (AGS) on the mechanical behavior of the polycrystalline Cu-Ta nanoalloy. The results indicate that the critical grain size of polycrystalline 9Cu-Ta is smaller than that of pure Cu. The critical grain size of polycrystalline Cu (6.86 nm) is reduced to 3.89 nm with the addition of approximately 10% Ta atoms. This reduction is attributed to the combined effects of dislocation slip and subgrain strengthening mechanisms. Furthermore, the investigation highlights the variation of mechanical properties with increasing temperature and the influence of temperature on the critical grain size. The analysis also reveals the existence of distinct plastic deformation mechanisms corresponding to the critical grain size in the polycrystalline Cu-Ta nanoalloy.

Methods

Molecular dynamic simulation has been carried out under a fixed strain rate of 1.0 × 1010 s−1 for specifically analyzing the effect of temperature and average grain size (AGS) of the polycrystalline nanoalloy using embedded atom method potential (EAM). The polycrystalline structures with different grain sizes were generated using the Voronoi construction method. Simulations were carried out to evaluate the effect of temperature and grain size on the deformation behavior. The obtained data were analyzed to determine the critical grain size, variation in mechanical properties, and the associated deformation mechanisms of the polycrystalline 9Cu-Ta alloy.

在目前的工作中,对多晶纳米材料的研究已经扩展到具有9:1原子浓度的铜和钽的特定纳米合金。研究温度和平均晶粒尺寸对多晶Cu-Ta纳米合金力学行为的影响。结果表明:多晶9Cu-Ta的临界晶粒尺寸小于纯Cu的临界晶粒尺寸;加入约10%的Ta原子后,多晶Cu的临界晶粒尺寸由6.86 nm降至3.89 nm。这是位错滑移和亚晶强化机制共同作用的结果。此外,研究还强调了力学性能随温度升高的变化以及温度对临界晶粒尺寸的影响。分析还揭示了多晶Cu-Ta纳米合金在临界晶粒尺寸下存在不同的塑性变形机制。方法采用嵌入原子法电势(EAM)对多晶纳米合金在1.0 × 1010 s−1固定应变速率下进行分子动力学模拟,具体分析温度和平均晶粒尺寸(AGS)对多晶纳米合金的影响。采用Voronoi构造法生成了不同晶粒尺寸的多晶结构。模拟了温度和晶粒尺寸对合金变形行为的影响。对所获得的数据进行分析,以确定多晶9Cu-Ta合金的临界晶粒尺寸、力学性能变化以及相关的变形机制。
{"title":"Effect of nanoalloying on dynamic thermophysical response of polycrystalline copper-tantalum","authors":"Mahesh Kumar Gupta,&nbsp;Santosh Kumar Rai,&nbsp;Vinay Panwar,&nbsp;R. P. Mahapatra,&nbsp;Abhishek Tevatia","doi":"10.1007/s00894-026-06630-w","DOIUrl":"10.1007/s00894-026-06630-w","url":null,"abstract":"<div><h3>Context</h3><p>In the present work, the investigation of polycrystalline nanomaterials has been extended to a specific nanoalloy of copper and tantalum having a 9:1 atomic concentration. The study aims to analyze the influence of temperature and average grain size (AGS) on the mechanical behavior of the polycrystalline Cu-Ta nanoalloy. The results indicate that the critical grain size of polycrystalline 9Cu-Ta is smaller than that of pure Cu. The critical grain size of polycrystalline Cu (6.86 nm) is reduced to 3.89 nm with the addition of approximately 10% Ta atoms. This reduction is attributed to the combined effects of dislocation slip and subgrain strengthening mechanisms. Furthermore, the investigation highlights the variation of mechanical properties with increasing temperature and the influence of temperature on the critical grain size. The analysis also reveals the existence of distinct plastic deformation mechanisms corresponding to the critical grain size in the polycrystalline Cu-Ta nanoalloy.</p><h3>Methods</h3><p>Molecular dynamic simulation has been carried out under a fixed strain rate of 1.0 × 10<sup>10</sup> s<sup>−1</sup> for specifically analyzing the effect of temperature and average grain size (AGS) of the polycrystalline nanoalloy using embedded atom method potential (EAM). The polycrystalline structures with different grain sizes were generated using the Voronoi construction method. Simulations were carried out to evaluate the effect of temperature and grain size on the deformation behavior. The obtained data were analyzed to determine the critical grain size, variation in mechanical properties, and the associated deformation mechanisms of the polycrystalline 9Cu-Ta alloy.\u0000</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099219","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}
引用次数: 0
First-principles insights into lead-free K2AgRhX6 (X = F, Cl, Br, I) halide double perovskites as stable platforms for next-generation optoelectronic and energy conversion devices 无铅K2AgRhX6 (X = F, Cl, Br, I)卤化物双钙钛矿作为下一代光电和能量转换器件稳定平台的第一性原理见解
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-025-06626-y
Farooq Afzaal, Rashid Jalil, Ibtsam Riaz, Nawaz Muhammad, G. Murtaza, Maha Naeem, Muhammad Moin, Hafiz Irfan Ali

Context

Double perovskites have emerged as promising candidates for renewable energy technologies due to their structural simplicity and thermodynamic stability. Among them, K2AgRhF6 is the most stable (-2.54 eV/atom), consistent with its highest bulk modulus (64.15 GPa), tolerance factor (0.85), and octahedral factor (0.86). Elastic analysis indicates ductile behavior for K2AgRhF6 (ν = 0.35, B/G = 3.13) and K2AgRhBr6 (ν = 0.28, B/G = 1.99), while K2AgRhCl6 (ν = 0.26, B/G = 1.75) and K2AgRhI6 (ν = 0.13, B/G = 2.37) lie near the brittle-ductile threshold. Band structure calculations reveal semiconducting gaps of 2.56 eV (F), 2.03 eV (Cl), 1.44 eV (Br), and 0.55 eV (I), with K2AgRhBr6 and K2AgRhI6 exhibiting strong optical absorption in the visible spectrum. Thermoelectric analysis yields figures of merit approaching 0.75 at room temperature across the series, highlighting their efficiency in energy conversion. Collectively, these findings position K2AgRhX6 halide double perovskites as robust, lead-free multifunctional materials with integrated structural stability, tunable optoelectronic response, and promising thermoelectric efficiency for next-generation optoelectronic devices.

Method

In this work, density functional theory (DFT) calculations within the WIEN2k framework were used to explore the structural, mechanical, electronic, optical, and thermoelectric properties of halide-based double perovskites K2AgRhX6 (X = F, Cl, Br, I). All compounds crystallize in the cubic Fm3ˉm (225) space group, and their negative formation energies confirm thermodynamic stability.

双钙钛矿因其结构简单和热力学稳定性而成为可再生能源技术的有前途的候选者。其中,K2AgRhF6最稳定(-2.54 eV/原子),具有最高的体积模量(64.15 GPa)、容差因子(0.85)和八面体因子(0.86)。弹性分析表明,K2AgRhF6 (ν = 0.35, B/G = 3.13)和K2AgRhBr6 (ν = 0.28, B/G = 1.99)具有延性,而K2AgRhCl6 (ν = 0.26, B/G = 1.75)和K2AgRhI6 (ν = 0.13, B/G = 2.37)处于脆性-延性阈值附近。能带结构计算表明,K2AgRhBr6和k2aghi6的半导体隙分别为2.56 eV (F)、2.03 eV (Cl)、1.44 eV (Br)和0.55 eV (I),在可见光谱中表现出较强的光吸收。热电分析结果表明,在室温下,整个系列的优点接近0.75,突出了它们在能量转换方面的效率。总的来说,这些发现将K2AgRhX6卤化物双钙钛矿定位为坚固的无铅多功能材料,具有集成的结构稳定性,可调谐的光电响应,以及下一代光电器件的热电效率。方法利用WIEN2k框架下的密度泛函理论(DFT)计算,研究了卤化物基双钙钛矿K2AgRhX6 (X = F, Cl, Br, I)的结构、力学、电子、光学和热电性质。所有化合物均在立方Fm3 - m(225)空间群中结晶,它们的负生成能证实了热力学稳定性。
{"title":"First-principles insights into lead-free K2AgRhX6 (X = F, Cl, Br, I) halide double perovskites as stable platforms for next-generation optoelectronic and energy conversion devices","authors":"Farooq Afzaal,&nbsp;Rashid Jalil,&nbsp;Ibtsam Riaz,&nbsp;Nawaz Muhammad,&nbsp;G. Murtaza,&nbsp;Maha Naeem,&nbsp;Muhammad Moin,&nbsp;Hafiz Irfan Ali","doi":"10.1007/s00894-025-06626-y","DOIUrl":"10.1007/s00894-025-06626-y","url":null,"abstract":"<div><h3>Context</h3><p>Double perovskites have emerged as promising candidates for renewable energy technologies due to their structural simplicity and thermodynamic stability. Among them, K<sub>2</sub>AgRhF<sub>6</sub> is the most stable (-2.54 eV/atom), consistent with its highest bulk modulus (64.15 GPa), tolerance factor (0.85), and octahedral factor (0.86). Elastic analysis indicates ductile behavior for K<sub>2</sub>AgRhF<sub>6</sub> (ν = 0.35, B/G = 3.13) and K<sub>2</sub>AgRhBr<sub>6</sub> (ν = 0.28, B/G = 1.99), while K<sub>2</sub>AgRhCl<sub>6</sub> (ν = 0.26, B/G = 1.75) and K<sub>2</sub>AgRhI<sub>6</sub> (ν = 0.13, B/G = 2.37) lie near the brittle-ductile threshold. Band structure calculations reveal semiconducting gaps of 2.56 eV (F), 2.03 eV (Cl), 1.44 eV (Br), and 0.55 eV (I), with K<sub>2</sub>AgRhBr<sub>6</sub> and K<sub>2</sub>AgRhI<sub>6</sub> exhibiting strong optical absorption in the visible spectrum. Thermoelectric analysis yields figures of merit approaching 0.75 at room temperature across the series, highlighting their efficiency in energy conversion. Collectively, these findings position K<sub>2</sub>AgRhX<sub>6</sub> halide double perovskites as robust, lead-free multifunctional materials with integrated structural stability, tunable optoelectronic response, and promising thermoelectric efficiency for next-generation optoelectronic devices.</p><h3>Method</h3><p>In this work, density functional theory (DFT) calculations within the WIEN2k framework were used to explore the structural, mechanical, electronic, optical, and thermoelectric properties of halide-based double perovskites K<sub>2</sub>AgRhX<sub>6</sub> (X = F, Cl, Br, I). All compounds crystallize in the cubic Fm3ˉm (225) space group, and their negative formation energies confirm thermodynamic stability.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099221","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}
引用次数: 0
Analysis of molecular conjugation influence on color characteristics of indigo compounds 分子偶联对靛蓝化合物颜色特性的影响分析
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-026-06638-2
Zengbo Ke, Jiahao Zhang, Hao Li, Xiaolong Li, Chengfang Qiao, Youying Di

Context

Indigo dyes are historically significant and possess a unique π-conjugated core, making them valuable for both traditional pigments and emerging applications in organic electronics. A fundamental challenge is understanding how subtle molecular modifications, particularly substituent effects, quantitatively influence their conjugation extent and resulting color properties. This study systematically investigates the parent indigo and three N, N′-substituted derivatives (phenyl, ethyl, and vinyl) to elucidate the precise relationship between molecular structure, electronic properties, and visible light absorption. Density functional theory (DFT) and Time-dependent DFT (TD-DFT) calculations reveal how substituents modulate the HOMO–LUMO gap and intramolecular interactions, directly correlating with calculated absorption wavelengths (628 to 807 nm) and predicted colors (medium blue to green-cyan) via complementary color theory. These findings provide a quantitative framework for designing indigo-based dyes with targeted optical properties, including near-infrared absorption.

Methods

Density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations were performed. Ground-state geometry optimizations and frequency calculations were carried out using the B3LYP-D3(BJ) functional with the def2TZVP basis set. Excited-state calculations employed the CAM-B3LYP-D3(BJ)/def2TZVP level, with the IEFPCM solvation model simulating aqueous conditions. The Multiwfn 3.8 and VMD 1.9.3 software packages were used for interaction region indicator (IRI) analysis, electrostatic potential (ESP) mapping, electron density difference (EDD) analysis, hole–electron analysis, and color prediction based on absorption spectra.

靛蓝染料具有重要的历史意义,具有独特的π共轭核,使其在传统颜料和有机电子领域的新兴应用中都有价值。一个基本的挑战是理解细微的分子修饰,特别是取代基效应,如何定量地影响它们的共轭程度和所产生的颜色性质。本研究系统地研究了母体靛蓝和三个N, N '取代衍生物(苯基,乙基和乙烯基),以阐明分子结构,电子性质和可见光吸收之间的精确关系。密度泛函理论(DFT)和时间依赖DFT (TD-DFT)计算揭示了取代基如何调节HOMO-LUMO间隙和分子内相互作用,直接与计算的吸收波长(628至807 nm)和通过互补色理论预测的颜色(中蓝至绿青色)相关。这些发现为设计具有目标光学特性(包括近红外吸收)的靛蓝基染料提供了定量框架。方法进行密度泛函理论(DFT)和时变DFT (TD-DFT)计算。利用def2TZVP基集的B3LYP-D3(BJ)泛函进行了基态几何优化和频率计算。激发态计算采用CAM-B3LYP-D3(BJ)/def2TZVP水平,IEFPCM溶剂化模型模拟水溶液条件。采用Multiwfn 3.8和VMD 1.9.3软件包进行相互作用区指示(IRI)分析、静电电位(ESP)作图、电子密度差(EDD)分析、空穴电子分析和基于吸收光谱的颜色预测。
{"title":"Analysis of molecular conjugation influence on color characteristics of indigo compounds","authors":"Zengbo Ke,&nbsp;Jiahao Zhang,&nbsp;Hao Li,&nbsp;Xiaolong Li,&nbsp;Chengfang Qiao,&nbsp;Youying Di","doi":"10.1007/s00894-026-06638-2","DOIUrl":"10.1007/s00894-026-06638-2","url":null,"abstract":"<div><h3>Context</h3><p>Indigo dyes are historically significant and possess a unique π-conjugated core, making them valuable for both traditional pigments and emerging applications in organic electronics. A fundamental challenge is understanding how subtle molecular modifications, particularly substituent effects, quantitatively influence their conjugation extent and resulting color properties. This study systematically investigates the parent indigo and three <i>N</i>, <i>N</i>′-substituted derivatives (phenyl, ethyl, and vinyl) to elucidate the precise relationship between molecular structure, electronic properties, and visible light absorption. Density functional theory (DFT) and Time-dependent DFT (TD-DFT) calculations reveal how substituents modulate the HOMO–LUMO gap and intramolecular interactions, directly correlating with calculated absorption wavelengths (628 to 807 nm) and predicted colors (medium blue to green-cyan) via complementary color theory. These findings provide a quantitative framework for designing indigo-based dyes with targeted optical properties, including near-infrared absorption.</p><h3>Methods</h3><p>Density functional theory (DFT) and time-dependent DFT (TD-DFT) calculations were performed. Ground-state geometry optimizations and frequency calculations were carried out using the B3LYP-D3(BJ) functional with the def2TZVP basis set. Excited-state calculations employed the CAM-B3LYP-D3(BJ)/def2TZVP level, with the IEFPCM solvation model simulating aqueous conditions. The Multiwfn 3.8 and VMD 1.9.3 software packages were used for interaction region indicator (IRI) analysis, electrostatic potential (ESP) mapping, electron density difference (EDD) analysis, hole–electron analysis, and color prediction based on absorption spectra.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099216","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}
引用次数: 0
Effect of single-solvent systems on the predicted crystal morphology of α-CL-20: a molecular dynamics and modified attachment energy study 单溶剂体系对α-CL-20预测晶体形态的影响:分子动力学和修饰的附着能研究
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-026-06643-5
Zhijiang Yue, Zhenfeng Jia, Qingying Duan, Shan Sha, Xianfeng Wei

A modified attachment energy model combined with classical molecular dynamics simulations is employed to study the solvent-dependent crystal growth of α-CL-20 in six solvent systems commonly used for the preparation of CL-20: maleic acid, pentaerythritol, glycine, citric acid, polyvinyl alcohol, and butyl carbamate. Surface roughness is quantified via the ratio of solvent-accessible area to geometric cross-sectional area, and electrostatic potential maps and radial distribution functions are used to characterize the polarity and non-bonded interactions at the interfaces. The results show that the interactions between α-CL-20 and solvent are dominated by hydrogen bonding, with some contributions from van der Waals contacts. The calculation results based on modified attachment energy theory show that in all six solvent systems, the growth of the (021) face is substantially inhibited, while the (020), (102), (111), and (002) faces tend to remain as persistent, moderately growing facets. The predicted morphologies in solutions evolve from the prismatic vacuum habit toward nearly polyhedral, spheroidal crystals, which suggest that strong, face-selective hydrogen bonding on rough, polar faces can serve as an effective microscopic mechanism for solvent-based morphology engineering of CL-20.

The Materials Studio software was used to calculate crystal growth, electrostatic potential, and radial distribution functions of α-CL-20. The VISTA software was used to predict morphologies of α-CL-20 in six solutions.

采用改进的附着能模型结合经典分子动力学模拟,研究了α-CL-20在马来酸、季戊四醇、甘氨酸、柠檬酸、聚乙烯醇和氨基甲酸丁酯等制备CL-20常用的六种溶剂体系中的溶剂依赖性晶体生长。表面粗糙度通过溶剂可及面积与几何横截面积的比值来量化,静电势图和径向分布函数用于表征界面上的极性和非键相互作用。结果表明:α-CL-20与溶剂的相互作用以氢键为主,范德华接触也有一定贡献;基于修正的附着能理论的计算结果表明,在所有6种溶剂体系中,(021)面生长受到抑制,而(020)、(102)、(111)和(002)面倾向于保持持久的、中等生长的面。预测溶液中的形态从棱柱形真空习惯演变为接近多面体的球形晶体,这表明在粗糙的极性表面上的强、面选择性氢键可以作为CL-20溶剂基形貌工程的有效微观机制。利用Materials Studio软件计算α-CL-20的晶体生长、静电势和径向分布函数。采用VISTA软件预测α-CL-20在6种溶液中的形态。
{"title":"Effect of single-solvent systems on the predicted crystal morphology of α-CL-20: a molecular dynamics and modified attachment energy study","authors":"Zhijiang Yue,&nbsp;Zhenfeng Jia,&nbsp;Qingying Duan,&nbsp;Shan Sha,&nbsp;Xianfeng Wei","doi":"10.1007/s00894-026-06643-5","DOIUrl":"10.1007/s00894-026-06643-5","url":null,"abstract":"<p>A modified attachment energy model combined with classical molecular dynamics simulations is employed to study the solvent-dependent crystal growth of α-CL-20 in six solvent systems commonly used for the preparation of CL-20: maleic acid, pentaerythritol, glycine, citric acid, polyvinyl alcohol, and butyl carbamate. Surface roughness is quantified via the ratio of solvent-accessible area to geometric cross-sectional area, and electrostatic potential maps and radial distribution functions are used to characterize the polarity and non-bonded interactions at the interfaces. The results show that the interactions between α-CL-20 and solvent are dominated by hydrogen bonding, with some contributions from van der Waals contacts. The calculation results based on modified attachment energy theory show that in all six solvent systems, the growth of the (021) face is substantially inhibited, while the (020), (102), (111), and (002) faces tend to remain as persistent, moderately growing facets. The predicted morphologies in solutions evolve from the prismatic vacuum habit toward nearly polyhedral, spheroidal crystals, which suggest that strong, face-selective hydrogen bonding on rough, polar faces can serve as an effective microscopic mechanism for solvent-based morphology engineering of CL-20.</p><p>The Materials Studio software was used to calculate crystal growth, electrostatic potential, and radial distribution functions of α-CL-20. The VISTA software was used to predict morphologies of α-CL-20 in six solutions.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099218","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}
引用次数: 0
An atomic-level investigation into the barrier properties of C3N5-modified epoxy composite coatings for concrete structures 混凝土结构用c3n5改性环氧复合涂料阻隔性能的原子水平研究
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-025-06628-w
Yongde Yao, Fujian Tang, Yufang He, Zhibin Lin

Context

Nanomaterial-modified epoxy coatings are considered an effective method for protecting concrete. Molecular dynamics (MD) simulations were used to study the effect of C3N5 and its derivative C3N5-OH on the properties of epoxy coating, as well as the interfacial properties between epoxy coating and concrete. Three MD models were established and their properties were investigated in both dry and sodium chloride solution. Results indicate that doping the epoxy coating with C3N5 and C3N5-OH reduces its free volume. The strength increases from 0.42 GPa to 0.70 GPa and 0.76 GPa after doping with C3N5 and C3N5-OH, respectively. Similarly, the elastic modulus increases from 3.09 GPa to 4.26 GPa and 4.95 GPa, respectively. The addition of C3N5 and C3N5-OH increases the interaction energy of epoxy from 99.3 kcal/mol to 132.2 kcal/mol and 167.2 kcal/mol, respectively. Correspondingly, the adsorption distance decreases. The doping of C3N5 effectively inhibits the diffusion of sodium chloride solution, mitigating its impact on the properties of the epoxy/concrete interface. This study demonstrates that C3N5-modified epoxy coatings are a promising solution for extending the service life of concrete structures.

Methods

The initial molecular model construction, cross-linked network generation, and the addition of NaCl solution were all completed using the Materials Studio software. MD simulations were conducted using LAMMPS software. The CVFF force field was selected for the epoxy resin composite layer, while the ClayFF force field was used for the CSH matrix.

纳米改性环氧涂料被认为是一种有效的混凝土防护方法。采用分子动力学(MD)模拟研究了C3N5及其衍生物C3N5- oh对环氧涂料性能的影响,以及环氧涂料与混凝土的界面性能。建立了三种MD模型,研究了它们在干溶液和氯化钠溶液中的性能。结果表明,C3N5和C3N5- oh的掺杂降低了环氧涂层的自由体积。C3N5和C3N5- oh掺杂后,强度分别从0.42 GPa提高到0.70 GPa和0.76 GPa。同样,弹性模量也从3.09 GPa增加到4.26 GPa和4.95 GPa。C3N5和C3N5- oh的加入使环氧树脂的相互作用能分别从99.3 kcal/mol提高到132.2 kcal/mol和167.2 kcal/mol。相应的,吸附距离减小。C3N5的掺杂有效抑制了氯化钠溶液的扩散,减轻了其对环氧/混凝土界面性能的影响。该研究表明,c3n5改性环氧涂料是一种很有前途的延长混凝土结构使用寿命的解决方案。方法采用Materials Studio软件完成初始分子模型构建、交联网络生成和NaCl溶液的加入。利用LAMMPS软件进行了MD模拟。环氧树脂复合层采用CVFF力场,CSH基体采用ClayFF力场。
{"title":"An atomic-level investigation into the barrier properties of C3N5-modified epoxy composite coatings for concrete structures","authors":"Yongde Yao,&nbsp;Fujian Tang,&nbsp;Yufang He,&nbsp;Zhibin Lin","doi":"10.1007/s00894-025-06628-w","DOIUrl":"10.1007/s00894-025-06628-w","url":null,"abstract":"<div><h3>Context</h3><p>Nanomaterial-modified epoxy coatings are considered an effective method for protecting concrete. Molecular dynamics (MD) simulations were used to study the effect of C<sub>3</sub>N<sub>5</sub> and its derivative C<sub>3</sub>N<sub>5</sub>-OH on the properties of epoxy coating, as well as the interfacial properties between epoxy coating and concrete. Three MD models were established and their properties were investigated in both dry and sodium chloride solution. Results indicate that doping the epoxy coating with C<sub>3</sub>N<sub>5</sub> and C<sub>3</sub>N<sub>5</sub>-OH reduces its free volume. The strength increases from 0.42 GPa to 0.70 GPa and 0.76 GPa after doping with C<sub>3</sub>N<sub>5</sub> and C<sub>3</sub>N<sub>5</sub>-OH, respectively. Similarly, the elastic modulus increases from 3.09 GPa to 4.26 GPa and 4.95 GPa, respectively. The addition of C<sub>3</sub>N<sub>5</sub> and C<sub>3</sub>N<sub>5</sub>-OH increases the interaction energy of epoxy from 99.3 kcal/mol to 132.2 kcal/mol and 167.2 kcal/mol, respectively. Correspondingly, the adsorption distance decreases. The doping of C<sub>3</sub>N<sub>5</sub> effectively inhibits the diffusion of sodium chloride solution, mitigating its impact on the properties of the epoxy/concrete interface. This study demonstrates that C<sub>3</sub>N<sub>5</sub>-modified epoxy coatings are a promising solution for extending the service life of concrete structures.</p><h3>Methods</h3><p>The initial molecular model construction, cross-linked network generation, and the addition of NaCl solution were all completed using the Materials Studio software. MD simulations were conducted using LAMMPS software. The CVFF force field was selected for the epoxy resin composite layer, while the ClayFF force field was used for the CSH matrix.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099220","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}
引用次数: 0
Theoretical study on the effects of intermolecular hydrogen bonding and solvent on the Raman spectroscopy of furfural in transformer oil 分子间氢键和溶剂对变压器油中糠醛拉曼光谱影响的理论研究
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-025-06622-2
Haiyang Shi, Guochao Gan, Shizheng Sun, Xue Han, Hao Ding, Shuang Liu, Yingzhou Huang, Kai Li, Weigen Chen

Context

The application of Raman spectroscopy for detecting furfural dissolved in transformer oil represents a highly promising approach for online monitoring for assessing the aging process of transformer insulating paper. Combining density functional theory (DFT) computational simulations with experimental measurements, the characteristic Raman spectra of the furfural molecule were analyzed, and the assignments of its key vibrational modes were determined. Furfural molecular clusters of varying sizes were constructed to investigate the effect of intermolecular hydrogen bonding forces on the Raman signals. The Raman spectra of furfural in solvents with different dielectric constants were measured. Comparison with simulation results indicates that enhanced solvent polarity induces a rearrangement of the electron cloud around the C=O bond, an increase in dipole moment, and a narrowing of the band gap. These changes thereby result in a red shift of the Raman peak position and an enhancement of scattering intensity. This work deepens the understanding, at the microscopic level, of the aggregation behavior of furfural in transformer oil and the influence of the solvent environment on its Raman characteristics.

Methods

All quantum chemical calculations were conducted using the Gaussian 16W program, and molecular structures were built using GaussView 6.0. The ground-state geometry of furfural was optimized using DFT with the B3LYP functional. The solvent effect caused by solvents with different dielectric constants was simulated and calculated using the polarized continuum model (PCM) of furfural. For all computations, the 6-311+G (2d, p) basis set was employed for C, H, and O atoms.

摘要利用拉曼光谱技术检测变压器油中溶解的糠醛,为在线监测变压器绝缘纸老化过程提供了一种很有前景的方法。结合密度泛函理论(DFT)计算模拟和实验测量,分析了糠醛分子的特征拉曼光谱,确定了其关键振动模式的归属。构建不同大小的糠醛分子簇,研究分子间氢键力对拉曼信号的影响。测定了糠醛在不同介电常数溶剂中的拉曼光谱。与模拟结果的比较表明,溶剂极性的增强引起了C=O键周围电子云的重排,偶极矩的增加和带隙的缩小。这些变化导致拉曼峰位置的红移和散射强度的增强。本研究在微观层面上加深了对糠醛在变压器油中的聚集行为以及溶剂环境对其拉曼特性的影响的认识。方法采用高斯16W程序进行量子化学计算,采用GaussView 6.0构建分子结构。利用B3LYP泛函对糠醛的基态几何结构进行了优化。利用糠醛的极化连续介质模型(PCM)对不同介电常数的溶剂效应进行了模拟和计算。对于所有的计算,C、H和O原子采用6-311+G (2d, p)基集。
{"title":"Theoretical study on the effects of intermolecular hydrogen bonding and solvent on the Raman spectroscopy of furfural in transformer oil","authors":"Haiyang Shi,&nbsp;Guochao Gan,&nbsp;Shizheng Sun,&nbsp;Xue Han,&nbsp;Hao Ding,&nbsp;Shuang Liu,&nbsp;Yingzhou Huang,&nbsp;Kai Li,&nbsp;Weigen Chen","doi":"10.1007/s00894-025-06622-2","DOIUrl":"10.1007/s00894-025-06622-2","url":null,"abstract":"<div><h3>Context</h3><p>The application of Raman spectroscopy for detecting furfural dissolved in transformer oil represents a highly promising approach for online monitoring for assessing the aging process of transformer insulating paper. Combining density functional theory (DFT) computational simulations with experimental measurements, the characteristic Raman spectra of the furfural molecule were analyzed, and the assignments of its key vibrational modes were determined. Furfural molecular clusters of varying sizes were constructed to investigate the effect of intermolecular hydrogen bonding forces on the Raman signals. The Raman spectra of furfural in solvents with different dielectric constants were measured. Comparison with simulation results indicates that enhanced solvent polarity induces a rearrangement of the electron cloud around the C=O bond, an increase in dipole moment, and a narrowing of the band gap. These changes thereby result in a red shift of the Raman peak position and an enhancement of scattering intensity. This work deepens the understanding, at the microscopic level, of the aggregation behavior of furfural in transformer oil and the influence of the solvent environment on its Raman characteristics.</p><h3>Methods</h3><p>All quantum chemical calculations were conducted using the Gaussian 16W program, and molecular structures were built using GaussView 6.0. The ground-state geometry of furfural was optimized using DFT with the B3LYP functional. The solvent effect caused by solvents with different dielectric constants was simulated and calculated using the polarized continuum model (PCM) of furfural. For all computations, the 6-311+G (2d, p) basis set was employed for C, H, and O atoms.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099215","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}
引用次数: 0
Molecular dynamics simulation of the influence of oxidized functional groups on asphalt-aggregate bonding properties 氧化官能团对沥青-骨料键合性能影响的分子动力学模拟
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-02-03 DOI: 10.1007/s00894-026-06644-4
Maoping Ran, Jiang Zuo, Xinglin Zhou, Jing Yuan, Xinxing Zhou
<div><h3>Context</h3><p>To elucidate the mechanism by which oxidized functional groups influence asphalt-aggregate adhesion behavior at the microscopic level, molecular dynamics (MD) simulations were employed based on the Derek four-component model—a model well-matched with real asphalt’s density and thermal expansion coefficient. Asphalt molecular models containing 12 characteristic molecules were constructed, with two stable oxidized functional groups (carbonyl, C=O; sulfoxyl, S=O) introduced to represent asphalt aging. Three concentration gradients (0%, 50%, 100%) were designed by adjusting the number of functional groups. Meanwhile, silica (SiO<sub>2</sub>) aggregate models were cleaved along the (0 0 1) crystal plane, transformed into an orthorhombic structure, and passivated with hydrogen atoms to simulate real aggregate surface properties. An asphalt-aggregate interface system with a “SiO<sub>2</sub> layer-10 Å vacuum layer-asphalt layer-50 Å vacuum layer” sandwich structure was established. The model rationality was verified by three key indicators: density (1.011~1.079 g/cm<sup>3</sup>), solubility parameter (18.604~19.237 (J cm<sup>−3</sup>)<sup>0.5</sup>), and glass transition temperature (237.97~247.61 K), ensuring consistency with published experimental and simulation data. The effects of functional group type and concentration on asphalt-aggregate adhesion were analyzed by evaluating interfacial energies (interaction energy, van der Waals energy, electrostatic energy). As the concentration of carbonyl (C=O) groups increased, the interaction energy, van der Waals energy, and electrostatic energy in the asphalt-aggregate system all showed an upward trend. The increase at 100% concentration was greater than that at 50% concentration, with electrostatic energy exhibiting the largest overall increase. As the sulfoxyl group (S=O) content increased, the interaction energy and van der Waals energy in the asphalt-aggregate system decreased, while the electrostatic energy increased. When both carbonyl (C=O) and sulfoxyl (S=O) groups coexist, the variation patterns of various interfacial energies resemble those observed under sulfoxide influence alone. This indicates that, under identical conditions, sulfoxyl (S=O) exerts a greater impact on interfacial energies. Oxidized functional groups significantly influence the adhesion properties between asphalt and aggregates.</p><h3>Methods</h3><p>To investigate the effects of oxygen-containing functional groups with different concentrations and types on the asphalt-aggregate adhesion behavior, asphalt molecular models were first constructed via the molecular assembly method using the Amorphous Cell module in Materials Studio 2023 software, referencing TLC-FID data of asphalt four components to determine the molecular quantity of each component. A total of 15 asphalt models were established, covering three functional group types (C=O, S=O, C=O + S=O) and three concentration gradients. SiO<sub>2</sub>
为了阐明氧化官能团在微观水平上影响沥青-骨料粘附行为的机理,基于Derek四组分模型进行了分子动力学(MD)模拟,该模型与真实沥青的密度和热膨胀系数很好地匹配。构建了包含12个特征分子的沥青分子模型,并引入两个稳定的氧化官能团(羰基,C=O;亚砜基,S=O)来表征沥青老化。通过调整官能团的数量,设计3个浓度梯度(0%、50%、100%)。同时,将二氧化硅(SiO2)聚集体模型沿(0 ~ 0 ~ 1)晶面切割,转变为正交结构,并进行氢原子钝化,模拟真实的聚集体表面性质。建立了“SiO2层-10 Å真空层-沥青层-50 Å真空层”夹层结构的沥青-骨料界面体系。通过密度(1.011~1.079 g/cm3)、溶解度参数(18.604~19.237 (J cm−3)0.5)和玻璃化转变温度(237.97~247.61 K)三个关键指标验证了模型的合理性,保证了与已发表的实验和仿真数据的一致性。通过评价界面能(相互作用能、范德华能、静电能),分析了官能团类型和浓度对沥青-骨料黏附性能的影响。随着羰基(C=O)基团浓度的增加,沥青-骨料体系的相互作用能、范德华能和静电能均呈上升趋势。100%浓度下的增幅大于50%浓度下的增幅,其中静电能的整体增幅最大。随着亚砜基(S=O)含量的增加,沥青-骨料体系的相互作用能和范德华能降低,而静电能增加。当羰基(C=O)和亚砜基(S=O)共存时,各种界面能的变化模式与亚砜单独影响下的变化模式相似。这说明在相同的条件下,亚砜基(S=O)对界面能的影响更大。氧化官能团对沥青与集料的粘附性能有显著影响。方法为了研究不同浓度和类型的含氧官能团对沥青-骨料黏附行为的影响,首先利用Materials Studio 2023软件中的Amorphous Cell模块,通过分子组装法构建沥青分子模型,参考沥青四组分的TLC-FID数据,确定各组分的分子量。共建立了15个沥青模型,涵盖了3种官能团类型(C=O、S=O、C=O + S=O)和3种浓度梯度。通过沿(0 ~ 0 ~ 1)方向切割SiO2晶体,插入10 Å真空层形成周期晶胞,并进行氢原子钝化,模拟真实团聚体表面性质,建立SiO2团聚体模型。随后,通过Forcite模块对沥青-骨料模型进行分子动力学模拟,其中采用COMPASS III力场来描述分子间键合和非键合相互作用。模拟过程包括三个关键步骤:(1)几何优化(中等精度,50,000次迭代)以消除不合理的原子构型;(2)在NVT集成(由Nose恒温器控制)下进行5次退火循环(298~498 K,共200 ps)以稳定系统;(3)在298 K和1.0 × 10−4 GPa(由Berendsen气压调节器控制)下进行500 ps NVT + 1000 ps NPT的平衡模拟。
{"title":"Molecular dynamics simulation of the influence of oxidized functional groups on asphalt-aggregate bonding properties","authors":"Maoping Ran,&nbsp;Jiang Zuo,&nbsp;Xinglin Zhou,&nbsp;Jing Yuan,&nbsp;Xinxing Zhou","doi":"10.1007/s00894-026-06644-4","DOIUrl":"10.1007/s00894-026-06644-4","url":null,"abstract":"&lt;div&gt;&lt;h3&gt;Context&lt;/h3&gt;&lt;p&gt;To elucidate the mechanism by which oxidized functional groups influence asphalt-aggregate adhesion behavior at the microscopic level, molecular dynamics (MD) simulations were employed based on the Derek four-component model—a model well-matched with real asphalt’s density and thermal expansion coefficient. Asphalt molecular models containing 12 characteristic molecules were constructed, with two stable oxidized functional groups (carbonyl, C=O; sulfoxyl, S=O) introduced to represent asphalt aging. Three concentration gradients (0%, 50%, 100%) were designed by adjusting the number of functional groups. Meanwhile, silica (SiO&lt;sub&gt;2&lt;/sub&gt;) aggregate models were cleaved along the (0 0 1) crystal plane, transformed into an orthorhombic structure, and passivated with hydrogen atoms to simulate real aggregate surface properties. An asphalt-aggregate interface system with a “SiO&lt;sub&gt;2&lt;/sub&gt; layer-10 Å vacuum layer-asphalt layer-50 Å vacuum layer” sandwich structure was established. The model rationality was verified by three key indicators: density (1.011~1.079 g/cm&lt;sup&gt;3&lt;/sup&gt;), solubility parameter (18.604~19.237 (J cm&lt;sup&gt;−3&lt;/sup&gt;)&lt;sup&gt;0.5&lt;/sup&gt;), and glass transition temperature (237.97~247.61 K), ensuring consistency with published experimental and simulation data. The effects of functional group type and concentration on asphalt-aggregate adhesion were analyzed by evaluating interfacial energies (interaction energy, van der Waals energy, electrostatic energy). As the concentration of carbonyl (C=O) groups increased, the interaction energy, van der Waals energy, and electrostatic energy in the asphalt-aggregate system all showed an upward trend. The increase at 100% concentration was greater than that at 50% concentration, with electrostatic energy exhibiting the largest overall increase. As the sulfoxyl group (S=O) content increased, the interaction energy and van der Waals energy in the asphalt-aggregate system decreased, while the electrostatic energy increased. When both carbonyl (C=O) and sulfoxyl (S=O) groups coexist, the variation patterns of various interfacial energies resemble those observed under sulfoxide influence alone. This indicates that, under identical conditions, sulfoxyl (S=O) exerts a greater impact on interfacial energies. Oxidized functional groups significantly influence the adhesion properties between asphalt and aggregates.&lt;/p&gt;&lt;h3&gt;Methods&lt;/h3&gt;&lt;p&gt;To investigate the effects of oxygen-containing functional groups with different concentrations and types on the asphalt-aggregate adhesion behavior, asphalt molecular models were first constructed via the molecular assembly method using the Amorphous Cell module in Materials Studio 2023 software, referencing TLC-FID data of asphalt four components to determine the molecular quantity of each component. A total of 15 asphalt models were established, covering three functional group types (C=O, S=O, C=O + S=O) and three concentration gradients. SiO&lt;sub&gt;2&lt;/sub&gt; ","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 3","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-02-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146099217","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}
引用次数: 0
A DFT-D2 study on adsorption of iodonitromethane on doped (B, N & Fe) and Fe-functionalized monovacancy graphene surfaces 掺杂(B, N, Fe)和Fe功能化单空位石墨烯表面吸附碘硝基甲烷的DFT-D2研究
IF 2.5 4区 化学 Q4 BIOCHEMISTRY & MOLECULAR BIOLOGY Pub Date : 2026-01-31 DOI: 10.1007/s00894-025-06614-2
Rajesh Kumar Chatra, Deepak Tomar, K. Dhana Sai Shree,  Jogender

Context

INM is a toxic halogen nitromethane and disinfection by-product that is very dangerous to the environment and human beings since it is highly cytotoxic and mutagenic. Here, we examine the adsorption of INM on pristine graphene (PG), monovacancy graphene (MVG), nitrogen-doped vacuity graphene (NVG), and boron-doped vacuity graphene (BVG) and their Fe-functionalized analogs (FeG, FNG, and FBG) using density functional theory (DFT). To explain the adsorption mechanism, adsorption energy calculations, Hirshfeld charge transfer analysis, and electronic structure evaluations, such as band gap, density of states (DOS), and partial DOS (PDOS), were used. INM shows poor physisorption on PG, MVG, and NVG, whereas BVG shows stronger chemisorption by direct bonding. Fe adsorption is very important in increasing the strength of adsorption and redistribution of charges which results in strong electronic structure alterations. Then Fe-doped vacancy graphene is explored with respect to band gap, DOS, and PDOS plots. INM is then adsorbed on the surface. Fe-doped vacancy graphene (FVG) exhibits the highest adsorption energy and the largest electronic modification which validates the fact that it is strongly chemisorptively interacting. These findings emphasize FVG as a promising and efficient material to remove toxic INM in the contaminated environment and to develop graphene-based adsorbents to clean up the environment.

Methods

All calculations were done in the DMol3 package of Materials Studio with a doubled numerical plus polarization (DNP) basis set and DFT semicore pseudopotentials (DSPP). The GGA-PBE functional was used with the DFT-D correction by Grimme which considered the effects of exchange-correlation and dispersion. Spin-unrestricted geometry optimizations were driven to strict energy (2.0 × 10−5 Ha) and force (4 × 10−3 Ha Å−1) convergence factors, on a grid of 6 × 6 × 1 Monkhorst Pack k-points. To assess charge transfer and bonding properties, Hirshfeld charge and Mayer bond order analyses were used.

inm是一种有毒的卤素硝基甲烷和消毒副产物,对环境和人类都非常危险,因为它具有高度的细胞毒性和诱变性。本研究利用密度泛函理论(DFT)研究了INM在原始石墨烯(PG)、单空位石墨烯(MVG)、氮掺杂真空石墨烯(NVG)、硼掺杂真空石墨烯(BVG)及其铁功能化类似物(FeG、FNG和FBG)上的吸附。为了解释吸附机理,使用了吸附能计算、Hirshfeld电荷转移分析和电子结构评价,如带隙、态密度(DOS)和部分态密度(PDOS)。INM在PG、MVG和NVG上表现出较差的物理吸附,而BVG通过直接键合表现出较强的化学吸附。铁的吸附对提高吸附强度和电荷的重新分配非常重要,从而导致强烈的电子结构改变。然后探讨了掺铁空位石墨烯的带隙、DOS和PDOS图。然后INM被吸附在表面。fe掺杂的空位石墨烯(FVG)表现出最高的吸附能和最大的电子修饰,这证实了它具有强烈的化学吸附相互作用。这些发现强调了FVG是一种有前途的高效材料,可以去除污染环境中的有毒INM,并开发基于石墨烯的吸附剂来清洁环境。方法采用双数值+极化(DNP)基集和DFT半核赝势(DSPP),在Materials Studio的DMol3软件包中进行计算。GGA-PBE泛函采用了考虑交换相关性和色散影响的grime DFT-D校正。在6 × 6 × 1个Monkhorst Pack k点的网格上,不受自旋限制的几何优化被驱动到严格的能量(2.0 × 10−5 Ha)和力(4 × 10−3 Ha Å−1)收敛因子。为了评估电荷转移和成键性质,使用了Hirshfeld电荷和Mayer键序分析。
{"title":"A DFT-D2 study on adsorption of iodonitromethane on doped (B, N & Fe) and Fe-functionalized monovacancy graphene surfaces","authors":"Rajesh Kumar Chatra,&nbsp;Deepak Tomar,&nbsp;K. Dhana Sai Shree,&nbsp; Jogender","doi":"10.1007/s00894-025-06614-2","DOIUrl":"10.1007/s00894-025-06614-2","url":null,"abstract":"<div><h3>Context</h3><p>INM is a toxic halogen nitromethane and disinfection by-product that is very dangerous to the environment and human beings since it is highly cytotoxic and mutagenic. Here, we examine the adsorption of INM on pristine graphene (PG), monovacancy graphene (MVG), nitrogen-doped vacuity graphene (NVG), and boron-doped vacuity graphene (BVG) and their Fe-functionalized analogs (FeG, FNG, and FBG) using density functional theory (DFT). To explain the adsorption mechanism, adsorption energy calculations, Hirshfeld charge transfer analysis, and electronic structure evaluations, such as band gap, density of states (DOS), and partial DOS (PDOS), were used. INM shows poor physisorption on PG, MVG, and NVG, whereas BVG shows stronger chemisorption by direct bonding. Fe adsorption is very important in increasing the strength of adsorption and redistribution of charges which results in strong electronic structure alterations. Then Fe-doped vacancy graphene is explored with respect to band gap, DOS, and PDOS plots. INM is then adsorbed on the surface. Fe-doped vacancy graphene (FVG) exhibits the highest adsorption energy and the largest electronic modification which validates the fact that it is strongly chemisorptively interacting. These findings emphasize FVG as a promising and efficient material to remove toxic INM in the contaminated environment and to develop graphene-based adsorbents to clean up the environment.</p><h3>Methods</h3><p>All calculations were done in the DMol<sup>3</sup> package of Materials Studio with a doubled numerical plus polarization (DNP) basis set and DFT semicore pseudopotentials (DSPP). The GGA-PBE functional was used with the DFT-D correction by Grimme which considered the effects of exchange-correlation and dispersion. Spin-unrestricted geometry optimizations were driven to strict energy (2.0 × 10<sup>−5</sup> Ha) and force (4 × 10<sup>−3</sup> Ha Å<sup>−1</sup>) convergence factors, on a grid of 6 × 6 × 1 Monkhorst Pack k-points. To assess charge transfer and bonding properties, Hirshfeld charge and Mayer bond order analyses were used.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"32 2","pages":""},"PeriodicalIF":2.5,"publicationDate":"2026-01-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146083068","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}
引用次数: 0
期刊
Journal of Molecular Modeling
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:604180095
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1