Valence-only model potential calculations on copper hydride molecule

Different sets of one-electron functions obtained according to the strong-orthogonal geminal theory (GEM) [1], the Generalized Molecular Orbital (GMO) method [2] and the exchange maximization between virtual and occupied orbitals (EVO) [3], are tested as basis for CI calculations. The efficiency of the three procedures is discussed investigating the electronic structure of the CuH molecule using an effective-core potential. The values computed for the bond length, the dissociation energy and the vibrational frequency of the ground electronic state are compared with the experimental ones. The charge distribution is examined to estimate the contribution of the d electrons to the Cu-H bond. Comparisons are made with the results obtained by other theoretical works in which the copper atom is treated as a one valence electron atom.     

The performance of energy extrapolation procedures in truncated averaged coupled-pair functionals

Summary Energy extrapolation techniques in conjunction with individual configuration selection are applied to averaged coupled-pair functional expansions. In order to test the quality of this approach, benchmark calculations have been performed for N₂, the open and ring forms of O₃, and for the ground and several excited states of CuH and PdH. Reliable energy estimates are obtained for N₂ and the two transition metal hydrides and spectroscopic properties are in close agreement with the values for the non-truncated expansions. In the case of O₃ the perturbation corrections substantially underestimate the complete singles and doubles results. These deviations cancel to a large extent, however, in the calculated isomerization energy. The accuracy of the one-particle density matrix is examined by computing dipole moments for several electronic states of CuH and PdH. Deviations are significant in some cases. For the evaluation of properties the current approach requires modifications.Dedicated to Prof. W. Kutzelnigg on the occasion of his sixtieth birthday     

Ground state potential energy curve and dissociation energy of MgH

