Vincent Pohl

ORBKIT: A modular python toolbox for cross‐platform postprocessing of quantum chemical wavefunction data

ORBKIT is a toolbox for postprocessing electronic structure calculations based on a highly modular and portable Python architecture. The program allows computing a multitude of electronic properties of molecular systems on arbitrary spatial grids from the basis set representation of its electronic wavefunction, as well as several grid‐independent properties. The required data can be extracted directly from the standard output of a large number of quantum chemistry programs. ORBKIT can be used as a standalone program to determine standard quantities, for example, the electron density, molecular orbitals, and derivatives thereof. The cornerstone of ORBKIT is its modular structure. The existing basic functions can be arranged in an individual way and can be easily extended by user‐written modules to determine any other derived quantity. ORBKIT offers multiple output formats that can be processed by common visualization tools (VMD, Molden, etc.). Additionally, ORBKIT possesses routines to order molecular orbitals computed at different nuclear configurations according to their electronic character and to interpolate the wavefunction between these configurations. The program is open‐source under GNU‐LGPLv3 license and freely available at https://github.com/orbkit/orbkit/. This article provides an overview of ORBKIT with particular focus on its capabilities and applicability, and includes several example calculations.

Multidirectional Angular Electronic Flux during Adiabatic Attosecond Charge Migration in Excited Benzene

Recently, adiabatic attosecond charge migration (AACM) has been monitored and simulated for the first time, with application to the oriented iodoacetylene cation where AACM starts from the initial superposition of the ground state (φ0) and an electronic excited state (φ1). Here, we develop the theory for electronic fluxes during AACM in ring-shaped molecules, with application to oriented benzene prepared in the superposition of the ground and first excited singlet states. The initial state and its time evolution are analogous to coherent tunneling where φ0 and φ1 have different meanings; however, they denote the wave functions of the lowest tunneling doublet. This analogy suggests to transfer the theory of electronic fluxes during coherent tunneling to AACM, with suitable modifications which account for (i) the different time scales and (ii) the different electronic states, and which make use of (iii) the preparation of the initial state for AACM by a linearly polarized laser pulse. Application to benzene yields the multidirectional angular electronic flux with a pincer-motion type pattern during AACM: this unequivocal result confirms a previous working hypothesis. Moreover, the theory of AACM allows quantification of the electronic flux; that is, the maximum number of electrons (out of 42) which flow concertedly during AACM in benzene is 6 × 0.08 = 0.48.

Computation of the Electronic Flux Density in the Born−Oppenheimer Approximation

A molecule in the electronic ground state described in the Born–Oppenheimer approximation (BOA) by the wave function ΨBO = Φ0χ0 (where Φ0 is the time-independent electronic energy eigenfunction and χ0 is a time-dependent nuclear wave packet) exhibits a nonzero nuclear flux density, whereas it always displays zero electronic flux density (EFD), because the electrons are in a stationary state. A hierarchical approach to the computation of the EFD within the context of the BOA, which utilizes only standard techniques of quantum chemistry (to obtain Φ0) and quantum dynamics (to describe the evolution of χ0 on the ground-state potential energy surface), provides a resolution of this puzzling, nonintuitive result. The procedure is applied to H2(+) oriented parallel with the z-axis and vibrating in the ground state (2)Σg(+). First, Φ0 and χ0 are combined by the coupled-channels technique to give the normally dominant z-component of the EFD. Imposition of the constraints of electronic continuity, cylindrical symmetry of Φ0 and two boundary conditions on the EFD through a scaling procedure yields an improved z-component, which is then used to compute the complementary orthogonal ρ-component. The resulting EFD agrees with its highly accurate counterpart furnished by a non-BOA treatment of the system.

Adiabatic electronic flux density: A Born-Oppenheimer broken-symmetry ansatz

The Born-Oppenheimer approximation leads to the counterintuitive result of a vanishing electronic flux density upon vibrational dynamics in the electronic ground state. To circumvent this long known issue, we propose using pairwise anti-symmetrically translated vibronic densities to generate a symmetric electronic density that can be forced to satisfy the continuity equation approximately. The so-called Born-Oppenheimer broken symmetry ansatz yields all components of the flux density simultaneously while requiring only knowledge about the nuclear quantum dynamics on the electronic adiabatic ground state potential energy surface. The underlying minimization procedure is transparent and computationally inexpensive, and the solution can be computed from the standard output of any quantum chemistry program. Taylor series expansion reveals that the implicit electron dynamics originates from non-adiabatic coupling to the explicit Born-Oppenheimer nuclear dynamics. The new approach is applied to the

An open-source framework for analyzing N-electron dynamics. I. Multideterminantal wave functions

{\rm H}_2^+

An open-source framework for analyzing N-electron dynamics. II. Hybrid density functional theory/configuration interaction methodology

molecular ion vibrating in its

Cyanographone and isocyanographone — Two asymmetrically functionalized graphene pseudohalides and their potential use in chemical sensing

{}^2\Sigma^+_g

Quantum control of electronic fluxes during adiabatic attosecond charge migration in degenerate superposition states of benzene

ground state. The electronic flux density is found to have the correct nodal structure and symmetry properties at all times.

Electronic and nuclear flux densities in the H 2 molecule

The aim of the present contribution is to provide a framework for analyzing and visualizing the correlated many‐electron dynamics of molecular systems, where an explicitly time‐dependent electronic wave packet is represented as a linear combination of N‐electron wave functions. The central quantity of interest is the electronic flux density, which contains all information about the transient electronic density, the associated phase, and their temporal evolution. It is computed from the associated one‐electron operator by reducing the multideterminantal, many‐electron wave packet using the Slater‐Condon rules. Here, we introduce a general tool for post‐processing multideterminant configuration‐interaction wave functions obtained at various levels of theory. It is tailored to extract directly the data from the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions. The procedure is implemented in the open‐source Python program detCI@ORBKIT, which shares and builds on the modular design of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). The new procedure is applied to ultrafast charge migration processes in different molecular systems, demonstrating its broad applicability. Convergence of the N‐electron dynamics with respect to the electronic structure theory level and basis set size is investigated. This provides an assessment of the robustness of qualitative and quantitative statements that can be made concerning dynamical features observed in charge migration simulations.

Attosecond angular flux of partial charges on the carbon atoms of benzene in non-aromatic excited state

In this contribution, we extend our framework for analyzing and visualizing correlated many‐electron dynamics to non‐variational, highly scalable electronic structure method. Specifically, an explicitly time‐dependent electronic wave packet is written as a linear combination of N‐electron wave functions at the configuration interaction singles (CIS) level, which are obtained from a reference time‐dependent density functional theory (TDDFT) calculation. The procedure is implemented in the open‐source Python program detCI@ORBKIT, which extends the capabilities of our recently published post‐processing toolbox (Hermann et al., J. Comput. Chem. 2016, 37, 1511). From the output of standard quantum chemistry packages using atom‐centered Gaussian‐type basis functions, the framework exploits the multideterminental structure of the hybrid TDDFT/CIS wave packet to compute fundamental one‐electron quantities such as difference electronic densities, transient electronic flux densities, and transition dipole moments. The hybrid scheme is benchmarked against wave function data for the laser‐driven state selective excitation in LiH. It is shown that all features of the electron dynamics are in good quantitative agreement with the higher‐level method provided a judicious choice of functional is made. Broadband excitation of a medium‐sized organic chromophore further demonstrates the scalability of the method. In addition, the time‐dependent flux densities unravel the mechanistic details of the simulated charge migration process at a glance.