The Fourier method for tri-atomic systems in the search for the optimal coordinate system

Dynamics of FHCl Including Two Bidirectional Dissociation Channels: Comparative Study Using Quantum Nuclear Wavepackets and Semiclassical Trajectory Surface Hopping

The dynamics starting from the charge transfer excited state of neutral FHCl includes two bidirectional dissociation channels, producing "FH" and "Cl" fragments or "F" and "HCl" fragments by electron-transfer (ET) or proton-transfer (PT) processes, respectively. The quantum features of its dynamics were studied using the quantum dynamics of nuclear wavepacket propagation and the semiclassical dynamics of trajectory surface hopping propagation. The MS-CASPT2(17,11)/aug-cc-pVTZ method was used for calculating the energies of electronic states. Two critical quantum features identified in previous studies were the dominance of PT over ET and the ∼80 fs time gap between the onset of the earlier PT and the later ET processes. These features, in contrast to classical anticipation, were confirmed again, and their classical interpretations were developed based on the results of semiclassical dynamics. The relative location of nonadiabatic regions with respect to the starting point of dynamics and the activation of angular motion altering the HFCl angle play crucial roles in governing the two quantum features.

Nuclear wave-packet-propagation-based study of the electron-coupled, proton-transfer process in the charge-transfer state of FHCl exhibiting three electronic states in full-dimensional space

The charge-transfer (CT) excited state of FHCl (F⁺H-Cl^-), generated by the photodetachment of an electron from its precursor anion (FHCl^-) by a photon energy of ∼9.5 eV, is a realistic prototype of two bidirectional-coupled reaction pathways, namely the proton-transfer (PT) and electron-transfer (ET) channels, that produce F + HCl and FH + Cl combinations, respectively. The early-time dynamics of the CT was studied via the time-dependent propagations of nuclear wave packets comprising three nonadiabatically coupled electronic states defined within a three-dimensional space. The detailed analyses of the early-time dynamics revealed an interesting phenomenon in which the onset of PT was ∼80 fs earlier than that of ET, indicating that PT dominated ET in this case. A more significant finding was that the proper adjustment of the electronic-charge distribution for the onset of ET was obtained ∼80 fs after the onset of PT; this adjustment was mediated by the initial movement of the H atom, i.e., the F-H vibration mode. To avail experimental observables, the branching ratio, χ = PT/(PT + ET), and absorption spectrum generating the neutral FHCl molecule from its precursor anion were also simulated. The results further demonstrated the dependences of the χs and spectrum on the change in the initial vibration level of the precursor anion, as well as the isotopic substitution of the connecting H atom with deuterium, tritium, and muonium.

Wavepacket propagations for the early time dynamics of proton-coupled electron transfer in the charge-transfer state of NH3Cl complex

A charge-transfer (CT) excited state of NH₃Cl, generated by photo-detachment of an electron from the anionic NH₃Cl^- precursor, can be represented as H₂N⁺-H-Cl^- and proceeds to two chemical reactions: one reaction generating NH₂ and HCl resulting from a proton transfer (PT) and the other reaction producing NH₃ and a Cl atom resulting from an electron transfer (ET); both are coupled to form a typical proton-coupled electron transfer (PCET) process. The early time dynamics of this CT were studied using time-dependent wavepacket propagation on three nonadiabatically coupled electronic states in a reduced three-dimensional space. The electronic states were treated using the XMS-CASPT2/aug-cc-pVTZ ab initio methodology. The population dynamics of the three coupled electronic states were analyzed in detail to reveal the initial stage of the PCET process up to ∼100 fs, while the branching ratio, χ = PT/(ET+PT), was determined after wavepacket propagations of up to 2000 fs. Another main result is the dependence of χ on the vibration levels of the initial precursor anion and the isotope substitution of the connecting H atom with deuterium and tritium. Our study reveals the detailed microscopic features of the PCET process embedded in the CT state of the NH₃Cl complex and certain systematic dependences of the branching ratio χ on the above factors.

Competitive bond breaking in floppy molecular trimers: HeNeH and HeNeH- calculations

Calculations for the J=0 bound states associated with the two title molecular systems are carried out using distributed Gaussian function approach and using the sum of two-body potential approximation to describe the overall interactions. The results yield one bound state for the neutral trimer and three bound states for the anionic triatom. The relative values of dissociation energies place the first dissociation threshold at the H emission for the neutral complex and at the He emission for the corresponding anion. The general spatial properties of the various bound states are analyzed on both systems.

Abstractive dissociation of oxygen over Al(111): a nonadiabatic quantum model

The dissociation of oxygen on a clean aluminum surface is studied theoretically. A nonadiabatic quantum dynamical model is used, based on four electronically distinct potential energy surfaces characterized by the extent of charge transfer from the metal to the adsorbate. A flat surface approximation is used to reduce the computation complexity. The conservation of the helicopter angular momentum allows Boltzmann averaging of the outcome of the propagation of a three degrees of freedom wave function. The dissociation event is simulated by solving the time-dependent Schrödinger equation for a period of 30 femtoseconds. As a function of incident kinetic energy, the dissociation yield follows the experimental trend. An attempt at simulation employing only the lowest adiabatic surface failed, qualitatively disagreeing with both experiment and nonadiabatic calculations. The final products, adsorptive dissociation and abstractive dissociation, are obtained by carrying out a semiclassical molecular dynamics simulation with surface hopping which describes the back charge transfer from an oxygen atom negative ion to the surface. The final adsorbed oxygen pair distribution compares well with experiment. By running the dynamical events backward in time, a correlation is established between the products and the initial conditions which lead to their production. Qualitative agreement is thus obtained with recent experiments that show suppression of abstraction by rotational excitation.

On the use of optimal internal vibrational coordinates for symmetrical bent triatomic molecules

The use of generalized internal coordinates for the variational calculation of excited vibrational states of symmetrical bent triatomic molecules is considered with applications to the SO2, O3, NO2, and H2O molecules. These coordinates depend on two external parameters which can be properly optimized. We propose a simple analytical method to determine the optimal internal coordinates for this kind of molecules based on the minimization with respect to the external parameters of the zero-point energy, assuming only quadratic terms in the Hamiltonian and no quadratic coupling between the optimal coordinates. The optimal values of the parameters thus obtained are shown to agree quite well with those that minimize the sum of a number of unconverged energies of the lowest vibrational states, computed variationally using a small basis function set. The unconverged variational calculation uses a basis set consisting of the eigenfunctions of the uncoupled anharmonic internal coordinate Hamiltonian. Variational calculations of the excited vibrational states for the four molecules considered carried out with an increasing number of basis functions, also evidence the excellent convergence properties of the optimal internal coordinates versus those provided by other normal and local coordinate systems.

Energy levels and wave functions of weakly-bound 4Hex 20NeyH (x+y=2) systems using Pekeris coordinates and a symmetry-adapted Lanczos approach

Energy levels and wave functions of floppy triatomic rare gas hydrides are calculated using a Pekeris coordinate system and the importance of various triangular configurations is assessed through the calculation of reduced distribution functions and relative weights. The calculations are performed using a symmetry-adapted Lanczos recursion within the discrete variable representation. For the 4He2H- anion, the present results are compared with those obtained from calculations based on other methods, and the accuracy of the present method is discussed. Calculations are also performed for the case of 4He2H and 20Ne2H, as well as for the mixed 4He 20NeH neutrals. Our results show that no bound states are found for 4He2H while only one bound state is found for both the 20Ne2H and 4He 20NeH complexes. Interestingly, a very important and common property of these systems is that there is a significant contribution from linear configurations to their bound states.