Jason N. Byrd

At What Chain Length Do Unbranched Alkanes Prefer Folded Conformations

Short unbranched alkanes are known to prefer linear conformations, whereas long unbranched alkanes are folded. It is not known with certainty at what chain length the linear conformation is no longer the global minimum. To clarify this point, we use ab initio and density functional methods to compute the relative energies of the linear and hairpin alkane conformers for increasing chain lengths. Extensive electronic structure calculations are performed to obtain optimized geometries, harmonic frequencies, and accurate single point energies for the selected alkane conformers from octane through octadecane. Benchmark CCSD(T)/cc-pVTZ single point calculations are performed for chains through tetradecane, whereas approximate methods are required for the longer chains up to octadecane. Using frozen natural orbitals to unambiguously truncate the virtual orbital space, we are able to compute composite CCSD FNO(T) single point energies for all the chain lengths. This approximate composite method has significant computational savings compared to full CCSD(T) while retaining ∼0.15 kcal/mol accuracy compared to the benchmark results. More approximate dual-basis resolution-of-the-identity double-hybrid DFT calculations are also performed and shown to have reasonable 0.2-0.4 kcal/mol errors compared with our benchmark values. After including contributions from temperature dependent internal energy shifts, we find the preference for folded conformations to lie between hexadecane and octadecane, in excellent agreement with recent experiments [ Lüttschwager , N. O. ; Wassermann , T. N. ; Mata , R. A. ; Suhm , M. A. Angew. Chem. Int. Ed. 2013 , 52 , 463 ].

Structure and thermochemistry of K 2 Rb, KRb 2 , and K 2 Rb 2

The formation and interaction of ultracold polar molecules is a topic of active research. Understanding possible reaction paths and molecular combinations requires accurate studies of the fragment and product energetics. We have calculated accurate gradient optimized ground-state structures and zero-point corrected atomization energies for the trimers and tetramers formed by the reaction of KRb with KRb and corresponding isolated atoms. The K{sub 2}Rb and KRb{sub 2} trimers are found to have global minima at the C{sub 2v} configuration with atomization energies of 6065 and 5931 cm{sup -1} while the tetramer is found to have two stable planar structures, of D{sub 2h} and C{sub s} symmetry, which have atomization energies of 11131 cm{sup -1} and 11133 cm{sup -1}, respectively. We have calculated the minimum energy reaction path for the reaction KRb + KRb to K{sub 2} + Rb{sub 2} and found it to be barrierless.

Approximating electronically excited states with equation-of-motion linear coupled-cluster theory

A new perturbative approach to canonical equation-of-motion coupled-cluster theory is presented using coupled-cluster perturbation theory. A second-order Møller-Plesset partitioning of the Hamiltonian is used to obtain the well known equation-of-motion many-body perturbation theory equations and two new equation-of-motion methods based on the linear coupled-cluster doubles and linear coupled-cluster singles and doubles wavefunctions. These new methods are benchmarked against very accurate theoretical and experimental spectra from 25 small organic molecules. It is found that the proposed methods have excellent agreement with canonical equation-of-motion coupled-cluster singles and doubles state for state orderings and relative excited state energies as well as acceptable quantitative agreement for absolute excitation energies compared with the best estimate theory and experimental spectra.

Long-range interactions between like homonuclear alkali metal diatoms

Long-range electrostatic and van der Waals coefficients up to terms of order R(-8) have been evaluated by the sum over states method using ab initio and time-dependent density functional theory. We employ several widely used density functionals and systematically investigate the convergence of the calculated results with basis set size. Static electric moments and polarizabilities up to octopole order are also calculated. We present values for Li(2) through K(2) which are in good agreement with existing values, in addition to new results for Rb(2) and Cs(2). Interaction potential curves calculated from these results are shown to agree well with high level ab initio theory. Preliminary results are reported that demonstrate the applicability of the method to larger alkali clusters.

Tuning ultracold chemical reactions via Rydberg-dressed interactions.

We show that ultracold chemical reactions with an activation barrier can be tuned using Rydberg-dressed interactions. Scattering in the ultracold regime is sensitive to long-range interactions, especially when weakly bound (or quasibound) states exist near the collision threshold. We investigate how, by Rydberg dressing a reactant, one enhances its polarizability and modifies the long-range van der Waals collision complex, which can alter chemical reaction rates by shifting the position of near-threshold bound states. We carry out a full quantum mechanical scattering calculation for the benchmark system H(2)+D, and show that resonances can be moved substantially and that rate coefficients at cold and ultracold temperatures can be increased by several orders of magnitude.

Long-range forces between polar alkali-metal diatoms aligned by external electric fields

Long range electrostatic, induction and dispersion coefficients including terms of order

Long range intermolecular interactions between the alkali diatomics Na2, K2, and NaK

R^{-8}

Molecular cluster perturbation theory. I. Formalism

have been calculated by the sum over states method using time dependent density functional theory. We also computed electrostatic moments and static polarizabilities of the individual diatoms up to the octopole order using coupled cluster and density functional theory. The laboratory-frame transformed electrostatic moments and van der Waals coefficients corresponding to the alignment of the diatomic molecules were found. We use this transformation to obtain the coupling induced by an external DC electric field, and present values for all XY combinations of like polar alkali diatomic molecules with atoms from Li to Cs. Analytic solutions to the dressed-state laboratory-frame electrostatic moments and long range intermolecular potentials are also given for the DC low-field limit.

Controllable binding of polar molecules and metastability of one-dimensional gases with attractive dipole forces.

Long range interactions between the ground state alkali diatomics Na(2)-Na(2), K(2)-K(2), Na(2)-K(2), and NaK-NaK are examined. Interaction energies are first determined from ab initio calculations at the coupled-cluster with singles, doubles, and perturbative triples [CCSD(T)] level of theory, including counterpoise corrections. Long range energies calculated from diatomic molecular properties (polarizabilities and dipole and quadrupole moments) are then compared with the ab initio energies. A simple asymptotic model potential E(LR)=E(elec)+E(disp)+E(ind) is shown to accurately represent the intermolecular interactions for these systems at long range.

Ab initio potential curves for the X2Σu+, A2Πu and B2Σg+ states of Ca2+

We present second-order molecular cluster perturbation theory (MCPT(2)), a linear scaling methodology to calculate arbitrarily large systems with explicit calculation of individual wave functions in a coupled-cluster framework. This new MCPT(2) framework uses coupled-cluster perturbation theory and an expansion in terms of molecular dimer interactions to obtain molecular wave functions that are infinite order in both the electronic fluctuation operator and all possible dimer (and products of dimers) interactions. The MCPT(2) framework has been implemented in the new SIA/Aces4 parallel architecture, making use of the advanced dynamic memory control and fine-grained parallelism to perform very large explicit molecular cluster calculations. To illustrate the power of this method, we have computed energy shifts, lattice site dipole moments, and harmonic vibrational frequencies via explicit calculation of the bulk system for the polar and non-polar polymorphs of solid hydrogen fluoride. The explicit lattice size (without using any periodic boundary conditions) was expanded up to 1000 HF molecules, with 32,000 basis functions and 10,000 electrons. Our obtained HF lattice site dipole moments and harmonic vibrational frequencies agree well with the existing literature.