Cluster isomerization induced by electron attachment

The mechanism of cluster isomerization induced by electron attachment to Na4Cl4 clusters was explored by the constant temperature quantum path integral molecular dynamics method over the temperature range 50°K–1200°K. The configurational modifications can be related to the role of an excess electron as a pseudonegative ion and to the effect of partial neutralization of the alkali ion(s) by the excess electron which induce isomerization at moderately low temperatures in comparison with the parent neutral Na4Cl4 cluster.KeywordsTemperature DomainExcess ElectronReorganization EnergyElectron AttachmentConfigurational ChangeThese keywords were added by machine and not by the authors. This process is experimental and the keywords may be updated as the learning algorithm improves.

We present the aspects of the path integral molecular dynamics (PIMD) method relative to (a) its theoretical fundamental principles, (b) its applicability to model quantum systems and (c) its implementation as a simulation tool. The PIMD method is based on the discretized path integral representation of quantum mechanics. In this representation, a quantum particle is isomorphic to a closed polymer chain. The problem of the indistinguishability of quantum particles is tackled with a non-local exchange potential. When the exact density matrix of the quantum particles is used, the exchange potential is exact. We use a high temperature approximation to the density matrix leading to an approximate form of the exchange potential. This quantum molecular dynamics method allows the simulation of collections of quantum particles at finite temperature. Our algorithm can be made to scale linearly with the number of quantum states on which the density matrix is projected. Therefore, it can be optimized to run efficiently on parallel computers. We apply the PIMD method to the electron plasma in 3-dimension. The kinetic and potential energies are calculated and compared with results for similar systems simulated with a variational Monte Carlo method. Both results show good agreements with each other at all the densities studied. The method is then use to model the thermodynamic behavior of a simple alkali metal. In these simulations, ions and valence electrons are treated as classical and quantum particles, respectively. The simple metal undergoes a phase transformation upon heating. Furthermore, to demonstrate the richness of behaviors that can be studied with the PIMD method, we also report on the metal to insulator transition in a hydrogenoid lattice. Finally, in previous studies of alkali metals, electrons interacted with ions via local pseudo-potentials, an extension of the method to modeling electrons in non-local pseudo-potentials is also presented with applications.

A variational path integral molecular dynamics study of a solid helium-4

A path integral centroid molecular dynamics method for Bose and Fermi statistics

Recent neutron Compton scattering experiments have detected the proton momentum distribution in water. The theoretical calculation of this property can be carried out via "open" path integral expressions. In this work, present an extension of the staging path integral molecular dynamics method, which is then employed to calculate the proton momentum distributions of water in the solid, liquid, and supercritical phases. We utilize a flexible, single point charge empirical force field to model the system's interactions. The calculated momentum distributions depict both agreement and discrepancies with experiment. The differences may be explained by the deviation of the force field from the true interactions. These distributions provide an abundance of information about the environment and interactions surrounding the proton.

Theoretical analysis of Muoniated compounds using path integral molecular dynamics method

D. Scharf, J. Jortner and U. Landman, atomic and molecular quantum mechanics by the path integral molecular dynamics method: Chem. Phys. Letters 130 (1986) 504

Efficient numerical algorithms are developed for use with two finite temperature semiclassical approximations to quantum dynamics both of which require trajectories generated on potentials of mean force derived from the path integral expression for the density matrix. The numerical algorithms are formed from the combination of a classical adiabatic relation similar to that used in the Car–Parrinello method and an efficient path integral molecular dynamics scheme. Results on model, an anharmonic oscillator and a realistic, fluid para-hydrogen, problem indicate that semiclassical dynamics can be obtained for virtually the same computational cost as structure and thermodynamics.

A local pseudopotential for the interaction of an electron with a water molecule in the electronic ground state is developed. The potential contains Coulomb, adiabatic polarization, exclusion, and exchange contributions. The potential is suitable for a description of excess electron states in water clusers, and for studies of electron solvation in water. Quantum path integral molecular dynamics simulations of electron localization in water clusters employing this potential yield results in agreement with available experimental data and all-electron quantum chemical calculations.

The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics simulations. The present results suggest that the inclusion of nuclear quantum effects into an empirical model for water enhances the ability of such model to faithfully represent experimental data, presumably through an increased ability of the model itself to capture realistic physical effects.

