A perturbation theory of isoelectronic molecules - Application to CO based on N sub 2.
Perturbation theory of heteronuclear diatomic molecules based on isoelectronic homonuclear molecules applied to carbon monoxide and nitrogen
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Perturbation theory of heteronuclear diatomic molecules based on isoelectronic homonuclear molecules applied to carbon monoxide and nitrogen
The universal energy relation, discovered for metallic and covalent solids as well as nuclear matter, is tested for diatomic molecules. It is found that it applies well to covalent diatomic bonds, but that ionic diatomic bonds are in a distinct class. A simple extension of the universal binding energy relation that includes the effects of ionicity ensues. It yields accurate prediction of spectroscopic data for both ionic and covalent bonds in 150 molecules. The form of the covalent part is given by the universal relation, suggesting an intimate relationship between the energetics of solids and diatomic molecules.
Wave functions and total energy of energy states of diatomic molecules - Hartree-Fock calculation
Cusp conditions to describe behavior of wave function at singularities of Coulomb potential corresponding to coalescence of two or more particles - diatomic molecules
Behavior of wave function for system of charged particles near coalescence of any two of them, noting electron-nucleus cusp conditions for diatomic molecules
Potential energy curves for bound states of diatomic molecules from spectroscopic constants, using Rydberg-Klein-Rees method
Transport coefficient expressions for collisions between rigid diatomic molecules, using perturbation method
Electron scattering from diatomic molecules in fixed nucleus approximation using polarized single-center orbitals method applied to positively ionized hydrogen
The main data need for stellar atmosphere and spectrum modeling remains atomic and molecular transition data, particularly energy levels and transition cross-sections. We emphasize that data is needed for bound-free (b - f) as well as bound-bound (b - b), and collisional as well as radiative transitions. Data is now needed for polyatomic molecules as well as atoms, ions, and diatomic molecules. In addition, data for the formation of, and extinction due to, liquid and solid phase dust grains is needed. A prioritization of species and data types is presented, and gives emphasis to Fe group elements, and elements important for the investigation of nucleosynthesis and Galactic chemical evolution, such as the -elements and n-capture elements. Special data needs for topical problems in the modeling of cool stars and brown dwarfs are described.
Many liquid phase studies assume that the potential energy surfaces of reacting molecules are the same as in the gas phase, neglecting complex solvent dynamics that can completely alter the nature of chemical reactivity. Even studies that include solvent effects typically only consider them in an average, equilibrium way as part of a potential of mean force (PMF). In this work, we use mixed quantum/classical simulations to compare how equilibrium and non-equilibrium solvent motions affect the photodissociation of a simple diatomic molecule, NaK + , in liquid tetrahydrofuran. A PMF analysis shows that as the excited-state molecule dissociates with the solvent at equilibrium, the bonding electron remains associated with K + at short bond distances but eventually localizes on Na + at the end of dissociation. When we examine non-equilibrium dynamical photodissociation trajectories, however, we find that they fall into three distinct categories: about a quarter of them have the bonding electron mainly associated with Na + , another quarter stay mainly associated with K + , and about half have the bonding electron shared roughly equally between the two ions. The results show that equilibrium PMFs cannot accurately describe the dynamics of bond-breaking chemical reactions in solution because there is insufficient time for the solvent to reach equilibrium on the time scale over which bond dissociation occurs. Furthermore, our analysis shows that the solvent coupling between the electronic energy surfaces is similar at and away from equilibrium, suggesting that other factors, such as solute velocity-driven solvent memory effects, play a more important role in explaining the failure of the equilibrium PMF to predict the non-equilibrium dynamics.
This software supports the paper, "Exploring the scaling limitations of the variational quantum eigensolver with the bond dissociation of hydride diatomic molecules," published in the International Journal of Quantum Chemistry, whose abstract is as follows: Materials simulations involving strongly correlated electrons pose fundamental challenges to state-of-the-art electronic structure methods but are hypothesized to be the ideal use case for quantum computing. To date, no quantum computer has simulated a molecule of a size and complexity relevant to real-world applications, despite the fact that the variational quantum eigensolver (VQE) algorithm can predict chemically accurate total energies. Nevertheless, because of the many applications of moderately-sized, strongly correlated systems, such as molecular catalysts, the successful use of the VQE stands as an important waypoint in the advancement toward useful chemical modeling on near-term quantum processors. In this paper, we take a significant step in this direction. We lay out the steps, write, and run parallel code for an (emulated) quantum computer to compute the bond dissociation curves of the TiH, LiH, NaH, and KH diatomic hydride molecules using VQE. TiH was chosen as a relatively simple chemical system that incorporates d orbitals and strong electron correlation. Because current VQE implementations on existing quantum hardware are limited by qubit error rates, the number of qubits available, and the allowable gate depth, recent studies have focused on chemical systems involving s and p block elements. Through VQE + UCCSD calculations of TiH, we evaluate the near-term feasibility of modeling a molecule with d-orbitals on real quantum hardware. We demonstrate that the inclusion of d-orbitals and the use of the UCCSD ansatz, which are both necessary to capture the correct TiH physics, dramatically increase the cost of this problem. We estimate the approximate error rates necessary to model TiH on current quantum computing hardware using VQE+UCCSD and show them to likely be prohibitive until significant improvements in hardware and error correction algorithms are available.
The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which have made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here, we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Furthermore, our algorithm extends the symmetry driven structure search method by adding constraints for the bond length and the number of atoms in a molecule while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range of 10 GPa–4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm3¯m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P42/mmc symmetry at 2500 GPa, and finally to its cubic analog form with Pm3¯n symmetry at 3000 GPa.
The adiabatic-nuclei approximation of vibrational-rotational excitation of homonuclear diatomic molecules can be simply augmented to describe the vibrational-rotational coupling by including the dependence of the vibrational wave function on j. Appropriate formulas are given, and the theory, is applied to e-H2 excitation, whereby it is shown that deviations from the simple Born-Oppenheimer approximation measured by Wong and Schultz can be explained. More important, it can be seen that the inclusion of the j-dependent centrifugal term is essential for transitions involving high-rotational quantum numbers.
The basis and techniques are presented for generating opacity probability distribution functions for the CN molecule (red and violet systems) and the C2 molecule (Swan, Phillips, Ballik-Ramsay systems), two of the more important diatomic molecules in the spectra of carbon stars, with a view to including these distribution functions in equilibrium model atmosphere calculations. Comparisons to the CO molecule are also shown. T he computation of the monochromatic absorption coefficient uses the most recent molecular data with revision of the oscillator strengths for some of the band systems. The total molecular stellar mass absorption coefficient is established through fifteen equations of molecular dissociation equilibrium to relate the distribution functions to each other on a per gram of stellar material basis.
The Langevin equation was used to explore an adsorbate desorption mechanism. Calculations were performed using iterative extended Hueckel on a silica model site with various small adsorbates, e.g., H, CH, OH, NO, CO. It was found that barriers to free traversal from one site to another are substantial (about 3-10 eV). A bootstrap desorption mechanism for some molecules in the process of forming at a site also became apparent from the calculations. The desorption mechanisms appear to be somewhat balanced by a counterforce - the attraction of sites for the newly desorbed molecule. The order of attraction to a silica grain site for the diatomic molecules considered was OH greater than CH greater than CO greater than NO, when these entities were sufficiently distant. The nature of the silica grain and that of the 'cold' desorption mechanism, when considered together, suggest that the abundance of very small grains might be less common than anticipated.
The Langevin equation was used to explore an adsorbate desorption mechanism. Calculations were performed using iterative extended Huckel on a silica model site with various small adsorbates, e.g., H, CH, OH, NO, CO. It was found that barriers to free traversal from one site to another are substantial (approximately 3 - 10 eV). A bootstrap desorption mechanism for some molecules in the process of forming at a site also became apparent from the calculations. The desorption mechanisms appear to be somewhat balanced by a counterforce--the attraction of sites for the newly desorbed molecule. The order of attraction to a silica grain site for the diatomic molecules considered was OH > CH > CO > NO, when these entities were sufficiently distant. The nature of the silica grain and that of the "cold" desorption mechanism, when considered together, suggest that the abundance of very small grains might be less common than anticipated.
A general approach is presented for calculating the reduction of the dissociation energy of diatomic molecules immersed in a dense (n = less than 10 exp 22/cu cm) gas of molecules and atoms. The dissociation energy of a molecule in a dense gas differs from that of the molecule in vacuum because the intermolecular forces change the intramolecular dynamics of the molecule, and, consequently, the energy of the molecular bond.
The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H – and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.