Thev1andv3bands of ND3

Wavepacket dynamical study of H-atom tunneling in catecholate monoanion: the role of intermode couplings and energy flow

We present a study of H-atom tunneling in catecholate monoanion through wavepacket dynamical simulations. In our earlier study of this symmetrical double-well system [Phys. Chem. Chem. Phys., 2022, 24, 10887], a limited number of transition state modes were identified as being important for the tunneling process. These include the imaginary frequency mode Q₁, the CO scissor mode Q₁₀, and the OHO bending mode Q₂₉. In this work, starting from non-stationary initial states prepared with excitations in these modes, we have carried out wavepacket dynamics in two and three dimensional spaces. We analyse the dynamical effects of the intermode couplings, in particular the role of energy flow between the studied modes on H-atom tunneling. We find that while Q₁₀ strongly modulates the donor-acceptor distance, it does not exchange energy with Q₁. However, excitation in Q₂₉ or Q₁ does lead to rapid energy exchange between these modes, which modifies the tunneling rate at early times.

Multipole-moment effects in ion-molecule reactions at low temperatures: part I - ion-dipole enhancement of the rate coefficients of the He+ + NH3 and He+ + ND3 reactions at collisional energies Ecoll/kB near 0 K

The energy dependence of the rates of the reactions between He+ and ammonia (NY3, Y = {H,D}), forming NY2+, Y and He as well as NY+, Y2 and He, and the corresponding product branching ratios have been measured at low collision energies Ecoll between 0 and kB·40 K using a recently developed merged-beam technique [Allmendinger et al., ChemPhysChem, 2016, 17, 3596]. To avoid heating of the ions by stray electric fields, the reactions are observed within the large orbit of a highly excited Rydberg electron. A beam of He Rydberg atoms was merged with a supersonic beam of ammonia using a curved surface-electrode Rydberg-Stark deflector, which is also used for adjusting the final velocity of the He Rydberg atoms, and thus the collision energy. A collision-energy resolution of about 200 mK was reached at the lowest Ecoll values. The reaction rate coefficients exhibit a sharp increase at collision energies below ∼kB·5 K and pronounced deviations from Langevin-capture behaviour. The experimental results are interpreted in terms of an adiabatic capture model describing the rotational-state-dependent orientation of the ammonia molecules by the electric field of the He+ atom. The model faithfully describes the experimental observations and enables the identification of three classes of |JKMp〉 rotational states of the ammonia molecules showing different low-energy capture behaviour: (A) high-field-seeking states with |KM| ≥ 1 correlating to the lower component of the umbrella-motion tunnelling doublet at low fields. These states undergo a negative linear Stark shift, which leads to strongly enhanced rate coefficients; (B) high-field-seeking states subject to a quadratic Stark shift at low fields and which exhibit only weak rate enhancements; and (C) low-field-seeking states with |KM| ≥ 1. These states exhibit a positive Stark shift at low fields, which completely suppresses the reactions at low collision energies. Marked differences in the low-energy reactivity of NH3 and ND3-the rate enhancements in ND3 are more pronounced than in NH3-are quantitatively explained by the model. They result from the reduced magnitudes of the tunnelling splitting and rotational intervals in ND3 and the different occupations of the rotational levels in the supersonic beam caused by the different nuclear-spin statistical weights. Thermal capture rate constants are derived from the model for the temperature range between 0 and 10 K relevant for astrochemistry. Comparison of the calculated thermal capture rate coefficients with the absolute reaction rates measured above 27 K by Marquette et al. (Chem. Phys. Lett., 1985, 122, 431) suggests that only 40% of the close collisions are reactive.

Isotope Labelling for Reaction Mechanism Analysis in DBD Plasma Processes

Dielectric barrier discharge (DBD) plasmas and plasma catalysis are becoming an alternative procedure to activate various gas phase reactions. A low-temperature and normal operating pressure are the main advantages of these processes, but a limited energy efficiency and little selectivity control hinder their practical implementation. In this work, we propose the use of isotope labelling to retrieve information about the intermediate reactions that may intervene during the DBD processes contributing to a decrease in their energy efficiency. The results are shown for the wet reforming reaction of methane, using D2O instead of H2O as reactant, and for the ammonia synthesis, using NH3/D2/N2 mixtures. In the two cases, it was found that a significant amount of outlet gas molecules, either reactants or products, have deuterium in their structure (e.g., HD for hydrogen, CDxHy for methane, or NDxHy for ammonia). From the analysis of the evolution of the labelled molecules as a function of power, useful information has been obtained about the exchange events of H by D atoms (or vice versa) between the plasma intermediate species. An evaluation of the number of these events revealed a significant progression with the plasma power, a tendency that is recognized to be detrimental for the energy efficiency of reactant to product transformation. The labelling technique is proposed as a useful approach for the analysis of plasma reaction mechanisms.

