Search

This paper demonstrates the performance of our previously suggested property-energy consistent method on the example of the generation of effective basis sets, pecS-1 and pecS-2, suited for the calculation of hydrogen, carbon, nitrogen, and oxygen chemical shifts. The new basis sets were successfully approbated in the GIAO-DFT calculations of the chemical shifts of 35 molecules using six different functionals. The pecS-1 basis set demonstrated very good accuracy, which makes this small basis set an effective means for the large-scale computations. At the same time, the pecS-2 basis set also gave very accurate results, thus putting it on a par with the other commensurate basis sets suited for the chemical shifts calculations. [ABSTRACT FROM AUTHOR]

The quest for "chemical accuracy" is becoming more and more demanded in the field of structure and kinetics of molecules at solid surfaces. In this paper, as an example, we focus on the barrier for hydrogen diffusion on a α-Al2O3(0001) surface, aiming for a couple cluster singles, doubles, and perturbative triples [CCSD(T)]-level benchmark. We employ the density functional theory (DFT) optimized minimum and transition state structures reported by Heiden, Usvyat, and Saalfrank [J. Phys. Chem. C 123, 6675 (2019)]. The barrier is first evaluated at the periodic Hartree–Fock and local Møller–Plesset second-order perturbation (MP2) level of theory. The possible sources of errors are then analyzed, which includes basis set incompleteness error, frozen core, density fitting, local approximation errors, as well as the MP2 method error. Using periodic and embedded fragment models, corrections to these errors are evaluated. In particular, two corrections are found to be non-negligible (both from the chemical accuracy perspective and at the scale of the barrier value of 0.72 eV): the correction to the frozen core-approximation of 0.06 eV and the CCSD(T) correction of 0.07 eV. Our correlated wave function results are compared to barriers obtained from DFT. Among the tested DFT functionals, the best performing for this barrier is B3LYP-D3. [ABSTRACT FROM AUTHOR]

In this paper, equilibrium properties of structure II hydrates of hydrogen were determined from Monte Carlo simulations in the isothermal-isobaric Gibbs ensemble. Water and hydrogen molecules are described by the TIP4P/Ice and Silvera-Goldman models, respectively. The use of the Gibbs ensemble has many key advantages for the simulation of hydrates. By the separation of hydrogen vapor and hydrate phases into their own domains, coupled with transfer moves of hydrogen molecules between domains, cage occupancies were determined. Furthermore, the choice of this ensemble also allows equilibrium lattice constants and guest molecule chemical potentials to be straightforwardly estimated. Results for hydrogen mass fractions indicate reasonable agreement with prior simulation data and theoretical models, while detailed analysis of cage occupancy distributions and neighboring cage pair occupancy combinations gives valuable insight into the behavior of this hydrate at the inter-cage scale. These results will aid in the construction of theoretical models, for which knowledge of the occupancy of neighboring cages is of great importance. In support of previous experimental and theoretical works, we also find evidence of double occupancy of a few small cages inside of the hydrate stability zone, albeit at very high pressures; approximately 0.1% of small cages are doubly occupied at 300 MPa, for temperatures of 225 K and 250 K. [ABSTRACT FROM AUTHOR]

The manipulation of spin states at the molecular scale is of fundamental importance for the development of molecular spintronic devices. One of the feasible approaches for the modification of a molecular spin state is through the adsorption of certain specific atoms or molecules including H, NO, CO, NH3, and O2. In this paper, we demonstrate that the local spin state of an individual iron phthalocyanine (FePc) molecule adsorbed on an Au(111) surface exhibits controllable switching by hydrogen adsorption, as evidenced by using first-principles calculations based on density functional theory. Our theoretical calculations indicate that different numbers of hydrogen adsorbed at the pyridinic N sites of the FePc molecule largely modify the structural and electronic properties of the FePc/Au(111) composite by forming extra N-H bonds. In particular, the adsorption of one or up to three hydrogen atoms induces a redistribution of charge (spin) density within the FePc molecule, and hence a switching to a lowspin state (S = 1/2) from an intermediate spin state (S = 1) is achieved, while the adsorption of four hydrogen atoms distorts the molecular conformation by increasing Fe-N bond lengths in FePc and thus breaks the ligand field exerted on the Fe 3d orbitals via stronger hybridization with the substrate, leading to an opposite switching to a high-spin state (S = 2). These findings obtained from the theoretical simulations could be useful for experimental manipulation or design of single-molecule spintronic devices. [ABSTRACT FROM AUTHOR]

We extend the statistical associating fluid theory of quantum corrected Mie potentials (SAFT-VRQ Mie), previously developed for pure fluids [Aasen et al., J. Chem. Phys. 151, 064508 (2019)], to fluid mixtures. In this model, particles interact via Mie potentials with Feynman–Hibbs quantum corrections of first order (Mie-FH1) or second order (Mie-FH2). This is done using a third-order Barker–Henderson expansion of the Helmholtz energy from a non-additive hard-sphere reference system. We survey existing experimental measurements and ab initio calculations of thermodynamic properties of mixtures of neon, helium, deuterium, and hydrogen and use them to optimize the Mie-FH1 and Mie-FH2 force fields for binary interactions. Simulations employing the optimized force fields are shown to follow the experimental results closely over the entire phase envelopes. SAFT-VRQ Mie reproduces results from simulations employing these force fields, with the exception of near-critical states for mixtures containing helium. This breakdown is explained in terms of the extremely low dispersive energy of helium and the challenges inherent in current implementations of the Barker–Henderson expansion for mixtures. The interaction parameters of two cubic equations of state (Soave–Redlich–Kwong and Peng–Robinson) are also fitted to experiments and used as performance benchmarks. There are large gaps in the ranges and properties that have been experimentally measured for these systems, making the force fields presented especially useful. [ABSTRACT FROM AUTHOR]

Chemical waves that arise in the H2 + O2 reaction on a bimetallic Rh(111)/Ni surface have been studied in the 10-6 and 10-5 mbar range at T = 773 K with photoelectron emission microscopy (PEEM), low-energy electron diffraction (LEED), and Auger electron spectroscopy (AES). Nickel coverages of 0.3, 0.6, and 1.0 monolayers were investigated. Coadsorbed with some oxygen, Ni starts to penetrate the Rh bulk region substantially only beyond 900 K. In the 10-5 mbar range, chemical waves are characterized by three distinct gray levels in PEEM. This number reduces to only two levels in the 10-6 mbar range. In situ LEED showed the periodic appearance of a (nx1) (n = 8, 10) pattern during chemical waves which was assigned to a 2D-Ni oxide. With in situ AES, one observes that the bright phase in PEEM correlates with a high Ni coverage and the dark phase with a low Ni coverage. [ABSTRACT FROM AUTHOR]

