Femtosecond two-photon photoassociation of hot magnesium atoms: A quantum dynamical study using thermal random phase wavefunctions

Diffuse Basis Functions for Relativistic s and d Block Gaussian Basis Sets

Diffuse s, p, and d functions have been optimized for use with previously reported relativistic basis sets for the s and d blocks of the periodic table. The functions were optimized on the 4:1 weighted average of the s² and p² configurations of the anion, with the d shell in the dⁿ⁺¹ configuration for the d blocks. Exponents were extrapolated for groups 2 and 12, which have unstable or weakly bound anions. The diffuse basis sets have been tested by application to calculations of electron affinities of the group 11 elements (Cu, Ag, and Au), double electron affinities of the group 11 monocations, and potential energy curves of Mg₂ and Ca₂ van der Waals dimers, as well as some response properties of the group 1 anions (Rb^-, Cs^-, and Fr^-), the group 2 elements (Sr, Ba, and Ra), and RbLi, CsLi, and FrLi molecules.

Investigation of photoassociation with full-dimensional thermal-random-phase wavefunctions

By taking the femtosecond two-photon photoassociation (PA) of magnesium atoms as an example, we propose a method to calculate the thermally averaged population, which is transferred from the ground X¹Σ_g ⁺ state to the target (1)¹Π_g state, based on the solution of full-dimensional time-dependent Schrödinger equation. In this method, named as method A, we use thermal-random-phase wavefunctions with the random phases expanded in both the vibrational and rotational degrees of freedom to model the thermal ensemble of the initial eigenstates. This method is compared with the other two methods (B and C) at different temperatures. Method B is also based on thermal-random-phase wavefunctions, except that the random-phase expansion is merely used for the vibrational degree of freedom. Method C is based on the independent propagation of every initial eigenstate, instead of the thermal-random-phase wavefunctions. Taking the (1)¹Π_g state as the target state, it is found that although these three methods can present the same population on the (1)¹Π_g state, the computation efficiency of method A increases dramatically with the increase in temperature. With this efficient method A, we find that the PA process at 1000 K can also induce rotational coherence, i.e., the molecular field-free alignment in the excited electronic states.

Quantum computation solves a half-century-old enigma: Elusive vibrational states of magnesium dimer found

The high-lying vibrational states of the magnesium dimer (Mg₂), which has been recognized as an important system in studies of ultracold and collisional phenomena, have eluded experimental characterization for half a century. Until now, only the first 14 vibrational states of Mg₂ have been experimentally resolved, although it has been suggested that the ground-state potential may support five additional levels. Here, we present highly accurate ab initio potential energy curves based on state-of-the-art coupled-cluster and full configuration interaction computations for the ground and excited electronic states involved in the experimental investigations of Mg₂. Our ground-state potential unambiguously confirms the existence of 19 vibrational levels, with ~1 cm^-1 root mean square deviation between the calculated rovibrational term values and the available experimental and experimentally derived data. Our computations reproduce the latest laser-induced fluorescence spectrum and provide guidance for the experimental detection of the previously unresolved vibrational levels.

Optimized sampling of mixed-state observables

Quantum dynamical simulations of statistical ensembles pose a significant computational challenge due to the fact that mixed states need to be represented. If the underlying dynamics is fully unitary, for example, in ultrafast coherent control at finite temperatures, then one approach to approximate time-dependent observables is to sample the density operator by solving the Schrödinger equation for a set of wave functions with randomized phases. We show that, on average, random-phase wave functions perform well for ensembles with high mixedness, whereas at higher purities a deterministic sampling of the energetically lowest-lying eigenstates becomes superior. We prove that minimization of the worst-case error for computing arbitrary observables is uniquely attained by eigenstate-based sampling. We show that this error can be used to form a qualitative estimate of the set of ensemble purities for which the sampling performance of the eigenstate-based approach is superior to random-phase wave functions. Furthermore, we present refinements to both schemes which remove redundant information from the sampling procedure to accelerate their convergence. Finally, we point out how the structure of low-rank observables can be exploited to further improve eigenstate-based sampling schemes.

