On the quantum dynamical treatment of surface vibrational modes for reactive scattering of H2 from Cu(111) at 925 K.

We construct the effective Hartree potential for H2 on Cu(111) as introduced in our earlier work [Dutta et al., J. Chem. Phys. 154, 104103 (2021), and Dutta et al., J. Chem. Phys. 157, 194112 (2022)] starting from the same gas-metal interaction potential obtained for 0 K. Unlike in that work, we now explicitly account for surface expansion at 925 K and investigate different models to describe the surface vibrational modes: (i) a cluster model yielding harmonic normal modes at 0 K and (ii) slab models resulting in phonons at 0 and 925 K according to the quasi-harmonic approximation-all consistently calculated at the density functional theory level with the same exchange-correlation potential. While performing dynamical calculations for the H2(v = 0, j = 0)-Cu(111) system employing Hartree potential constructed with 925 K phonons and surface temperature, (i) the calculated chemisorption probabilities are the highest compared to the other approaches over the energy domain and (ii) the threshold for the reaction probability is the lowest, in close agreement with the experiment. Although the survival probabilities (v' = 0) depict the expected trend (lower in magnitude), the excitation probabilities (v' = 1) display a higher magnitude since the 925 K phonons and surface temperature are more effective for the excitation process compared to the phonons/normal modes obtained from the other approaches investigated to describe the surface.

Effect of surface temperature on quantum dynamics of D2 on Cu(111) using a chemically accurate potential energy surface.

The effect of surface mode vibrations on the reactive scattering of D2, initialized in the ground rovibrational state (v = 0, j = 0), from a Cu(111) surface is investigated for different surface temperature situations. We adopt a time and temperature dependent effective Hamiltonian [Dutta et al., J. Chem. Phys. 154, 104103 (2021)] constructed by combining the linearly coupled many oscillator model [Sahoo et al., J. Chem. Phys. 136, 084306 (2012)] and the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)] potential within the mean-field approach. Such an effective Hamiltonian is employed for six-dimensional quantum dynamical calculations to obtain temperature dependent reaction and state-to-state scattering probability profiles as a function of incidence energy of colliding D2 molecules. As reported in the experimental studies, the movements of surface atoms modify the dissociative scattering dynamics at higher surface temperature by exhibiting vibrational quantum and surface atoms' recoil effects in the low and high collision energy domains, respectively. Finally, we compare our present theoretical results with the experimental and other theoretical outcomes, as well as discuss the novelty of our findings.

A Novel Apparatus for Simultaneous Laser-Light-Sheet Optical Particle Counting and Video Recording in the Same Measurement Chamber at High Temperature.

A novel apparatus was developed, to investigate the detachment of particle structures consisting of soot and ash from a single fibre or a fibre array in hot gas flow. Key features of the novel apparatus are operation at high temperatures while two different measurement techniques are applied simultaneously in the same measurement chamber to observe particle structure detachment from a loaded fibre array. A heated inlet can heat the air stream at the position of the fibre array up to 470 °C, allowing detachment investigations at temperatures relevant for the operation of, e.g., soot particle filters. The first measurement technique integrated in the setup is video recording of the fibre array, which gives qualitative information on the rearrangement or detachment of particulate matter on the fibre. Because it is often difficult to distinguish rearrangement and detachment from pure visual observations, a second measurement technique is applied. This technique is a laser-light-sheet optical particle counter, which can detect detached particle structures and determine their size. The measurable size range is 257 to 1523 µm for glass spheres. This paper presents and discusses the novel apparatus, its calibration and first detachment measurement results.

IR spectroscopic characterization of the co-adsorption of CO2 and H2 onto cationic Cun+ clusters.

To understand elementary reaction steps in the hydrogenation of CO2 over copper-based catalysts, we experimentally study the adsorption of CO2 and H2 onto cationic Cun+ clusters. For this, we react Cun+ clusters formed by laser ablation with a mixture of H2 and CO2 in a flow tube-type reaction channel and characterize the products formed by IR multiple-photon dissociation spectroscopy employing the IR free-electron laser FELICE. We analyze the spectra by comparing them to literature spectra of Cun+ clusters reacted with H2 and with new spectra of Cun+ clusters reacted with CO2. The latter indicate that CO2 is physisorbed in an end-on configuration when reacted with the clusters alone. Although the spectra for the co-adsorption products evidence H2 dissociation, no signs for CO2 activation or reduction are observed. This lack of reactivity for CO2 is rationalized by density functional theory calculations, which indicate that CO2 dissociation is hindered by a large reaction barrier. CO2 reduction to formate should energetically be possible, but the lack of formate observation is attributed to kinetic hindering.