Chemical waves in the H2 + O2 reaction on a Rh(111) surface alloyed with Ni [ ΘNi < 1.5 monolayers (ML)] have been investigated in the 10-7 and 10-6 mbar range at T = 773 K using scanning photoelectron microscopy and x-ray photoelectron spectroscopy as in situ methods. The local intensity variations of the O 1s and the Ni 2p signal display an anticorrelated behavior. The coincidence of a high oxygen signal with a low Ni 2p intensity, which seemingly contradicts the chemical attraction between O and Ni, has been explained with a phase separation of the oxygen covered Rh(111)/Ni surface into a 3D-Ni oxide and into a Ni poor metallic phase. Macroscopic NiO islands (≈1 μm size) formed under reaction conditions have been identified as 2D-Ni oxide. Titration experiments of the oxygen covered Rh(111)/Ni surface with H2 demonstrated that the reactivity of oxygen is decreased by an order of magnitude through the addition of 0.6 ML Ni. An excitation mechanism is proposed in which the periodic formation and reduction of NiO modulate the catalytic activity. [ABSTRACT FROM AUTHOR]

The hydrogen transfer reaction catalysed by soybean lipoxygenase (SLO) has been the focus of intense study following observations of a high kinetic isotope effect (KIE). Today high KIEs are generally thought to indicate departure from classical rate theory and are seen as a strong signature of tunnelling of the transferring particle, hydrogen or one of its isotopes, through the reaction energy barrier. In this paper, we build a qualitative quantum rate model with few free parameters that describes the dynamics of the transferring particle when it is exposed to energetic potentials exerted by the donor and the acceptor. The enzyme's impact on the dynamics is modelled by an additional energetic term, an oscillatory contribution known as "gating." By varying two key parameters, the gating frequency and the mean donor-acceptor separation, the model is able to reproduce well the KIE data for SLO wild-type and a variety of SLO mutants over the experimentally accessible temperature range. While SLO-specific constants have been considered here, it is possible to adapt these for other enzymes.

There are many atomistic simulation methods with very different costs, accuracies, transferabilities, and numbers of empirical parameters. I show how statistical model selection can compare these methods fairly, even when they are very different. These comparisons are also useful for developing new methods that balance cost and accuracy. As an example, I build a semiempirical model for hydrogen clusters. [ABSTRACT FROM AUTHOR]

We present a classical density functional theory (DFT) for fluid mixtures that is based on a third-order thermodynamic perturbation theory of Feynman-Hibbs-corrected Mie potentials. The DFT is developed to study the interfacial properties of hydrogen, helium, neon, deuterium, and their mixtures, i.e., fluids that are strongly influenced by quantum effects at low temperatures. White Bear fundamental measure theory is used for the hard-sphere contribution of the Helmholtz energy functional, and a weighted density approximation is used for the dispersion contribution. For mixtures, a contribution is included to account for non-additivity in the Lorentz–Berthelot combination rule. Predictions of the radial distribution function from DFT are in excellent agreement with results from molecular simulations, both for pure components and mixtures. Above the normal boiling point and 5% below the critical temperature, the DFT yields surface tensions of neon, hydrogen, and deuterium with average deviations from experiments of 7.5%, 4.4%, and 1.8%, respectively. The surface tensions of hydrogen/deuterium, para-hydrogen/helium, deuterium/helium, and hydrogen/neon mixtures are reproduced with a mean absolute error of 5.4%, 8.1%, 1.3%, and 7.5%, respectively. The surface tensions are predicted with an excellent accuracy at temperatures above 20 K. The poor accuracy below 20 K is due to the inability of Feynman–Hibbs-corrected Mie potentials to represent the real fluid behavior at these conditions, motivating the development of new intermolecular potentials. This DFT can be leveraged in the future to study confined fluids and assess the performance of porous materials for hydrogen storage and transport. [ABSTRACT FROM AUTHOR]

An alternative approach to density functional theory based on self-consistent field theory for ring polymers is applied to neutral atoms hydrogen to neon in their ground-states. The spontaneous emergence of an atomic shell structure and spherical symmetry-breaking of the total electron density are predicted by the model using the ideas of polymer excluded-volume between pairs of electrons to enforce the Pauli-exclusion principle and an exact electron self-interaction correction. The Pauli potential is approximated by neglecting inter-atomic correlations along with other types of correlations, and comparisons to Hartree–Fock theory are made, which also ignores correlations. The model shows excellent agreement with Hartree–Fock theory to within the standards of orbital-free density functional theory for the atomic binding energies and density profiles of the first six elements, providing exact matches for the elements hydrogen and helium. The predicted shell structure starts to deviate significantly past the element neon, and spherical symmetry-breaking is first predicted to occur at carbon instead of boron. The self-consistent field theory energy functional that describes the model is decomposed into thermodynamic components to trace the origin of spherical symmetry-breaking. It is found to arise from the electron density approaching closer to the nucleus in non-spherical distributions, which lowers the energy despite resulting in frustration between the quantum kinetic energy, electron–electron interaction, and the Pauli exclusion interaction. The symmetry-breaking effect is found to have a minimal impact on the binding energies, which suggests that the spherical-averaging approximation used in previous work is physically reasonable when investigating atomic systems. The pair density contour plots display behavior similar to polymer macro-phase separation, where individual electron pairs occupy single lobe structures that together form a dumbbell shape analogous to the 2p orbital shape. It is further shown that the predicted densities satisfy known constraints and produce the same total electronic density profile that is predicted by other formulations of quantum mechanics. [ABSTRACT FROM AUTHOR]

Non-Maxwellian evaporation of light atoms and molecules (particles) such as He and H2 from liquids has been observed experimentally. In this work, we use simulations to study systematically the evaporation of Lennard-Jones particles from liquid water. We find instances of sub- and super-Maxwellian evaporation, depending on the mass of the particle and the particle-water interaction strength. The observed trends are in qualitative agreement with experiment. We interpret these trends in terms of the potential of mean force and the effectiveness and frequency of collisions during the evaporation process. The angular distribution of evaporating particles is also analyzed, and it is shown that trends in the energy from velocity components tangential and normal to the liquid surface must be understood separately in order to interpret properly the angular distributions. [ABSTRACT FROM AUTHOR]