In view of the important implications of excess electrons (EEs) interacting with CO2-H2O clusters in many fields, using ab initio molecular dynamics simulation technique, we reveal the structures and dynamics of an EE associated with its localization and subsequent time evolution in heterogeneous CO2-H2O mixed media. Our results indicate that although hydration can increase the electron-binding ability of a CO2 molecule, it only plays an assisting role. Instead, it is the bending vibrations that play the major role in localizing the EE. Due to enhanced attraction of CO2, an EE can stably reside in the empty, low-lying π(*) orbital of a CO2 molecule via a localization process arising from its initial binding state. The localization is completed within a few tens of femtoseconds. After EE trapping, the ∠OCO angle of the core CO2 (-) oscillates in the range of 127°∼142°, with an oscillation period of about 48 fs. The corresponding vertical detachment energy of the EE is about 4.0 eV, which indicates extreme stability of such a CO2-bound solvated EE in [CO2(H2O)n](-) systems. Interestingly, hydration occurs not only on the O atoms of the core CO2 (-) through formation of O⋯H-O H-bond(s), but also on the C atom, through formation of a C⋯H-O H-bond. In the latter binding mode, the EE cloud exhibits considerable penetration to the solvent water molecules, and its IR characteristic peak is relatively red-shifted compared with the former. Hydration on the C site can increase the EE distribution at the C atom and thus reduce the C⋯H distance in the C⋯H-O H-bonds, and vice versa. The number of water molecules associated with the CO2 (-) anion in the first hydration shell is about 4∼7. No dimer-core (C2O4 (-)) and core-switching were observed in the double CO2 aqueous media. This work provides molecular dynamics insights into the localization and time evolution dynamics of an EE in heterogeneous CO2-H2O media.

Nonphotochemical quenching (NPQ) is the fundamental process by which plants exposed to high light intensities dissipate the potentially harmful excess energy as heat. Recently, it has been shown that efficient energy dissipation can be induced in the major light-harvesting complexes of photosystem II (LHCII) in the absence of protein-protein interactions. Spectroscopic measurements on these samples (LHCII gels) in the quenched state revealed specific alterations in the absorption and circular dichroism bands assigned to neoxanthin and lutein 1 molecules. In this work, we investigate the changes in conformation of the pigments involved in NPQ using resonance Raman spectroscopy. By selective excitation we show that, as well as the twisting of neoxanthin that has been reported previously, the lutein 1 pigment also undergoes a significant change in conformation when LHCII switches to the energy dissipative state. Selective two-photon excitation of carotenoid (Car) dark states (Car S(1)) performed on LHCII gels shows that the extent of electronic interactions between Car S(1) and chlorophyll states correlates linearly with chlorophyll fluorescence quenching, as observed previously for isolated LHCII (aggregated versus trimeric) and whole plants (with versus without NPQ).

Localization of excess electrons in sodium chloride clusters is studied using first-principles calculations based on the local-density approximation. In clusters with a rocksaltlike structure containing an anion vacancy, an excess electron was found to localize on an anion vacancy site, in analogy to the bulk F center, stabilizing the structure. Approximately the same degree of localization was found for surface, corner, and interior anion vacancies. For clusters without an anion vacancy, such as ${\mathrm{Na}}_{14}$${\mathrm{Cl}}_{13}$, the excess electron is found to be delocalized over the cluster surface. In contrast to previous theoretical predictions, no tendency is found for the excess electron to localize at a single sodium site, and localization-induced dissociation of a neutral alkali-metal atom is not expected to occur for an isolated cluster.

The solvation dynamics and reactivity of localized excess electrons in aqueous environments have attracted great attention in many areas of physics, chemistry, and biology. This manifold attraction results from the importance of water as a solvent in nature as well as from the key role of low-energy electrons in many chemical reactions. One prominent example is the electron-induced dissociation of chlorofluorocarbons (CFCs). Low-energy electrons are also critical in the radiation chemistry that occurs in nuclear reactors. Excess electrons in an aqueous environment are localized and stabilized by the local rearrangement of the surrounding water dipoles. Such solvated or hydrated electrons are known to play an important role in systems such as biochemical reactions and atmospheric chemistry. Despite numerous studies over many years, little is known about the microscopic details of these electron-induced chemical processes, and interest in the fundamental processes involved in the reactivity of trapped electrons continues. In this Account, we present a surface science study of the dynamics and reactivity of such localized low-energy electrons at D(2)O crystallites that are supported by a Ru(001) single crystal metal surface. This approach enables us to investigate the generation and relaxation dynamics as well as dissociative electron attachment (DEA) reaction of excess electrons under well-defined conditions. They are generated by photoexcitation in the metal template and transferred to trapping sites at the vacuum interface of crystalline D(2)O islands. In these traps, the electrons are effectively decoupled from the electronic states of the metal template, leading to extraordinarily long excited state lifetimes on the order of minutes. Using these long-lived, low-energy electrons, we study the DEA to CFCl(3) that is coadsorbed at very low concentrations (∼10(12) cm(-2)). Using rate equations and direct measurement of the change of surface dipole moment, we estimated the electron surface density for DEA, yielding cross sections that are orders of magnitude higher than the electron density measured in the gas phase.

Excess electron capture by hydrated histidine side-chain group