Infrared Chemiluminescence Study of the Reaction of Hydroxyl Radical with Formamide and the Secondary Unimolecular Reaction of Chemically Activated Carbamic Acid

Reactions of OH and OD radicals with NH₂CHO and ND₂CHO were studied by Fourier transform infrared emission spectroscopy of the product molecules from a fast-flow reactor at 298 K. Vibrational distributions of the HOD and H₂O molecules from the primary reactions with the C-H bond were obtained by computer simulation of the emission spectra. The vibrational distributions resemble those for other direct H atom abstraction reactions, such as with acetaldehyde. The highest observed level gives an estimate of the C-H bond dissociation energy in formamide of 90.5 ± 1.3 kcal mol^-1. Observation of CO₂, ammonia, and secondary water chemiluminescence gave evidence that recombination of OH and NH₂CO forms carbamic acid (NH₂COOH) with excitation energy of 103 kcal mol^-1, which decomposes through two pathways forming either NH₃ + CO₂ or H₂O + HNCO. The branching fraction for ammonia formation was estimated to be 2-3 times larger than formation of water. This observation was confirmed by RRKM calculation of the decomposition rate constants. A new simulation method was developed to analyze infrared emission from NH₃, NH₂D, ND₂H, and ND₃. Dynamical aspects of the primary and secondary reactions are discussed based on the vibrational distributions of CO₂ and those of H/D isotopes of water and ammonia.

Promoting and inhibiting tunneling via nuclear motions

Accurate, experimental rotational-vibrational energy levels determined via the MARVEL (Measured Active Rotational-Vibrational Energy Levels) algorithm and published recently for the symmetric-top (14)NH3 molecule in J. Quant. Spectrosc. Radiat. Transfer, 2015, 116, 117-130 are analyzed to unravel the promoting and inhibiting effects of vibrations and rotations on the tunneling splittings of the corresponding symmetric (s) and antisymmetric (a) rovibrational energy level pairs. The experimental transition data useful from the point of view of the present analysis cover the range 0.7-7000 cm(-1), sufficiently detailed rovibrational energy sets worth analyzing are available for 20 vibrational bands. The highest J value, where J stands for the rotational quantum number, within the experimental dataset employed is 30. Coupling of the "umbrella" motion of (14)NH3 with other vibrational degrees of freedom has only a minor effect on the a-s tunneling splitting characterizing the ground vibrational state, 0.79436(70) cm(-1). In the majority of the cases rotation around the C3 axis increases, while rotation around the two perpendicular axes decreases the tunneling splittings. For example, for the pair of vibrational ground states, 0(+) and 0(-), the tunneling splitting basically disappears at around J = 25 for the (J,K) = (J,1) states, where K = |k| is the usual quantum number characterizing the projection of the rotational angular momentum on the principal axis. The tunneling splittings, defined as energy differences E(a) - E(s) of corresponding energy level pairs, as a function of J and K show a very regular behavior for the ground state (GS) and the nν2 bands. For the other bands investigated exceptions from a regular behavior do occur, especially for bands characterized by degenerate vibrations, and occasionally the data available are not sufficient to arrive at definitive conclusions. The most irregular behavior is observed for rotational states characterized by the k - l = 3n rule (l is the vibrational angular momentum quantum number), with n = 0, 1, 2,… High-quality, variationally computed rovibrational data support all the conclusions of this study based on experimental energy levels.

Paramagnetic Nanoparticles Leave Their Mark on Nuclear Spins of Transiently Adsorbed Proteins