We performed fixed-node diffusion Monte Carlo (DMC) calculations to investigate structural and energetic properties of graphenylene (GPNL), a two-dimensional network of sp2-bonded carbon atoms with large near-circular pores, and its H2 separation performance for gas mixtures. We have found that the energetic stability of a GPNL monolayer is comparable to that of γ-graphyne, as evidenced by its large cohesive energy of 6.755(3) eV/atom. Diffusion barriers of several gas molecules, including hydrogen, through a GPNL membrane were determined from the analysis of their adsorption energies depending on the adsorption distance, which led to our estimation for hydrogen selectivity with respect to other target molecules. DMC hydrogen selectivity of a GPNL monolayer was found to be exceptionally high at 300 K, as high as 1010–1011 against CO and N2 gases. This, along with high hydrogen permeance due to its generic pore structure, leads us to conclude that GPNL is a promising membrane to be used as a high-performance hydrogen separator from gas mixtures. We find that when compared to our DMC results, DFT calculations tend to overestimate H2 selectivity, which is mostly due to their inaccurate description of short-range repulsive interactions. [ABSTRACT FROM AUTHOR]

We describe a simple scheme to perform phonon calculations with quantum Monte Carlo (QMC) methods and demonstrate it on metallic hydrogen. Because of the energy and length scales of metallic hydrogen and the statistical noise inherent to QMC methods, the conventional manner of calculating force constants is prohibitively expensive. We show that our alternate approach is nearly 100 times more efficient in resolving the force constants needed to calculate the phonon spectrum in the harmonic approximation. This requires only the calculation of atomic forces, as in the conventional approach, and otherwise little or no programmatic modification. [ABSTRACT FROM AUTHOR]

The activity of Pt(111) electrodes for the hydrogen evolution reaction (HER) in 0.5M H2SO4 solution is found to increase with continuous potential cycling in the HER potential region. In addition, the basic cyclic voltammograms obtained in 0.5M H2SO4 saturated with N2 after HER show several characteristic changes: the current waves for hydrogen adsorption in the region of0.2 < E < 0.35 V and for sulfate adsorption at 0.35 < E < 0.5 V decrease and the current spike at 0.44 V for the phase transition of the sulfate adlayer gradually disappears. We suggest that these changes are caused by the absorption of a small amount of hydrogen in the subsurface layer and propose a mechanism by which this enhances hydrogen evolution. [ABSTRACT FROM AUTHOR]

We introduce hybrid gausslet/Gaussian basis sets, where a standard Gaussian basis is added to a gausslet basis in order to increase accuracy near the nuclei while keeping the spacing of the grid of gausslets relatively large. The Gaussians are orthogonalized to the gausslets, which are already orthonormal, and approximations are introduced to maintain the diagonal property of the two electron part of the Hamiltonian so that it continues to scale as the second power of the number of basis functions rather than the fourth. We introduce several corrections to the Hamiltonian designed to enforce certain exact properties, such as the values of certain two-electron integrals. We also introduce a simple universal energy correction that compensates for the incompleteness of the basis stemming from the electron–electron cusps based on the measured double occupancy of each basis function. We perform a number of Hartree Fock and full configuration interaction (full-CI) test calculations on two electron systems and Hartree Fock on a ten-atom hydrogen chain to benchmark these techniques. The inclusion of the cusp correction allows us to obtain complete basis set full-CI results for the two electron cases at the level of several micro-Hartrees, and we see similar apparent accuracy for Hartree Fock on the ten-atom hydrogen chain. [ABSTRACT FROM AUTHOR]

In this work, we present an in-depth study of the role of mass transport conditions in tuning the hydrogen evolution kinetics on gold by means of rotation rate control. Interestingly, we find that the hydrogen evolution reaction (HER) activity decreases with the increasing rotation rate of the electrode. As we increase the rotation (mass transport) rate, the locally generated hydroxyl ions (2H2O +2e− → H2 + 2OH−) are transported away from the electrode surface at an accelerated rate. This results in decreasing local pH and, because of the need to satisfy local electroneutrality, decreasing near-surface cation concentration. This decrease in the near-surface cation concentration results in the suppression of HER. This is because the cations near the surface play a central role in stabilizing the transition state for the rate determining Volmer step (*H–OHδ−–cat+). Furthermore, we present a detailed analytical model that qualitatively captures the observed mass transport dependence of HER solely based on the principle of electroneutrality. Finally, we also correlate the cation identity dependence of HER on gold (Li+ < Na+ < K+) to the changes in the effective concentration of the cations in the double layer with the changes in their solvation energy. [ABSTRACT FROM AUTHOR]

The catalytic reduction in carbon dioxide is a crucial step in many chemical industrial reactions, such as methanol synthesis, the reverse water–gas shift reaction, and formic acid synthesis. Here, we investigate the role of bulk hydrogen, where hydrogen atoms are found deep inside a metal surface as opposed to subsurface ones, upon CO2 reduction over a Ni(110) surface using density functional theory and ab initio molecular dynamics simulations. While it has previously been shown that subsurface hydrogen stabilizes CO2 and can aid in overcoming reaction barriers, the role of bulk hydrogen is less studied and thus unknown with regard to CO2 reduction. We find that the presence of bulk hydrogen can significantly alter the electronic structure of the Ni(110) surface, particularly the work function and d-band center, such that CO2 adsorbs more strongly to the surface and is more easily reduced. Our results show an enhanced CO2 dissociation in the presence of bulk hydrogen, shedding light on a hitherto underappreciated mechanistic pathway for CO2 reduction on metal surfaces. [ABSTRACT FROM AUTHOR]

Phosphorus is of particular interest in astrochemistry because it is a biogenic element together with hydrogen, carbon, nitrogen, oxygen, and sulfur. However, the chemical evolution of such element in the interstellar medium (ISM) is still far from an accurate characterization, with the chemistry of P-bearing molecules being poorly understood. To provide a contribution in this direction, we have carried out an accurate investigation of the potential energy surface for the reaction between the CP radical and methanimine (CH2NH), two species already detected in the ISM. In analogy to similar systems, i.e., CH2NH + X, with X = OH, CN, and CCH, this reaction can occur—from an energetic point of view—under the harsh conditions of the ISM. Furthermore, since the major products of the aforementioned reaction, namely, E- and Z-2-phosphanylidyneethan-1-imine (HN=CHCP) and N-(phosphaneylidynemethyl)methanimine (H2C=NCP), have not been spectroscopically characterized yet, some effort has been made for filling this gap by means of accurate computational approaches. [ABSTRACT FROM AUTHOR]

