-
2500 vibronic eigenstates of the NO$_3$ radical
Authors:
Henrik R. Larsson,
Alexandra Viel
Abstract:
The nitrate radical NO$_3$ plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic $\tilde X ^2A_2'$ state. Using recently developed tensor network state…
▽ More
The nitrate radical NO$_3$ plays an important role in atmospheric chemistry, yet many aspects of its coupled and anharmonic vibronic structure remain elusive. Here, using an accurate, coupled full-dimensional diabatic potential that includes five electronic states, we revisit the vibronic spectrum associated with the electronic $\tilde X ^2A_2'$ state. Using recently developed tensor network state methods, we are able to compute more than 2500 vibronic states, thereby increasing the number of computed full-dimensional states by a factor of 50, compared to previous work. While we obtain good agreement with experiment for most of the assigned vibronic levels, for several others, we observe striking disagreement. Further, for the antisymmetric bending motion we find remarkably large symmetry-induced level splittings that are larger than the zero-order reference. We discuss non-negligible nonadiabatic effects and show that the Born-Oppenheimer approximation leads to significant errors in the spectrum.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
A tensor network view of multilayer multiconfiguration time-dependent Hartree methods
Authors:
Henrik R. Larsson
Abstract:
The multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization group (DMRG) are powerful workhorses applied mostly in different scientific fields. Although both methods are based on tensor network states, very different mathematical languages are used for describing them. This severely limits knowledge transfer and sometimes leads to re-inventions…
▽ More
The multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method and the density matrix renormalization group (DMRG) are powerful workhorses applied mostly in different scientific fields. Although both methods are based on tensor network states, very different mathematical languages are used for describing them. This severely limits knowledge transfer and sometimes leads to re-inventions of ideas well known in the other field. Here, we review ML-MCTDH and DMRG theory using both MCTDH expressions and tensor network diagrams. We derive the ML-MCTDH equations of motions using diagrams and compare them with time-dependent and time-independent DMRG algorithms. We further review two selected recent advancements. The first advancement is related to optimizing unoccupied single-particle functions in MCTDH, which corresponds to subspace enrichment in the DMRG. The second advancement is related to finding optimal tree structures and on highlighting similarities and differences of tensor networks used in MCTDH and DMRG theories. We hope that this contribution will foster more fruitful cross-fertilization of ideas between ML-MCTDH and DMRG.
△ Less
Submitted 12 February, 2024; v1 submitted 20 November, 2023;
originally announced November 2023.
-
Block2: a comprehensive open source framework to develop and apply state-of-the-art DMRG algorithms in electronic structure and beyond
Authors:
Huanchen Zhai,
Henrik R. Larsson,
Seunghoon Lee,
Zhi-Hao Cui,
Tianyu Zhu,
Chong Sun,
Linqing Peng,
Ruojing Peng,
Ke Liao,
Johannes Tölle,
Junjie Yang,
Shuoxue Li,
Garnet Kin-Lic Chan
Abstract:
Block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matr…
▽ More
Block2 is an open source framework to implement and perform density matrix renormalization group and matrix product state algorithms. Out-of-the-box it supports the eigenstate, time-dependent, response, and finite-temperature algorithms. In addition, it carries special optimizations for ab initio electronic structure Hamiltonians and implements many quantum chemistry extensions to the density matrix renormalization group, such as dynamical correlation theories. The code is designed with an emphasis on flexibility, extensibility, and efficiency, and to support integration with external numerical packages. Here we explain the design principles and currently supported features and present numerical examples in a range of applications.
△ Less
Submitted 24 November, 2023; v1 submitted 5 October, 2023;
originally announced October 2023.
-
State-resolved infrared spectrum of the protonated water dimer: Revisiting the characteristic proton transfer doublet peak
Authors:
Henrik R. Larsson,
Markus Schröder,
Richard Beckmann,
Fabien Brieuc,
Christoph Schran,
Dominik Marx,
Oriol Vendrell
Abstract:
The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water di…
▽ More
The infrared (IR) spectra of protonated water clusters encode precise information on the dynamics and structure of the hydrated proton. However, the strong anharmonic coupling and quantum effects of these elusive species remain puzzling up to the present day. Here, we report unequivocal evidence that the interplay between the proton transfer and the water wagging motions in the protonated water dimer (Zundel ion) giving rise to the characteristic doublet peak is both more complex and more sensitive to subtle energetic changes than previously thought. In particular, hitherto overlooked low-intensity satellite peaks in the experimental spectrum are now unveiled and mechanistically assigned. Our findings rely on the comparison of IR spectra obtained using two highly accurate potential energy surfaces in conjunction with highly accurate state-resolved quantum simulations. We demonstrate that these high-accuracy simulations are important for providing definite assignments of the complex IR signals of fluxional molecules.
△ Less
Submitted 5 September, 2022; v1 submitted 23 June, 2022;
originally announced June 2022.
-
The chromium dimer: closing a chapter of quantum chemistry
Authors:
Henrik R. Larsson,
Huanchen Zhai,
C. J. Umrigar,
Garnet Kin-Lic Chan
Abstract:
The complex electronic structure and unusual potential energy curve of the chromium dimer have fascinated scientists for decades, with agreement between theory and experiment so far elusive. Here, we present a new ab initio simulation of the potential energy curve and vibrational spectrum that significantly improves on all earlier estimates. Our data support a shift in earlier experimental assignm…
▽ More
The complex electronic structure and unusual potential energy curve of the chromium dimer have fascinated scientists for decades, with agreement between theory and experiment so far elusive. Here, we present a new ab initio simulation of the potential energy curve and vibrational spectrum that significantly improves on all earlier estimates. Our data support a shift in earlier experimental assignments of a cluster of vibrational frequencies by one quantum number. The new vibrational assignment yields an experimentally derived potential energy curve in quantitative agreement with theory across all bond lengths and across all measured frequencies. By solving this long-standing problem, our results raise the possibility of quantitative quantum chemical modeling of transition metal clusters with spectroscopic accuracy.
△ Less
Submitted 24 August, 2022; v1 submitted 21 June, 2022;
originally announced June 2022.
-
Matrix product states with large sites
Authors:
Henrik R. Larsson,
Huanchen Zhai,
Klaas Gunst,
Garnet Kin-Lic Chan
Abstract:
We explore various ways to group orbitals into clusters in a matrix product state (MPS). We explain how a generic cluster MPS can often lead to an increase in computational cost and instead propose a special cluster structure, involving only the first and last orbitals/sites, with a wider scope for computational advantage. This structure is a natural formalism to describe correlated multireference…
▽ More
We explore various ways to group orbitals into clusters in a matrix product state (MPS). We explain how a generic cluster MPS can often lead to an increase in computational cost and instead propose a special cluster structure, involving only the first and last orbitals/sites, with a wider scope for computational advantage. This structure is a natural formalism to describe correlated multireference (MR) theories. We demonstrate the flexibility and usefulness of this approach by implementing various uncontracted MR configuration interaction, perturbation and linearized coupled cluster theories using an MPS with large cluster sites. Applications to the nitrogen dimer, the chromium dimer, and benzene, including up to triple excitations in the external space, demonstrate the utility of an MPS with up to two large sites. We use our results to analyze the quality of different multireference approximations.
△ Less
Submitted 21 January, 2022; v1 submitted 22 September, 2021;
originally announced September 2021.
-
Control of concerted back-to-back double ionization dynamics in helium
Authors:
Henrik R. Larsson,
David J. Tannor
Abstract:
Double ionization (DI) is a fundamental process that despite its apparent simplicity provides rich opportunities for probing and controlling the electronic motion. Even for the simplest multielectron atom, helium, new DI mechanisms are still being found. To first order in the field strength, a strong external field doubly ionizes the electrons in helium such that they are ejected into the same dir…
▽ More
Double ionization (DI) is a fundamental process that despite its apparent simplicity provides rich opportunities for probing and controlling the electronic motion. Even for the simplest multielectron atom, helium, new DI mechanisms are still being found. To first order in the field strength, a strong external field doubly ionizes the electrons in helium such that they are ejected into the same direction (front-to-back motion). The ejection into opposite directions (back-to-back motion) cannot be described to first order, making it a challenging target for control. Here, we address this challenge and optimize the field with the objective of back-to-back double ionization using a (1 + 1)-dimensional model. The optimization is performed using four different control procedures: (1) short-time control, (2) derivative-free optimization of basis expansions of the field, (3) the Krotov method, and (4) control of the classical equations of motion. All four procedures lead to fields with dominant back-to-back motion. All the fields obtained exploit essentially the same two-step mechanism leading to back-to-back motion: first, the electrons are displaced by the field into the same direction. Second, after the field turns off, the nuclear attraction and the electron-electron repulsion combine to generate the final motion into opposite directions for each electron. By performing quasi-classical calculations, we confirm that this mechanism is essentially classical.
△ Less
Submitted 11 October, 2021; v1 submitted 21 May, 2020;
originally announced May 2020.
-
Minimal matrix product states and generalizations of mean-field and geminal wavefunctions
Authors:
Henrik R. Larsson,
Carlos A. Jiménez-Hoyos,
Garnet Kin-Lic Chan
Abstract:
Simple wavefunctions of low computational cost but which can achieve qualitative accuracy across the whole potential energy surface (PES) are of relevance to many areas of electronic structure as well as to applications to dynamics. Here, we explore a class of simple wavefunctions, the minimal matrix product state (MMPS), that generalizes many simple wavefunctions in common use, such as projected…
▽ More
Simple wavefunctions of low computational cost but which can achieve qualitative accuracy across the whole potential energy surface (PES) are of relevance to many areas of electronic structure as well as to applications to dynamics. Here, we explore a class of simple wavefunctions, the minimal matrix product state (MMPS), that generalizes many simple wavefunctions in common use, such as projected mean-field wavefunctions, geminal wavefunctions, and generalized valence bond states. By examining the performance of MMPSs for PESs of some prototypical systems, we find that they yield good qualitative behavior across the whole PES, often significantly improving on the aforementioned ansätze.
△ Less
Submitted 26 June, 2020; v1 submitted 7 May, 2020;
originally announced May 2020.
-
Computing vibrational eigenstates with tree tensor network states (TTNS)
Authors:
Henrik R. Larsson
Abstract:
We present how to compute vibrational eigenstates with tree tensor network states (TTNSs), the underlying ansatz behind the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. The eigenstates are computed with an algorithm that is based on the density matrix renormalization group (DMRG). We apply this to compute the vibrational spectrum of acetonitrile (CH$_3$CN) to high accura…
▽ More
We present how to compute vibrational eigenstates with tree tensor network states (TTNSs), the underlying ansatz behind the multilayer multiconfiguration time-dependent Hartree (ML-MCTDH) method. The eigenstates are computed with an algorithm that is based on the density matrix renormalization group (DMRG). We apply this to compute the vibrational spectrum of acetonitrile (CH$_3$CN) to high accuracy and compare TTNSs with matrix product states (MPSs), the ansatz behind the DMRG. The presented optimization scheme converges much faster than ML-MCTDH-based optimization. For this particular system, we found no major advantage of the more general TTNS over MPS. We highlight that for both TTNS and MPS, the usage of an adaptive bond dimension significantly reduces the amount of required parameters. We furthermore propose a procedure to find good trees.
△ Less
Submitted 26 November, 2019; v1 submitted 30 September, 2019;
originally announced September 2019.
-
Global optimization of parameters in the reactive force field ReaxFF for SiOH
Authors:
H. R. Larsson,
A. C. T. van Duin,
B. Hartke
Abstract:
We have used unbiased global optimization to fit a reactive force field to a given set of reference data. Specifically, we have employed genetic algorithms (GA) to fit ReaxFF to SiOH data, using an in-house GA code that is parallelized across reference data items via the message-passing interface (MPI). Details of GA tuning turn out to be far less important for global optimization efficiency than…
▽ More
We have used unbiased global optimization to fit a reactive force field to a given set of reference data. Specifically, we have employed genetic algorithms (GA) to fit ReaxFF to SiOH data, using an in-house GA code that is parallelized across reference data items via the message-passing interface (MPI). Details of GA tuning turn out to be far less important for global optimization efficiency than using suitable ranges within which the parameters are varied. To establish these ranges, either prior knowledge can be used or successive stages of GA optimizations, each building upon the best parameter vectors and ranges found in the previous stage. We finally arrive at optimized force fields with smaller error measures than those published previously. Hence, this optimization approach will contribute to converting force-field fitting from a specialist task to an everyday commodity, even for the more difficult case of reactive force fields.
△ Less
Submitted 15 September, 2019;
originally announced September 2019.
-
Resonance dynamics of DCO ($\widetilde{X}\,{}^2A'$) simulated with the dynamically pruned discrete variable representation (DP-DVR)
Authors:
Henrik R. Larsson,
Jens Riedel,
Jie Wei,
Friedrich Temps,
Bernd Hartke
Abstract:
Selected resonance states of the deuterated formyl radical in the electronic ground state ($\widetilde{X}\,{}^2A'$) are computed using our recently introduced dynamically pruned discrete variable representation (DP-DVR) [H. R. Larsson, B. Hartke and D. J. Tannor, J. Chem. Phys., 145, 204108 (2016)]. Their decay and asymptotic distributions are analyzed and, for selected resonances, compared to exp…
▽ More
Selected resonance states of the deuterated formyl radical in the electronic ground state ($\widetilde{X}\,{}^2A'$) are computed using our recently introduced dynamically pruned discrete variable representation (DP-DVR) [H. R. Larsson, B. Hartke and D. J. Tannor, J. Chem. Phys., 145, 204108 (2016)]. Their decay and asymptotic distributions are analyzed and, for selected resonances, compared to experimental results obtained by a combination of stimulated emission pumping (SEP) and velocity-map imaging of the product D atoms. The theoretical results show good agreement with the experimental kinetic energy distributions. The intramolecular vibrational energy redistribution (IVR) is analyzed and compared with previous results from an effective polyad Hamiltonian. Specifically, we analyzed the part of the wavefunction that remains in the interaction region during the decay. The results from the polyad Hamiltonian could mainly be confirmed. The C=O stretch quantum number is typically conserved, while the D-C=O bend quantum number decreases. Differences are due to strong anharmonic coupling such that all resonances have major contributions from several zero-order states. For some of the resonances, the coupling is so strong that no further zero-order states appear during the dynamics in the interaction region, even after propagating for 300 ps.
△ Less
Submitted 8 May, 2018; v1 submitted 20 February, 2018;
originally announced February 2018.
-
Dynamical pruning of the multiconfiguration time-dependent Hartree method (DP-MCTDH): An efficient approach for multidimensional quantum dynamics
Authors:
Henrik R. Larsson,
David J. Tannor
Abstract:
We present two strategies for combining dynamical pruning with the multiconfiguration time-dependent Hartree method (DP-MCTDH), where dynamical pruning means on-the-fly selection of relevant basis functions. The first strategy prunes the primitive basis that represents the single-particle functions (SPFs). This is useful for smaller systems that require many primitive basis functions per degree of…
▽ More
We present two strategies for combining dynamical pruning with the multiconfiguration time-dependent Hartree method (DP-MCTDH), where dynamical pruning means on-the-fly selection of relevant basis functions. The first strategy prunes the primitive basis that represents the single-particle functions (SPFs). This is useful for smaller systems that require many primitive basis functions per degree of freedom, as we will illustrate for NO$_2$. Furthermore, this allows for higher-dimensional mode combination and partially lifts the sum-of-product-form requirement onto the structure of the Hamiltonian, as we illustrate for nonadiabatic 24-dimensional pyrazine. The second strategy prunes the set of configurations of SPF at each time step. We show that this strategy yields significant speed-ups with factors between 5 and 50 in computing time, making it competitive with the multilayer MCTDH method.
△ Less
Submitted 17 July, 2017; v1 submitted 8 May, 2017;
originally announced May 2017.
-
Efficient molecular quantum dynamics in coordinate and phase space using pruned bases
Authors:
Henrik R. Larsson,
Bernd Hartke,
David J. Tannor
Abstract:
We present an efficient implementation of dynamically pruned quantum dynamics, both in coordinate space and in phase space. We combine the ideas behind the biorthogonal von Neumann basis (PvB) with the orthogonalized momentum-symmetrized Gaussians (Weylets) to create a new basis, projected Weylets, that takes the best from both methods. We benchmark pruned dynamics using phase-space-localized PvB,…
▽ More
We present an efficient implementation of dynamically pruned quantum dynamics, both in coordinate space and in phase space. We combine the ideas behind the biorthogonal von Neumann basis (PvB) with the orthogonalized momentum-symmetrized Gaussians (Weylets) to create a new basis, projected Weylets, that takes the best from both methods. We benchmark pruned dynamics using phase-space-localized PvB, projected Weylets, and coordinate-space-localized DVR bases, with real-world examples in up to six dimensions. We show that coordinate-space localization is most important for efficient pruning and that pruned dynamics is much faster compared to unpruned, exact dynamics. Phase-space localization is useful for more demanding dynamics where many basis functions are required. There, projected Weylets offer a more compact representation than pruned DVR bases.
△ Less
Submitted 15 November, 2016; v1 submitted 13 June, 2016;
originally announced June 2016.
-
Correlation effects in strong-field ionization of heteronuclear diatomic molecules
Authors:
Henrik R. Larsson,
Sebastian Bauch,
Lasse Kragh Sørensen,
Michael Bonitz
Abstract:
We develop a time-dependent theory to investigate electron dynamics and photoionization processes of diatomic molecules interacting with strong laser fields including electron-electron correlation effects. We combine the recently formulated time-dependent generalized-active-space configuration interaction theory [D. Hochstuhl and M. Bonitz, Phys. Rev. A 86, 053424 (2012); S. Bauch, et al., Phys. R…
▽ More
We develop a time-dependent theory to investigate electron dynamics and photoionization processes of diatomic molecules interacting with strong laser fields including electron-electron correlation effects. We combine the recently formulated time-dependent generalized-active-space configuration interaction theory [D. Hochstuhl and M. Bonitz, Phys. Rev. A 86, 053424 (2012); S. Bauch, et al., Phys. Rev. A 90, 062508 (2014)] with a prolate spheroidal basis set including localized orbitals and continuum states to describe the bound electrons and the outgoing photoelectron. As an example, we study the strong-field ionization of the two-center four-electron lithium hydride molecule in different intensity regimes. By using single-cycle pulses, two orientations of the asymmetric heteronuclear molecule are investigated: Li-H, with the electrical field pointing from H to Li, and the opposite case of H-Li. The preferred orientation for ionization is determined and we find a transition from H-Li, for low intensity, to Li-H, for high intensity. The influence of electron correlations is studied at different levels of approximation, and we find a significant change in the preferred orientation. For certain intensity regimes, even an interchange of the preferred configuration is observed, relative to the uncorrelated simulations. Further insight is provided by detailed comparisons of photoelectron angular distributions with and without correlation effects taken into account.
△ Less
Submitted 1 February, 2016; v1 submitted 15 July, 2015;
originally announced July 2015.