Source author record

Jason N. Byrd

Jason N. Byrd appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2020arXiv

Leveraging the Super Instruction Architecture to Develop Massively Parallel Computational and Environmental Chemistry Applications

The task of developing high performing parallel software must be made easier and more cost effective in order to fully exploit existing and emerging large scale computer systems for the advancement of science. The Super Instruction Architecture is a parallel programming platform geared towards applications that need to manage large amounts of data stored in potentially sparse multidimensional arrays during calculations. The SIA platform was originally designed for the Quantum Chemistry software package ACESIII. More recently, the SIA was reimplemented to overcome limitations in the original ACESIII program. It has now been successfully employed in the new Aces4 Quantum Chemistry software package and the development of the atmospheric transport application MATLOC, thus demonstrating the versatility of the SIA approach. MATLOC calculates transport and dispersion of mass over regions in the range of 100-1000s of square kilometers and is a significant improvement over existing community codes. This paper describes results from both the transport and dispersion application as well as some difficult Quantum Chemistry open shell coupled cluster benchmark calculations using Aces4.

preprint2016arXiv

Predictive coupled-cluster isomer orderings for some Si${}_n$C${}_m$ ($m, n\le 12$) clusters; A pragmatic comparison between DFT and complete basis limit coupled-cluster benchmarks

The accurate determination of the preferred ${\rm Si}_{12}{\rm C}_{12}$ isomer is important to guide experimental efforts directed towards synthesizing SiC nano-wires and related polymer structures which are anticipated to be highly efficient exciton materials for opto-electronic devices. In order to definitively identify preferred isomeric structures for silicon carbon nano-clusters, highly accurate geometries, energies and harmonic zero point energies have been computed using coupled-cluster theory with systematic extrapolation to the complete basis limit for set of silicon carbon clusters ranging in size from SiC$_3$ to ${\rm Si}_{12}{\rm C}_{12}$. It is found that post-MBPT(2) correlation energy plays a significant role in obtaining converged relative isomer energies, suggesting that predictions using low rung density functional methods will not have adequate accuracy. Utilizing the best composite coupled-cluster energy that is still computationally feasible, entailing a 3-4 SCF and CCSD extrapolation with triple-$ζ$ (T) correlation, the {\it closo} ${\rm Si}_{12}{\rm C}_{12}$ isomer is identified to be the preferred isomer in support of previous calculations [J. Chem. Phys. 2015, 142, 034303]. Additionally we have investigated more pragmatic approaches to obtaining accurate silicon carbide isomer energies, including the use of frozen natural orbital coupled-cluster theory and several rungs of standard and double-hybrid density functional theory. Frozen natural orbitals as a way to compute post MBPT(2) correlation energy is found to be an excellent balance between efficiency and accuracy.

preprint2015arXiv

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 (EOM-MBPT(2)) equations and two new equation-of-motion methods based on the linear coupled-cluster doubles (EOM-LCCD) and linear coupled-cluster singles and doubles (EOM-LCCSD) wavefunctions. This is achieved by performing a short-circuiting procedure on the MBPT(2) similarity transformed Hamiltonian. 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 EOM-CCSD 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.

preprint2015arXiv

Computationally efficient double hybrid density functional theory using dual basis methods

We examine the application of the recently developed dual basis methods of Head-Gordon and co-workers to double hybrid density functional computations. Using the B2-PLYP, B2GP-PLYP, DSD-BLYP and DSD-PBEP86 density functionals, we assess the performance of dual basis methods for the calculation of conformational energy changes in C$_4$-C$_7$ alkanes and for the S22 set of noncovalent interaction energies. The dual basis methods, combined with resolution-of-the-identity second-order Møller-Plesset theory, are shown to give results in excellent agreement with conventional methods at a much reduced computational cost.

preprint2015arXiv

Molecular Cluster Perturbation Theory. I. Formalism

We present second-order molecular cluster perturbation theory (MCPT(2)), a linear scaling methodology to calculate arbitrarily large systems with explicit calculation of individual wavefunctions 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 wavefunctions 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 1,000 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.

preprint2014arXiv

Ab initio potential curves for the X $^2Σ_u^+$, A $^2Π_u$ and B $^2Σ_g^+$ states of Ca$_{2}^+$

We report \textit{ab initio} calculations of the X $^2 Σ_{u}^+$, A $^2Π_u$ and B $^2 Σ_{g}^+$ states of the Ca$_{2}^+$ dimer. All electron CAS+MRCI calculations are performed for the X $^2 Σ_{u}^+$ and B $^2 Σ_{g}^+$ states, while valence CAS+MRCI calculations using an effective core potential are used to describe the A $^2Π_u$ state. A double well is found in the B $^2 Σ_{g}^+$ state. Spectroscopic constants, vibrational levels, transition moments and radiative lifetimes are calculated for the most abundant isotope of calcium ($^{40}$Ca). The static dipole and quadrupole polarizabilities, and the leading order van der Waals coefficients are also calculated for all three states.

preprint2014arXiv

At what chain length do unbranched alkanes prefer folded conformations?

Short unbranched alkanes are known to prefer linear conformations, while 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 {\it 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, while 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 $\sim0.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. {\it Angew. Chem. Int. Ed.} 2013, 52, 463].

preprint2014arXiv

Controlling ultracold chemical reactions via Rydberg-dressed interactions

We show that ultracold chemical reactions can be manipulated and controlled by using Rydberg-dressed interactions. Scattering in the ultracold regime is sensitive to long-range interactions, especially when weakly bound (or quasi-bound) 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.

preprint2014arXiv

Correlation correction to configuration interaction singles from coupled cluster perturbation theory