An exhaustive investigation of state-to-state H+ + H2(v, j) → H+ + H2(v′, j′) transitions for rovibrational levels of molecular hydrogen below 1.3 eV from the bottom of the H2 well is carried out by means of a statistical quantum method, which assumes the complex-forming nature of the process. Integral cross sections for transitions involving states H2(v = 0, j = 0–12), H2(v = 1, j = 0–8), and H2(v = 2, j = 0–3) are obtained for collision energies within a range of Emin = 10−5 eV and Emax = 2 eV. Rate constants are then calculated between T = 5 K and 3000 K, and they are compared, when possible, with previous values reported in the literature. As a first application, the cooling rate coefficient of H2 excited by protons is determined and compared with a recent estimate. [ABSTRACT FROM AUTHOR]

The corrosion and oxidation of actinide metals, leading to the formation of metal-oxide surface layers with the catalytic evolution of hydrogen, impacts the management of nuclear materials. Here, the interaction of hydrogen with actinide dioxide (AnO2, An = U, Np, or Pu) (011) surfaces by Hubbard corrected density functional theory (PBEsol+U) has been studied, including spin–orbit interactions and non-collinear 3k anti-ferromagnetic behavior. The actinide dioxides crystalize in the fluorite-type structure, and although the (111) surface dominates the crystal morphology, the (011) surface energetics may lead to more significant interaction with hydrogen. The dissociative adsorption of hydrogen on the UO2 (0.44 eV), NpO2 (−0.47 eV), and PuO2 (−1.71 eV) (011) surfaces has been calculated. It is found that hydrogen dissociates on the PuO2 (011) surface; however, UO2 (011) and NpO2 (011) surfaces are relatively inert. Recombination of hydrogen ions is likely to occur on the UO2 (011) and NpO2 (011) surfaces, whereas hydroxide formation is shown to occur on the PuO2 (011) surface, which distorts the surface structure. [ABSTRACT FROM AUTHOR]

Nickel-based materials have been used in electrochemical water splitting and as a co-catalyst to effectively improve photocatalytic activity. However, the mechanism of the catalytic effect on hydrogen evolution of NiO and Ni(OH)2 needs further investigation because of the high risk in self-reduction with photo-induced electrons. In this work, the self-reduction of Ni2+ was confirmed during the photocatalytic H2 evolution in the CdS/Ni(OH)2 hybrid materials. Meanwhile, in situ formed metallic Ni plays an important role in the increase in the catalytic activity. Fortunately, only a few photo-induced electrons participate in the reduction of Ni2+ in CdS/Ni(OH)2 hybrid materials. The existence of metallic Ni would prevent more Ni2+ from self-reduction. The synergistic effect of Ni(OH)2 and metallic Ni contributes to the improvement of H2 evolution of CdS nanorods. [ABSTRACT FROM AUTHOR]

Multicomponent density functional theory (DFT) has many practical advantages for incorporating nuclear quantum effects into quantum chemistry calculations. Within the nuclear-electronic orbital (NEO) framework, specified nuclei, typically protons, are treated quantum mechanically on the same level as the electrons. Previously, electron-proton correlation functionals based on the local density approximation (LDA), denoted epc17 and epc18, were developed and shown to provide more accurate proton densities and energies compared to the neglect of electron-proton correlation, but a quantitatively accurate description of both densities and energies simultaneously has remained elusive. Herein, an electron-proton correlation functional that depends on the electron and proton density gradients, as well as the densities, is derived and implemented. Compared to the LDA functionals, the resulting generalized gradient approximation functional, denoted epc19, is able to simultaneously provide accurate proton densities and energies, as well as reproduce the impact of nuclear quantum effects on optimized geometries. In addition, without further parameterization, the NEO-DFT/epc19 method provides accurate densities and energies for deuterium as well as hydrogen. These results demonstrate that the form of the epc19 functional is able to capture the essential aspects of electron-proton correlation and highlights the importance of including gradient terms. This approach will enable the exploration of nuclear quantum effects and isotope effects in a wide range of systems. [ABSTRACT FROM AUTHOR]

We present a perturbation theory that combines the use of a third-order Barker–Henderson expansion of the Helmholtz energy with Mie-potentials that include first- (Mie-FH1) and second-order (Mie-FH2) Feynman–Hibbs quantum corrections. The resulting equation of state, the statistical associating fluid theory for Mie potentials of variable range corrected for quantum effects (SAFT-VRQ-Mie), is compared to molecular simulations and is seen to reproduce the thermodynamic properties of generic Mie-FH1 and Mie-FH2 fluids accurately. SAFT-VRQ Mie is exploited to obtain optimal parameters for the intermolecular potentials of neon, helium, deuterium, ortho-, para-, and normal-hydrogen for the Mie-FH1 and Mie-FH2 formulations. For helium, hydrogen, and deuterium, the use of either the first- or second-order corrections yields significantly higher accuracy in the representation of supercritical densities, heat capacities, and speed of sounds when compared to classical Mie fluids, although the Mie-FH2 is slightly more accurate than Mie-FH1 for supercritical properties. The Mie-FH1 potential is recommended for most of the fluids since it yields a more accurate representation of the pure-component phase equilibria and extrapolates better to low temperatures. Notwithstanding, for helium, where the quantum effects are largest, we find that none of the potentials give an accurate representation of the entire phase envelope, and its thermodynamic properties are represented accurately only at temperatures above 20 K. Overall, supercritical heat capacities are well represented, with some deviations from experiments seen in the liquid phase region for helium and hydrogen. [ABSTRACT FROM AUTHOR]

The Diffusion Monte Carlo (DMC) method was applied to anionic hydrogen clusters H − ( H 2 ) n (n = 1–16, 32) and their deuterated analogs using a polarizable all-atom potential energy surface (PES) developed by Calvo and Yurtsever. For the hydrogen clusters, the binding energy ΔEn appears to be a smooth function of the cluster size n, thus contradicting the previous claim that n = 12 is a "magic number" cluster. The structures of the low energy minima of the PES for these clusters belong to the icosahedral motif with the H2 molecules aligned toward the central H− ion. However, their ground state wavefunctions are highly delocalized and resemble neither the structures of the global nor local minima. Moreover, the strong nuclear quantum effects result in a nearly complete orientational disordering of the H2 molecules. For the deuterium clusters, the ground state wavefunctions are localized and the D2 molecules are aligned toward the central D− ion. However, their structures are still characterized as disordered and, as such, do not display size sensitivity. In addition, DMC simulations were performed on the mixed H − ( H 2 ) n ( D 2 ) p clusters with (n, p) = (6, 6) and (16, 16). Again, in contradiction to the previous claim, we found that the "more quantum" H2 molecules prefer to reside farther from the central H− ion than the D2 molecules. [ABSTRACT FROM AUTHOR]

