We compute the entanglement between the electronic and vibrational motions in the simplest molecular system, the hydrogen molecular ion, considering the molecule as a bipartite system, electron and vibrational motion. For that purpose we compute an accurate total non-Born-Oppenheimer wave function in terms of a huge expansion using nonorthogonal B-spline basis sets that expand separately the electronic and nuclear wave functions. According to the Schmidt decomposition theorem for bipartite systems, widely used in quantum-information theory, it is possible to find a much shorter but equivalent expansion in terms of the natural orbitals or Schmidt bases for the electronic and nuclear half spaces. Here we extend the Schmidt decomposition theorem to the case in which nonorthogonal bases are used to span the partitioned Hilbert spaces. This extension is first illustrated with two simple coupled systems, the former without an exact solution and the latter exactly solvable. In these model systems of distinguishable coupled particles it is shown that the entanglement content does not increase monotonically with the excitation energy, but only within the manifold of states that belong to an existing excitation mode, if any. In the hydrogen molecular ion the entanglement content for each non-Born-Oppenheimer vibronic state is quantified through the von Neumann and linear entropies and we show that entanglement serves as a witness to distinguish vibronic states related to different Born-Oppenheimer molecular energy curves or electronic excitation modes.
We compute the entanglement between the electronic and vibrational motions in the simplest molecular system, the hydrogen molecular ion, considering the molecule as a bipartite system, electron and vibrational motion. For that purpose we compute an accurate total non-Born-Oppenheimer wave function in terms of a huge expansion using nonorthogonal B-spline basis sets that expand separately the electronic and nuclear wave functions. According to the Schmidt decomposition theorem for bipartite systems, widely used in quantum-information theory, it is possible to find a much shorter but equivalent expansion in terms of the natural orbitals or Schmidt bases for the electronic and nuclear half spaces. Here we extend the Schmidt decomposition theorem to the case in which nonorthogonal bases are used to span the partitioned Hilbert spaces. This extension is first illustrated with two simple coupled systems, the former without an exact solution and the latter exactly solvable. In these model systems of distinguishable coupled particles it is shown that the entanglement content does not increase monotonically with the excitation energy, but only within the manifold of states that belong to an existing excitation mode, if any. In the hydrogen molecular ion the entanglement content for each non-Born-Oppenheimer vibronic state is quantified through the von Neumann and linear entropies and we show that entanglement serves as a witness to distinguish vibronic states related to different Born-Oppenheimer molecular energy curves or electronic excitation modes.
ABSTRACT: In the present work, we have computed the entanglement between the electronic and nuclear motions in two molecular model systems: the one-dimensional hydrogen molecular ion (H2+) and the Shin-Metiu model, considering the molecules as a bipartite systems: electron and nuclear motion. For that purpose, we have computed the Born-Oppenheimer and non-Born-Oppenheimer (Born-Huang) wave function in terms of the Fourier Grid Hamiltonian basis that expands both the electronic and nuclear wave functions. Also, according to the Schmidt decomposition theorem for bipartite systems, widely used in quantum-information theory, there is a much shorter but equivalent expansion in terms of the Schmidt bases for the electronic and nuclear sub-spaces. In these models of distinguishable coupled particles we have shown that the entanglements contents do not increase monotonically with the excitation energy. In the hydrogen molecular ion and in the ShinMetiu model, the entanglements contents for each Born-Oppenheimer electro-nuclear state is quantified through the von-Neumann and linear entropies and we have shown that entanglement serves as a witness of distinguishability of nuclear states related to different Born-Oppenheimer molecular energy curves or electronic excitation modes.
We analyze the entanglement between electronic and nuclear motions in molecular wave functions widely used by theoretical chemists, namely, (i) Born-Oppenheimer factorization in the adiabatic picture, (ii) the transformation into a diabatic picture, (iii) the use of a Born-Huang expansion, and (iv) the eigenfunction of the full molecular Hamiltonian. Our showcase is based on two one-electron one-dimensional molecular Hamiltonians (H2+ and the Shin-Metiu model). We find that within the Born-Oppenheimer approximation, any molecular state (although factorizable) is always entangled, and its entanglement content may be assessed by the variation of the electronic wave function along the different nuclear geometries, with the nuclear wave function indeed playing the role of a tester. The presence of avoided crossings among the adiabatic potential energy curves brings about dramatic changes in the entanglement content of the wave function: sharp avoided crossings favor a diabatic picture (real crossings between potential energy curves), while in broad avoided crossings, the adiabatic picture prevails. The total eigenfunction of the molecular Shin-Metiu Hamiltonian indicates that nuclear densities accommodate well within the diabatic curves for strong adiabatic couplings but within adiabatic curves for weak ones. Consequently, we find that the electron-nuclei entanglement content is a valid witness to unveil strong or weak nonadiabatic couplings in molecules. In terms of entanglement, we also find that the Born-Huang expansion, based on Born-Oppenheimer adiabatic electronic states, does not provide a correct trend of entanglement compared with that of the total molecular eigenfunction, thus indicating a very slow convergence of this expansion.
We investigate the entanglement between electronic and vibrational degrees of freedom produced by a vibronic coupling in a molecular system described in the Born-Oppenheimer approximation. Entanglement in a pure state of the Hilbert space H=Hel⨂Hvib is quantified using the von Neumann entropy of the reduced density matrix and the reduced linear entropy. Expressions for these entanglement measures are derived for the $2\ifmmode\times\else\texttimes\fi{}{N}_{v}$ and $3\ifmmode\times\else\texttimes\fi{}{N}_{v}$ cases of the bipartite entanglement, where 2 and 3 are the dimensions of the electronic Hilbert space Hel, and Nv is the dimension of Hvib. We study the entanglement dynamics for two electronic states coupled by a laser pulse (a $2\ifmmode\times\else\texttimes\fi{}{N}_{v}$ case), taking as an example a coupling between the ${a}^{3}{{\ensuremath{\Sigma}}_{u}}^{+}(6s,6s)$ and 1g(6s,6p3/2) states of the Cs2 molecule. The reduced linear entropy expression obtained for the $3\ifmmode\times\else\texttimes\fi{}{N}_{v}$ case is used to follow the entanglement evolution in a scheme proposed for the control of the vibronic dynamics in a Cs2 cold molecule, implying the ${a}^{3}{{\ensuremath{\Sigma}}_{u}}^{+}(6s,6s)$, ${0}_{g}{}^{\ensuremath{-}}(6s,6{p}_{3/2})$, and ${0}_{g}{}^{\ensuremath{-}}(6s,5d)$ electronic states, which are coupled by a nonadiabatic radial coupling and a sequence of chirped laser pulses.
We characterize both entanglement and quantum coherence in a molecular system by connecting the linear entropy of electronic-nuclear entanglement with Wigner-Yanase skew information measuring vibronic coherence and local quantum uncertainty on electronic energy. Linear entropy of entanglement and quantifiers of quantum coherence are derived for a molecular system described in a bipartite Hilbert space H=Hel⨂Hvib of finite dimension ${N}_{el}\ifmmode\times\else\texttimes\fi{}{N}_{v}$, and relations between them are established. For the specific case of the electronic-vibrational entanglement, we find the linear entropy of entanglement as having a more complex informational content than the von Neumann entropy. By keeping the information carried by the vibronic coherences in a molecule, linear entropy seizes vibrational motion in the electronic potentials as entanglement dynamics. We analyze entanglement oscillations in an isolated molecule, and show examples for the control of entanglement dynamics in a molecule through the creation of coherent vibrational wave packets in several electronic potentials by using chirped laser pulses.
Nuclear quantum phenomena beyond the Born-Oppenheimer approximation are known to play an important role in a growing number of chemical and biological processes. While there exists no unique consensus on a rigorous and efficient implementation of coupled electron-nuclear quantum dynamics, it is recognized that these problems scale exponentially with system size on classical processors and, therefore, may benefit from quantum computing implementations. Here, we introduce a methodology for the efficient quantum treatment of the electron-nuclear problem on near-term quantum computers, based upon the Nuclear-Electronic Orbital (NEO) approach. We generalize the electronic two-qubit tapering scheme to include nuclei by exploiting symmetries inherent in the NEO framework, thereby reducing the Hamiltonian dimension, number of qubits, gates, and measurements needed for calculations. We also develop parameter transfer and initialization techniques, which improve convergence behavior relative to conventional initialization. These techniques are applied to H2 and malonaldehyde for which results agree with NEO full configuration interaction and NEO complete active space configuration interaction benchmarks for ground state energy to within 10-6 hartree and entanglement entropy to within 10-4. These implementations therefore significantly reduce resource requirements for full quantum simulations of molecules on near-term quantum devices while maintaining high accuracy.
Broad in energy optical pulses induce ultrafast molecular dynamics where nuclear degrees of freedom are entangled with electronic ones. We discuss a matrix representation of wave functions of such entangled systems. Singular Value Decomposition (SVD) of this matrix provides a representation as a sum of separable terms. Their weights can be arranged in decreasing order. The representation provided by the SVD is equivalent to a Schmidt decomposition. If there is only one term or if one term is already a good approximation, the system is not entangled. The SVD also provides either an exact or a few term approximation for the partial traces. A simple example, the dynamics of LiH upon ultrafast excitation to several non-adiabatically coupled electronic states, is provided. The major contribution to the entanglement is created during the exit from the Franck Condon region. An additional contribution is the entanglement due to the nuclear motion induced non-adiabatic transitions.