Abstract:Developing full-dimensional machine-learned potentials with the current gold-standard coupled-cluster (CC) level is a challenging already for medium-sized molecules due to the high computational cost. Consequently, researchers are often bound to use lower-level electronic structure methods such as density functional theory or second-order Moller-Plesset perturbation theory (MP2). Here, we demonstrate on a representative example that gold-standard potentials can now be effectively constructed for molecules of 15 atoms using off-the-shelf hardware. This is achieved by accelerating the CCSD(T) computations via the accurate and cost-effective frozen natural orbital (FNO) approach. The Delta-machine learning (Delta-ML) approach is employed with the use of permutationally invariant polynomials to fit a full-dimensional PES of the acetylacetone molecule, but any other effective descriptor and ML approach can similarly benefit from the accelerated data generation proposed here. Our benchmarks for the global minima, H-transfer TS, and many high-lying configurations show the excellent agreement of FNO-CCSD(T) results with conventional CCSD(T) while achieving a significant time advantage of about a factor of 30-40. The obtained Delta-ML PES shows high fidelity from multiple perspectives including energetic, structural, and vibrational properties. We obtain the symmetric double well H-transfer barrier of 3.15 kcal/mol in excellent agreement with the direct FNO-CCSD(T) barrier of 3.11 kcal/mol as well as with the benchmark CCSD(F12*)(T+)/CBS value of 3.21 kcal/mol. Furthermore, the tunneling splitting due to H-atom transfer is calculated using a 1D double-well potential, providing improved estimates over previous ones obtained using an MP2-based PES. The methodology introduced here represents a significant advancement in the efficient and precise construction of potentials at the CCSD(T) level for molecules above the current limit of 15 atoms.

Eigenstate calculation in the state-averaged (multi-layer) multi-configurational time-dependent Hartree approach

A non-hierarchical multi-layer multi-configurational time-dependent Hartree approach for quantum dynamics on general potential energy surfaces

Computing vibrational eigenstates with tree tensor network states (TTNS)

A multi-layer multi-configurational time-dependent Hartree approach to lattice models beyond one dimension

A New Fragment-Based Approach For Calculating Electronic Excitation Energies Of Large Systems

Calculations of Excited Electronic States by Converging on Saddle Points Using Generalized Mode Following

An Efficient Step-Forward Way To Solve The Schrodinger Eigenvalue Equation In Self-Consistent Calculations

Using a pruned basis and a sparse collocation grid with more points than basis functions to do efficient and accurate MCTDH calculations with general potential energy surfaces

Optimization of highly excited matrix product states with an application to vibrational spectroscopy

Addressing anharmonic effects with density-fitted multicomponent density functional theory

Combining state-of-the-art quantum chemistry and machine learning make gold standard potential energy surfaces accessible for medium-sized molecules

Turbo charging time-dependent density-functional theory with Lanczos chains

Molecule-optimized Basis Sets and Hamiltonians for Accelerated Electronic Structure Calculations of Atoms and Molecules

Efficient and systematic calculation of arbitrary observables for the matrix product state excitation ansatz

Efficient Implementation of Local Excitation Approximation for Treating Excited States of Molecules in Condensed Phase

Linear Scaling Calculations of Excitation Energies with Active-Space Particle-Particle Random Phase Approximation

Improved on-the-Fly MCTDH Simulations with Many-Body-Potential Tensor Decomposition and Projection Diabatization

A trust-region augmented Hessian implementation for state-specific and state-averaged CASSCF wave functions

Excited States via Coupled Cluster Theory without Equation-of-Motion Methods: Seeking Higher Roots with Application to Doubly Excited States and Double Core Hole States