New analytical gradients of the electronic energy of a confined molecular system within the extreme-pressure continuum model are presented and applied to the study of the equilibrium geometries of linear chains of hydrogen molecules nH2 under pressures. The decrease in inter- and intramolecular H–H distances with the increase in the pressure has been studied up to 80 GPa. We have also shown that the compression of the bond-lengths can be interpreted in terms of the effect of the confining potential of the electron density of the molecular systems. [ABSTRACT FROM AUTHOR]

The potential energy surface describing the interaction of the SH+ ion in its ground X3Σ− electronic state with molecular hydrogen has been computed through restricted coupled cluster calculations including single, double, and (perturbative) triple excitations [RCCSD(T)], with the assumption of fixed bond lengths. The computed points were fit to an analytical form suitable for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate coefficients. Representative sets of energy dependent state-to-state rotationally inelastic cross sections and rate coefficients are presented and discussed. A propensity to conserve the fine-structure label is observed, as expected from a formal analysis of the scattering equations. [ABSTRACT FROM AUTHOR]

We report state-to-state quasi-classical trajectory rate constants for the D + H2 reactive collision, using the accurate H3 global potential energy surface of Mielke et al. [J. Chem. Phys. 116, 4142 (2002)]. High relative collision energies (up to ≈56 000 K) and high rovibrational levels of H2 (up to ≈50 000 K), relevant to various astrophysical media, are considered. The HD product molecule is formed in highly excited rovibrational states, over a wide collision energy range. The collision-induced dissociation channel (often overlooked in fully quantum reaction dynamics calculations) is found to be significantly populated, even at collision energies as low as 1500 K. [ABSTRACT FROM AUTHOR]

The free-complement (FC) theory for solving the Schrödinger equation (SE) was applied to calculate the potential energy curves of the ground and excited states of the hydrogen molecule (H2) with the 1Σg+, 1Σu+, 3Σg+, 3Σu+, 1Πg, 1Πu, 3Πg, 3Πu, 1Δg, 1Δu, 3Δg, 3Δu, 1Φg, 1Φu, 3Φg, and 3Φu symmetries (in total, 54 states). The initial functions of the FC theory were formulated based on the atomic states of the hydrogen atom and its positive and negative ions at the dissociation limits. The local Schrödinger equation (LSE) method, which is a simple sampling-type integral-free methodology, was employed instead of the ordinary variational method and highly accurate results were obtained stably and smoothly along the potential energy curves. Thus, with the FC-LSE method, we succeeded to perform the comprehensive studies of the H2 molecule from the ground to excited states belonging up to higher angular momentum symmetries and from equilibriums to dissociation limits with almost satisfying spectroscopic accuracy, i.e., 10−6 hartree order around 1 cm−1, as absolute solutions of the SE by moderately small calculations. [ABSTRACT FROM AUTHOR]

The PM6 implementation in the GAMESS program is extended to elements requiring d-integrals and interfaced with the conducter-like polarized continuum model of solvation, including gradients. The accuracy of aqueous solvation energies computed using AM1, PM3, PM6, and DFT tight binding (DFTB) and the Solvation Model Density (SMD) continuum solvation model is tested using the Minnesota Solvation Database data set. The errors in SMD solvation energies predicted using Neglect of Diatomic Differential Overlap (NDDO)-based methods are considerably larger than when using density functional theory (DFT) and HF, with root mean square error (RMSE) values of 3.4-5.9 (neutrals) and 6-15 kcal/mol (ions) compared to 2.4 and ~5 kcal/mol for HF/6-31G(d). For the NDDO-based methods, the errors are especially large for cations and considerably higher than the corresponding conductor-like screening model results, which suggests that the NDDO/SMD results can be improved by re-parameterizing the SMD parameters focusing on ions. We found that the best results are obtained by changing only the radii for hydrogen, carbon, oxygen, nitrogen, and sulfur, and this leads to RMSE values for PM3 (neutrals: 2.8/ions: ~5 kcal/mol), PM6 (4.7/~5 kcal/mol), and DFTB (3.9/~5 kcal/mol) that are more comparable to HF/6-31G(d) (2.4/~5 kcal/mol). Although the radii are optimized to reproduce aqueous solvation energies, they also lead more accurate predictions for other polar solvents such as dimethyl sulfoxide, acetonitrile, and methanol, while the improvements for non-polar solvents are negligible. [ABSTRACT FROM AUTHOR]

The Electric Stark effect of a Kramers-Henneberger (KH) state of hydrogen atoms in both linearly and circularly polarized laser fields is studied. For the ground KH state of H atoms with a small quiver amplitude, the quadratic Stark effect is observed. For a large quiver amplitude, the Stark effect is quadratic only in a weak electric field and quickly changes to linear as the electric field increases. The atomic structure of the KH state is very sensitive to the electric field and can be easily polarized. The huge polarizability and induced dipole moment are comparable to those of Rydberg atoms. [ABSTRACT FROM AUTHOR]

Progress in the hydrogen fuel field requires a clear understanding and characterization of how materials of interest interact with hydrogen. Due to the inherently quantum mechanical nature of hydrogen nuclei, any theoretical studies of these systems must be treated quantum dynamically. One class of material that has been examined in this context are dihydrogen complexes. Since their discovery by Kubas in 1984, many such complexes have been studied both experimentally and theoretically. This particular study examines the rotational dynamics of the dihydrogen ligand in the Fe(H)2(H2)(PEtPh2)3 complex, allowing for full motion in both the rotational degrees of freedom and treating the quantum dynamics (QD) explicitly. A "gas-phase" global potential energy surface is first constructed using density functional theory with the Becke, 3-parameter, Lee-Yang-Parr functional; this is followed by an exact QD calculation of the corresponding rotation/libration states. The results provide insight into the dynamical correlation of the two rotation angles as well as a comprehensive analysis of both ground- and excited-state librational tunneling splittings. The latter was computed to be 6.914 cm-1--in excellent agreement with the experimental value of 6.4 cm-1. This work represents the first full-dimensional ab initio exact QD calculation ever performed for dihydrogen ligand rotation in a coordination complex. [ABSTRACT FROM AUTHOR]