New high-resolution visible emission spectra of the MgH molecule have been recorded with high signal-to-noise ratios using a Fourier transform spectrometer. Many bands of the A 2Pi-->X 2Sigma+ and B' 2Sigma+-->X 2Sigma+ electronic transitions of 24MgH were analyzed; the new data span the v' = 0-3 levels of the A 2Pi and B'2Sigma+ excited states and the v'=0-11 levels of the X 2Sigma+ ground electronic state. The vibration-rotation energy levels of the perturbed A 2Pi and B' 2Sigma+ states were fitted as individual term values, while those of the X 2Sigma+ ground state were fitted using the direct-potential-fit approach. A new analytic potential energy function that imposes the theoretically correct attractive potential at long-range, and a radial Hamiltonian that includes the spin-rotation interaction were employed, and a significantly improved value for the ground state dissociation energy of MgH was obtained. The v'=11 level of the X 2Sigma+ ground electronic state was found to be the highest bound vibrational level of 24MgH, lying only about 13 cm(-1) below the dissociation asymptote. The equilibrium dissociation energy for the X 2Sigma+ ground state of 24MgH has been determined to be De=11104.7+/-0.5 cm(-1) (1.37681+/-0.00006 eV), whereas the zero-point energy (v'=0) is 739.11+/-0.01 cm(-1). The zero-point dissociation energy is therefore D0=10365.6+/-0.5 cm(-1) (1.28517+/-0.00006 eV). The uncertainty in the new experimental dissociation energy of MgH is more than 2 orders of magnitude smaller than that for the best value available in the literature. MgH is now the only hydride molecule other than H2 itself for which all bound vibrational levels of the ground electronic state are observed experimentally and for which the dissociation energy is determined with subwavenumber accuracy.     

Evaluating the variation of ion energy under different parameter settings in traveling wave ion mobility mass spectrometry

Ion mobility mass spectrometry (IM-MS) can be used to differentiate and identify isobaric ions. To improve IM-MS resolution, the second generation of traveling wave ion mobility (TWIM) technology was launched. There were reports showing ions were heated up by TWIM. With higher ion energy, it could alter the conformation of larger ions or MS/MS experiments. To monitor the energy exchange relating to the TWIM process, the combined use of thermometer ions with unique molecular structure and theoretical calculations to determine the effective temperature of ions had been explored. In this report, the use of a simple experimental approach to estimate the variation on the ion energy that result from changing a TWIM parameter setting is demonstrated. The approach aims to achieve the same percentage of ion dissociation in a collision cell, which is part of the original instrument and located at the exit of TWIM cell. Similar to the traditional MS/MS experiments, the same level of ion dissociation could be achieved by adjusting the electrical potential that was applied to the collision cell. The higher the ion energy after the TWIM separation, the lower the electrical potential was required to achieve the same level of ion dissociation. Together with the information on the number of electrical charge in the selected ion, the difference in the required electrical potentials could be converted into electron volt of ion energy that resulted from changing the TWIM parameter setting. The results showed ion energy could be changed 1–9 eV when the parameter of TWIM was adjusted.     

Single-Crystal to Single-Crystal Transformations: Stepwise CO2 Insertions into Bridging Hydrides of [(NHC)CuH]2 Complexes

Mechanistic studies of substrate insertion into dimeric [(NHC)CuH]₂ (NHC=N-heterocyclic carbene) complexes with two bridging hydrides have been shown to require dimer dissociation to generate transient, highly reactive (NHC)Cu−H monomers in solution. Using single-crystal to single-crystal (SC-SC) transformations, we discovered a new pathway of stepwise insertion of CO₂ into [(NHC)CuH]₂ without complete dissociation of the dimer. The first CO₂ insertion into dimeric [(IPr*OMe)CuH]₂ (IPr*OMe=N,N′-bis(2,6-bis(diphenylmethyl)-4-methoxy-phenyl)imidazole-2-ylidene) produced a dicopper formate hydride [(IPr*OMe)Cu]₂(μ-1,3-O₂CH)(μ-H). A second CO₂ insertion produced a dicopper bis(formate), [(IPr*OMe)Cu]₂(μ-1,3-O₂CH)(μ-1,1-O₂CH), containing two different bonding modes of the bridging formate. These dicopper formate complexes are inaccessible from solution reactions since the dicopper core cleanly ruptures to monomeric complexes when dissolved in a solvent.     

Translational to vibrational energy conversion during surface-induced dissociation of n-butylbenzene molecular ions colliding at self-assembled monolayer surfaces

Translational to vibrational (T-->V) energy conversion in the course of inelastic collisions of n-butylbenzene molecular ions with thiolate self-assembled monolayer (SAM) gold surfaces is studied to better understand internal energy uptake by the hyperthermal projectile ions. The projectile ion is selected by a mass spectrometer of BE configuration and product ions are analyzed using a quadrupole mass analyzer after kinetic energy selection with an electric sector. The branching ratio for formation of the fragment ions m/z 91 and m/z 92, measured over a range of collision energies, is used to estimate the average internal energy with the aid of calculations based on unimolecular dissociation kinetics [Rice-Ramsperger-Kassel-Marcus (RRKM) theory]. The measured T-->V conversion efficiencies (the fraction of the laboratory kinetic energy converted into internal energy) are 11 approximately 12% for dodecanethiolate SAM (H-SAM) and 19 approximately 20% for 2-perfluorooctylethanethiolate SAM (F-SAM), respectively, over ranges of a few 10s of eV. The values are similar to those reported earlier for other thermometer molecules undergoing surface collisions. Chemical sputtering leading to ionization of the surface is a prominent feature of the surface-induced dissociation (SID) spectra of n-butylbenzene acquired using the H-SAM surface but not the F-SAM surface because of the lower ionization energy of the former.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

Efficient calculation of low energy statistical rates for gas phase dissociation using umbrella sampling

Monte Carlo (MC) simulations can be used to compute microcanonical statistical rates of gas phase dissociation reactions. Unfortunately, the MC approach may suffer from a slow convergence and large statistical errors for energies just above the dissociation threshold. In this work, umbrella sampling is proposed as a device to reduce the statistical error of MC rate constants. The method is tested by computing the classical dissociation rate for the reaction [H5O2+]* --> H2O + H3O(+) over the range of internal energy 38 < E < or = 100 kcal/mol. Comparing with other literature methods, it is found that umbrella sampling reduces the computational effort by up to two orders of magnitude when used in conjunction with a careful choice of sampling distributions. The comparison between MC rate constants and classical Rice-Ramsperberg-Kassel harmonic theory shows that anharmonicity plays an important role in the dissociation process of the Zundel cation (H5O2+) at all energies.     

Use of a single trajectory to study product energy partitioning in unimolecular dissociation: mass effects for halogenated alkanes

A single trajectory (ST) direct dynamics approach is compared with quasiclassical trajectory (QCT) direct dynamics calculations for determining product energy partitioning in unimolecular dissociation. Three comparisons are made by simulating C(2)H(5)F-->HF + C(2)H(4) product energy partitioning for the MP26-31G(*) and MP26-311 + + G(**) potential energy surfaces (PESs) and using the MP26-31G(*) PES for C(2)H(5)F dissociation as a model to simulate CHCl(2)CCl(3)-->HCl + C(2)Cl(4) dissociation and its product energy partitioning. The trajectories are initiated at the transition state with fixed energy in reaction-coordinate translation E(t) (double dagger). The QCT simulations have zero-point energy (ZPE) in the vibrational modes orthogonal to the reaction coordinate, while there is no ZPE for the STs. A semiquantitative agreement is obtained between the ST and QCT average percent product energy partitionings. The ST approach is used to study mass effects for product energy partitioning in HX(X = F or Cl) elimination from halogenated alkanes by using the MP26-31G(*) PES for C(2)H(5)F dissociation and varying the masses of the C, H, and F atoms. There is, at most, only a small mass effect for partitioning of energy to HX vibration and rotation. In contrast, there are substantial mass effects for partitioning to relative translation and the polyatomic product's vibration and rotation. If the center of mass of the polyatomic product is located away from the C atom from which HX recoils, the polyatomic has substantial rotation energy. Polyatomic products, with heavy atoms such as Cl atoms replacing the H atoms, receive substantial vibration energy that is primarily transferred to the wag-bend motions. For E(t) (double dagger) of 1.0 kcalmol, the ST calculations give average percent partitionings to relative translation, polyatomic vibration, polyatomic rotation, HX vibration, and HX rotation of 74.9%, 6.8%, 1.5%, 14.4%, and 2.4% for C(2)H(5)F dissociation and 39.7%, 38.1%, 0.2%, 16.1%, and 5.9% for a model of CHCl(2)CCl(3) dissociation.     

Origins of enantioselectivity in the chiral diphosphine-ligated CuH-catalyzed asymmetric hydrosilylation of ketones

Computational investigations on the asymmetric hydrosilylation of acetophenone over ligated CuH catalysts were performed with the DFT method. The calculations predict that the catalytic reaction involves two steps: (1) CuH addition to the carbonyl group via a four-membered transition state (TS) with the formation of copper-alkoxide intermediates; (2) regeneration of the ligated CuH catalyst by an external SiH(4) through a metathesis process to yield the corresponding silyl ether. The calculations in the chiral diphosphine-ligated CuH systems suggest that the metathesis process is the rate-determining step (RDS). The CuH addition step is vital for the distribution of the racemic products and therefore represents the stereo-controlling step (SCT). In this step, the greater steric hindrance between the aromatic rings of the ligands and the substrate is identified as the major factor for enantioselectivity. The corresponding TS in the face-to-face mode, suffering less steric hindrance, is more stable than its analogue in the edge-to-face mode. The enantioselectivities are calculated to be related not only to the P-Cu-P bite angles in the stereo-controlling TSs, but also to the substituents at the P-aryl rings of the chiral ligands. In short, a larger P-Cu-P bite angle and suitably modified P-aryl rings together are necessary to achieve excellent ee values.     

A comparison of the on-top dissociation of H2 on Ni(100) and Cu(100)

The on-top dissociations of H₂ on Ni(100) and Cu(100) are studied using a cluster approach. Correlation effects are accounted for through the use of CASSCF and CCI methods. The central metal atom is treated with all its electrons whereas the other cluster atoms are described by recently developed one electron ECP's. A molecular chemisorbed H₂ state on nickel, similar to that recently observed experimentally, was identified in the cluster calculations and also for the triatomic NiH₂. No such state was found on copper. The large differences found for the on top dissociation of H₂ on nickel and copper are attributed solely to the difference in 3d orbital occupation. The parallel between the on top dissociation reaction on the cluster and the dissociation on a single atom is also studied. While the neutral triatomic NiH₂ represents a qualitatively correct model in the nickel case, the negatively charged CuH ₂ ^– is required as a model in the copper case.     

Recoil energy distributions for dissociation of the van der Waals molecule p-difluorobenzene-Ar with 450-3000 cm(-1) excess energy

Velocity map imaging has been used to measure the distributions of translational energy released in the dissociation of p-difluorobenzene-Ar van der Waals complexes from the 5(1), 3(1), 5(2), 3(1)5(1), 5(3), 3(2), and 3(2)5(1) states. These states span 818-3317 cm(-1) of vibrational energy and correspond to a range of energies above dissociation of 451-2950 cm(-1). The translational energy release (recoil energy) distributions are remarkably similar, peaking at very low energy (10-20 cm(-1)) and decaying in an exponential fashion to approach zero near 300 cm(-1). The average translational energy released is small, shows no dependence on the initial vibrational energy, and spans the range 58-72 cm(-1) for the vibrational levels probed. The average value for the seven levels studied is 63 cm(-1). The low fraction of transfer to translation is qualitatively in accord with Ewing's momentum gap model [G. E. Ewing, Faraday Discuss. 73, 325 (1982)]. No evidence is found in the distributions for a high energy tail, although it is likely that the experiment is not sufficiently sensitive to detect a low fraction of transfer at high translational energies. The average translational energy released is lower than has been seen in comparable systems dissociating from triplet and cation states.     

An application of the novel quantum mechanical/molecular mechanical method combined with the theory of energy representation: An ionic dissociation of a water molecule in the supercritical water

A novel quantum chemical approach recently developed has been applied to an ionic dissociation of a water molecule (2H(2)O-->H(3)O(+)+OH(-)) in ambient and supercritical water. The method is based on the quantum mechanical/molecular mechanical (QM/MM) simulations combined with the theory of energy representation (QM/MM-ER), where the energy distribution function of MM solvent molecules around a QM solute serves as a fundamental variable to determine the hydration free energy of the solute according to the rigorous framework of the theory of energy representation. The density dependence of the dissociation free energy in the supercritical water has been investigated for the density range from 0.1 to 0.6 g/cm(3) with the temperature fixed at a constant. It has been found that the product ionic species significantly stabilizes in the high density region as compared with the low density. Consequently, the dissociation free energy decreases monotonically as the density increases. The decomposition of the hydration free energy has revealed that the entropic term (-TDeltaS) strongly depends on the density of the solution and dominates the behavior of the dissociation free energy with respect to the variation of the density. The increase in the entropic term in the low density region can be attributed to the decrease in the translational degrees of freedom brought about by the aggregation of solvent water molecules around the ionic solute.     

Activation of Peptide ions by blackbody radiation: factors that lead to dissociation kinetics in the rapid energy exchange limit

Unimolecular rate constants for blackbody infrared radiative dissociation (BIRD) were calculated for the model protonated peptide (AlaGly)(n) (n = 2-32) using a variety of dissociation parameters. Combinations of dissociation threshold energies ranging from 0.8 to 1.7 eV and transition entropies corresponding to Arrhenius preexponential factors ranging from very "tight" (A(infinity) = 10(9.9) s(-1)) to "loose" (A(infinity) = 10(16.8) s(-1)) were selected to represent dissociation parameters within the experimental temperature range (300-520 K) and kinetic window (k(uni) = 0.001-0.20 s(-1)) typically used in the BIRD experiment. Arrhenius parameters were determined from the temperature dependence of these values and compared to those in the rapid energy exchange (REX) limit. In this limit, the internal energy of a population of ions is given by a Boltzmann distribution, and kinetics are the same as those in the traditional high-pressure limit. For a dissociation process to be in this limit, the rate of photon exchange between an ion and the vacuum chamber walls must be significantly greater than the dissociation rate. Kinetics rapidly approach the REX limit either as the molecular size or threshold dissociation energy increases or as the transition-state entropy or experimental temperature decreases. Under typical experimental conditions, peptide ions larger than 1.6 kDa should be in the REX limit. Smaller ions may also be in the REX limit depending on the value of the threshold dissociation energy and transition-state entropy. Either modeling or information about the dissociation mechanism must be known in order to confirm REX limit kinetics for these smaller ions. Three principal factors that lead to the size dependence of REX limit kinetics are identified. With increasing molecular size, rates of radiative absorption and emission increase, internal energy distributions become relatively narrower, and the microcanonical dissociation rate constants increase more slowly over the energy distribution of ions. Guidelines established here should make BIRD an even more reliable method to obtain information about dissociation energetics and mechanisms for intermediate size molecules.     

The nascent translational energy of difluorocarbene produced in IR MPD of the pentafluoroethyl radical

We report upon the direct detection of difluorocarbene following infrared multiphoton photolysis of pentafluoroethyl iodide using well-defined (SLM, TEM₀₀, 80 ns pulse width) TEA CO₂-laser pulses. The rate of appearance of CF₂ at 1 mTorr pressure and RRKM modelling of the unimolecular dissociation of C₂F₃I and C₂F₅ using reasonable input parameters are presented. These support a mechanism whereby CF₂ is produced by secondary photolysis of pentafluoroethyl radicals. Measurements of the velocity of CF₂ by the transient diffusion technique lead to an estimate of 2.6 kcal/mol for its average translational energy acquired from the homolytic cleavage of the CI and CC bonds. This value is higher than that predicted from the models using reasonable spontaneous dissociation rates ( = 10⁹ s^?1). An inherent assumption of the models is that the excess energy of dissociation is distributed statistically among the vibrational modes of the reaction complex and that there are no small barriers in the exit channel.     

Measuring internal energy deposition in collisional activation using hydrated ion nanocalorimetry to obtain peptide dissociation energies and entropies

The internal energy deposited in both on- and off-resonance collisional activation in Fourier transform ion cyclotron resonance mass spectrometry is measured with ion nanocalorimetry and is used to obtain information about the dissociation energy and entropy of a protonated peptide. Activation of Na⁺(H₂O)₃₀ results in sequential loss of water molecules, and the internal energy of the activated ion can be obtained from the abundances of the product ions. Information about internal energy deposition in on-resonance collisional activation of protonated peptides is inferred from dissociation data obtained under identical conditions for hydrated ions that have similar m/z and degrees-of-freedom. From experimental internal energy deposition curves and Rice-Ramsperger-Kassel-Marcus (RRKM) theory, dissociation data as a function of collision energy for protonated leucine enkephalin, which has a comparable m/z and degrees-of-freedom as Na⁺(H₂O)₃₀, are modeled. The threshold dissociation energies and entropies are correlated for data acquired at a single time point, resulting in a relatively wide range of threshold dissociation energies (1.1 to 1.7 eV) that can fit these data. However, this range of values could be significantly reduced by fitting data acquired at different dissociation times. By measuring the internal energy of an activated ion, the number of fitting parameters necessary to obtain information about the dissociation parameters by modeling these data is reduced and could result in improved accuracy for such methods.     

Analysis of the HO2 vibrational spectrum on an accurate ab initio potential energy surface

The complete vibrational spectrum of the HO2(X(2)A' ') radical, up to the H + O2 dissociation limit, has been determined quantum mechanically on an accurate potential energy surface (PES), based on approximately 15000 ab initio points at the icMRCI+Q/aug-cc-pVQZ level of theory. The vibrational states are found to be assignable at low energies but become more irregular as the energy approaches the dissociation limit. However, even at very high energies, regularity still exists, in sharp contrast to earlier results based on the double many-body expansion (DMBE) IV potential. Several Fermi resonances have been identified, and the spectrum is fit with a spectroscopic Hamiltonian. In addition, the vibrational dynamics is analyzed using a periodic orbit approach.     

An ab initio potential energy surface and vibrational energy levels of HXeBr

A three-dimensional global potential energy surface for the electronic ground state of HXeBr molecule is constructed from more than 4200 ab initio points. These points are generated using an internally contracted multi-reference configuration interaction method with the Davidson correction （icMRCI ＋ Q） and large basis sets. The stabilities and dissociation barriers are identified from the potential energy surfaces. The three-body dissociation channel is found to be the dominate dissociation channel for HXeBr. Based on the obtained potentials, low-lying vibrational energy levels of HXeBr calculated using the Lanczos algorithm is found to be in good agreement with the available experimental band origins.