Journal of Molecular Modeling最新文献

{"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":"https://doi.org/10.1007/s00894-025-06531-4","url":null,"abstract":"<p><strong>Context: </strong>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><p><strong>Methods: </strong>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>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":"307"},"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}

{"title":"DFT and TDDFT studies on π-conjugated ligands for copper sensing: analyzing electronic structures and optical performance.","authors":"Susheel, Abdullah Saad Alsubaie, Mukhtar Ahmed, Sumit Sahil Malhotra, Asha, Manoj Kumar Gupta, Azaj Ansari","doi":"10.1007/s00894-025-06543-0","DOIUrl":"https://doi.org/10.1007/s00894-025-06543-0","url":null,"abstract":"<p><strong>Context: </strong>Structural, bonding aspects and optical characteristics of a set of π-conjugated quinoline-based ligands (L1-L4) and coordinated copper (Cu²⁺) ion were investigated by using density functional theory (DFT) and time dependent DFT methods. DFT results showed that L3 exhibits the lowest HOMO-LUMO energy gap (3.05 eV) indicating high reactivity and strong charge transfer ability while its copper complex further reduces the gap to 2.52 eV. Electrostatic potential maps highlighted a negative potential region around nitrogen and carbonyl oxygen sites confirming their role in copper coordination. Natural bond orbital analysis of the L3 complex revealed the highest stabilization energy of 79.15 kcal/mol indicating substantial donor-acceptor interactions. ELF and LOL plots further supported efficient π-delocalization in L3-Cu²⁺ while NCI analysis further confirmed reduced steric repulsion around the Cu²⁺ coordination sphere compared to other complexes which support its favorable geometry and stability. QTAIM analysis indicated a mixed electrostatic covalent character of Cu-N/O bonds. TDDFT results showed strong ligand to metal charge transfer bands in the visible spectrum for L3-Cu²⁺ at 452 and 667 nm which lend credence to a mechanism of chelation-enhanced charge transfer. Non-linear optical analysis revealed enhanced first hyperpolarizability upon complexation particularly for L1-Cu²⁺ (β = 9.46 × 10^-30 esu) and L4-Cu²⁺ (β = 9.12 × 10^-30 esu). These observations provide a useful layout for generating metal ion sensors in the future with improved optical response and selectivity. These theoretical findings agree with the coordination behavior seen in experiments supporting L3-based systems in the copper ion detection applications.</p><p><strong>Methods: </strong>Geometry optimization and frequency analyses were performed using DFT at the B3LYP/6-311G(d,p) level for non-metal atoms and LANL2DZ basis set for copper. The polarized continuum model was used for the solvation as implemented in Gaussian 16. The NBO6.0 program was utilized to investigate the bonding nature and stabilization energies of the complexes. The ORCA4.2 program was used to simulate the absorption spectrum. The Multiwfn and VMD programs were used for the topological analysis.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":"306"},"PeriodicalIF":2.5,"publicationDate":"2025-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145342745","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}