We employ the on-the-fly surface hopping classical trajectory algorithm to study the electron detachment process in low-energy H- + H2 collisions. The ground-state and the first-excited-state Local Complex Potentials (LCPs) calculated by the generalized diatomics-in-molecule method are used for the full three-dimensional nonadiabatic nuclear dynamics. Two kinds of nonadiabatic effects are taken into account: discrete-discrete transitions and discrete-continuous transitions. Discrete-discrete nonadiabatic transition probabilities are calculated by means of the adiabatic-potential-based formula within the Landau-Zener model for each individual trajectory computed along real parts of the LCPs. Discrete-continuous (electron detachment) nonadiabatic transition probabilities are calculated via quasi-stationary widths which are related to the imaginary parts of the LCPs of both the electronic states of the H- 3 anion. Two mechanisms of the electron detachment process are treated and discussed: the direct mechanism based on quasi-stationarity of the ground state and the indirect mechanism based on both nonadiabatic transitions from the ground state to the first excited state and quasi-stationarity of the excited state. It is shown that the direct mechanism prevails at low collision energies, while the indirect mechanism makes a substantial contribution at relatively high collision energies, roughly higher than 5 eV. At collision energies higher than 2 eV, the electron detachment probability has rather high values and this affects noticeably other inelastic processes in these collisions. The electron detachment cross section in H- + H2 collisions is calculated for the collision energy range from 1 to 100 eV and a reasonable agreement with available experimental data is obtained. [ABSTRACT FROM AUTHOR]

We solve the time-dependent Schrödinger equation describing a water molecule driven by a superposition of the extreme ultraviolet and IR pulses typical for a reconstruction of attosecond beating by interference of two-photon transitions experiment. This solution is obtained by a combination of the time-dependent coordinate scaling and the density functional theory with self-interaction correction. Results of this solution are used to determine the time delay in photoionization of the water and hydrogen molecules. [ABSTRACT FROM AUTHOR]

The investigation of partially deuterated methoxy radicals is important because the symmetry lowering from C3v to Cs provides new insights into the couplings between rovibronic states via Jahn-Teller and spin-orbit interactions. The vibrational spectrum of the partially deuterated methoxy radical CH2DO in a matrix of p-H2 has been recorded. This species was prepared by irradiating a p-H2 matrix containing deuterated d1-nitritomethane (CH2DONO) at 3.3 K with laser light at 355 nm. The identification of the radical is based on the photochemical behavior of the precursor and comparison of observed vibrational wavenumbers and infrared (IR) intensities with those predicted from a refined quartic, curvilinear, internal coordinate force field calculated with the coupled-cluster singles and doubles with perturbative triples/cc-pVTZ method. CH2DO reacts with H2 with a rate coefficient (3.5 - 1.0) x 10-3 s-1. Predominantly c-CHDOH and a negligibly small amount of t- CHDOH were produced. This stereoselectivity results from the reaction H + Cs-CH2DOH, which was demonstrated by an additional experiment on irradiation of a CH2DOH/Cl2/p-H2 matrix with ultraviolet and IR light to induce the H + CH2DOH reaction; only c-CHDOH was observed from this experiment. Even though the energies of transition states and products for the formation of c-CHDOH and t-CHDOH differ by only ~10 cm-1, the selective formation of c-CHDOH can be explained by tunneling of the hydrogen atom via an optimal tunneling path. Similarly, the vibronic spectrum for the partially deuterated specie d2-methoxy radical (CHD2O) was obtained upon irradiation of d2-nitritomethane (CHD2ONO) at 355 nm. Lines associated with the fundamental vibrational modes were observed and assigned; line positions agree with theoretically predicted vibrational wavenumbers. CHD2O reacts with H2 with a rate coefficient (6.0 ± 1.4) x 10-3 s-1; CD2OH was produced as a major product because the barrier for the formation of CHDOH from H + CHD2OH is greater by ***400 cm-1. Rate coefficients of the decays of CH3O, CH2DO, CHD2O, and CD3O and their corresponding potential energy surfaces are compared. [ABSTRACT FROM AUTHOR]

Hydrogen abstraction from NH3 by OH to produce H2O and NH2--an important reaction in combustion of NH3 fuel--was studied with a theoretical approach that combines high level quantum chemistry and advanced chemical kinetics methods. Thermal rate constants calculated from first principles agree well (within 5%-20%) with available experimental data over a temperature range that extends from 200 to 2500 K. Quantum mechanical tunneling effects were found to be important; they lead to a decided curvature and non-Arrhenius behavior for the rate constant. [ABSTRACT FROM AUTHOR]

Graphene-based derivatives with covalent functionalization and well-defined stoichiometry are highly desirable in viewof their application as functional surfaces. Here, we have evaluated by ab initio calculations the energy of formation and the phase diagram of hydroxylated graphane structures, i.e., fully functionalized graphene derivatives coordinated with -H and -OH groups. We compared these structures to different hydrogenated and non-hydrogenated graphene oxide derivatives, with high level of epoxide and hydroxyl groups functionalization. Based on our calculations, stable phases of hydroxylated graphane with lowand high contents of hydrogen are demonstrated for high oxygen and hydrogen partial pressure, respectively. Stable phases of graphene oxide with a mixed carbon hybridization are also found. Notably, the synthesis of hydroxylated graphane has been recently reported in the literature. [ABSTRACT FROM AUTHOR]

Quantum reactive scattering calculations are reported for the ultracold hydrogen-exchange reaction and its non-reactive atom-exchange isotopic counterparts, proceeding from excited rotational states. It is shown that while the geometric phase (GP) does not necessarily control the reaction to all final states, one can always find final states where it does. For the isotopic counterpart reactions, these states can be used to make a measurement of the GP effect by separately measuring the even and odd symmetry contributions, which experimentally requires nuclear-spin final-state resolution. This follows from symmetry considerations that make the even and odd identical-particle exchange symmetrywavefunctions which include the GP locally equivalent to the opposite symmetrywavefunctions which do not. It is shown how this equivalence can be used to define a constant which quantifies the GP effect and can be obtained solely from experimentally observable rates. This equivalence reflects the important role that discrete symmetries play in ultracold chemistry and highlights the key role that ultracold reactions can play in understanding fundamental aspects of chemical reactivity more generally. [ABSTRACT FROM AUTHOR]

