Origin of the complex dynamics in structural isomerization of small clusters: The effects of potential topography

Energy natural orbital characterization of nonadiabatic electron wavepackets in the densely quasi-degenerate electronic state manifold

Dynamics and energetic structure of largely fluctuating nonadiabatic electron wavepackets are studied in terms of Energy Natural Orbitals (ENOs) [K. Takatsuka and Y. Arasaki, J. Chem. Phys. 154, 094103 (2021)]. Such huge fluctuating states are sampled from the highly excited states of clusters of 12 boron atoms (B₁₂), which have densely quasi-degenerate electronic excited-state manifold, each adiabatic state of which gets promptly mixed with other states through the frequent and enduring nonadiabatic interactions within the manifold. Yet, the wavepacket states are expected to be of very long lifetimes. This excited-state electronic wavepacket dynamics is extremely interesting but very hard to analyze since they are usually represented in large time-dependent configuration interaction wavefunctions and/or in some other complicated forms. We have found that ENO gives an invariant energy orbital picture to characterize not only the static highly correlated electronic wavefunctions but also those time-dependent electronic wavefunctions. Hence, we first demonstrate how the ENO representation works for some general cases, choosing proton transfer in water dimer and electron-deficient multicenter chemical bonding in diborane in the ground state. We then penetrate with ENO deep into the analysis of the essential nature of nonadiabatic electron wavepacket dynamics in the excited states and show the mechanism of the coexistence of huge electronic fluctuation and rather strong chemical bonds under very random electron flows within the molecule. To quantify the intra-molecular energy flow associated with the huge electronic-state fluctuation, we define and numerically demonstrate what we call the electronic energy flux.

Quantum Chaos in the Dynamics of Molecules

Quantum chaos is reviewed from the viewpoint of "what is molecule?", particularly placing emphasis on their dynamics. Molecules are composed of heavy nuclei and light electrons, and thereby the very basic molecular theory due to Born and Oppenheimer gives a view that quantum electronic states provide potential functions working on nuclei, which in turn are often treated classically or semiclassically. Therefore, the classic study of chaos in molecular science began with those nuclear dynamics particularly about the vibrational energy randomization within a molecule. Statistical laws in probabilities and rates of chemical reactions even for small molecules of several atoms are among the chemical phenomena requiring the notion of chaos. Particularly the dynamics behind unimolecular decomposition are referred to as Intra-molecular Vibrational energy Redistribution (IVR). Semiclassical mechanics is also one of the main research fields of quantum chaos. We herein demonstrate chaos that appears only in semiclassical and full quantum dynamics. A fundamental phenomenon possibly giving birth to quantum chaos is "bifurcation and merging" of quantum wavepackets, rather than "stretching and folding" of the baker's transformation and the horseshoe map as a geometrical foundation of classical chaos. Such wavepacket bifurcation and merging are indeed experimentally measurable as we showed before in the series of studies on real-time probing of nonadiabatic chemical reactions. After tracking these aspects of molecular chaos, we will explore quantum chaos found in nonadiabatic electron wavepacket dynamics, which emerges in the realm far beyond the Born-Oppenheimer paradigm. In this class of chaos, we propose a notion of Intra-molecular Nonadiabatic Electronic Energy Redistribution (INEER), which is a consequence of the chaotic fluxes of electrons and energy within a molecule.

Local integrals and their globally connected invariant structure in phase space giving rise to a promoting mode of chemical reaction

We develop a method to extract local integrals, that is, integrals defined locally in the linear regime of an arbitrary point in phase space. The individual integral represents a vibrational mode. We also propose an index that quantifies the extent of connection between neighboring local integrals. Those pieces that are smoothly connected over a wide range represent a global structure of phase space. With a six-atomic Lennard-Jones cluster, we show that it is possible to identify which vibrational mode in the potential basin correlates smoothly to that in the area of transition state, which is nothing but a reactive mode. As an application of the method, we attempt to enhance the structural transition by exciting the reactive mode thus found. This method works successfully as shown in numerical calculations.

Nonempirical Statistical Theory for Atomic Evaporation from Nonrigid Clusters: Applications to the Absolute Rate Constant and Kinetic Energy Release

A high energy atomic cluster undergoing frequent structural isomerization behaves like a liquid droplet, from which atoms or molecules can be emitted. Even after evaporation, the daughter cluster may still keep changing its structure. We study the dynamics of such an evaporation process of atomic evaporation. To do so, we develop a statistical rate theory for dissociation of highly nonrigid molecules and propose a simple method to calculate the absolute value of classical phase-space volume for a potential function that has many locally stable basins. The statistical prediction of the final distribution of the released kinetic energy is also developed. A direct application of the Rice-Ramsperger-Kassed-Marcus (RRKM) theory to this kind of multichannel chemical reaction is prohibitively difficult, unless further modeling and/or assumptions are made. We carry out a completely nonempirical statistical calculation for these dynamical quantities, in that nothing empirical is introduced like remodeling (or reparametrization) of artificial potential energy functions or recalibration of the phase-space volume referring to other "empirical" values such as those estimated with the molecular dynamics method. The so-called dividing surface is determined variationally, at which the flux is calculated in a consistent manner with the estimate of the phase-space volume in the initial state. Also, for the correct treatment of a highly nonrigid cluster, the phase-space volume and flux are estimated without the separation of vibrational and rotational motions. Both the microcanonical reaction rate and the final kinetic energy distribution thus obtained have quite accurately reproduced the corresponding quantities given by molecular dynamics calculations. This establishes the validity of the statistical arguments, which in turn brings about the deeper physical insight about the evaporation dynamics.

Semiclassical quantization of strongly chaotic vibrations in an M7-like cluster

We report semiclassical energy spectra of vibrational state of a cluster composed of seven identical atoms like Ar7 in terms of our previously developed semiclassical wave function, which we call the action-decomposed function. The classical dynamics of this vibrational state is strongly chaotic and undergoes a large amplitude motion due to structural isomerization, which demands a long run of trajectory calculation. Permutation of identical particles should also be taken into account as a quantum effect, since a single molecular shape can be shared by many permutational isomers. Furthermore, chaos causes a spurious divergence in the amplitude factor of a correlation function in the initial value representation, which arises from the amplitude factor (prefactor) of a semiclassical wave function, while the final-state representation is suffered from the well-known divergence arising at caustics. Both approaches therefore face tremendous difficulty in a long-time calculation of the correlation functions. We challenge to extract some limited number of spectral lines from such chaotic dynamics. We further apply a correlation function that is free of such a troublesome amplitude factor. Numerical results from all these schemes are reported.