Metal-water interface formation: Thermodynamics from ab initio molecular dynamics simulations.

Metal-water interfaces are central to many electrochemical, (electro)catalytic, and materials science processes and systems. However, our current understanding of their thermodynamic properties is limited by the scarcity of accurate experimental and computational data and procedures. In this work, thermodynamic quantities for metal-water interface formation are computed for a range of FCC(111) surfaces (Pd, Pt, Au, Ag, Rh, and PdAu) through extensive density functional theory based molecular dynamics and the two-phase entropy model. We show that metal-water interface formation is thermodynamically favorable and that most metal surfaces studied in this work are completely wettable, i.e., have contact angles of zero. Interfacial water has higher entropy than bulk water due to the increased population of low-frequency translational modes. The entropic contributions also correlate with the orientational water density, and the highest solvation entropies are observed for interfaces with a moderately ordered first water layer; the entropic contributions account for up to ∼25% of the formation free energy. Water adsorption energy correlates with the water orientation and structure and is found to be a good descriptor of the internal energy part of the interface formation free energy, but it alone cannot satisfactorily explain the interfacial thermodynamics; the interface formation is driven by the competition between energetic and entropic contributions. The obtained results and insight can be used to develop, parameterize, and benchmark theoretical and computational methods for studying metal-water interfaces. Overall, our study yields benchmark-quality data and fundamental insight into the thermodynamic forces driving metal-water interface formation.

Addressing Dynamics at Catalytic Heterogeneous Interfaces with DFT-MD: Anomalous Temperature Distributions from Commonly Used Thermostats.

Density functional theory-based molecular dynamics (DFT-MD) has been widely used for studying the chemistry of heterogeneous interfacial systems under operational conditions. We report frequently overlooked errors in thermostated or constant-temperature DFT-MD simulations applied to study (electro)catalytic chemistry. Our results demonstrate that commonly used thermostats such as Nosé-Hoover, Berendsen, and simple velocity-rescaling methods fail to provide a reliable temperature description for systems considered. Instead, nonconstant temperatures and large temperature gradients within the different parts of the system are observed. The errors are not a "feature" of any particular code but are present in several ab initio molecular dynamics implementations. This uneven temperature distribution, due to inadequate thermostatting, is well-known in the classical MD community, where it is ascribed to the failure in kinetic energy equipartition among different degrees of freedom in heterogeneous systems (Harvey et al. J. Comput. Chem. 1998, 726-740) and termed the flying ice cube effect. We provide tantamount evidence that interfacial systems are susceptible to substantial flying ice cube effects and demonstrate that the traditional Nosé-Hoover and Berendsen thermostats should be applied with care when simulating, for example, catalytic properties or structures of solvated interfaces and supported clusters. We conclude that the flying ice cube effect in these systems can be conveniently avoided using Langevin dynamics.

Influence of a Cu-zirconia interface structure on CO2 adsorption and activation.

CO2 adsorption and activation on a catalyst are key elementary steps for CO2 conversion to various valuable products. In the present computational study, we screened different Cu-ZrO2 interface structures and analyzed the influence of the interface structure on CO2 binding strength using density functional theory calculations. Our results demonstrate that a Cu nanorod favors one position on both tetragonal and monoclinic ZrO2 surfaces, where the bottom Cu atoms are placed close to the lattice oxygens. In agreement with previous calculations, we find that CO2 prefers a bent bidentate configuration at the Cu-ZrO2 interface and the molecule is clearly activated being negatively charged. Straining of the Cu nanorod influences CO2 adsorption energy but does not change the preferred nanorod position on zirconia. Altogether, our results highlight that CO2 adsorption and activation depend sensitively on the chemical composition and atomic structure of the interface used in the calculations. This structure sensitivity may potentially impact further catalytic steps and the overall computed reactivity profile.

A fugacity model assessment of ibuprofen, diclofenac, carbamazepine, and their transformation product concentrations in an aquatic environment.

An updated version of FATEMOD, a multimedia fugacity model for environmental fate of organic chemicals, was set up to assess environmental behaviour of three pharmaceuticals in northern Lake Päijänne, Finland. Concentrations of ibuprofen, diclofenac, and carbamazepine were estimated at various depths at two sites: near a wastewater treatment plant and 3.5 km downstream the plant. When compared with environmental sampling data from corresponding depths and sites, the predicted concentrations, ranging from nanograms to hundreds of nanograms per litre, were found to be in good agreement. Weather data were utilised with the model to rationalise the effects of various environmental parameters on the sampling results, and, e.g. the roles of various properties of lake dynamics and photodegradation were identified. The new model also enables simultaneous assessment of transformation products. Environmentally formed transformation product concentrations were estimated to be at highest an order of magnitude lower than those of the parent compounds, and unlikely to reach a detectable level. However, a possibility that conjugates of ibuprofen are present at higher levels than the parent compound was identified. Simulation results suggest that environmental degradation half-lives of the inspected contaminants under stratified lake conditions are in the range of some weeks to months.