An effective six-dimensional ab initio potential energy surface (PES) for H2-OCS which explicitly includes the intramolecular stretch normal modes of carbonyl sulfide (OCS) is presented. The electronic structure computations are carried out using the explicitly correlated coupled cluster [CCSD(T)-F12] method with the augmented correlation-consistent aug-cc-pVTZ basis set, and the accuracy is critically tested by performing a series of benchmark calculations. Analytic four-dimensional PESs are obtained by least-squares fitting vibrationally averaged interaction energies to the Morse/long-range potential model. These fits to 13 485 points have a root-mean-square deviation (RMSD) of 0.16 cm-1. The combined radial discrete variable representation/angular finite basis representation method and the Lanczos algorithm were employed to evaluate the rovibrational energy levels for five isotopic species of the OCS-hydrogen complexes. The predicted transition frequencies and intensities based on the resulting vibrationally averaged PESs are in good agreement with the available experimental values, whose RMSDs are smaller than 0.004 cm-1 for five different species of OCS-hydrogen complexes. The calculated infrared band origin shifts for all five species of OCS-hydrogen complexes are only 0.03 cm-1 smaller than the corresponding experimental values. These validate the high quality of our PESs which can be used for modeling OCS doped in hydrogen clusters to further study quantum solution and microscopic superfluidity. In addition, the analytic coordinate transformation functions between isotopologues are also derived due to the center of mass shifting of different isotope substitutes. [ABSTRACT FROM AUTHOR]

The bleaching process in the C-F stretching mode (v3 band) of CH3F-(ortho-H2)n [n = 0 and 1] clusters in solid para-H2 was monitored using pump and probe laser spectroscopy on the C-H stretching mode (v1 and 2v5 bands). From an analysis of the depleted spectral profiles, the transition frequency and linewidth of each cluster were directly determined. The results agree with the values previously derived from a deconvolution analysis of the broadened v1/2v5 spectrum observed by FTIR spectroscopy. The complementary increase and decrease between the n = 0 and 1 components were also verified through monitoring the v1 and 2v5 bands, which suggests a closed system among the CH3F-(ortho-H2)n clusters. These observations provide experimental verification of the CH3F-(ortho-H2)n cluster model. On the other hand, a trial to observe the bleaching process by pumping the C-H stretching mode was not successful. This result may be important for understanding the dynamics of vibrational relaxation processes in CH3F-(ortho-H2)n in solid para-H2. [ABSTRACT FROM AUTHOR]

Our first-principles calculations reveal that there exist contrasting modes of initial nucleations of graphene on transition metal surfaces, in which hydrogen plays the role. On Cu(100) and Cu(111) surfaces, an sp2-type network of carbons can be automatically formed with the help of hydrogen under very low carbon coverages. Thus, by tuning the chemical potential of hydrogen, both of the nucleation process and the following growth can be finely controlled. In contrast, on the Ni(111) surface, instead of hydrogen, the carbon coverage is the critical factor for the nucleation and growth. These findings serve as new insights for further improving the poor quality of the grown graphene on transition metal substrates. [ABSTRACT FROM AUTHOR]

Silicon 1s Near Edge X-ray Absorption Fine Structure (NEXAFS) spectra of silicon nanocrystals have been examined as a function of nanocrystal size (3-100 nm), varying surface functionalization (hydrogen or 1-pentyl termination), or embedded in oxide. The NEXAFS spectra are characterized as a function of nanocrystal size and surface functionalization. Clear spectroscopic evidence for long range order is observed silicon nanocrystals that are 5-8 nm in diameter or larger. Energy shifts in the silicon 1s NEXAFS spectra of covalently functionalized silicon nanocrystals with changing size are attributed to surface chemical shifts and not to quantum confinement effects. [ABSTRACT FROM AUTHOR]

We have performed calculations on the dissociative chemisorption of H2 on un-reconstructed and reconstructed Au(111) with density functional theory, and dynamics calculations on this process on un-reconstructed Au(111). Due to a very late barrier for dissociation, H2 + Au(111) is a candidate H2-metal system for which the dissociative chemisorption could be considerably affected by the energy transfer to electron-hole pairs. Minimum barrier geometries and potential energy surfaces were computed for six density functionals. The functionals tested yield minimum barrier heights in the range of 1.15-1.6 eV, and barriers that are even later than found for the similar H2 + Cu(111) system. The potential energy surfaces have been used in quasi-classical trajectory calculations of the initial (v,J) state resolved reaction probability for several vibrational states v and rotational states J of H2 and D2. Our calculations may serve as predictions for state-resolved associative desorption experiments, from which initial state-resolved dissociative chemisorption probabilities can be extracted by invoking detailed balance. The vibrational efficacy ηv=0→1 reported for D2 dissociating on un-reconstructed Au(111) (about 0.9) is similar to that found in earlier quantum dynamics calculations on H2 + Ag(111), but larger than found for D2 + Cu(111). With the two functionals tested most extensively, the reactivity of H2 and D2 exhibits an almost monotonic increase with increasing rotational quantum number J. Test calculations suggest that, for chemical accuracy (1 kcal/mol), the herringbone reconstruction of Au(111) should be modeled. [ABSTRACT FROM AUTHOR]

The addition reactions of chlorine atom with isobutene (i-C4H8) in solid para-hydrogen (p-H2) were investigated with infrared (IR) absorption spectra. When a p-H2 matrix containing Cl2 and isobutene was irradiated with ultraviolet light at 365 nm, intense lines in a set at 534.5, 1001.0, 1212.9, 1366.0, 2961.6, and 2934.7 cm-1, and several weaker others due to the 1-chloro-2-methyl-2-propyl radical, ·C(CH3)2CH2Cl, and those in a second set including intense ones at 642.7, 799.2, 1098.2, 1371.8, and 3027.3 cm-1 due to the 2-chloro-2-methylpropyl radical, ·CH2C(CH3)2Cl, appeared; the ratio of ·C(CH3)2CH2Cl to ·CH2C(CH3)2Cl was approximately (3 ± 1):1. The observed wavenumbers and relative intensities agree with the vibrational wavenumbers and IR intensities predicted with the B3PW91/aug-cc-pVTZ method. That the Cl atom adds to both carbons of the C=C bond of isobutene with the terminal site slightly favored is consistent with the energies of products predicted theoretically, but is in contrast to the reaction of Cl + propene in solid p-H2 in which the addition of Cl to mainly the central C atom was previously reported. The role of the p-H2 matrix in affecting the reaction paths is discussed. Absorption lines of the complex i-C4H8·Cl2 and the dichloro-product anti-1,2-dichloro-2-methylpropane, a-CH2ClCCl(CH3)2, are also characterized. [ABSTRACT FROM AUTHOR]

