首页 > 最新文献

Journal of Chemical Theory and Computation最新文献

英文 中文
Performance of Universal Machine-Learned Potentials with Explicit Long-Range Interactions in Biomolecular Simulations. 生物分子模拟中具有显式远程相互作用的通用机器学习势的性能。
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-05 DOI: 10.1021/acs.jctc.5c01828
Viktor Zaverkin, Matheus Ferraz, Francesco Alesiani, Mathias Niepert

Universal machine-learned potentials promise transferable accuracy across compositional and vibrational degrees of freedom, yet their application to biomolecular simulations remains underexplored. This work systematically evaluates equivariant message-passing architectures trained on the SPICE-v2 data set with and without explicit long-range dispersion and electrostatics. We assess the impact of model size, training data composition, and electrostatic treatment across in- and out-of-distribution benchmark data sets, as well as molecular simulations of bulk liquid water, aqueous NaCl solutions, and biomolecules, including alanine tripeptide, the mini-protein Trp-cage, and Crambin. While larger models improve accuracy on benchmark data sets, this trend does not consistently extend to properties obtained from simulations. Predicted properties also depend on the composition of the training data set. Long-range electrostatics show no systematic impact across systems. However, for Trp-cage, their inclusion yields increased conformational variability. Our results suggest that imbalanced data sets and immature evaluation practices currently challenge the applicability of universal machine-learned potentials to biomolecular simulations.

