Journal of Molecular Modeling最新文献

{"title":"Theoretical study of potential energetic material CL-20/DNAN eutectic explosive based on molecular dynamics method.","authors":"Jihang Du, Baoguo Wang, Yafang Chen, Xinyi Li, Chunguang Wang","doi":"10.1007/s00894-024-06109-6","DOIUrl":"https://doi.org/10.1007/s00894-024-06109-6","url":null,"abstract":"<p><strong>Context: </strong>The exploration of CL-20 eutectic has been a subject of fervent interest within the realm of high-energy material modification. Through the utilization of density functional and molecular dynamics methods, an investigation into the characteristics of hexanitrohexaazaisowurtzitane (CL-20)/2,4-dinitroanisole (DNAN) within the molar ratio range of 9:1-1:9 was conducted. This inquiry encompassed the scrutiny of molecular interaction pathway, attachment force, initiating molecular distance, unified energy concentration, and physical characteristics. Furthermore, EXPLO-5 was harnessed to prognosticate the explosion features and byproducts of unadulterated CL-20, DNAN, and CL-20/DNAN frameworks. The findings delineate a substantial differentiation in the electrostatic charge distribution on the surface between CL-20 and DNAN particles, signifying the preeminence of intermolecular interactions between disparate entities over those within similar entities, thus intimating the plausibility of eutectic constitution. Remarkably, the identification of maximal attachment force at a molar ratio of 4:6 suggests the heightened likelihood of eutectic formation, propelled primarily by electrostatic and van der Waals forces. The resultant eutectic explosive evinces intermediate reactivity and exemplary mechanical attributes. Moreover, the detonation achievement of the eutectic with a molar proportion of 4:6 straddles that of CL-20 and DNAN, representing a new type of insensitive high-energy material.</p><p><strong>Methods: </strong>The testing method employs the Materials Studio software and utilizes the molecular dynamics (MD) method to predict the properties of CL-20/DNAN co-crystals with different ratios and crystal faces. The MD simulation time step is set to 1 fs, and the total MD simulation time is 2 ns. An isothermal-isobaric (NPT) ensemble is used for the 2-ns MD simulation. The COMPASS force field is employed, with the temperature set to 295 K. The prediction of detonation characteristics and products is conducted using the EXPLO-5 software.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141999208","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":"<ArticleTitle xmlns:ns0=\"http://www.w3.org/1998/Math/MathML\">New global minimum conformers for the Pt <ns0:math><ns0:mmultiscripts><ns0:mrow /> <ns0:mn>19</ns0:mn> <ns0:mrow /></ns0:mmultiscripts> </ns0:math> and Pt <ns0:math><ns0:mmultiscripts><ns0:mrow /> <ns0:mn>20</ns0:mn> <ns0:mrow /></ns0:mmultiscripts> </ns0:math> clusters: low symmetric species featuring different active sites.","authors":"José Manuel Guevara-Vela, Miguel Gallegos, Tomás Rocha-Rinza, Álvaro Muñoz-Castro, Peter L Rodríguez Kessler, Ángel Martín Pendás","doi":"10.1007/s00894-024-06099-5","DOIUrl":"10.1007/s00894-024-06099-5","url":null,"abstract":"<p><strong>Context: </strong>The study of platinum (Pt) clusters and nanoparticles is essential due to their extensive range of potential technological applications, particularly in catalysis. The electronic properties that yield optimal catalytic performance at the nanoscale are significantly influenced by the size and structure of Pt clusters. This research aimed to identify the lowest-energy conformers for Pt <math><mmultiscripts><mrow></mrow> <mn>18</mn> <mrow></mrow></mmultiscripts> </math> , Pt <math><mmultiscripts><mrow></mrow> <mn>19</mn> <mrow></mrow></mmultiscripts> </math> , and Pt <math><mmultiscripts><mrow></mrow> <mn>20</mn> <mrow></mrow></mmultiscripts> </math> species using Density Functional Theory (DFT). We discovered new low-symmetry conformers for Pt <math><mmultiscripts><mrow></mrow> <mn>19</mn> <mrow></mrow></mmultiscripts> </math> and Pt <math><mmultiscripts><mrow></mrow> <mn>20</mn> <mrow></mrow></mmultiscripts> </math> , which are 3.0 and 1.0 kcal/mol more stable, respectively, than previously reported structures. Our study highlights the importance of using density functional approximations that incorporate moderate levels of exact Hartree-Fock exchange, alongside basis sets of at least quadruple-zeta quality. The resulting structures are asymmetric with varying active sites, as evidenced by sigma hole analysis on the electrostatic potential surface. This suggests a potential correlation between electronic structure and catalytic properties, warranting further investigation.</p><p><strong>Methods: </strong>An equivariant graph neural network interatomic potential (NequIP) within the Atomic Simulation Environment suite (ASE) was used to provide initial geometries of the aggregates under study. DFT calculations were performed with the ORCA 5 package, using functional approximations that included Generalized Gradient Approximation (PBE), meta-GGA (TPSS, M06-L), hybrid (PBE0, PBEh), meta-GGA hybrid (TPSSh), and range-separated hybrid ( <math><mi>ω</mi></math> B97x) functionals. Def2-TZVP and Def2-QZVP as well as members of the cc-pwCVXZ-PP family to check basis set convergence were used. QTAIM calculations were performed using the AIMAll suite. Structures were visualized with the AVOGADRO code.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11330413/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141995072","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}

{"title":"Charge transfer interaction revisited by a Fermi-Dirac derived approach.","authors":"Oscar Gayraud, Frédéric Guégan","doi":"10.1007/s00894-024-06097-7","DOIUrl":"https://doi.org/10.1007/s00894-024-06097-7","url":null,"abstract":"<p><strong>Context: </strong>In this article, we adapt a recent proposition to use a Fermi-Dirac-type population scheme on Kohn-Sham molecular orbitals to the case of an interaction with a thermalised electrode. This allows to derive a fundamental non-linear equation linking the chemical potential of the electrode and the amount of charge transferred to the system under study, hence allows to quantify the propensity to charge transfer (philicity). This methodology is applied to a large set of common electrophiles and nucleophiles, showing decent relation with more standard philicity descriptors. Chemical hardness is also revisited by this approach.</p><p><strong>Methods: </strong>All calculations were performed using the Gaussian 16 software package at the M062X/aug-cc-pvtz level of theory. Data analysis was then performed through a Python3 dedicated program (relying on the fsolve numerical solver from the SciPy package), using Gaussian output files, and available as supplementary material.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141970381","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":"Enhanced thermoelectric performance of Hf-doped ZrNiSn: a first principle study.","authors":"Di Cao, Jiannong Cao","doi":"10.1007/s00894-024-06102-z","DOIUrl":"https://doi.org/10.1007/s00894-024-06102-z","url":null,"abstract":"<p><strong>Context and results: </strong>In this work, we perform a systematic study on the thermoelectric properties of Zr_1-xNiSnHf_x using first-principles calculations combined with Boltzmann transport equations. The power factor of Zr_1-xNiSnHf_x increases as the temperature increases from 300 to 1200 K, because the increase in electrical conductivity is greater than the decrease in the Seebeck coefficient. The power factor of Zr_7/8NiSnHf_1/8 is larger than that of other Zr_1-xNiSnHf_x thermoelectric materials, but the thermoelectric figure of merit (ZT) is similar to that of others materials. This is due to the higher electronic thermal conductivity of Zr_7/8NiSnHf_1/8 compared to other materials. The maximum ZT of p-type (n-type) Zr_1-xNiSnHf_x is 0.98 (0.97), 0.9 (0.89), 0.83 (0.80), and 0.72 (0.73) at 300 K, 600 K, 900 K, and 1200 K, respectively, which are greater than those of the pure ZrNiSn. In conclusion, Hf-doped ZrNiSn can enhance the thermoelectric performance and are promising candidates for thermoelectric materials.</p><p><strong>Computational method: </strong>This paper uses FP-LAPW implemented in the WIEN2K code. The thermoelectric performance is calculated based on the semi-classical Boltzmann theory implanted using the BoltzTraP code. The electronic thermal conductivity (κ_e) and the carrier concentration (n) have been calculated using the density functional theory.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141974825","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":"Density functional theory guide for copolymerization mechanism between allyl radical with radicalophile: photo-driven radical mediated [3 + 2] cyclization.","authors":"Ou Liu, Piaoyi Chen, Qinglin Xiao, Chengfeng Yue, Yugang Huang, Guodong Ye","doi":"10.1007/s00894-024-06104-x","DOIUrl":"https://doi.org/10.1007/s00894-024-06104-x","url":null,"abstract":"<p><strong>Context: </strong>The challenge of activating inert allyl monomers for polymerization has persisted, prompting our proposal of the photo-driven radical mediated [3 + 2] cyclization reaction (PRMC). This innovative approach significantly expedites the homopolymerization of multi-allyl monomers, enabling the synthesis of embolic microspheres for hepatocellular carcinoma interventions. PRMC involves allyl monomers to form allylic radicals and then radicals participating in a cycloaddition reaction with unsaturated olefins as radicalophiles to form cyclopentane-based radical products. While extensively studied in the theoretical and experimental homopolymerization, PRMC's application in copolymerization remains unexplored. To address this knowledge gap, we explored the elementary reaction, selecting allyl methyl ether radicals (AMER) and α,β-unsaturated ketones as radicalophiles for copolymerization investigations by density functional theory (DFT) analysis. We quantified energy differences between ground and excited states of reactants, elucidated frontier molecular orbitals, and assessed thermodynamic data for copolymerization feasibility. We also evaluated the electronic properties of reactants, predicting the reactivity of radicalophiles and the interactions of intermolecular reactions. Additionally, we applied transition state theory and interaction/deformation models and conducted a local orbital analysis to comprehensively study excess electron distribution and gyration radius of cyclic radical product. Our findings offer vital insights into PRMC's potential in copolymerization. This research provides a robust theoretical foundation for practical application, enhancing the polymerization field.</p><p><strong>Methods: </strong>Based on density functional theory (DFT), the calculations were performed at the M06-2X/6-311 +  + G(d,p) level in/by Gaussian 16 package. Subsequently, our analytical results apply time-dependent density-functional theory (TD-DFT) and solvent modeling (SMD). Single-point energy calculations determine the driving force behind the radicals' reaction with radicalophiles. Furthermore, we assessed the electrostatic potential (ESP) of the reactants. The results of the calculations were visualized by the Multiwfn 3.6 and VMD 1.9 programs.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141970382","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":"Mechanical strain effect on the optoelectronic properties and photocatalysis applications of layered AlN/GaN nanoheterostructure.","authors":"Nitika, Sandeep Arora, Dharamvir Singh Ahlawat","doi":"10.1007/s00894-024-06103-y","DOIUrl":"https://doi.org/10.1007/s00894-024-06103-y","url":null,"abstract":"<p><strong>Context: </strong>The aim of this work is to use first principles calculations to examine the effects of different mechanical strains on the optoelectronic and photocatalytic capabilities of the 2D/2D nanoheterostructure of AlN/GaN. By utilizing the lmBJ (Meta-GGA) and PBEsol (GGA) functional, the bandgap of the nanoheterostructure is calculated and found to be 4.89 eV and 3.24 eV. Simulated 2D AlN/GaN nanoheterostructure exhibits exceptional optical and electronic characteristics under applied biaxial tensile and compressive strains. The band gap changes from 4.89 to 3.77 eV, while the energy gap nature transitions from direct to indirect during tensile strain fluctuations of 0% to 8%. Strain is also found to have a significant effect on the optical absorption peaks. And a 0-8% rise in tensile strain causes the initial absorption peak of the 2D AlN/GaN nanoheterostructure to shift from 4.88 to 4.20 eV, which results in a 14% red shift in photon energy for every 2% change in strain. Furthermore, the optimum bandgap and band edge positions of the 2D AlN/GaN nanoheterostructure enable the water redox process to produce hydrogen and oxygen for wide range of pH. Thus, modification via strain may be an effective method for altering the optical as well as electronic characteristics of a 2D AlN/GaN nanoheterostructure, and this study may pave the way for new applications of this material in optoelectronic devices in the future.</p><p><strong>Methods: </strong>In the current work, density functional theory is used to explore every attribute of the 2D AlN/GaN nanoheterostructure. To characterize the electronic exchange-correlation, we used the PBEsol functional. In order to prevent any interlayer contact between periodicity of images, a vacuum is produced along the z-direction of approximately 10 Å. To increase the precision of bandgap prediction, the electronic and optical characteristics were computed using the meta-GGA lmBJ functional. To account for interlayer van der Waals interactions, nanoheterostructure computations were performed using the DFT-D3 functional.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141974826","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":"Structural, wide band gap half-metallic, and pressure-dependent thermodynamic predictions of Li₂TMMgO₆ (TM = V, Nb, and Ta) double perovskites.","authors":"Wisam Ayad Ahmed Ahmed, Evren Görkem Özdemir, H A Rahnamaye Aliabad","doi":"10.1007/s00894-024-06107-8","DOIUrl":"https://doi.org/10.1007/s00894-024-06107-8","url":null,"abstract":"<p><strong>Context: </strong>Li₂VMgO₆, Li₂NbMgO₆, and Li₂TaMgO₆ double perovskite compounds were energetically the most stable in the FM phase. The lattice constants were 7.63 Å, 7.94 Å, and 7.95 Å, and the Curie temperatures were 910.451 K, 930.739 K, and 1258.821 K, respectively. The wide bandgap semiconductor characters were provided in the GGA-PBE methods as 2.139 eV, 4.209 eV, and 5.007 eV, respectively. This wide band gap semiconductor state in the majority carriers and the metallic state in the minority states made these double perovskites true half-metallic ferromagnetics. The bulk modulus obtained in the ground state calculations and the values obtained from thermodynamic calculations were relatively close. Debye temperatures in the initial state conditions were 747 K, 685.13 K, and 587.77 K, respectively. The total magnetic moment values were calculated as 3.00 µ_B/f.u. The most significant contribution to this value came from oxygen atoms.</p><p><strong>Methods: </strong>The theoretical calculations of Li₂VMgO₆, Li₂NbMgO₆, and Li₂TaMgO₆ double perovskite alloys were performed using the WIEN2k program developed by Blaha et al. The electronic calculations were made with GGA-PBE, GGA + mBJ, and GGA + U approximations in the space number 225 and the Fm-3 m symmetry group. The thermodynamic calculations were performed using Gibbs2. In thermodynamic calculations, temperature increases were determined as 100 K and temperature values were increased from 0 to 1200 K.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141970383","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":"First-principle calculations of the electronic, vibrational, and thermodynamic properties of nitrogen-rich salt of 3,6-dinitramino-1,2,4,5-tetrazine [(NH₄)₂(DNAT)].","authors":"Si-Jia Lei, Qi-Jun Liu, Fu-Sheng Liu, Zheng-Tang Liu, Wen-Shuo Yuan","doi":"10.1007/s00894-024-06098-6","DOIUrl":"https://doi.org/10.1007/s00894-024-06098-6","url":null,"abstract":"<p><strong>Context: </strong>Energy-containing materials such as explosives have attracted considerable interest recently. In the field of high-energy materials, tetrazine and its derivatives can largely meet the requirements of high nitrogen content and oxygen balance. Nitrogen-rich energetic salts are important research subjects. Nitrogen-rich salt of 3,6-dinitramino-1,2,4,5-tetrazine is a high-energy nitrogen-rich material, but there are few related studies. This paper systematically studies the crystal structure and electronic, vibrational, and thermodynamic properties of (NH₄)₂(DNAT). The lattice parameters of (NH₄)₂(DNAT) are observed to align well with the experimental values. The properties of electrons are analyzed by band structure and density of states (DOS). The phonon dispersion curves indicate that the compound is dynamically stable. The vibrational modes of bonds and chemical groups are described in detail, and the peaks in the Raman and infrared spectra are assigned to different vibration modes. Based on the vibration characteristics, thermodynamic properties such as enthalpy (H), Helmholtz free energy (F), entropy (S), Gibbs free energy (G), constant volume heat capacity (C_V), and Debye temperature (Θ) are analyzed. This article can pave the way for subsequent work or provide data support to other researchers, promoting further research.</p><p><strong>Methods: </strong>In this study, we utilized the density functional theory (DFT) for our calculations. The exchange-correlation potential and van der Waals interactions were characterized based on the GGA-PBE + G function calculation. We obtained Brillouin zone integrals using Monkhorst-Pack k-point grids, with the k-point of the Brillouin zone set to a 2 × 2 × 2 grid. During the self-consistent field operation, we set the total energy convergence tolerance to 5 × 10^-6 eV per atom. The cut-off energy for the calculation was established at 830 eV. Additionally, the states of H (1s¹), C (2s² 2p²), N (2s² 2p³), and O (2s² 2p⁴) were treated as valence electrons in our study.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141905430","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":"Molecular dynamics simulation of DNAN/DNB cocrystal PBXs.","authors":"Xin-Yi Li, Bao-Guo Wang, Ya-Fang Chen, Jian-Bo Fu, Ji-Hang Du, Chun-Guang Wang","doi":"10.1007/s00894-024-06096-8","DOIUrl":"https://doi.org/10.1007/s00894-024-06096-8","url":null,"abstract":"<p><strong>Context: </strong>The DNAN/DNB eutectic is a high-energy explosive eutectic with superior safety and thermal stability compared to traditional melt-cast explosives. However, the addition of polymer binders can effectively enhance its mechanical properties, allowing for continued production demands without the need for changes to existing factory equipment. In this paper, a model of the DNAN/DNB eutectic explosive was established, and five different types of polymers-cis-1,4-polybutadiene (BR), ethylene-vinyl acetate copolymer (EVA), polyethylene glycol (PEG), fluorinated polymer (F2603), and polyvinylidene fluoride (PVDF)-were added to the (1 0 - 1), (1 0 1), and (0 1 1) cleavage planes, respectively, to form polymer-bonded explosives (PBXs). The stability, trigger bond length, mechanical properties, and detonation performance of the various polymer-bound PBXs were predicted retrogressively. Among the five PBX models, the DNAN/DNB/PEG model exhibited the highest binding energy and the shortest trigger bond length, indicating a significant improvement in stability, compatibility, and sensitivity compared to the original eutectic. Additionally, although the detonation performance of DNAN/DNB decreased after the addition of binders, the final results were still satisfactory. Overall, the DNAN/DNB/PEG model demonstrated excellent comprehensive performance, proving that among the many polymer binders, PEG is the optimal choice for DNAN/DNB.</p><p><strong>Methods: </strong>Within the Materials Studio software, molecular dynamics (MD) simulations were employed to predict the properties of the DNAN/DNB eutectic PBX. The MD simulation timestep was set to 1 fs, with a cumulative simulation duration of 2 ns. A 2 ns MD simulation was conducted using the isothermal-isobaric ensemble (NPT). The COMPASS force field was applied, and the temperature was fixed at 295 K.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141900549","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":"Exchange-correlation kernel for perturbation dependent auxiliary functions in auxiliary density perturbation theory.","authors":"Luis I Hernández-Segura, Flor A Olvera-Rubalcava, Roberto Flores-Moreno, Patrizia Calaminici, Andreas M Köster","doi":"10.1007/s00894-024-06091-z","DOIUrl":"10.1007/s00894-024-06091-z","url":null,"abstract":"<p><strong>Context: </strong>Analytic exchange-correlation kernel formulations are of the outermost importance for density functional theory (DFT) perturbation calculations. In this paper, the working equation for the exchange-correlation kernel of the generalized gradient approximation (GGA) for perturbation dependent auxiliary functions is derived and discussed in the framework of auxiliary density functional theory (ADFT). The presented new formulation is extended to the unrestricted approach, too. A comprehensive discussion of the implementation of the GGA ADFT kernel, using either the native exchange-correlation functional implementations in deMon2k or the ones from the LibXC library, is given. Calculations with analytic exchange-correlation kernels are compared to their finite difference counterparts. The obtained results are in quantitative agreement. Nevertheless, analytic GGA ADFT kernel implementations show substantial improvement in the computational performance. Similar results are reported for analytic second derivatives of effective core potential (ECP) and model core potential (MCP) matrix elements when compared to their finite difference counterparts in molecular frequency analyses.</p><p><strong>Method: </strong>All calculations are performed in the framework of ADFT as implemented in deMon2k. In the ADFT analytic frequency calculations, auxiliary density perturbation theory was used. The underlying two-center exchange-correlation kernel matrix elements are calculated by numerical integration either with analytic or finite difference kernel expressions. Validation calculations are performed with the VWN and PBE functionals employing DFT-optimized DZVP basis sets in conjunction with automatically generated GEN-A2 auxiliary density function sets. In the (Pt₃Cu)_n cluster benchmark calculations, the RPBE functional was used. For Pt atoms, the quasi-relativistic LANL2DZ effective core potential with the corresponding valence basis set was employed, whereas for Cu atoms, the all-electron DFT-optimized TZVP basis was applied. The auxiliary density was expanded by the automatically generated GEN-A2* auxiliary function set. We run all benchmark calculations in parallel on 24 cores.</p>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":null,"pages":null},"PeriodicalIF":2.1,"publicationDate":"2024-08-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC11310252/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141900548","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}