The time-dependent wave packet method has been employed to calculate state-to-state differential cross sections for the title reaction in full dimensions. It is found that the majority of H2O is produced in the first stretching excited states, with a large fraction of available energy for the reaction ending up as product internal motion. The differential cross sections for collision energy up to 0.4 eV are all peaked in the backward direction, but the width of the angular distribution increases considerably as the increase of collision energy. The isotope effect was also examined by comparing the scattering angular distribution for the title reaction with those for the HD + OH and D2 + OH reactions obtained in our previous work. [ABSTRACT FROM AUTHOR]

We have experimentally performed the coherent control of delocalized ro-vibrational wave packets (RVWs) of solid para-hydrogen (p-H2) by the wave packet interferometry (WPI) combined with coherent anti-Stokes Raman scattering (CARS). RVWs of solid p-H2 are delocalized in the crystal, and the wave function with wave vector k ~ 0 is selectively excited via the stimulated Raman process. We have excited the RVW twice by a pair of femtosecond laser pulses with delay controlled by a stabilized Michelson interferometer. Using a broad-band laser pulse, multiple ro-vibrational states can be excited simultaneously. We have observed the time-dependent Ramsey fringe spectra as a function of the inter-pulse delay by a spectrally resolved CARS technique using a narrowband probe pulse, resolving the different intermediate states. Due to the different fringe oscillation periods among those intermediate states, we can manipulate their amplitude ratio by tuning the inter-pulse delay on the sub-femtosecond time scale. The state-selective manipulation and detection of the CARS signal combined with the WPI is a general and efficient protocol for the control of the interference of multiple quantum states in various quantum systems. [ABSTRACT FROM AUTHOR]

We have determined four-dimensional ab initio quasi-diabatic potential energy surfaces describing the interaction of CH(X2Π) with H2, under the assumption of fixed CH and H2 internuclear separations. These calculations employed the multi-reference configuration interaction method [MRCISD+Q(Davidson)]. The computed points were fit to an analytical form suitable for time-independent quantum scattering calculations of rotationally inelastic cross sections and rate constants. [ABSTRACT FROM AUTHOR]

Fundamental details concerning the interaction between H2 and CH3-Si(111) have been elucidated by the combination of diffractive scattering experiments and electronic structure and scattering calculations. Rotationally inelastic diffraction (RID) of H2 and D2 from this model hydrocarbon-decorated semiconductor interface has been confirmed for the first time via both time-of-flight and diffraction measurements, with modest j = 0 → 2 RID intensities for H2 compared to the strong RID features observed for D2 over a large range of kinematic scattering conditions along two high-symmetry azimuthal directions. The Debye-Waller model was applied to the thermal attenuation of diffraction peaks, allowing for precise determination of the RID probabilities by accounting for incoherent motion of the CH3-Si(111) surface atoms. The probabilities of rotationally inelastic diffraction of H2 and D2 have been quantitatively evaluated as a function of beam energy and scattering angle, and have been compared with complementary electronic structure and scattering calculations to provide insight into the interaction potential between H2 (D2) and hence the surface charge density distribution. Specifically, a six-dimensional potential energy surface (PES), describing the electronic structure of the H2(D2)/CH3-Si(111) system, has been computed based on interpolation of density functional theory energies. Quantum and classical dynamics simulations have allowed for an assessment of the accuracy of the PES, and subsequently for identification of the features of the PES that serve as classical turning points. A close scrutiny of the PES reveals the highly anisotropic character of the interaction potential at these turning points. This combination of experiment and theory provides new and important details about the interaction of H2 with a hybrid organic-semiconductor interface, which can be used to further investigate energy flow in technologically relevant systems. [ABSTRACT FROM AUTHOR]

We report an investigation of the translation-rotation (TR) level structure of H2 entrapped in C60, in the rigid-monomer approximation, by means of a low-order perturbation theory (PT). We focus in particular on the degree to which PT can accurately account for that level structure, by comparison with the variational quantum five-dimensional calculations. To apply PT to the system, the interaction potential of H2@C60 is decomposed into a sum over bipolar spherical tensors. A zeroth-order Hamiltonian, Ħ 0, is then constructed as the sum of the TR kinetic-energy operator and the one term in the tensor decomposition of the potential that depends solely on the radial displacement of the H2 center of mass (c.m.) from the cage center. The remaining terms in the potential are treated as perturbations. The eigenstates of Ħ0, constructed to also account for the coupling of the angular momentum of the H2 c.m. about the cage center with the rotational angular momentum of the H2 about the c.m., are taken as the PT zeroth-order states. This zeroth-order level structure is shown to be an excellent approximation to the true one except for two types of TR-level splittings present in the latter. We then show that first-order PT accounts very well for these splittings, with respect to both their patterns and magnitudes. This allows one to connect specific features of the level structure with specific features of the potential-energy surface, and provides important new physical insight into the characteristics of the TR level structure. [ABSTRACT FROM AUTHOR]

The hydrogens in protonated acetylene are very mobile and can easily migrate around the C2 core by moving between classical and non-classical structures of the cation. The lowest energy structure is the T-shaped, non-classical cation with a hydrogen bridging the two carbons. Conversion to the classical H2CCH+ ion requires only 4 kcal/mol. The effect of circularly polarized light on the migration of hydrogens in oriented C2H3+ has been simulated by Born-Oppenheimer molecular dynamics. Classical trajectory calculations were carried out with the M062X/6-311+G(3df,2pd) level of theory using linearly and circularly polarized 32 cycle 7 µm cosine squared pulses with peak intensity of 5.6 × 1013 W/cm² and 3.15 × 1013 W/cm², respectively. These linearly and circularly polarized pulses transfer similar amounts of energy and total angular momentum to C2H3+. The average angular momentum vectors of the three hydrogens show opposite directions of rotation for right and left circularly polarized light, but no directional preference for linearly polarized light. This difference results in an appreciable amount of angular displacement of the three hydrogens relative to the C2 core for circularly polarized light, but only an insignificant amount for linearly polarized light. Over the course of the simulation with circularly polarized light, this corresponds to a propeller-like motion of the three hydrogens around the C2 core of protonated acetylene. [ABSTRACT FROM AUTHOR]