{"title":"Performance of Universal Machine-Learned Potentials with Explicit Long-Range Interactions in Biomolecular Simulations.","authors":"Viktor Zaverkin, Matheus Ferraz, Francesco Alesiani, Mathias Niepert","doi":"10.1021/acs.jctc.5c01828","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01828","url":null,"abstract":"<p><p>Universal machine-learned potentials promise transferable accuracy across compositional and vibrational degrees of freedom, yet their application to biomolecular simulations remains underexplored. This work systematically evaluates equivariant message-passing architectures trained on the SPICE-v2 data set with and without explicit long-range dispersion and electrostatics. We assess the impact of model size, training data composition, and electrostatic treatment across in- and out-of-distribution benchmark data sets, as well as molecular simulations of bulk liquid water, aqueous NaCl solutions, and biomolecules, including alanine tripeptide, the mini-protein Trp-cage, and Crambin. While larger models improve accuracy on benchmark data sets, this trend does not consistently extend to properties obtained from simulations. Predicted properties also depend on the composition of the training data set. Long-range electrostatics show no systematic impact across systems. However, for Trp-cage, their inclusion yields increased conformational variability. Our results suggest that imbalanced data sets and immature evaluation practices currently challenge the applicability of universal machine-learned potentials to biomolecular simulations.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146123205","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}
引用次数: 0
Reduced Density Matrix and Cumulant Approximations of Quantum Linear Response. 量子线性响应的约化密度矩阵和累积量近似。
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-05 DOI: 10.1021/acs.jctc.5c01353
Theo Juncker von Buchwald, Erik Rosendahl Kjellgren, Jacob Kongsted, Stephan P A Sauer, Sonia Coriani, Karl Michael Ziems

Linear response (LR) is an important tool in the computational chemist's toolbox. It is therefore unsurprising that the emergence of quantum computers has led to a quantum counterpart known as quantum LR (qLR). However, the current quantum era of near-term intermediate-scale quantum (NISQ) computers is dominated by noise, short decoherence times, and slow measurement speeds. It is therefore of interest to find approximations that can greatly reduce the quantum workload while only slightly impacting the quality of a method. In an effort to achieve this, we approximate the naive qLR with the singles and doubles (qLRSD) method, by either directly approximating the reduced density matrices (RDMs) or indirectly through their respective reduced density cumulants (RDCs). We present an analysis of the measurement costs associated with qLR using RDMs and report qLR results for model hydrogen ladder systems; for varying active space sizes in OCS, SeH2, and H2S; and for symmetrically stretched H2O and BeH2. Discouragingly, while approximations to the 4-body RDMs and RDCs seem to produce good results for systems at the equilibrium geometry and for some types of core excitations, they both tend to fail when the system exhibits strong correlation. All approximations to the 3-body RDMs and/or RDCs severely affect the results and cannot be applied.

线性响应(LR)是计算化学家工具箱中的一个重要工具。因此,量子计算机的出现导致了量子对应物量子LR (qLR)的出现也就不足为奇了。然而,当前的量子时代,短期的中等规模量子(NISQ)计算机被噪声、短退相干时间和缓慢的测量速度所主导。因此,找到可以大大减少量子工作负载,同时只略微影响方法质量的近似值是很有意义的。为了实现这一目标,我们用单双元(qLRSD)方法近似朴素qLR,通过直接近似约简密度矩阵(rdm)或间接通过其各自的约简密度累积量(rdc)。我们使用rdm分析了与qLR相关的测量成本,并报告了模型氢阶梯系统的qLR结果;用于OCS、SeH2和H2S中不同的活性空间大小;以及对称拉伸的H2O和BeH2。令人沮丧的是,虽然对四体rdm和rdc的近似似乎对平衡几何和某些类型的核心激励的系统产生良好的结果,但当系统表现出强相关性时,它们往往都失败。所有对三体rdm和/或rdc的近似都严重影响结果,不能应用。
{"title":"Reduced Density Matrix and Cumulant Approximations of Quantum Linear Response.","authors":"Theo Juncker von Buchwald, Erik Rosendahl Kjellgren, Jacob Kongsted, Stephan P A Sauer, Sonia Coriani, Karl Michael Ziems","doi":"10.1021/acs.jctc.5c01353","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01353","url":null,"abstract":"<p><p>Linear response (LR) is an important tool in the computational chemist's toolbox. It is therefore unsurprising that the emergence of quantum computers has led to a quantum counterpart known as quantum LR (qLR). However, the current quantum era of near-term intermediate-scale quantum (NISQ) computers is dominated by noise, short decoherence times, and slow measurement speeds. It is therefore of interest to find approximations that can greatly reduce the quantum workload while only slightly impacting the quality of a method. In an effort to achieve this, we approximate the naive qLR with the singles and doubles (qLRSD) method, by either directly approximating the reduced density matrices (RDMs) or indirectly through their respective reduced density cumulants (RDCs). We present an analysis of the measurement costs associated with qLR using RDMs and report qLR results for model hydrogen ladder systems; for varying active space sizes in OCS, SeH<sub>2</sub>, and H<sub>2</sub>S; and for symmetrically stretched H<sub>2</sub>O and BeH<sub>2</sub>. Discouragingly, while approximations to the 4-body RDMs and RDCs seem to produce good results for systems at the equilibrium geometry and for some types of core excitations, they both tend to fail when the system exhibits strong correlation. All approximations to the 3-body RDMs and/or RDCs severely affect the results and cannot be applied.</p>","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":" ","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146117165","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}
引用次数: 0
Transition State Theory for Dissociation of Dynamic Bonding Networks 动态键网络解离的过渡态理论
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01797
Eric V. Anslyn,Dmitrii E. Makarov
Cell adhesion, molecular recognition, biomolecular folding and unfolding, dynamic cross-linking in soft materials, and many other phenomena involve formation or dissociation of multiple chemical bonds. Here, we study the overall time scale required to break or form N bonds. Strictly speaking, this time scale depends on the initial conditions, e.g., the number and which bonds are formed/broken, and its estimation requires kinetic details about forming and breaking of each individual bond influenced by the larger network of other bonds. We show, however, that a simple estimate, analogous to transition-state theory in chemical kinetics, accurately predicts the mean first passage time to form or break all the bonds in terms of single-bond properties and thermodynamic properties of the network. As the thermodynamics of bond networks can often be described by well-studied statistical-mechanical models, such as the Ising model and its extensions, our theory provides a link between the global dynamics and thermodynamics of multibond arrays and networks.
细胞粘附,分子识别,生物分子折叠和展开,软材料中的动态交联以及许多其他现象涉及多个化学键的形成或解离。在这里,我们研究了打破或形成N键所需的总体时间尺度。严格地说,这个时间尺度取决于初始条件,例如,形成/断裂的键的数量和哪些键,它的估计需要受其他更大的键网络影响的每个单独键的形成和断裂的动力学细节。然而,我们表明,一个简单的估计,类似于化学动力学中的过渡状态理论,准确地预测了在单键性质和网络的热力学性质方面形成或打破所有键的平均首次通过时间。由于键网络的热力学通常可以通过充分研究的统计力学模型来描述,例如Ising模型及其扩展,我们的理论提供了多键阵列和网络的全局动力学和热力学之间的联系。
{"title":"Transition State Theory for Dissociation of Dynamic Bonding Networks","authors":"Eric V. Anslyn,Dmitrii E. Makarov","doi":"10.1021/acs.jctc.5c01797","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01797","url":null,"abstract":"Cell adhesion, molecular recognition, biomolecular folding and unfolding, dynamic cross-linking in soft materials, and many other phenomena involve formation or dissociation of multiple chemical bonds. Here, we study the overall time scale required to break or form N bonds. Strictly speaking, this time scale depends on the initial conditions, e.g., the number and which bonds are formed/broken, and its estimation requires kinetic details about forming and breaking of each individual bond influenced by the larger network of other bonds. We show, however, that a simple estimate, analogous to transition-state theory in chemical kinetics, accurately predicts the mean first passage time to form or break all the bonds in terms of single-bond properties and thermodynamic properties of the network. As the thermodynamics of bond networks can often be described by well-studied statistical-mechanical models, such as the Ising model and its extensions, our theory provides a link between the global dynamics and thermodynamics of multibond arrays and networks.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"88 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111077","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}
引用次数: 0
Kinetic Monte Carlo Framework for Coupled Degradation and Dehydration of Anion Exchange Membranes 阴离子交换膜耦合降解脱水的动力学蒙特卡罗框架
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c02063
Esteban D. Gadea,Shakkira Erimban,Ignacio J. Bombau,Damian A. Scherlis,John J. Karnes,Valeria Molinero
Kinetic Monte Carlo (kMC) simulations, augmented with temporal-acceleration schemes, can efficiently handle stiff reaction-transport networks when fast processes rapidly relax to quasi-equilibrium on a fixed lattice. However, in glassy anion-exchange membranes (AEM), rare and irreversible chemical degradation events continuously reshape the nanoscale morphology, and the associated hydration and transport degrees of freedom remain far from a well-defined local equilibrium. This combination of evolving state space and nonequilibrated fast dynamics lies outside the scope of existing kMC acceleration frameworks. To address this challenge, we introduce an auxiliary-particle kinetic Monte Carlo (AP-kMC) scheme. In AP-kMC, short-lived mobile particles spawned at degradation sites execute hop, water-elimination, and decay moves, enforcing rapid local relaxation of the hydration structure while preserving the stochastic rules of kMC. Parameterized with molecular-dynamics morphologies and experimental solution degradation kinetics, AP-kMC reproduces the evolution of ion-exchange capacity, water uptake, and conductivity, and reveals a feedback loop in which poorly hydrated sites degrade first and each degradation event induces further local dehydration. The resulting thinning and fragmentation of water channels cause loss of hydrophilic percolation and abrupt conductivity collapse well before complete charge loss. AP-kMC thus reframes AEM durability as a coupled degradation–drying–percolation problem and provides a transferable strategy to simulate reactive, out-of-equilibrium polymer electrolytes where local solvation controls reactivity.
动力学蒙特卡罗(kMC)模拟,增加了时间加速方案,可以有效地处理刚性反应-输运网络,当快速过程在固定晶格上迅速放松到准平衡状态时。然而,在玻璃阴离子交换膜(AEM)中,罕见且不可逆的化学降解事件不断重塑纳米级形态,相关的水化和运输自由度仍远未达到明确的局部平衡。这种演化状态空间和非平衡快速动力学的结合超出了现有kMC加速框架的范围。为了解决这一挑战,我们引入了辅助粒子动力学蒙特卡罗(AP-kMC)方案。在AP-kMC中,在降解位点产生的短寿命移动颗粒执行跳跃、水消除和衰变运动,在保持kMC随机规则的同时,强制水化结构的局部快速松弛。通过分子动力学形态学和实验溶液降解动力学参数化,AP-kMC再现了离子交换能力、吸水率和电导率的演变,并揭示了一个反馈回路,在这个反馈回路中,水化程度较低的部位首先降解,每次降解事件都会引起进一步的局部脱水。由此导致的水通道变薄和破碎导致亲水性渗透的损失和电导率的突然崩溃,甚至在完全电荷损失之前。因此,AP-kMC将AEM耐久性重新定义为一个耦合的降解-干燥-渗透问题,并提供了一种可转移的策略来模拟反应性,非平衡聚合物电解质,其中局部溶剂化控制反应性。
{"title":"Kinetic Monte Carlo Framework for Coupled Degradation and Dehydration of Anion Exchange Membranes","authors":"Esteban D. Gadea,Shakkira Erimban,Ignacio J. Bombau,Damian A. Scherlis,John J. Karnes,Valeria Molinero","doi":"10.1021/acs.jctc.5c02063","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c02063","url":null,"abstract":"Kinetic Monte Carlo (kMC) simulations, augmented with temporal-acceleration schemes, can efficiently handle stiff reaction-transport networks when fast processes rapidly relax to quasi-equilibrium on a fixed lattice. However, in glassy anion-exchange membranes (AEM), rare and irreversible chemical degradation events continuously reshape the nanoscale morphology, and the associated hydration and transport degrees of freedom remain far from a well-defined local equilibrium. This combination of evolving state space and nonequilibrated fast dynamics lies outside the scope of existing kMC acceleration frameworks. To address this challenge, we introduce an auxiliary-particle kinetic Monte Carlo (AP-kMC) scheme. In AP-kMC, short-lived mobile particles spawned at degradation sites execute hop, water-elimination, and decay moves, enforcing rapid local relaxation of the hydration structure while preserving the stochastic rules of kMC. Parameterized with molecular-dynamics morphologies and experimental solution degradation kinetics, AP-kMC reproduces the evolution of ion-exchange capacity, water uptake, and conductivity, and reveals a feedback loop in which poorly hydrated sites degrade first and each degradation event induces further local dehydration. The resulting thinning and fragmentation of water channels cause loss of hydrophilic percolation and abrupt conductivity collapse well before complete charge loss. AP-kMC thus reframes AEM durability as a coupled degradation–drying–percolation problem and provides a transferable strategy to simulate reactive, out-of-equilibrium polymer electrolytes where local solvation controls reactivity.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"8 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111136","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}
引用次数: 0
A Hardware-Feasible Quantum Machine Learning Framework for Structure-Based Virtual Screening 一种硬件可行的基于结构的虚拟筛选量子机器学习框架
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01872
Pei-Kun Yang
In structure-based virtual screening, evaluating the binding free energy of protein–ligand complexes requires accounting for both molecular conformations and spatial transformations, such as shifts and rotations, which can lead to an exponential increase in possible configurations. Classical computing approaches are limited in handling this combinatorial explosion, whereas quantum computing offers a promising alternative due to its inherent parallelism. In this study, we propose a quantum machine learning framework that encodes molecular features into quantum states and processes them through parametrized quantum gates, with all architectural and representational choices deliberately guided by near-term hardware feasibility and an explicit focus on minimal qubit counts, shallow circuit depth, and compact input representations. The model is implemented and optimized in PyTorch, and its predictive performance is examined under three conditions: ideal simulation, finite-shot sampling, and quantum-noise simulation. With six quantum circuit units, the model achieves a root-mean-square deviation of 2.37 kcal/mol and a Pearson correlation coefficient of 0.650. The predictions remain stable with 100,000 measurement shots, demonstrating compatibility with near-term quantum hardware. Although the introduction of noise slightly reduces absolute accuracy, the Pearson correlation coefficient remains stable, indicating that the ranking of ligand affinities is preserved. These results highlight a practical, scalable quantum approach that balances predictive power and robustness, providing a feasible pathway to accelerate virtual screening using moderately deep quantum circuits.
在基于结构的虚拟筛选中,评估蛋白质-配体复合物的结合自由能需要考虑分子构象和空间转换,如移位和旋转,这可能导致可能的构型呈指数增长。经典计算方法在处理这种组合爆炸方面受到限制,而量子计算由于其固有的并行性提供了一个有前途的替代方案。在本研究中,我们提出了一个量子机器学习框架,该框架将分子特征编码为量子态,并通过参数化量子门对其进行处理,所有架构和表征选择都刻意以近期硬件可行性为指导,并明确关注最小量子位计数、浅电路深度和紧凑的输入表示。在PyTorch中对该模型进行了实现和优化,并在理想模拟、有限采样和量子噪声模拟三种条件下对其预测性能进行了测试。采用6个量子电路单元,模型的均方根偏差为2.37 kcal/mol, Pearson相关系数为0.650。这些预测在10万次测量中保持稳定,证明了与近期量子硬件的兼容性。虽然噪声的引入略微降低了绝对精度,但Pearson相关系数保持稳定,表明配体亲和度的排序保持不变。这些结果突出了一种实用的、可扩展的量子方法,它平衡了预测能力和鲁棒性,为使用适度深度量子电路加速虚拟筛选提供了可行的途径。
{"title":"A Hardware-Feasible Quantum Machine Learning Framework for Structure-Based Virtual Screening","authors":"Pei-Kun Yang","doi":"10.1021/acs.jctc.5c01872","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01872","url":null,"abstract":"In structure-based virtual screening, evaluating the binding free energy of protein–ligand complexes requires accounting for both molecular conformations and spatial transformations, such as shifts and rotations, which can lead to an exponential increase in possible configurations. Classical computing approaches are limited in handling this combinatorial explosion, whereas quantum computing offers a promising alternative due to its inherent parallelism. In this study, we propose a quantum machine learning framework that encodes molecular features into quantum states and processes them through parametrized quantum gates, with all architectural and representational choices deliberately guided by near-term hardware feasibility and an explicit focus on minimal qubit counts, shallow circuit depth, and compact input representations. The model is implemented and optimized in PyTorch, and its predictive performance is examined under three conditions: ideal simulation, finite-shot sampling, and quantum-noise simulation. With six quantum circuit units, the model achieves a root-mean-square deviation of 2.37 kcal/mol and a Pearson correlation coefficient of 0.650. The predictions remain stable with 100,000 measurement shots, demonstrating compatibility with near-term quantum hardware. Although the introduction of noise slightly reduces absolute accuracy, the Pearson correlation coefficient remains stable, indicating that the ranking of ligand affinities is preserved. These results highlight a practical, scalable quantum approach that balances predictive power and robustness, providing a feasible pathway to accelerate virtual screening using moderately deep quantum circuits.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"1 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111083","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}
引用次数: 0
Multistate Iterative Qubit Coupled Cluster (MS-iQCC): A Quantum-Inspired, State-Averaged Approach to Ground- And Excited-State Energies 多态迭代量子比特耦合簇(MS-iQCC):一种量子启发的、状态平均的基态和激发态能量方法
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01849
Robert A. Lang,Shashank G. Mehendale,Ilya G. Ryabinkin,Artur F. Izmaylov
We introduce the multistate iterative qubit coupled cluster (MS-iQCC) method, a quantum-inspired algorithm that runs efficiently on classical hardware and is designed to predict both ground and excited electronic states of molecules. Accurate excited-state energetics are essential for interpreting spectroscopy and chemical reactivity, but standard electronic structure methods are either too computationally expensive for larger systems or lose reliability in the presence of strong electron correlation. MS-iQCC addresses this challenge by simultaneously optimizing multiple electronic states in a single, state-averaged procedure that treats ground and excited states on equal footing. This removes the energetic bias that is introduced when excited states are computed one at a time and constrained to remain orthogonal to previously optimized states. The approach supports multireference electronic structure by allowing multideterminantal initial guesses and by adaptively building a compact exponential ansatz from a pool of qubit excitation generators. We apply MS-iQCC to H4, H2O, N2, and C2, including strongly correlated geometries, and observe robust convergence of all targeted state energies to chemically meaningful accuracy across their potential energy surfaces.
我们介绍了多态迭代量子比特耦合簇(MS-iQCC)方法,这是一种在经典硬件上高效运行的量子启发算法,旨在预测分子的基态和激发态。准确的激发态能量学对于解释光谱和化学反应性至关重要,但是标准的电子结构方法对于较大的系统来说计算成本太高,或者在存在强电子相关性的情况下失去可靠性。MS-iQCC解决了这一挑战,同时优化多个电子状态在一个单一的,状态平均的程序,对待基态和激发态在平等的基础上。这消除了一次计算一个激发态时引入的能量偏差,并且约束保持与先前优化的状态正交。该方法支持多参考电子结构,允许多确定性初始猜测,并自适应地从量子位激发发生器池中构建紧凑的指数ansatz。我们将MS-iQCC应用于H4, H2O, N2和C2,包括强相关的几何结构,并观察到所有目标状态能量在其势能表面上的鲁棒收敛,具有化学意义的准确性。
{"title":"Multistate Iterative Qubit Coupled Cluster (MS-iQCC): A Quantum-Inspired, State-Averaged Approach to Ground- And Excited-State Energies","authors":"Robert A. Lang,Shashank G. Mehendale,Ilya G. Ryabinkin,Artur F. Izmaylov","doi":"10.1021/acs.jctc.5c01849","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01849","url":null,"abstract":"We introduce the multistate iterative qubit coupled cluster (MS-iQCC) method, a quantum-inspired algorithm that runs efficiently on classical hardware and is designed to predict both ground and excited electronic states of molecules. Accurate excited-state energetics are essential for interpreting spectroscopy and chemical reactivity, but standard electronic structure methods are either too computationally expensive for larger systems or lose reliability in the presence of strong electron correlation. MS-iQCC addresses this challenge by simultaneously optimizing multiple electronic states in a single, state-averaged procedure that treats ground and excited states on equal footing. This removes the energetic bias that is introduced when excited states are computed one at a time and constrained to remain orthogonal to previously optimized states. The approach supports multireference electronic structure by allowing multideterminantal initial guesses and by adaptively building a compact exponential ansatz from a pool of qubit excitation generators. We apply MS-iQCC to H4, H2O, N2, and C2, including strongly correlated geometries, and observe robust convergence of all targeted state energies to chemically meaningful accuracy across their potential energy surfaces.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"295 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111078","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}
引用次数: 0
Braess’ Paradox in Enzyme Kinetics: Asymmetry from Population Balance without Direct Cooperativity 酶动力学中的Braess悖论:非直接协同性种群平衡的不对称性
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01269
Malte Schäffner,Colin A. Smith,Robert Tampé,Helmut Grubmüller
The ATPase ABCE1, a member of the ubiquitous ATP-Binding Cassette protein superfamily, is essential in eukaryotic and archaeal ribosome recycling. It comprises a pair of homologous nucleotide-binding domains (NBDs), each containing a consensus nucleotide-binding site (NBS), where ATP hydrolysis takes place. Each of these sites can be in either an open or closed conformation. Despite the near symmetry of the two NBDs, and quite unexpectedly, their hydrolysis kinetics are highly asymmetric. While substitution of the catalytic glutamate (E238Q) in NBSI reduced the overall turnover rate of the ATPase by a factor of 2, as one might expect, the corresponding substitution in NBSII (E485Q) shows a so far unexplained 10-fold increase. To address this issue, we used Markov models to study how such a drastic asymmetry can arise. Specifically, we asked whether this observation can be explained without previously proposed direct allosteric interactions, such as electrostatic interactions, between the two NBSs. Indeed, using a Bayesian approach, we found Markov models that quantitatively predict the experimentally observed kinetics, as well as additional steady-state ATP occupancy data, both without such direct allosteric interaction. In particular, our results show that the observed remarkable asymmetry is fully explained by the structure-induced property that opening and closing always involves both NBSs. These models can explain the unexpected fast kinetics of the mutant of NBSII in terms of a drastic population shift due to the mutation, which circumvents a kinetic trap state that slows wild-type kinetics. Our Bayesian Markov approach may help to quantitatively explain similar nonintuitive Braess-type kinetics also in other enzymes where chemical/conformation coupling is essential.
atp酶ABCE1是普遍存在的atp结合盒蛋白超家族的成员,在真核生物和古细菌核糖体循环中起着至关重要的作用。它包括一对同源核苷酸结合域(nbd),每个都包含一个一致的核苷酸结合位点(NBS), ATP水解发生在这里。每个位点都可以是开放构象或封闭构象。出乎意料的是,尽管这两种nbd接近对称,但它们的水解动力学却是高度不对称的。虽然NBSI中催化谷氨酸(E238Q)的取代使atp酶的总周转率降低了2倍,正如人们所预料的那样,NBSII中相应的取代(E485Q)显示出迄今尚未解释的10倍的增加。为了解决这个问题,我们使用马尔可夫模型来研究这种剧烈的不对称是如何产生的。具体来说,我们想知道,如果没有先前提出的两种nbs之间的直接变构相互作用(如静电相互作用),是否可以解释这一观察结果。事实上,使用贝叶斯方法,我们发现马尔可夫模型可以定量预测实验观察到的动力学,以及额外的稳态ATP占用数据,两者都没有这种直接的变构相互作用。特别是,我们的研究结果表明,观察到的显著不对称性完全可以用结构诱导的性质来解释,即打开和关闭总是涉及两个nbs。这些模型可以解释NBSII突变体意想不到的快速动力学,因为突变引起了急剧的种群转移,从而绕过了减缓野生型动力学的动力学陷阱状态。我们的贝叶斯马尔可夫方法可能有助于定量解释类似的非直观的braess型动力学,也在其他酶的化学/构象耦合是必不可少的。
{"title":"Braess’ Paradox in Enzyme Kinetics: Asymmetry from Population Balance without Direct Cooperativity","authors":"Malte Schäffner,Colin A. Smith,Robert Tampé,Helmut Grubmüller","doi":"10.1021/acs.jctc.5c01269","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01269","url":null,"abstract":"The ATPase ABCE1, a member of the ubiquitous ATP-Binding Cassette protein superfamily, is essential in eukaryotic and archaeal ribosome recycling. It comprises a pair of homologous nucleotide-binding domains (NBDs), each containing a consensus nucleotide-binding site (NBS), where ATP hydrolysis takes place. Each of these sites can be in either an open or closed conformation. Despite the near symmetry of the two NBDs, and quite unexpectedly, their hydrolysis kinetics are highly asymmetric. While substitution of the catalytic glutamate (E238Q) in NBSI reduced the overall turnover rate of the ATPase by a factor of 2, as one might expect, the corresponding substitution in NBSII (E485Q) shows a so far unexplained 10-fold increase. To address this issue, we used Markov models to study how such a drastic asymmetry can arise. Specifically, we asked whether this observation can be explained without previously proposed direct allosteric interactions, such as electrostatic interactions, between the two NBSs. Indeed, using a Bayesian approach, we found Markov models that quantitatively predict the experimentally observed kinetics, as well as additional steady-state ATP occupancy data, both without such direct allosteric interaction. In particular, our results show that the observed remarkable asymmetry is fully explained by the structure-induced property that opening and closing always involves both NBSs. These models can explain the unexpected fast kinetics of the mutant of NBSII in terms of a drastic population shift due to the mutation, which circumvents a kinetic trap state that slows wild-type kinetics. Our Bayesian Markov approach may help to quantitatively explain similar nonintuitive Braess-type kinetics also in other enzymes where chemical/conformation coupling is essential.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"9 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111082","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}
引用次数: 0
Machine-Learning Ice Spectra: From 1 to 256 Features 机器学习冰光谱:从1到256个特征
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01413
Shokirbek Shermukhamedov,Jolla Kullgren,Daniel Sethio,Kersti Hermansson
The study explores how well machine learning and structural fingerprints can predict spectroscopic properties of ice (OH vibrational frequencies and 1H chemical shifts). A large theoretical data set (55 ice polymorphs, 1010 DFT data points both for the vibrations and for the NMR shifts) and a smaller cross-validation set are employed. The Message Passing Atomic Cluster Expansion (MACE) model performs the best, with high accuracy (root-mean-square deviation, RMSD, of 0.06 ppm for chemical shifts and ∼10 cm–1 for vibrational frequencies). Simpler descriptors like ACSF and SOAP, when paired with suitable regressors, nearly match MACE’s performance. At the other end of the complexity scale, it is found that using the simplest possible physics-based descriptor of the environment (a single H-bond distance) yields RMSD values three times as large for the vibrations and four times as large for the proton chemical shift compared to the MACE model. Depending on the context, those RMSD values may still be considered modest and useful, considering the gain in simplicity and transparency.
该研究探索了机器学习和结构指纹如何很好地预测冰的光谱特性(OH振动频率和1H化学位移)。一个大的理论数据集(55冰多态,1010 DFT数据点的振动和核磁共振位移)和一个较小的交叉验证集被采用。消息传递原子簇扩展(Message Passing Atomic Cluster Expansion, MACE)模型表现最好,精度高(化学位移的均方根偏差RMSD为0.06 ppm,振动频率为~ 10 cm-1)。更简单的描述符,如ACSF和SOAP,当与合适的回归器配对时,几乎可以匹配MACE的性能。在复杂性尺度的另一端,研究人员发现,与MACE模型相比,使用最简单的基于物理的环境描述符(单个氢键距离)产生的RMSD值是振动的三倍,是质子化学位移的四倍。根据上下文,考虑到简单性和透明度的增加,这些RMSD值可能仍然被认为是适度和有用的。
{"title":"Machine-Learning Ice Spectra: From 1 to 256 Features","authors":"Shokirbek Shermukhamedov,Jolla Kullgren,Daniel Sethio,Kersti Hermansson","doi":"10.1021/acs.jctc.5c01413","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01413","url":null,"abstract":"The study explores how well machine learning and structural fingerprints can predict spectroscopic properties of ice (OH vibrational frequencies and 1H chemical shifts). A large theoretical data set (55 ice polymorphs, 1010 DFT data points both for the vibrations and for the NMR shifts) and a smaller cross-validation set are employed. The Message Passing Atomic Cluster Expansion (MACE) model performs the best, with high accuracy (root-mean-square deviation, RMSD, of 0.06 ppm for chemical shifts and ∼10 cm–1 for vibrational frequencies). Simpler descriptors like ACSF and SOAP, when paired with suitable regressors, nearly match MACE’s performance. At the other end of the complexity scale, it is found that using the simplest possible physics-based descriptor of the environment (a single H-bond distance) yields RMSD values three times as large for the vibrations and four times as large for the proton chemical shift compared to the MACE model. Depending on the context, those RMSD values may still be considered modest and useful, considering the gain in simplicity and transparency.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"17 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111079","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}
引用次数: 0
Seniority-Zero Canonical Transformation Theory: Error Reduction via Late Truncation 资历零正则变换理论:通过后期截断减少误差
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01892
Daniel F. Calero-Osorio,Paul W. Ayers
We show how to add the effects of residual electron correlation to a reference seniority-zero wave function by transforming the true electronic Hamiltonian into seniority-zero form. The transformation is treated via the Baker–Campbell–Hausdorff (BCH) expansion, and the seniority-zero structure of the reference is exploited to evaluate the first three commutators exactly; the remaining contributions are handled with a recursive commutator approximation, as is typical in canonical transformation methods. By choosing a seniority-zero reference and using parallel computation, this method is practical for small- to medium-sized systems. Numerical tests show high accuracy, with errors ∼10–4 Hartree.
我们展示了如何通过将真正的电子哈密顿量转换为优先级零形式,将剩余电子相关效应添加到参考优先级零波函数中。通过Baker-Campbell-Hausdorff (BCH)展开对变换进行处理,并利用参考的资历零结构精确计算前三个换向子;其余的贡献用递归换向子近似处理,这是典型的规范转换方法。该方法通过选择一个优先级为零的参考点并采用并行计算的方法,适用于中小型系统。数值测试显示出很高的精度,误差为10-4哈特里。
{"title":"Seniority-Zero Canonical Transformation Theory: Error Reduction via Late Truncation","authors":"Daniel F. Calero-Osorio,Paul W. Ayers","doi":"10.1021/acs.jctc.5c01892","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01892","url":null,"abstract":"We show how to add the effects of residual electron correlation to a reference seniority-zero wave function by transforming the true electronic Hamiltonian into seniority-zero form. The transformation is treated via the Baker–Campbell–Hausdorff (BCH) expansion, and the seniority-zero structure of the reference is exploited to evaluate the first three commutators exactly; the remaining contributions are handled with a recursive commutator approximation, as is typical in canonical transformation methods. By choosing a seniority-zero reference and using parallel computation, this method is practical for small- to medium-sized systems. Numerical tests show high accuracy, with errors ∼10–4 Hartree.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"8 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111081","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}
引用次数: 0
Computing Exchange Coupling Constants in Transition Metal Complexes with Tensor Product Selected Configuration Interaction 用张量积选择构型相互作用计算过渡金属配合物中的交换耦合常数
IF 5.5 1区 化学 Q2 CHEMISTRY, PHYSICAL Pub Date : 2026-02-04 DOI: 10.1021/acs.jctc.5c01817
Arnab Bachhar,Nicholas J. Mayhall
Transition metal complexes present significant challenges for electronic structure theory due to strong electron correlation arising from partially filled d-orbitals. We compare our recently developed Tensor Product Selected Configuration Interaction (TPSCI) with Density Matrix Renormalization Group (DMRG) for computing exchange coupling constants in six transition metal systems, including dinuclear Cr, Fe, and Mn complexes and a tetranuclear Ni-cubane. TPSCI uses a locally correlated tensor product state basis to capture electronic structure efficiently while maintaining interpretability. From calculations on active spaces ranging from (22e,29o) to (42e,49o), we find that TPSCI consistently yields higher variational energies than DMRG due to truncation of local cluster states, but provides magnetic exchange coupling constants (J) generally within 10–30 cm–1 of DMRG results. Key advantages include natural multistate capability enabling direct J extrapolation with smaller statistical errors, and computational efficiency for challenging systems. However, cluster state truncation represents a fundamental limitation requiring careful convergence testing, particularly for large local cluster dimensions. We identify specific failure cases where current truncation schemes break down, highlighting the need for improved cluster state selection methods and distributed memory implementations to realize TPSCI’s full potential for strongly correlated systems.
过渡金属配合物由于部分填满d轨道而产生的强电子相关性对电子结构理论提出了重大挑战。我们将最近开发的张量积选择构型相互作用(TPSCI)与密度矩阵重正化群(DMRG)进行了比较,用于计算六种过渡金属体系的交换耦合常数,包括双核Cr, Fe和Mn配合物和四核Ni-cubane。TPSCI使用局部相关张量积状态基来有效捕获电子结构,同时保持可解释性。从(22e, 290)到(42e, 490)的有效空间计算中,我们发现由于局部簇态的截断,TPSCI始终比DMRG产生更高的变分能量,但提供的磁交换耦合常数(J)通常在DMRG结果的10-30 cm-1范围内。主要优势包括自然的多状态能力,支持以较小的统计误差进行直接J外推,以及具有挑战性系统的计算效率。然而,集群状态截断代表了一个基本限制,需要仔细的收敛测试,特别是对于大的局部集群维度。我们确定了当前截断方案失效的具体故障案例,强调了改进集群状态选择方法和分布式内存实现的必要性,以实现TPSCI在强相关系统中的全部潜力。
{"title":"Computing Exchange Coupling Constants in Transition Metal Complexes with Tensor Product Selected Configuration Interaction","authors":"Arnab Bachhar,Nicholas J. Mayhall","doi":"10.1021/acs.jctc.5c01817","DOIUrl":"https://doi.org/10.1021/acs.jctc.5c01817","url":null,"abstract":"Transition metal complexes present significant challenges for electronic structure theory due to strong electron correlation arising from partially filled d-orbitals. We compare our recently developed Tensor Product Selected Configuration Interaction (TPSCI) with Density Matrix Renormalization Group (DMRG) for computing exchange coupling constants in six transition metal systems, including dinuclear Cr, Fe, and Mn complexes and a tetranuclear Ni-cubane. TPSCI uses a locally correlated tensor product state basis to capture electronic structure efficiently while maintaining interpretability. From calculations on active spaces ranging from (22e,29o) to (42e,49o), we find that TPSCI consistently yields higher variational energies than DMRG due to truncation of local cluster states, but provides magnetic exchange coupling constants (J) generally within 10–30 cm–1 of DMRG results. Key advantages include natural multistate capability enabling direct J extrapolation with smaller statistical errors, and computational efficiency for challenging systems. However, cluster state truncation represents a fundamental limitation requiring careful convergence testing, particularly for large local cluster dimensions. We identify specific failure cases where current truncation schemes break down, highlighting the need for improved cluster state selection methods and distributed memory implementations to realize TPSCI’s full potential for strongly correlated systems.","PeriodicalId":45,"journal":{"name":"Journal of Chemical Theory and Computation","volume":"91 1","pages":""},"PeriodicalIF":5.5,"publicationDate":"2026-02-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146111080","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}
引用次数: 0
期刊
Journal of Chemical Theory and Computation
全部 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