libterachem.py: Modular development of electronic structure methods with graphical processing units
Abstract
We present libterachem.py, a Python interface and method framework that facilitates development and prototyping of electronic structure methods while leveraging TeraChem's optimized graphical processing unit (GPU) operations. Direct access to these building blocks and their intermediate quantities allows users to assemble and modify methods within Python while reusing the underlying GPU implementations. A layered architecture, programmable at both the operation level and the method level, supports numerical experimentation, composition of complete electronic structure methods with interchangeable components, and integration with external scientific libraries. We illustrate operation-level use with grid pruning for least-squares tensor hypercontraction and with exchange--correlation functionals supplied through LibXC and by the machine-learned Skala functional. We illustrate method-level use with dispersion corrections, periodic calculations, and geometry optimization with an external optimizer. Python-driven self-consistent-field calculations using libterachem.py incur no significant computational overhead compared to native TeraChem and preserve its multi-GPU scaling. A photoactive yellow protein QM/MM molecular dynamics application further demonstrates throughput comparable to native TeraChem in a workflow composed with external simulation libraries. These results establish libterachem.py as a platform for method development and interoperable molecular simulation that retains the performance of TeraChem's GPU operations.
Related
- VIBE-QC: an open-source, general-purpose electronic-structure code for molecules and solids developed with generative AI
- Electronic Descriptors Lower Oxygen-Vacancy Formation-Energy Error Relative to Crystal Structure Alone in Non-Magnetic Oxides
- Electronic Structure, Bonding Topology, and Photophysical Properties of the Low-Valent Beryllium Cluster CpBe2NO: A Multireference and Real-Space Wavefunction Analysis
- Development of surface modified Mg-2Gd-2Yb alloy through friction stir processing: corrosion behaviour and short-term anodic performance
- Synthetic Mechanism of a Fe(II) N -Heterocyclic Carbene Bidentate Complex Revealed by Electronic Structure Methods
BibTeX
@misc{zhang2026libterachempy,
title = {{libterachem.py: Modular development of electronic structure methods with graphical processing units}},
author = {Juncheng Harry Zhang and Ruiyan Wang and Andreas Erbs Hillers-Bendtsen and Todd J. Martínez},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.02540},
eprint = {2610.02540},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.02540},
}
RIS
TY - UNPB
TI - libterachem.py: Modular development of electronic structure methods with graphical processing units
AU - Juncheng Harry Zhang
AU - Ruiyan Wang
AU - Andreas Erbs Hillers-Bendtsen
AU - Todd J. Martínez
PY - 2026
DA - 2026/10/01
PB - arXiv
DO - 10.48550/arxiv.2610.02540
UR - https://arxiv.org/abs/2610.02540
AB - We present libterachem.py, a Python interface and method framework that facilitates development and prototyping of electronic structure methods while leveraging TeraChem's optimized graphical processing unit (GPU) operations. Direct access to these building blocks and their intermediate quantities allows users to assemble and modify methods within Python while reusing the underlying GPU implementations. A layered architecture, programmable at both the operation level and the method level, supports numerical experimentation, composition of complete electronic structure methods with interchangeable components, and integration with external scientific libraries. We illustrate operation-level use with grid pruning for least-squares tensor hypercontraction and with exchange--correlation functionals supplied through LibXC and by the machine-learned Skala functional. We illustrate method-level use with dispersion corrections, periodic calculations, and geometry optimization with an external optimizer. Python-driven self-consistent-field calculations using libterachem.py incur no significant computational overhead compared to native TeraChem and preserve its multi-GPU scaling. A photoactive yellow protein QM/MM molecular dynamics application further demonstrates throughput comparable to native TeraChem in a workflow composed with external simulation libraries. These results establish libterachem.py as a platform for method development and interoperable molecular simulation that retains the performance of TeraChem's GPU operations.
ER -