Dynamics behind the long-lived coherences of I2 in solid Xe.

The absorption spectrum of I2 in solid Xe shows resolved zero-phonon lines and phonon side bands near the origin of the B←X transition (550-625 nm). The long-lived |B⟩⟨X| coherence in this energy range (T2 = 600 fs on average) emerges as vibrationally unrelaxed fluorescence in resonance Raman (RR) spectra. Upon excitation in the structureless continuum at 532 nm, the oscillatory RR progression exhibits electronic dephasing time of T2 = 150 fs. Two RR progressions with markedly different vibrational coherence on the X-state are observed. The main progression of sharp overtones (T2 > 21 ps) is assigned to molecules trapped in double-substitution sites. The minor progression, which shows dephasing times T2 = 6-0.6 ps for v = 1-8, is assigned to molecules in triple-substitution sites. The line progressions allow a detailed characterization of the solvated B- and X-state potentials. Time-resolved coherent anti-Stokes Raman scattering is used to probe selected vibrational coherences on the X-state. Assignments are obtained through molecular dynamics simulations, which reproduce the relative dephasing rates between the two sites, clarify the role of rotation-translation dynamics, and enable quantum dynamics simulations of the spectra by the potentials of mean force that accurately describe the molecule-surrounding interactions.

Raman spectroscopy and crystal-field split rotational states of photoproducts CO and H2 after dissociation of formaldehyde in solid argon.

Raman signal is monitored after 248 nm photodissociation of formaldehyde in solid Ar at temperatures of 9-30 K. Rotational transitions J = 2 ← 0 for para-H(2) fragments and J = 3 ← 1 for ortho-H(2) are observed as sharp peaks at 347.2 cm(-1) and 578.3 cm(-1), respectively, which both are accompanied by a broader shoulder band that shows a split structure. The rovibrational spectrum of CO fragments has transitions at 2136.5 cm(-1), 2138.3 cm(-1), 2139.9 cm(-1), and 2149 cm(-1). To explain the observations, we performed adiabatic rotational potential calculations to simulate the Raman spectrum. The simulations indicate that the splitting of rotational transitions is a site effect, where H(2) molecules can reside in a substitution site, in addition to an interstitial site. In the former site, rotational motion is unperturbed by the electrostatic field of the host atoms, while the latter site splits the excited rotational manifolds, J = 2 and 3, into doublet and triplet structures, respectively. For CO, the spectrum can be ascribed to monomeric species in single- and double-substitution sites, to a dimeric species (CO)(2), and to a CO-H(2)O complex. The simulations show that a nearest-neighbor molecular complex CO-H(2) is not responsible for any of the observed spectral fingerprints. The cause of the exit of the molecular hydrogen from the initial cage can be traced to high translational energy of the fragment after the photodissociation. After the matrix has reached a thermal equilibrium, a diffusion driven formation of the complex is possibly hindered by the high rotational zero-point energy developed upon complexation.

Long-Lived Electronic Coherence of Iodine in the Condensed Phase: Sharp Zero-Phonon Lines in the B↔X Absorption and Emission of I2 in Solid Xe.

Our study of B←X absorption of molecular iodine (I2) isolated in a low-temperature crystalline xenon has revealed an exceptionally long-lived electronic coherence in condensed phase conditions. The visible absorption spectrum shows prominent vibronic structure in the form of zero-phonon lines (ZPLs) and phonon side bands (PSBs). The resolved spectrum implies weak interaction of the chromophore to the lattice degrees of freedom. The coherence extends past the vibrational period of the excited state molecule, unlike that observed in any condensed phase environment for I2 so far. The ZP transitions from the relaxing B-state populations were resolved in the hot luminescence when the 532 nm laser was used for excitation.

Rotational coherence imaging and control for CN molecules through time-frequency resolved coherent anti-Stokes Raman scattering.

Numerical wave packet simulations are performed for studying coherent anti-Stokes Raman scattering (CARS) for CN radicals. Electronic coherence is created by femtosecond laser pulses between the X(2)Σ and B(2)Σ states. Due to the large energy separation of vibrational states, the wave packets are superpositions of rotational states only. This allows for a specially detailed inspection of the second- and third-order coherences by a two-dimensional imaging approach. We present the time-frequency domain images to illustrate the intra- and intermolecular interferences, and discuss the procedure to rationally control and experimentally detect the interferograms in solid Xe environment.

