Krystyna Syty, Grzegorz Czekało, Khanh Ngoc Pham, Marcin Modrzejewski
{"title":"Multi-Level Coupled-Cluster Description of Crystal Lattice Energies.","authors":"Krystyna Syty, Grzegorz Czekało, Khanh Ngoc Pham, Marcin Modrzejewski","doi":"10.1021/acs.jctc.5c00428","DOIUrl":"10.1021/acs.jctc.5c00428","url":null,"abstract":"<p><p>The many-body expansion (MBE) of the lattice energy enables an ab initio description of molecular solids using correlated wave function approximations. However, the practical application of MBE requires computing the large number of <i>n</i>-body contributions efficiently. To this end, we employ a multi-level coupled-cluster approach which adapts the approximation level based on interaction type and intermolecular distance. The high-level method, including connected triple excitations, is applied only to monomer relaxation and dimer interactions roughly within the first and second coordination shells. Long-range dimers and trimers are treated using a simplified coupled-cluster description based on the random-phase approximation (RPA). A key feature is an energy correction which mitigates the underbinding error of the base RPA. Convergence to the bulk limit is accelerated by the addition of the periodic Hartree-Fock correction. The proposed approach is validated against recent diffusion Monte Carlo reference data for the X23 data set, achieving a mean absolute error of 3.1 kJ/mol, i.e., chemical accuracy for absolute lattice energies.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5533-5544"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12159999/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144172215","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Driven Similarity Renormalization Group with a Large Active Space: Applications to Oligoacenes, Zeaxanthin, and Chromium Dimer.","authors":"Chenyang Li, Xiaoxue Wang, Huanchen Zhai, Wei-Hai Fang","doi":"10.1021/acs.jctc.5c00322","DOIUrl":"10.1021/acs.jctc.5c00322","url":null,"abstract":"<p><p>We present a new implementation of the driven similarity renormalization group (DSRG) based on a density matrix renormalization group (DMRG) reference. The explicit build of high-order reduced density matrices is avoided by forming matrix-product-state compressed intermediates. This algorithm facilitates the application of DSRG second- and third-order perturbation theories to dodecacene with an active space of 50 electrons in 50 orbitals. This active space appears the largest employed to date within the framework of internally contracted multireference formalism. The DMRG-DSRG approach is applied to several challenging systems, including the singlet-triplet gaps (Δ<sub>ST</sub>) of oligoacenes ranging from naphthalene to dodecacene, the vertical excitation energies of zeaxanthin, and the ground-state potential energy curve (PEC) of Cr<sub>2</sub> molecule. Our best estimate for the vertical Δ<sub>ST</sub> of dodecacene is 0.22 eV, showing an excellent agreement with that of the linearized adiabatic connection method (0.24 eV). For zeaxanthin, all DSRG schemes suggest the order of <math><mn>2</mn><mspace></mspace><mmultiscripts><msubsup><mi>A</mi><mi>g</mi><mo>-</mo></msubsup><mprescripts></mprescripts><none></none><mn>1</mn></mmultiscripts><mo><</mo><mn>1</mn><mspace></mspace><mmultiscripts><msubsup><mi>B</mi><mi>u</mi><mo>+</mo></msubsup><mprescripts></mprescripts><none></none><mn>1</mn></mmultiscripts><mo><</mo><mn>1</mn><mspace></mspace><mmultiscripts><msubsup><mi>B</mi><mi>u</mi><mo>-</mo></msubsup><mprescripts></mprescripts><none></none><mn>1</mn></mmultiscripts></math> for excited states. Both the equilibrium and the shoulder regions of the Cr<sub>2</sub> PEC are reasonably reproduced by the linearized DSRG with one- and two-body operators.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5445-5456"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144126131","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Analytical Gradient Theory for Density-Fitted Exact Two-Component Hartree-Fock, State-Specific Complete Active Space Self-Consistent Field, and Second-Order Møller-Plesset Perturbation Theories.","authors":"Jae Woo Park","doi":"10.1021/acs.jctc.5c00405","DOIUrl":"10.1021/acs.jctc.5c00405","url":null,"abstract":"<p><p>The exact two-component (X2C) relativistic quantum chemistry calculations can be used to describe scalar relativistic effects and spin-orbit couplings at reasonable computational cost. However, they have limited applicability to wave function-based quantum chemistry methods, particularly geometric optimizations and dynamics simulations, owing to the high computational demands of these methods in sizable molecular systems. In this work, we report our implementation of an analytical gradient algorithm with a density-fitting approximation for Hartree-Fock, state-specific complete active space self-consistent field (CASSCF), and second-order Møller-Plesset perturbation theory (MP2) calculations with the X2C one-electron Hamiltonian. This implementation uses a second-order orbital optimization scheme to facilitate convergence in X2C-CASSCF calculations, as well as a response (<i>Z</i>-vector) equation for evaluation of the X2C-MP2 nuclear gradient. We demonstrate the applicability of the algorithm for optimization of the geometry of Ir(ppy)<sub>2</sub>(bpy)<sup>+</sup> and evaluate its computational cost and parallelization (multithreading) efficiency.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5523-5532"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144179690","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Matthew Mostrom, Adam Hogan, Logan Ritter, William Morris, Brian Space
{"title":"PHAST-MBD: Implementing Many-Body Dispersion in the PHAST 2.0 Potential, Results for Noble Gases.","authors":"Matthew Mostrom, Adam Hogan, Logan Ritter, William Morris, Brian Space","doi":"10.1021/acs.jctc.5c00448","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c00448","url":null,"abstract":"<p><p>A recently published empirical force field (herein PHAST or PHAST 2.0) is employed in its many-body dispersion-corrected form (PHAST-MBD) to examine the effects of collective dispersion interactions. Rare gases are used as a systematic way to test increasing importance of van der Waals attractions in systems dominated by repulsion-dispersion that are a challenge to extant force fields. The effects of many-body dispersion were studied for liquid and supercritical fluid regime for the series Neon, Argon, Krypton and Xenon. The PHAST force field is a condensed phase atomistic molecular modeling potential that includes contributions from repulsion-dispersion, permanent electrostatics, and many-body polarization. Each of these pieces is physics based and seeks to mimic their constituent first-principles counterparts with as few fitting parameters as possible. Critically, it is built to reproduce accurate gas phase pair interactions. This facilitates the efficacy of mixing rules for unlike interactions while many-body effects are added via explicit polarization and dispersion models. The effectiveness of PHAST-MBD is demonstrated calculating rare gas densities as compared to experiment over a wide pressure range. Pair potentials fail systematically at high pressure and density as dispersion grows while PHAST-MBD reproduces experiment in all regimes. This is strong evidence in favor of the PHAST 2.0 paradigm of physically motivated empirical potentials that reproduce gas phase interactions and facilitate accurate mixing rules with many-body effects included explicitly. This work suggests a hybrid future approach that will be adopted in PHAST-MBD that keeps the accurate PHAST pair interactions and only includes many-body terms via the coupled dipole method (CDM); such an approach avoids the issues identified here that the CDM many body van der Waals (MBVDWs) formalism has reasonable but nonoptimal implicit mixing rules and can alter pair potentials.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":""},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144256714","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Xiao Wang, Kewei Sun, Haomin Xiao, Yang Zhao and Maxim F. Gelin*,
{"title":"","authors":"Xiao Wang, Kewei Sun, Haomin Xiao, Yang Zhao and Maxim F. Gelin*, ","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"21 11","pages":"XXX-XXX XXX-XXX"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://pubs.acs.org/doi/pdf/10.1021/acs.jctc.5c00153","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144354534","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Approximation to Second Order N-Electron Valence State Perturbation Theory: Limiting the Wave Function within Singles.","authors":"Yang Guo, Katarzyna Pernal","doi":"10.1021/acs.jctc.5c00582","DOIUrl":"10.1021/acs.jctc.5c00582","url":null,"abstract":"<p><p>Inspired by the linearized adiabatic connection (AC0) theory, an approximation to second-order N-electron valence state perturbation theory (NEVPT2) has been developed, termed NEVPT within singles (NEVPTS). This approach utilizes amplitudes derived from approximate single-excitation wave functions, requiring only 3rd-order reduced density matrices (RDMs). Consequently, it avoids the computational bottleneck associated with the construction of 4th-order RDMs in NEVPT2. The NEVPTS method demonstrates comparable performance to NEVPT2 in describing potential energy curves for diatomic molecules and singlet-triplet gaps in biradicals, while achieving superior accuracy to AC0 in these applications. For excitation energies of organic molecules, NEVPTS is less accurate than NEVPT2. The overall performance and computational costs of the NEVPTS method lie between those of NEVPT2 and AC0.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5545-5558"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144172211","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Balint Dudas, Christina Athanasiou, Juan Carlos Mobarec, Edina Rosta
{"title":"Quantifying Cooperativity through Binding Free Energies in Molecular Glue Degraders.","authors":"Balint Dudas, Christina Athanasiou, Juan Carlos Mobarec, Edina Rosta","doi":"10.1021/acs.jctc.5c00064","DOIUrl":"10.1021/acs.jctc.5c00064","url":null,"abstract":"<p><p>Molecular glues represent a novel therapeutic modality facilitating the stabilization of protein-protein interactions (PPIs), thus enabling the targeting of previously \"undruggable\" proteins. We develop a computational procedure to screen for molecular glues using a pathway-independent free energy calculation method for accurately assessing the cooperativity. We employ a combined ligand and protein free energy perturbation (FEP) method to calculate the cooperative effect of a ligand for ternary binding. We study the cooperative binding mechanisms of molecular glue degraders, specifically cereblon (CRBN) modulators targeting Ikaros family zinc finger 2 (IKZF2), a transcription factor implicated in cancer immunotherapy. We present a comprehensive computational protocol for screening large molecular libraries to identify potent molecular glues. By leveraging cooperative binding principles in ternary complex formation, our approach effectively predicts ligand-induced PPIs and their degradation potential. Benchmarking against experimental data for CRBN-Ikaros complexes, our protocol demonstrates high accuracy in identifying superior molecular glues, highlighting L4 and L5 as top performers. Furthermore, our high-throughput screening identified novel candidates from extensive chemical libraries, validated through advanced FEP+ simulations. This study not only underscores the transformative potential of molecular glues in targeted protein degradation but also sets the stage for their broader application across diverse protein targets, paving the way for innovative therapeutic strategies in drug discovery.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5712-5723"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.ncbi.nlm.nih.gov/pmc/articles/PMC12159975/pdf/","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143952402","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Amir Ayati, Hugh G A Burton, Stijn De Baerdemacker
{"title":"Spin-Symmetry Projected Constrained Unrestricted Hartree-Fock.","authors":"Amir Ayati, Hugh G A Burton, Stijn De Baerdemacker","doi":"10.1021/acs.jctc.5c00068","DOIUrl":"10.1021/acs.jctc.5c00068","url":null,"abstract":"<p><p>We introduce an electronic structure approach for spin-symmetry breaking and restoration at the mean-field level. The spin-projected constrained-unrestricted Hartree-Fock (SPcUHF) method restores the broken spin symmetry inherent in spin-constrained-UHF determinants by employing a nonorthogonal configuration interaction (NOCI) projection method. This method includes all possible configurations in spin space compatible with a Clebsch-Gordon recoupling scheme in a NOCI calculation. The tunable one-pair-at-a-time characteristics of the symmetry-breaking process in c-UHF allow us to reduce the computational costs of full projection. SPcUHF is tested on 4-, 6-, and 8-electron systems that exhibit dominant static and/or dynamic correlations.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5409-5424"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144118402","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Lucas Lang, Letizia Fiorucci, Giacomo Parigi, Claudio Luchinat, Enrico Ravera
{"title":"Theory of Field-Dependent NMR Shifts in Paramagnetic Molecules.","authors":"Lucas Lang, Letizia Fiorucci, Giacomo Parigi, Claudio Luchinat, Enrico Ravera","doi":"10.1021/acs.jctc.5c00433","DOIUrl":"10.1021/acs.jctc.5c00433","url":null,"abstract":"<p><p>NMR chemical shifts depend on the applied magnetic flux density, and this becomes more and more important as stronger and stronger magnetic fields are becoming available. Herein we develop a theory of the field dependence of NMR shifts of paramagnetic molecules in solution. Our derivation leads to two distinct approaches: a finite-field approach that describes the shift up to infinite order in the applied field <i>B</i><sub>0</sub> but requires numerical integration for the orientational average, and a second-order approach that is valid up to second order in <i>B</i><sub>0</sub>. In this latter approach, the orientational average can be performed analytically, and the field dependence cleanly separates into two additive terms: the well-known \"indirect\" field dependence due to incomplete averaging in solution and the \"direct\" field dependence due to the nonlinear response to the external field. In analogy to the diamagnetic case, the direct field dependence involves a fourth-order tensor <b>τ</b> whose elements are fourth derivatives of the electronic Helmholtz free energy. Generalizing the Van den Heuvel-Soncini equation, we provide analytical sum-over-states equations for these higher-order derivatives. Using the NiSAL-HDPT complex as an example, we demonstrate the applicability of the second-order approach at room temperature and the highest commercially available field strength and show that it agrees well with the field dependence measured experimentally.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5642-5660"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144155242","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Silvie Illésová, Martin Beseda, Saad Yalouz, Benjamin Lasorne, Bruno Senjean
{"title":"Transformation-Free Generation of a Quasi-Diabatic Representation from the State-Average Orbital-Optimized Variational Quantum Eigensolver.","authors":"Silvie Illésová, Martin Beseda, Saad Yalouz, Benjamin Lasorne, Bruno Senjean","doi":"10.1021/acs.jctc.5c00327","DOIUrl":"10.1021/acs.jctc.5c00327","url":null,"abstract":"<p><p>In the present work, we examine how the recent quantum-computing algorithm known as the state-average orbital-optimized variational quantum eigensolver (SA-OO-VQE), viewed within the context of quantum chemistry as a type of multiconfiguration self-consistent field (MCSCF) electronic-structure approach, exhibits a propensity to produce an <i>ab initio</i> quasi-diabatic representation \"for free\" if considered as a least-transformed block-diagonalization procedure, as alluded to in our previous work (Yalouz, S. et al. <i>J. Chem. Theory Comput.</i> <b>2022</b>, <i>18</i>, 776-794) and thoroughly assessed herein. To this end, we introduce intrinsic and residual descriptors of diabaticity and re-explore the definition and linear-algebra properties, as well as their consequences on the vibronic nonadiabatic couplings, of an optimal diabatic representation within this context and how much one may deviate from it. Such considerations are illustrated numerically in the prototypical case of formaldimine, which presents a well-known conical intersection between its ground and first-excited singlet electronic states.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":"5457-5480"},"PeriodicalIF":5.7,"publicationDate":"2025-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144131808","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":1,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}