Quantum thermodynamics and open-systems modeling

A comprehensive approach to modeling open quantum systems consistent with thermodynamics is presented. The theory of open quantum systems is employed to define system bath partitions. The Markovian master equation defines an isothermal partition between the system and bath. Two methods to derive the quantum master equation are described: the weak coupling limit and the repeated collision model. The role of the eigenoperators of the free system dynamics is highlighted, in particular, for driven systems. The thermodynamical relations are pointed out. Models that lead to loss of coherence, i.e., dephasing are described. The implication of the laws of thermodynamics to simulating transport and spectroscopy is described. The indications for self-averaging in large quantum systems and thus its importance in modeling are described. Basic modeling by the surrogate Hamiltonian is described, as well as thermal boundary conditions using the repeated collision model and their use in the stochastic surrogate Hamiltonian. The problem of modeling with explicitly time dependent driving is analyzed. Finally, the use of the stochastic surrogate Hamiltonian for modeling ultrafast spectroscopy and quantum control is reviewed.

Transition moments between excited electronic states from the Hermitian formulation of the coupled cluster quadratic response function

We introduce a new method for the computation of the transition moments between the excited electronic states based on the expectation value formalism of the coupled cluster theory [B. Jeziorski and R. Moszynski, Int. J. Quantum Chem. 48, 161 (1993)]. The working expressions of the new method solely employ the coupled cluster operator T and an auxiliary operator S that is expressed as a finite commutator expansion in terms of T and T^†. In the approximation adopted in the present paper, the cluster expansion is limited to single, double, and linear triple excitations. The computed dipole transition probabilities for the singlet-singlet and triplet-triplet transitions in alkali earth atoms agree well with the available theoretical and experimental data. In contrast to the existing coupled cluster response theory, the matrix elements obtained by using our approach satisfy the Hermitian symmetry even if the excitations in the cluster operator are truncated, but the operator S is exact. The Hermitian symmetry is slightly broken if the commutator series for the operator S are truncated. As a part of the numerical evidence for the new method, we report calculations of the transition moments between the excited triplet states which have not yet been reported in the literature within the coupled cluster theory. Slater-type basis sets constructed according to the correlation-consistency principle are used in our calculations.

Rotational Control of Asymmetric Molecules: Dipole- versus Polarizability-Driven Rotational Dynamics

We experimentally study the optical- and terahertz-induced rotational dynamics of asymmetric molecules in the gas phase. Terahertz and optical fields are identified as two distinct control handles over asymmetric molecules, as they couple to the rotational degrees of freedom via the molecular dipole and polarizability selectively. The distinction between those two rotational handles is highlighted by different types of quantum revivals observed in long-duration (>100  ps) field-free rotational evolution. The experimental results are in excellent agreement with random phase wave function (RPWF) simulations [Phys. Rev. A 91, 063420 (2015)] and provide verification of the RPWF as an efficient method for calculating asymmetric molecular dynamics at ambient temperatures, where exact calculation methods are practically not feasible. Our observations and analysis pave the way for orchestrated excitations by both optical and terahertz fields as complementary rotational handles that enable a plethora of new possibilities in three-dimensional rotational control of asymmetric molecules.

Analytical optimal pulse shapes obtained with the aid of genetic algorithms

We propose a methodology to design optimal pulses for achieving quantum optimal control on molecular systems. Our approach constrains pulse shapes to linear combinations of a fixed number of experimentally relevant pulse functions. Quantum optimal control is obtained by maximizing a multi-target fitness function using genetic algorithms. As a first application of the methodology, we generated an optimal pulse that successfully maximized the yield on a selected dissociation channel of a diatomic molecule. Our pulse is obtained as a linear combination of linearly chirped pulse functions. Data recorded along the evolution of the genetic algorithm contained important information regarding the interplay between radiative and diabatic processes. We performed a principal component analysis on these data to retrieve the most relevant processes along the optimal path. Our proposed methodology could be useful for performing quantum optimal control on more complex systems by employing a wider variety of pulse shape functions.