IR Spectroscopic Characterization of H2 Adsorption on Cationic Cun+ (n = 4-7) Clusters.

IR spectra of cationic copper clusters Cun+ (n = 4-7) complexed with hydrogen molecules are recorded via IR multiple-photon dissociation (IRMPD) spectroscopy. To this end, the copper clusters are generated via laser ablation and reacted with H2 and D2 in a flow-tube-type reaction channel. The complexes formed are irradiated using IR light provided by the free-electron laser for intracavity experiments (FELICE). The spectra are interpreted by making use of isotope-induced shifts of the vibrational bands and by comparing them to density functional theory calculated spectra for candidate structures. The structural candidates have been obtained from global sampling with the minima hopping method, and spectra are calculated at the semilocal (PBE) and hybrid (PBE0) functional level. The highest-quality spectra have been recorded for [5Cu, 2H/2D]+, and we find that the semilocal functional provides better agreement for the lowest-energy isomers. The interaction of hydrogen with the copper clusters strongly depends on their size. Binding energies are largest for Cu5+, which goes hand in hand with the observed predominantly dissociative adsorption. Due to smaller binding energies for dissociated H2 and D2 for Cu4+, also a significant amount of molecular adsorption is observed as to be expected according to the Evans-Polanyi principle. This is confirmed by transition-state calculations for Cu4+ and Cu5+, which show that hydrogen dissociation is not hindered by an endothermic reaction barrier for Cu5+ and by a slightly endothermic barrier for Cu4+. For Cu6+ and Cu7+, it was difficult to draw clear conclusions because the IR spectra could not be unambiguously assigned to structures.

Effect of surface temperature on quantum dynamics of H2 on Cu(111) using a chemically accurate potential energy surface.

The effect of surface atom vibrations on H2 scattering from a Cu(111) surface at different temperatures is being investigated for hydrogen molecules in their rovibrational ground state (v = 0, j = 0). We assume weakly correlated interactions between molecular degrees of freedom and surface modes through a Hartree product type wavefunction. While constructing the six-dimensional effective Hamiltonian, we employ (a) a chemically accurate potential energy surface according to the static corrugation model [M. Wijzenbroek and M. F. Somers, J. Chem. Phys. 137, 054703 (2012)]; (b) normal mode frequencies and displacement vectors calculated with different surface atom interaction potentials within a cluster approximation; and (c) initial state distributions for the vibrational modes according to Bose-Einstein probability factors. We carry out 6D quantum dynamics with the so-constructed effective Hamiltonian and analyze sticking and state-to-state scattering probabilities. The surface atom vibrations affect the chemisorption dynamics. The results show physically meaningful trends for both reaction and scattering probabilities compared to experimental and other theoretical results.

Importance of zero-point energy for crystalline ice phases: A comparison of force fields and density functional theory.

Density functional theory (DFT) including van der Waals (vdW) interactions and accounting for zero-point energy (ZPE) is believed to provide a good description of crystalline ice phases [B. Pamuk et al., Phys. Rev. Lett. 108, 193003 (2012)]. Given the computational cost of DFT, it is not surprising that extensive phonon calculations, which yield the ZPE, have only been done for a limited amount of ice structures. Computationally convenient force fields on the other hand are the method of choice for large systems and/or dynamical simulations, e.g., of supercooled water. Here, we present a systematic comparison for seven hydrogen-ordered crystalline ice phases (Ih, IX, II, XIII, XIV, XV, and VIII) between many commonly used nonpolarizable force fields and density functionals, including some recently developed meta-GGA functionals and accounting for vdW interactions. Starting from the experimentally determined crystal structures, we perform space-group-constrained structural relaxations. These provide the starting point for highly accurate phonon calculations that yield effectively volume-dependent ZPEs within the quasiharmonic approximation. In particular, when including ZPE, the force fields show a remarkably good performance for equilibrium volumes and cohesive energies superior to many density functionals. A decomposition of the cohesive energies into intramolecular deformation, electrostatic, and vdW contributions quantifies the differences between force fields and DFT. Results for the equilibrium volumes and phase transition pressures for all studied force fields are much more strongly affected by ZPE than all studied density functionals. We track this down to significantly smaller shifts of the O-H-stretch modes and compare with experimental data from Raman spectroscopy.