Rotation of methyl radicals in molecular solids.

Electron spin resonance (ESR) measurements were carried out to study the rotation of methyl radicals (CH(3)) in solid carbon monoxide, carbon dioxide, and nitrogen matrices. The radicals were produced by dissociating methane by plasma bursts generated by a focused 193 nm ArF excimer laser radiation during the gas condensation on the substrate. The ESR spectra exhibit anisotropic features that persist over the temperature range examined, and in most cases this indicates a restriction of rotation about the C(2) symmetry axis. A nonrotating CH(3) was also observed in a CO(2) matrix. The intensity ratio between the symmetric (A) and antisymmetric (E) nuclear spin states was recorded as a function of temperature for each molecular matrix. The rotational energy levels are modified from their gas phase structure with increasing crystal field strength. An anomalous situation was observed where the A/E ratio extended below the high temperature limit of 1/2.

Iodine-benzene complex as a candidate for a real-time control of a bimolecular reaction. Spectroscopic studies of the properties of the 1:1 complex isolated in solid krypton.

The properties of the 1:1 iodine-benzene complex isolated in a solid Kr matrix at low temperatures have been studied using UV-vis absorption, FTIR, resonance Raman, and femtosecond coherent anti-Stokes Raman spectroscopy (fs-CARS). The use of all these techniques on similar samples provides a wide view on the spectroscopic properties of the complex and allows comparison and combination of the results from different methods. The results for the complex cover its structure, the changes in the iodine molecule's vibrational frequencies and electronic absorption spectrum upon complexation, and the dynamics of the complexed I(2) molecule on both ground and excited electronic states. In addition, polarization beats between uncomplexed benzene and iodine molecules are detected in the fs-CARS spectra, showing an amplification of an electronically nonresonant CARS signal by the resonant iodine signal. The possibility of controlling the charge-transfer reaction of the I(2)-Bz complex using the excitation of a well-defined ground-state vibrational wavepacket, according to the Tannor-Rice-Kosloff scheme, is discussed on the basis of the experimental findings.

Rotation of methyl radicals in a solid krypton matrix.

Electron spin resonance (ESR) measurements were carried out to study the rotation of methyl radicals (CH(3)) in a solid krypton matrix at 17-31 K temperature range. The radicals were produced by dissociating methane by plasma bursts generated by a focused 193 nm excimer laser radiation during the krypton gas condensation on the substrate. The ESR spectrum exhibits only isotropic features at the temperature range examined, and the intensity ratio between the symmetric (A) and antisymmetric (E) spin state lines exhibits weaker temperature dependence than in a solid argon matrix. However, the general appearance of the methyl radical spectrum depends strongly on temperature due to the pronounced temperature dependency of the E state linewidths. The rotational energy level populations are analyzed based on the static crystal field model, pseudorotating cage model, and quantum chemical calculations for an axially symmetric, planar rotor. Crystal field strength parameter values of -140 cm(-1) in Ar and -240 cm(-1) in Kr match most closely the experimentally observed rotational energy level shifts from the gas phase value. In the alternative model, considering the lattice atom movement in a pseudorotating cage, the effective lowering of the rotational constants B and C to 80%-90% leads to similar effects.

From monomer to bulk: appearance of the structural motif of solid iodine in small clusters.

Formation of iodine clusters in a solid krypton matrix was studied using resonance Raman spectroscopy with a 1 cm(-1) resolution. The clusters were produced by annealing of the solid and recognized by appearance of additional spectral transitions. Two distinct regions, red-shifted from the fundamental vibrational wavenumber of the isolated I(2) at 211 cm(-1), were observed in the signal. The intermediate region spans the range 196-208 cm(-1), and the ultimate region consists of two peaks at 181 and 190 cm(-1) nearly identical to crystalline I(2). The experimental results were compared to DFT-D level electronic structure calculations of planar (I(2))(n) clusters (n = 1-7). The dimer, trimer, and tetramer structures, where the I(2) molecule is complexed from one end, were found to exhibit vibrational shifts corresponding to the intermediate size clusters. The larger, bulklike shift appears when the iodine molecule is coordinated from two opposite directions as in the case of a pentamer and higher clusters. Starting from the pentamer, the structural motif of crystalline iodine is clearly recognized in the clusters.

Vibrational characterization of the 1:1 iodine-benzene complex isolated in solid krypton.