The successful application of nanomaterials in biosciences necessitates an in-depth understanding of how they interface with biomolecules. Transient associations of proteins with nanoparticles (NPs) are accessible by solution NMR spectroscopy, albeit with some limitations. The incorporation of paramagnetic centers into NPs offers new opportunities to explore bio-nano interfaces. We propose NMR paramagnetic relaxation enhancement as a new tool to detect NP-binding surfaces on proteins with increased sensitivity, also extending the applicability of NMR investigations to heterogeneous biomolecular mixtures. The adsorption of ubiquitin on gadolinium-doped fluoride-based NPs produced residue-specific NMR line-broadening effects mapping to a contiguous area on the surface of the protein. Importantly, an identical paramagnetic fingerprint was observed in the presence of a competing protein-protein association equilibrium, exemplifying possible interactions taking place in crowded biological media. The interaction was further characterized using isothermal titration calorimetry and upconversion emission measurements. The data indicate that the used fluoride-based NPs are not biologically inert but rather are capable of biomolecular recognition.

Collision dynamics of symmetric top molecules: a comparison of the rotationally inelastic scattering of CD3 and ND3 with He

We compare rotationally inelastic scattering of deuterated methyl radicals (CD3) and ammonia (ND3) in collisions with helium using close-coupling quantum-mechanical scattering calculations performed with ab initio potential energy surfaces (PESs). The theoretical methods have been rigorously tested against angle-resolved experimental measurements obtained using crossed molecular beam apparatuses in combination with velocity map imaging [O. Tkáč, A. G. Sage, S. J. Greaves, A. J. Orr-Ewing, P. J. Dagdigian, Q. Ma, and M. H. Alexander, Chem. Sci. 4, 4199 (2013); O. Tkáč, A. K. Saha, J. Onvlee, C.-H. Yang, G. Sarma, C. K. Bishwakarma, S. Y. T. van de Meerakker, A. van der Avoird, D. H. Parker, and A. J. Orr-Ewing, Phys. Chem. Chem. Phys. 16, 477 (2014)]. Common features of the scattering dynamics of these two symmetric top molecules, one closed-shell and the other an open-shell radical, are identified and discussed. Two types of anisotropies in the PES influence the interaction of an atom with a nonlinear polyatomic molecule. The effects of these anisotropies can be clearly seen in the state-to-state integral cross sections out of the lowest CD3 rotational levels of each nuclear spin symmetry at a collision energy of 440 cm(-1). Similarities and differences in the differential cross sections for the ND3-He and CD3-He systems can be linked to the coupling terms derived from the PESs which govern particular initial to final rotational level transitions.

Matrix isolation spectroscopy and nuclear spin conversion of NH3 and ND3 in solid parahydrogen

We present matrix isolation infrared absorption spectra of NH3 and ND3 trapped in solid parahydrogen (pH2) at temperatures around 1.8 K. We used the relatively slow nuclear spin conversion (NSC) of NH3 and ND3 in freshly deposited pH2 samples as a tool to assign the sparse vibration-inversion-rotation (VIR) spectra of NH3 in the regions of the ν2, ν4, 2ν4, ν1, and ν3 bands and ND3 in the regions of the ν2, ν4, ν1, and ν3 fundamentals. Partial assignments are also presented for various combination bands of NH3. Detailed analysis of the ν2 bands of NH3 and ND3 indicates that both isotopomers are nearly free rotors; that the vibrational energy is blue-shifted by 1-2%; and that the rotational constants and inversion tunneling splitting are 91-94% and 67-75%, respectively, of the gas-phase values. The line shapes of the VIR absorptions are narrow (0.2-0.4 cm(-1)) for upper states that cannot rotationally relax and broad (>1 cm(-1)) for upper states that can rotationally relax. We report and assign a number of NH3-induced infrared absorption features of the pH2 host near 4150 cm(-1), along with a cooperative transition that involves simultaneous vibrational excitation of a pH2 molecule and rotation-inversion excitation of NH3. The NSCs of NH3 and ND3 were found to follow first-order kinetics with rate constants at 1.8 K of k = 1.88(16) × 10(-3) s(-1) and k = 1.08(8) × 10(-3) s(-1), respectively. These measured rate constants are compared to previous measurements for NH3 in an Ar matrix and with the rate constants measured for other dopant molecules isolated in solid pH2.

Ab initio-based double many-body expansion potential energy surface for the first excited triplet state of the ammonia molecule