A new state specific correlation correction to configuration interaction singles (CIS) excitation energies is preseted using coupled cluster perturbation theory (CCPT). General expressions for CIS-CCPT are derived and expanded explicitly to first order in the wavefunction and second order in the energy. By virtue of the nature of CCPT this method is {\it a priori} size extensive and incorporates infinite order effects into the wavefunction. This results in a balanced singles space excited state theory that at second order is an improvement over the ubiquitous CIS(D) method and comparable in quality to equation of motion coupled cluster (EOM-CC). A modest test set composed of the first four excited states from nine small organic molecules was used to quantify the accuracy and consistency of the CIS-CCPT2 excitation energies and density of states. We find that CIS-CCPT2 has a standard deviation error of 0.18 eV for excitation energies and 0.14 eV for density of states compared to EOM-CC, a factor of two better than CIS(D) with a significant reduction in the maximum deviation as well

preprint2013arXiv

Associative detachment of rubidium hydroxide

We performed calculations of the optimized structure, harmonic vibrational frequencies and dissociation energies of RbOH and its anion, and investigate the interactions between Rb and OH$^-$ leading to possible associative detachment pathways. The electron affinity of RbOH was computed to be 0.2890 eV, with a bond energy of Rb+OH$^-$ half that of Rb+OH. To determine other possible charge loss pathways, the Rb+OH and Rb+OH$^-$ dissociation curves were computed using couple cluster methods along all possible collisional angles. An adiabatic curve crossing between the neutral and charged molecule was found at the inner wall of the molecular potential curve for linear geometries. Associative detachment rates were estimated using the Langevin ion capture cross-section for hydroxide. We find for $v\ge 2$ an associative detachment rate of $>2\times 10^{-9}$ cm$^3$s$^{-1}$, while for $v=0$ and 1 no appreciable rate exists. This strong dependence on vibrational level suggests the ability to control the associative detachment rate directly.

preprint2012arXiv

Controllable binding of polar molecules and meta-stability of 1-D gases with attractive dipole forces

We explore one-dimensional (1-D) samples of ultracold polar molecules with attractive dipole-dipole interactions and show the existence of a repulsive barrier due to a strong quadrupole interaction between molecules. This barrier can stabilize a gas of ultracold KRb molecules and even lead to long-range wells supporting bound states between molecules. The properties of these wells can be controlled by external electric fields, allowing the formation of long polymer-like chains of KRb, and studies of quantum phase transitions by varying the effective interaction between molecules. We discuss the generalization of those results to other systems.

preprint2012arXiv

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

Long range electrostatic, induction and dispersion coefficients including terms of order $R^{-8}$ 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.

preprint2012arXiv

Long-range three-body atom-diatom potential for doublet Li${}_3$

An accurate long-range {\em ab initio} potential energy surface has been calculated for the ground state ${}^2A'$ lithium trimer in the frozen diatom approximation using all electron RCCSD(T). The {\em ab initio} energies are corrected for basis set superposition error and extrapolated to the complete basis limit. Molecular van der Waals dispersion coefficients and three-body dispersion damping terms for the atom-diatomic dissociation limit are presented from a linear least squares fit and shown to be an essentially exact representation of the {\em ab initio} surface at large range.

preprint2011arXiv

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 Li2 through K2 which are in good agreement with existing values, in addition to new results for Rb2 and Cs2. Interaction potential curves calculated from these results are shown to agree well with high level it ab initio theory.

preprint2010arXiv

Ab initio potential curves for the X $^2 Σ_{u}^+$ and B $^2 Σ_{g}^+$ states of Be$_{2}^+$: Existence of a double minimum

We report ab initio calculations of the X $^2 Σ_{u}^+ $ and B $^2 Σ_{g}^+$ states of the Be$_{2}^+$ dimer. Full valence configuration interaction calculations were performed using the aug-cc-pVnZ basis sets and the results were extrapolated to the CBS limit. Core-core, core-valence effects are included at the CCSDT/MTsmall level of theory. Two local minima, separated by a large barrier, are found in the expected repulsive B $^2 Σ_{g}^+$ state. Spectroscopic constants have been calculated and good agreement is found with the recent measurements of Merritt et al. Bound vibrational levels, transition moments and lifetimes have also been calculated.

preprint2010arXiv

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${}_2$Rb andKRb${}_2$ trimers are found to have global minima at the $C_{2v}$ configurationwith atomization energies of 6065 and 5931 cm${}^{-1}$ while the tetramer is found to have two stable planar structures, of $D_{2h}$ and $C_s$ symmetry, which have atomization energies of 11131 cm${}^{-1}$ an 11133 cm${}^{-1}$, respectively. We have calculated the minimum energy reaction path for the reaction KRb+KRb to K$_2+$Rb$_2$ and found it to be barrierless.

preprint2009arXiv

Potential energy surface of the 2A' Li2+Li doublet ground state

The lowest doublet electronic state for the lithium trimer (2A') is calculated for use in three-body scattering calculations using the valence electron FCI method with atomic cores represented using an effective core potential. It is shown that an accurate description of core-valence correlation is necessary for accurate calculations of molecular bond lengths, frequencies and dissociation energies. Interpolation between 2A' ab initio surface data points in a sparse grid is done using the global interpolant moving least squares method with a smooth radial data cutoff function included in the fitting weights and bivariate polynomials as a basis set. The Jahn-Teller splitting of the 2E' surface into the 2A1 and 2B2 states is investigated using a combination of both CASSCF and FCI levels of theory. Additionally, preliminary calculations of the 2A'' surface are also presented using second order spin restricted open-shell Moller-Plesset perturbation theory.