Nonadiabatic Vibrational Damping of Molecular Adsorbates: Insights into Electronic Friction and the Role of Electronic Coherence.

We present a perturbation approach rooted in time-dependent density-functional theory to calculate electron-hole (e-h) pair excitation spectra during the nonadiabatic vibrational damping of adsorbates on metal surfaces. Our analysis for the benchmark systems CO on Cu(100) and Pt(111) elucidates the surprisingly strong influence of rather short electronic coherence times. We demonstrate how in the limit of short electronic coherence times, as implicitly assumed in prevalent quantum nuclear theories for the vibrational lifetimes as well as electronic friction, band structure effects are washed out. Our results suggest that more accurate lifetime or chemicurrentlike experimental measurements could characterize the electronic coherence.

Design and validation of a parallelized micro-photobioreactor enabling phototrophic bioprocess development at elevated throughput.

Microalgae offer great potential for the industrial production of numerous compounds, but most of the currently available processes fail on economic aspects. Due to the lack of appropriate microcultivation systems, especially screening and early stage laboratory process characterization limit throughput in process development. Consequently, a demand for high throughput photobioreactors has recently been identified upon which some prototype systems emerged. However, compared to microbial microbioreactors, the systems so far introduced suffer from at least one of several drawbacks, that is, inhomogeneous conditions, poor mixing or excessive evaporation. In this context, a microtiter plate based micro-photobioreactor was developed enabling 48-fold parallelized cultivation. Strict control of the process conditions enabled a high comparability between the distinct wells of one plate (±5% fluctuation in biomass formation). The small scale, resulting in a beneficial surface to volume ratio, as well as the fast mixing due to rigorous orbital shaking, ensured an excellent light supply of the cultures. Moreover, non-invasive online biomass quantification was implemented via a scattered light analyzer that is capable of biomass measurements during continuous illumination of the cultures. The system was shown to be especially qualified for parallelized laboratory screening applications like for instance media optimization. Easy automation via integration into a liquid handling platform is given by design. Thereby, the presented micro-photobioreactor system significantly contributes to improving the time efficiency during the development of phototrophic bioprocesses. Biotechnol. Bioeng. 2017;114: 122-131. © 2016 Wiley Periodicals, Inc.

Dynamics of H2 dissociation on the close-packed (111) surface of the noblest metal: H2 + Au(111).

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.

Electronic Friction-Based Vibrational Lifetimes of Molecular Adsorbates: Beyond the Independent-Atom Approximation.

We assess the accuracy of vibrational damping rates of diatomic adsorbates on metal surfaces as calculated within the local-density friction approximation (LDFA). An atoms-in-molecules (AIM) type charge partitioning scheme accounts for intramolecular contributions and overcomes the systematic underestimation of the nonadiabatic losses obtained within the prevalent independent-atom approximation. The quantitative agreement obtained with theoretical and experimental benchmark data suggests the LDFA-AIM scheme as an efficient and reliable approach to account for electronic dissipation in ab initio molecular dynamics simulations of surface chemical reactions.

Fingerprints of energy dissipation for exothermic surface chemical reactions: O2 on Pd(100).

We present first-principles calculations of the sticking coefficient of O2 at Pd(100) to assess the effect of phononic energy dissipation on this kinetic parameter. For this, we augment dynamical simulations on six-dimensional potential energy surfaces (PESs) representing the molecular degrees of freedom with various effective accounts of surface mobility. In comparison to the prevalent frozen-surface approach, energy dissipation is found to qualitatively affect the calculated sticking curves. At the level of a generalized Langevin oscillator model, we achieve good agreement with experimental data. The agreement is similarly reached for PESs based on two different semi-local density-functional theory functionals. This robustness of the simulated sticking curve does not extend to the underlying adsorption mechanism, which is predominantly directly dissociative for one functional or molecularly trapped for the other. Completely different adsorption mechanisms therewith lead to rather similar sticking curves that agree equally well with the experimental data. This highlights the danger of the prevalent practice to extract corresponding mechanistic details from simple fingerprints of measured sticking data for such exothermic surface reactions.

Cytoskeletal transition in patterned cells correlates with interfacial energy model.