The structure and properties of a 1:1 iodine-benzene complex isolated in an inert krypton matrix at low temperature have been studied with infrared and resonance Raman spectroscopy and with MP2 calculations. The structure of the ground-state complex is found to be unsymmetric, and the I-I vibrational frequency is found to be red-shifted by 3.94 cm(-1) upon complexation. The experimental data agree well with computational results, leading to the conclusion that the I2-Bz complex structure is not axial but of above-bond type, identically with other halogen-benzene complexes.

Rotation of methyl radicals in a solid argon matrix.

Electron spin resonance (ESR) measurements were carried out to study the rotation of methyl radicals (CH3) in a solid argon matrix at 14-35 K temperatures. The radicals were produced by dissociating methane by plasma bursts generated either by a focused 193 nm laser radiation or a radio frequency discharge device during the gas condensation on the substrate. The ESR spectrum exhibits axial symmetry at the lowest temperature and is ascribed to ground state molecules with symmetric total nuclear spin function I=3/2. The hyperfine anisotropy (Aparallel)-Aperpendicular) was found to be -0.01 mT, whereas that of the g value was 2.5x10(-5). The anisotropy is observed for the first time in Ar and is manifested by the splitting of the low-field transition. Elevation of temperature leads reversibly to the appearance of excited state contribution having antisymmetric I=1/2. As a function of the sample temperature, the relative intensities of symmetric and antisymmetric spin states corresponding to ground and excited rotor states, respectively, proton hyperfine and electron g-tensor components, and spin-lattice relaxation rates were determined by a numerical fitting procedure. The experimental observations were interpreted in terms of a free rotation about the C3 axis and a thermal activation of the C2-type rotations above 15 K. The ground and excited rotational state energy levels were found to be separated by 11.2 cm-1 and to exhibit significantly different spin-lattice coupling. A crystal field model has been applied to evaluate the energy levels of the hindered rotor in the matrix, and crystal field parameter varepsilon4=-200 cm-1, corresponding to a 60 cm-1 effective potential barrier for rotation of the C3 axis, was obtained.

Femtosecond coherent anti-Stokes Raman-scattering polarization beat spectroscopy of I2-Xe complex in solid krypton.

Time-resolved coherent anti-Stokes Raman-scattering (CARS) measurements are carried out to study the interaction between xenon atom and iodine molecule in a solid krypton matrix. Interference between the CARS polarizations of the "free" and complexed iodine molecules is observed, while the quantum beats of the complex are not detected due to low concentration. Vibrational analysis based on the polarization beats yields accurate molecular constants for the I2-Xe complex. The harmonic frequency of the I2-Xe complex is found to be redshifted by 0.90 cm-1 when compared to the free I2, whereas the anharmonicity is approximately the same. The dephasing rate of the complex is found to be somewhat higher than that of the free iodine molecule in solid Kr, showing that the complexation affects dephasing, although not dramatically. Molecular dynamics simulations are carried out to find the conformation of the complex, and wave packet simulations are used to reproduce the CARS signal to confirm the assignments of the observed beatings as quantum and polarization beats. The results show that the polarization beats are a useful tool for investigating weak interactions in condensed phase.

Time-dependent alignment of molecules trapped in octahedral crystal fields.

The hindered rotational states of molecules confined in crystal fields of octahedral symmetry, and their time-dependent alignment obtained by pulsed nonresonant laser fields, are studied computationally. The control over the molecular axis direction is discussed based on the evolution of the rotational wave packet generated in the cubic crystal-field potential. The alignment degree obtained in a cooperative case, where the alignment field is applied in a favorable crystal-field direction, or in a competitive direction, where the crystal field has a saddle point, is presented. The investigation is divided into two time regimes where the pulse duration is either ultrashort, leading to nonadiabatic dynamics, or long with respect to period of molecular libration, which leads to synchronous alignment due to nearly adiabatic following. The results are contrasted to existing gas phase studies. In particular, the irregularity of the crystal-field energies leads to persistent interference patterns in the alignment signals. The use of nonadiabatic alignment for interrogation of crystal-field energetics and the use of adiabatic alignment for directional control of molecular dynamics in solids are proposed as practical applications.

Intense-field alignment of molecules confined in octahedral fields.

The combined effect of static octahedral potential and anisotropic interactions with intense linearly polarized light on the rotational motion of linear molecules is investigated. Avoided crossings of quantized energy levels corresponding to pendular states with different degrees of alignment are found by varying the strength parameters for the light and static fields. High alignment is achieved for both cooperative and competitive choices of the relative directionality of the two fields, thus presenting means for controlling the dynamics of impurity centers with respect to the surrounding media.