A global single-sheeted double many-body expansion potential energy surface is reported for the first excited triplet state of NH(3). It employs an approximate cluster expansion of the molecular potential that utilizes previously reported functions of the same family for the triatomic fragments. Four-body energy terms have been calibrated from extensive accurate ab initio data so as to reproduce the main features of the title system. A new switching function formalism has been reported to approximate the true multisheeted nature of NH(3)((3)A(2) ('')) potential energy surface, thus allowing the correct behavior at the NH(2)((2)A(")) + H((2)S) and NH(2)((4)A(")) + H((2)S) dissociation limits. The resulting fully six-dimensional potential energy function reproduces the correct symmetry under the permutation of identical atoms, and predicts the correct behavior at all dissociation channels while providing a realistic representation at all interatomic separations. The major attributes of the NH(3) double many-body expansion potential energy surface have also been characterized, and found to be in good agreement, both with the calculated ones from the raw ab initio energies and the theoretical results available in the literature.

Mode selective tunneling dynamics observed by high resolution spectroscopy of the bending fundamentals of N14H2D and N14D2H

High resolution (0.004 and 0.01 cm(-1) instrumental bandwidth) interferometric Fourier transform infrared spectra of (14)NH2D and (14)ND2H were measured on a Bomem DA002 spectrometer in a supersonic jet expansion and at room temperature. We report the analysis of the bending fundamentals of (14)NH2D with term values Tv(s)=1389.9063(2) cm(-1) and Tv(a)=1390.4953(2) cm(-1) for the nu(4b) fundamental and Tv(s)=1605.6404(7) cm(-1) and Tv(a)=1591.0019(7) cm(-1) for the nu(4a) fundamental, and of (14)ND2H with term values of Tv(s)=1233.3740(2) cm(-1) and Tv(a)=1235.8904(2) cm(-1) for the nu(4a) fundamental and Tv(s)=1461.7941(9) cm(-1) and Tv(a)=1461.9918(19) cm(-1) for the nu(4b) fundamental. In all cases Tv(s) gives the position of the symmetric inversion sublevel (with positive parity) and Tv(a) the position of the antisymmetric inversion sublevel (with negative parity). The notation for the fundamentals nu(4a) and nu(4b) is chosen by correlation with the degenerate nu(4) mode in the C(3v) symmetric molecules NH3 and ND3. The degeneracy is lifted in Cs symmetry and a indicates the symmetric, b the antisymmetric normal mode with respect to the Cs symmetry plane in NH2D and ND2H. Assignments were established with certainty by means of ground state combination differences. About 20 molecular parameters of the effective S-reduced Hamiltonian could be determined accurately for each fundamental. In particular, the effect of Fermi resonances of the 2nu(2) overtone with the nu(4a) bending mode was observed, leading to an increased inversion splitting in the case of ND2H and to a strongly increased inversion splitting and an inverted order of the two inversion levels in NH2D. Rotational perturbations observed with the nu(4b) bending fundamentals are probably due to Coriolis interactions with the inversion overtone 2nu(2). The results are important for understanding isotope effects on the inversion in ammonia as well as its selective catalysis and inhibition by excitation of different vibrational modes, as treated by quantum dynamics on high dimensional potential hypersurfaces of this molecule.

Global Analytical Potential Energy Surface for Large Amplitude Nuclear Motions in Ammonia

An analytical, full-dimensional, and global representation of the potential energy surface of NH(3) in the lowest adiabatic electronic state is developed, and parameters are determined by adjustment to ab initio data and thermochemical data for several low-lying dissociation channels. The electronic structure is calculated at the CASPT2 level within an [8,7] active space. The representation is compared to other recently published potential energy surfaces for this molecule. The present representation is distinguished by giving a qualitatively correct description of the potential energy for very large amplitude displacements of the nuclei from equilibrium. Other characteristic features of the present surface are the equilibrium geometries r(eq)(NH(3)) approximately 101.24 pm, r(eq)(NH(2)) approximately 102.60 pm, alpha(eq)(NH(3)) approximately 106.67 degrees, and the inversion barrier at r(inv)(NH(3)) approximately 99.80 pm and 1781 cm(-1) above the NH(3) minimum. The barrier to linearity in NH(2) is 11,914 cm(-1) above the NH(2)((2)B(1)) minimum. While the quartic force field for NH(3) from the present representation is significantly different from that of the other potential energy surfaces, the vibrational structures obtained from perturbation theory are quite similar for all representations studied here.