A cell's morphology is intricately regulated by microenvironmental cues and intracellular feedback signals. Besides biochemical factors, cell fate can be influenced by the mechanics and geometry of the surrounding matrix. The latter point was addressed herein, by studying cell adhesion on two-dimensional micropatterns. Endothelial cells were grown on maleic acid copolymer surfaces structured with stripes of fibronectin by microcontact printing. Experiments showed a biphasic behaviour of actin stress fibre spacing in dependence on the stripe width with a critical size of approx. 15 µm. In a concurrent modelling effort, cells on stripes were simulated as droplet-like structures, including variations of interfacial energy, total volume and dimensions of the nucleus. A biphasic behaviour with regard to cell morphology and area was found, triggered by the minimum of interfacial energy, with the phase transition occurring at a critical stripe width close to the critical stripe width found in the cell experiment. The correlation of experiment and simulation suggests a possible mechanism of the cytoskeletal rearrangements based on interfacial energy arguments.

Modeling heat dissipation at the nanoscale: an embedding approach for chemical reaction dynamics on metal surfaces.

We present an embedding technique for metallic systems that makes it possible to model energy dissipation into substrate phonons during surface chemical reactions from first principles. The separation of chemical and elastic contributions to the interaction potential provides a quantitative description of both electronic and phononic band structure. Application to the dissociation of O2 at Pd(100) predicts translationally "hot" oxygen adsorbates as a consequence of the released adsorption energy (ca. 2.6 eV). This finding questions the instant thermalization of reaction enthalpies generally assumed in models of heterogeneous catalysis.

From Jekyll to Hyde and Beyond: Hydrogen's Multifaceted Role in Passivation, H-Induced Breakdown, and Charging of Amorphous Silicon Nitride.

In semiconductor devices, hydrogen has traditionally been viewed as a panacea for defects, being adept at neutralizing dangling bonds and consequently purging the related states from the band gap. With amorphous silicon nitride (a-Si3N4)─a material critical for electronic, optical, and mechanical applications─this belief holds true as hydrogen passivates both silicon and nitrogen dangling bonds. However, there is more to the story. Our density functional theory calculations unveil hydrogen's multifaceted role upon incorporation in a-Si3N4. On the "Jekyll" side, hydrogen atoms are indeed restorative, healing coordination defects in a-Si3N4. However, "Hyde" emerges as hydrogen induces Si-N bond breaking, particularly in strained regions of the amorphous network. Beyond these dual roles, our study reveals an intricate balance between hydrogen defect centers and intrinsic charge traps that already exist in pristine a-Si3N4: the excess charges provided by the H atoms result in charging of the a-Si3N4 dielectric layer.

Best-of-Both-Worlds Predictive Approach to Dissociative Chemisorption on Metals.

Predictive capability, accuracy, and affordability are essential features of a theory that is capable of describing dissociative chemisorption on a metal surface. This type of reaction is important for heterogeneous catalysis. Here we present an approach in which we use diffusion Monte Carlo (DMC) to pin the minimum barrier height and construct a density functional that reproduces this value. This predictive approach allows the construction of a potential energy surface at the cost of density functional theory while retaining near DMC accuracy. Scrutinizing effects of energy dissipation and quantum tunneling, dynamics calculations suggest the approach to be of near chemical accuracy, reproducing molecular beam sticking experiments for the showcase H2 + Al(110) system to ∼1.4 kcal/mol.

Adsorption structure determination of a large polyaromatic trithiolate on Cu(111): combination of LEED-I(V) and DFT-vdW.

The adsorption geometry of 1,3,5-tris(4-mercaptophenyl)benzene (TMB) on Cu(111) is determined with high precision using two independent methods, experimentally by quantitative low energy electron diffraction (LEED-I(V)) and theoretically by dispersion corrected density functional theory (DFT-vdW). Structural refinement using both methods consistently results in similar adsorption sites and geometries. Thereby a level of confidence is reached that allows deduction of subtle structural details such as molecular deformations or relaxations of copper substrate atoms.

Experimental investigation on the change of pull-off force between bulk particulate material and an elastic polymeric filter fiber.