{"title":"From nonpolar to polar: the reconstruction mechanism of electronic structure in alkane C–C bonds induced by fluorine substitution","authors":"Chuan Wang,&nbsp;Ziqiu Wang,&nbsp;Kun Yang,&nbsp;Guangxi Zhang,&nbsp;Nayang Li","doi":"10.1007/s00894-025-06542-1","DOIUrl":"10.1007/s00894-025-06542-1","url":null,"abstract":"<div><h3>Context</h3><p>The C–C bonds in alkanes are generally considered nonpolar covalent bonds, but the mechanism by which fluorine substitution affects the polarity and cleavage mode of C–C bonds in alkanes remains unclear. In this study, using 1,1,1-trifluoroethane (CH₃-CF₃) as a model, we systematically investigated the reconstruction mechanism of the electronic structure of C–C bonds induced by fluorine substitution through density functional theory calculations, combined with bond critical point (BCP) topological analysis, atomic charge calculation, and flexible scan simulation. The results demonstrate that fluorine substitution transforms the C–C bond into a polar covalent bond via a strong electron-withdrawing inductive effect, leading to the separation of positive and negative charges on the two carbon atoms and breaking the nonpolar symmetric distribution. Flexible scan simulations reveal that the C–C bond cleavage in CH₃-CH₃ follows a typical covalent bond homolysis, whereas that in CH₃-CF₃ exhibits ionic bond characteristics. This study uncovers the reconstruction mechanism of C–C bonds from nonpolar to polar induced by fluorine substitution and its disruptive impact on the bond cleavage pathway, providing a theoretical basis for understanding the structure–property relationships of complex fluoroalkanes.</p><h3>Methods</h3><p>All quantum chemical calculations adopted the M06-2X functional paired with the Def2-TZVP basis set. Initial molecular structures were built in GaussView, with geometry optimization performed via Gaussian 16 to obtain optimized configurations and wavefunction files. Wavefunction analyses (including BCP topology, Hirshfeld atomic charge, Mayer bond order, and spin population/density calculations) were conducted using Multiwfn.</p><h3>Graphic Abstract</h3>\u0000<div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145306831","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}

{"title":"Influence of auxiliary acceptor substitution at D-A1-π-A2 structured highly efficient organic molecules for dye-sensitized solar cells using computational study","authors":"A. Arunkumar","doi":"10.1007/s00894-025-06529-y","DOIUrl":"10.1007/s00894-025-06529-y","url":null,"abstract":"<div><h3>Context</h3><p>This study used the PTRA-based molecule to design and investigate seven new donor-acceptor1-π-acceptor2 (D-A1-π-A2) organic dyes (PTRA1-PTRA7) for dye-sensitized solar cells (DSSCs). All dyes had an electron-A2 group made of rhodanine-3-acetic acid and a thiophene group as a spacer, while the electron-A1 unit varied from D-A1-π-A2. Optoelectronics was explored in relation to a structure and its influences. For the dyes PTRA1–PTRA7, the computational analysis of density functional theory (DFT) and its extended time-dependent DFT (TD-DFT) approaches was performed. The driving force of electron injection (<span>(Delta G_{inject})</span>), dye regeneration (<span>(Delta G_{reg})</span>), exciton binding energy (<i>E</i>_b), molecular orbitals (MOs) energy levels, optical ultraviolet-visible (UV-Vis) spectra, and electronic properties were all thoroughly discussed. The findings show that PTRA1–PTRA7 have smaller energy gaps (<i>E</i>_g) and higher absorption wavelength (<i>λ</i>_max) than PTRA. Out of all of them, PTRA6 has the lowest <i>E</i>_g (2.15 eV) and the red-shifted <i>λ</i>_max (559 nm). It has been determined that dyes PTRA1–PTRA7 are the most promising option for having highly efficient DSSCs, especially PTRA6. Additionally, these molecules serve as the most promising functional group in D-A1-π-A2 dyes for auxiliary electron-A. Due to its excellent electronic, optical, and photovoltaic (PV) characteristics, it could be utilized as a potential sensitizer for DSSCs.</p><h3>Methods</h3><p>Formalisms of the conductor-like polarizable continuum model (CPCM) have been used to study solvent effects. The results of the CPCM/TD-DFT demonstrate that precise absorption energies can only be obtained when the solvent effect is taken into account in the geometries of the excited states. The Gaussian 09w software package is used for related calculations.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145306791","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}

{"title":"Interface properties and nucleation thermodynamics of Al2Cu precipitated phase by first-principles calculations","authors":"Yan Huang,&nbsp;Jiyi Li,&nbsp;Cuifan Chen,&nbsp;Liang Zhang","doi":"10.1007/s00894-025-06534-1","DOIUrl":"10.1007/s00894-025-06534-1","url":null,"abstract":"<div><h3>Context</h3><p>The Al₂Cu precipitated phase plays a critical role in governing the stability and mechanical properties of age-hardened Al–Cu alloys, which are widely used in the aerospace and automotive industries due to their high strength-to-weight ratio. However, the influence of alloying elements (e.g., Fe, Mn, Mg, Sc, and Zr) on the interfacial stability and nucleation behavior of Al₂Cu remains insufficiently understood, limiting the rational design of advanced Al–Cu alloys with optimized performance. This study addresses this gap by investigating how elemental segregation affects the structural stability and nucleation thermodynamics of Al₂Cu interfaces, offering insights into strategies for enhancing the mechanical properties and thermal stability of these alloys.</p><h3>Methods</h3><p>First-principles calculations based on density functional theory (DFT) were employed to evaluate the effects of segregated elements (Fe, Mn, Mg, Sc, and Zr) on the coherent strain energy and interface energy of Al/Al₂Cu interfaces. Electronic structure calculations were performed using the Vienna ab-initio simulation package (VASP), employing the Perdew–Burke–Ernzerhof (PBE) functional within the generalized gradient approximation (GGA) to account for exchange–correlation effects. To further elucidate the bonding mechanisms, interface stability was analyzed through detailed electronic structure investigations. Classical nucleation theory was applied to predict phase nucleation thermodynamics at aging temperatures, examining the formation of core–shell precipitates (Al₂Cu cores with solute-rich shells) as a function of precipitate size.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145298025","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}

{"title":"Computational mechanistic study on N-nitrosation reaction of secondary amines","authors":"Meire Y. Kawamura,&nbsp;David J. Ponting,&nbsp;Chris G. Barber,&nbsp;Michael J. Burns","doi":"10.1007/s00894-025-06520-7","DOIUrl":"10.1007/s00894-025-06520-7","url":null,"abstract":"<div><h3>Context</h3><p>The presence of potentially carcinogenic nitrosamines in drugs has been a worldwide concern, driving strategies to control or mitigate their formation to protect patient health. Understanding the critical factors for <i>N</i>-nitrosation, such as mechanisms and energy barriers, enhances the risk assessment process to understand potential nitrosamine formation. Evaluation of the structural impact of amines on the <i>N</i>-nitrosation rate in the presence of nitrites and acidic media is of great interest to pharmaceutical companies assessing the risk of nitrosamine drug substance–related impurities. A range of secondary amines was explored using DFT calculations to assess the impact of electronic and steric effects on activation energy. <i>Asym</i>-N₂O₃ was selected as the nitrosating agent since its reaction was shown to be favorable following screening of pathways employing nitrosyl chloride, nitrous acid, <i>asym</i>-N₂O₃, <i>sym</i>-N₂O₃, and <i>trans-cis</i>-N₂O₃. The relatively low activation energies obtained for all amines indicate the reaction is very likely to occur if the reactive components encounter, even for amines with sterically hindered and electron-withdrawing groups. Understanding the interaction between the amine and nitrosating agent is therefore the defining factor in the risk of formation of more complex nitrosamines within drugs.</p><h3>Methods</h3><p>Calculations were performed using the Gaussian-16 program. The B3LYP-D3/def2-TZVP level of theory was applied for structure optimizations. The IEF-PCM implicit model was used for the solvent effect. Intrinsic reaction coordinate calculations were carried out to connect the transition state with the associated minimum.\u0000</p><h3>Graphical Abstract</h3>\u0000<div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145297978","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}

{"title":"DFT-based evaluation of covalent organic frameworks for adsorption, optoelectronic, clean energy storage, and gas sensor applications","authors":"Abhay P. Srivastava,&nbsp;Brijesh K. Pandey","doi":"10.1007/s00894-025-06535-0","DOIUrl":"10.1007/s00894-025-06535-0","url":null,"abstract":"<div><h3>Context</h3><p>Covalent Organic Frameworks (COFs), which are frameworks composed of light atoms held together by strong covalent bonds, are generating interest as potential materials for applications such as renewable energy and gas capture. We employed Density Functional Theory (DFT) calculations, as implemented in the VASP code, to look at both 2D and 3D COFs. We systematically analyzed various properties, including structural stability, phonon dispersion, electronic structures, density of states, adsorption behavior, and mechanical properties. To get better accuracy, we took into account van der Waals interactions and even used hybrid functionals. What we found was that 3D COFs generally exhibit greater mechanical strength and, in most cases, better gas adsorption, which seems to come from their interconnected pore structures. On the other hand, 2D COFs exhibit enhanced π-electron delocalization and direct band gaps of approximately 2.5 eV, which may be helpful in sensors and optoelectronics. Phonon analyses verified the dynamical stability of both structures. Ultimately, these results underscore the importance of dimensionality in tailoring COF properties for energy and electronic applications.</p><h3>Method</h3><p>First-principles simulations were performed using Density Functional Theory (DFT) within the Vienna Ab initio Simulation Package (VASP). To account for exchange–correlation effects, we employed the Generalised Gradient Approximation (GGA) in the Perdew–Burke–Ernzerhof (PBE) formulation, and we also utilised projector-augmented wave (PAW) pseudopotentials. Hybrid functional (HSE06) and DFT-D3 van der Waals corrections were introduced to improve the accuracy of our band gap prediction. The plane-wave cutoffs were set at 500 eV for the calculations, and Monkhorst–Pack k-point meshes were used with a 3 × 3 × 1 (2D) and 2 × 2 × 2 (3D) grid. Evaluated were structural optimisations, band structures, total and projected DOS, adsorption energies, and charge transfer (using Bader analysis). Assessment of bonding features utilised the Electron Localisation Function (ELF) and charge density difference (Δρ) visualisations. The Phonopy package was used to calculate Phonon dispersions and thus confirm the dynamic stability of the COFs.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145298043","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}

{"title":"Comparative quantum chemical analysis of dexamethasone and hydrocortisone: electronic structure, and reactivity indices using DFT","authors":"Masoumeh Eskandari-Nasab,&nbsp;Zainab Moosavi-Tekyeh,&nbsp;Mansoureh Zahedi-Tabrizi","doi":"10.1007/s00894-025-06522-5","DOIUrl":"10.1007/s00894-025-06522-5","url":null,"abstract":"<div><h3>Context</h3><p>The biological activity of steroidal compounds such as dexamethasone (DX) and hydrocortisone (HC) is closely linked to subtle variations in their molecular structure and electronic properties. This study provides a comparative quantum chemical analysis of DX and HC to clarify how these differences influence hydrogen bonding strength, reactivity, and their potential interactions with the glucocorticoid receptor (GR). Optimized geometries, natural bond orbital (NBO) analyses, frontier molecular orbitals (FMO), global reactivity descriptors, and average local ionization energy (ALIE) calculations demonstrate that DX exhibits greater polarity and electrophilic character compared to HC. These differences help explain the stronger receptor binding affinity observed for DX. Indeed, notably, despite the inherent limitations of gas-phase DFT calculations compared to experimental X-ray data, the theoretical results exhibit good agreement with experimental observations, suggesting the reliability of the computational approach in predicting molecular interactions within the GR active site.</p><h3>Methods</h3><p>All quantum chemical calculations were performed using density functional theory (DFT) with the B3LYP functional and 6-311++G(d,p) basis set. Structural optimization, FMO analysis, global reactivity descriptors, and dipole moment evaluations were carried out in Gaussian 09. NBO analysis was performed with NBO 5.0. Average local ionization energy (ALIE) surfaces were generated using Multiwfn 3.8, and molecular visualizations were produced with GaussView 5.0.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145298005","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}

{"title":"Antioxidant and chelating properties of phenolic compounds of agro-industrial waste of Carapa guianensis: theoretical insights for food and pharmaceutical applications","authors":"Neidy S. S. dos Santos,&nbsp;Loubenky Surfin,&nbsp;William Silva,&nbsp;Antonio Rodrigues da Cunha,&nbsp;Herbert Georg,&nbsp;Kaline Coutinho,&nbsp;Sylvio Canuto,&nbsp;Rodrigo Gester","doi":"10.1007/s00894-025-06517-2","DOIUrl":"10.1007/s00894-025-06517-2","url":null,"abstract":"<div><h3>Context</h3><p>This study investigates the antioxidant behavior and chelating properties of thirteen phenolic compounds identified in <i>Carapa guianensis</i> agro-industrial waste. Using quantum mechanical methods, we analyzed the molecular structure and solvent effects, revealing potential applications in the food and pharmaceutical industries. Approximately 85% of these compounds demonstrated superior antioxidant performance compared to the phenol backbone, with three compounds rivaling quercetin and ascorbic acid in all tested environments. Solvent polarity significantly influenced the antioxidant mechanism: while hydrogen atom transfer (HAT) dominated in the gas phase, sequential proton loss electron transfer (SPLET) became prevalent in polar solvents. Hydrogen abstraction occurred primarily at the <i>meta</i> position, though polar solvents increased activity at the <i>para</i> site. These predictions are confirmed by simulating the chemical reactions between the aromatic compounds and free OH radicals. The analysis of the hydrogen abstraction reactions indicates that the inclusion of Hartree-Fock exchange-correlation and dispersion corrections is essential to describe hydrogen abstraction by explicitly free radicals. These findings not only underscore the commercial potential of <i>Carapa guianensis</i> waste but also provide a comprehensive understanding of its antioxidant mechanisms, contributing to the sustainable use of natural resources for health and environmental benefits.</p><h3>Methods</h3><p>This is a two-step methodology. First, we select the best antioxidant molecules through thermochemical analysis. Then, we explore the hydrogen scavenging mechanism using transition state theory. All quantum mechanics and thermodynamic analyses were carried out within the framework of the density functional theory using the M06-2X, <span>(omega )</span>B97, <span>(omega )</span>B97XD, and CAM-B3LYP variants associated with the Pople’s 6-311++G(<i>d</i>, <i>p</i>) basis set. The effects of the solvent are also systematically investigated by considering the Solvent Density Model for different polar and nonpolar environments. Naturally, the unrestricted wave-function formalism is accounted for once the thermodynamic description of the antioxidant mechanisms involves optimizing open-shell structures. Additionally, the chemical reaction of the hydrogen abstraction due to OH radicals is ensured by applying the Synchronous Transit-Guided Quasi-Newton Method, which allows us to estimate the transition state structures and the energy barrier of the reaction. At this stage, DFT methods like CAM-B3LYP and <span>(omega )</span>B97XD were applied in association with the 6-31+G(<i>d</i>). All calculations were carried out taking advantage of the Gaussian 16 program.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145298050","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}

{"title":"Machine learning-driven prediction of organic solar cell performance: a data-centric approach to molecular design","authors":"Victor dos Reis Rodrigues,&nbsp;Víctor de Souza Assumção Bonfim,&nbsp;Demétrio Antônio da Silva Filho","doi":"10.1007/s00894-025-06514-5","DOIUrl":"10.1007/s00894-025-06514-5","url":null,"abstract":"<div><h3>Context</h3><p>Organic solar cells (OSCs) offer a promising route toward flexible and sustainable energy technologies, yet predictive modeling of device parameters remains challenging due to the chemical diversity of donor–acceptor systems and morphology-dependent effects. In this work, we present the first systematic demonstration of using autoencoder-compressed molecular fingerprints with tree-based machine learning models to predict key OSC performance metrics—power conversion efficiency (PCE), open-circuit voltage (V_oc), short-circuit current (J_sc), and fill factor (FF)—from a broad experimental dataset of  2500 donor–acceptor pairs, including both fullerene and non-fullerene acceptors. These compact models, trained on compressed descriptors of only 32 dimensions, achieved strong predictive accuracy (Pearson <span>(r &gt; 0.95)</span>, <span>(MAE &lt; 0.4)</span>, <span>(RMSE &lt; 0.95)</span>) while remaining lightweight enough to run on standard computing hardware. As a complementary result, some <i>k</i>-nearest neighbor models achieved near-perfect correlations (<span>(r sim 0.99)</span>) and quite small errors (<span>(MAE &lt; 0.044)</span> and <span>(RMSE&lt;0.4)</span>) in general, demonstrating the surprising strength of simple, instance-based learners when sufficient descriptive features are available. Supporting analyses reveal that fullerene datasets are more easily modeled than chemically diverse non-fullerene sets, that fingerprints encode substantial structural information, and that kernel density analyses identify critical ranges of molecular weight and energy offsets for high-efficiency devices. Collectively, this study establishes compressed fingerprint descriptors as a powerful, computationally inexpensive foundation for predictive modeling in OSCs, while also showcasing the unexpected efficacy of k-NN models trained on conventional descriptors. Together, these approaches provide a scalable path toward high-throughput prediction and guided molecular design of next-generation organic photovoltaic materials.</p><h3>Methods</h3><p>The dataset used in this work comprises approximately 2500 experimentally characterized donor–acceptor pairs from bulk heterojunction OSCs. These include both fullerene and non-fullerene acceptor systems. For each pair, the database provides electronic descriptors, polymerization-related metrics, and the SMILES representations of the donor and acceptor molecules. Molecular fingerprints were computed from SMILES codes using the RDKit and CDK cheminformatics toolkits. A variety of machine learning models were explored, including feedforward neural networks, autoencoders for feature compression, tree-based ensemble methods, and kernel-based regression algorithms. Hyperparameter tuning was carried out using the Optuna and BayesSearchCV libraries to ensure optimal model performance.</p></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 11","pages":""},"PeriodicalIF":2.5,"publicationDate":"2025-10-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145290539","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}