HYPOTHESIS: Adhesion between particles and a filter fiber is an important process of the filtration as it dictates the process of separation and in the following the detachment process of particles during filter regeneration. In addition to the shear stress that a new polymeric stretchable filter fiber implements into the particulate structure, the elongation of the substrate (fiber) is also expected to cause a structural change in the surface of the polymer. Thus, the changed contact area and surface energy could affect the adhesion force between particles and fibers. EXPERIMENTS: Systematic measurements of adhesion forces between a single particle and the stretchable substrate were performed using Atomic Force Microscope (AFM). The substrate surface characteristics (roughness) was changed directly beneath the modified measurement head using piezo-motors to achieve stepless elongation state. Polystyrene particles and particles made of Spheriglass were applied. FINDINGS: In the experiments, a reduced adhesion force between the particles and the filter fiber was found for a new high range of substrate roughness and peak-to-peak distance, in which the Rabinovich model has not been used before [1]. Further, the influence of high and low energy surface particulate material was evaluated to understand the detachment process in the new real adaptive filter and in DEM-simulation.

Floating in Space: How to Treat the Weak Interaction between CO Molecules in Interstellar Ices.

In the interstellar medium, six molecules have been conclusively detected in the solid state in interstellar ices, and a few dozen have been hypothesized and modeled to be present in the solid state as well. The icy mantles covering micrometer-sized dust grains are, in fact, thought to be at the core of complex molecule formation as a consequence of the local high density of molecules that are simultaneously adsorbed. From a structural perspective, the icy mantle is considered to be layered, with an amorphous water-rich inner layer surrounding the dust grain, covered by an amorphous CO-rich outer layer. Moreover, recent studies have suggested that the CO-rich layer might be crystalline and possibly even be segregated as a single crystal atop the ice mantle. If so, there are far-reaching consequences for the formation of more complex organic molecules, such as methanol and sugars, that use CO as a backbone. Validation of these claims requires further investigation, in particular on acquiring atomistic insight into surface processes, such as adsorption, diffusion, and reactivity on CO ices. Here, we present the first detailed computational study toward treating the weak interaction of (pure) CO ices. We provide a benchmark of the performance of various density functional theory methods in treating the binding of pure CO ices. Furthermore, we perform an atomistic and in-depth study of the binding energy of CO on amorphous and crystalline CO ices using a pair-potential-based force field. We find that CO adsorption is represented by a large distribution of binding energies (200-1600 K) on amorphous CO, including a significant amount of weak binding sites (<350 K). Increasing both the cluster size and the number of neighbors increases the mean of the observed binding energy distribution. Finally, we find that CO binding energies are dominated by dispersion and, as such, exchange-correlation functionals need to include a treatment of dispersion to accurately simulate surface processes on CO ices. In particular, we find the ωB97M-V functional to be a strong candidate for such simulations.

Characterization of CO Adsorbed to Clean and Partially Oxidized Cu(211) and Cu(111).

Copper-based catalysts gain activity through the presence of poorly coordinated Cu atoms and incomplete oxidation at the surface. The catalytic mechanisms can in principle be observed by controlled dosing of reactants to single-crystal substrates. However, the interconnected influences of surface defects, partial oxidation, and adsorbate coverage present a large matrix of conditions that have not been fully explored in the literature. We recently characterized oxygen and carbon monoxide coadsorption on Cu(111), a nominally defect-free surface, and now extend our study to the stepped surface Cu(211). Temperature-programmed desorption of CO adsorbed to bare metal surfaces confirms that two sites dominate desorption from a saturated layer: atop terrace atoms of local (111) character and atop step edge atoms with CO bound more strongly to the latter. At low coverage, discrete CO resonances in reflection adsorption infrared spectra can be assigned to these sites: 2077 cm-1 for extended (111) terraces, 2093 cm-1 for step sites, and additional kink-adsorbed molecules at 2110 cm-1. With increasing coverage, in contrast to Cu(111), the infrared spectral features on Cu(211) evolve and shift as a consequence of dipole-dipole coupling between differentially occupied types of sites. Auger electron spectroscopy shows that exposure to background O2 oxidizes the (211) surface at a rate nearly 1 order of magnitude greater than (111); we argue that the resulting surface is stoichiometric Cu2O, as previously found for Cu(111). This oxide binds CO less strongly than the bare metal and the underlying crystal cut continues to influence the adsorption sites available to CO. On oxidized (111) terraces, broad absorption peaks at 2115-2120 cm-1; on oxidized Cu(211), CO adsorbed to step sites appears as a resolved secondary peak at 2144 cm-1. This suite of spectroscopic signatures, obtained under carefully controlled conditions, will help to determine the origin and fate of adsorbed species in future studies of reaction mechanisms on copper.