@article{FischerAndersSaalfrank2022, author = {Fischer, Eric Wolfgang and Anders, Janet and Saalfrank, Peter}, title = {Cavity-altered thermal isomerization rates and dynamical resonant localization in vibro-polaritonic chemistry}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {156}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {15}, publisher = {American Institute of Physics}, address = {Melville, NY}, issn = {0021-9606}, doi = {10.1063/5.0076434}, pages = {16}, year = {2022}, abstract = {It has been experimentally demonstrated that reaction rates for molecules embedded in microfluidic optical cavities are altered when compared to rates observed under "ordinary" reaction conditions. However, precise mechanisms of how strong coupling of an optical cavity mode to molecular vibrations affects the reactivity and how resonance behavior emerges are still under dispute. In the present work, we approach these mechanistic issues from the perspective of a thermal model reaction, the inversion of ammonia along the umbrella mode, in the presence of a single-cavity mode of varying frequency and coupling strength. A topological analysis of the related cavity Born-Oppenheimer potential energy surface in combination with quantum mechanical and transition state theory rate calculations reveals two quantum effects, leading to decelerated reaction rates in qualitative agreement with experiments: the stiffening of quantized modes perpendicular to the reaction path at the transition state, which reduces the number of thermally accessible reaction channels, and the broadening of the barrier region, which attenuates tunneling. We find these two effects to be very robust in a fluctuating environment, causing statistical variations of potential parameters, such as the barrier height. Furthermore, by solving the time-dependent Schrodinger equation in the vibrational strong coupling regime, we identify a resonance behavior, in qualitative agreement with experimental and earlier theoretical work. The latter manifests as reduced reaction probability when the cavity frequency omega(c) is tuned resonant to a molecular reactant frequency. We find this effect to be based on the dynamical localization of the vibro-polaritonic wavepacket in the reactant well.}, language = {en} } @article{ScholzFlossSaalfranketal.2016, author = {Scholz, Robert and Floss, Gereon and Saalfrank, Peter and F{\"u}chsel, Gernot and Loncaric, Ivor and Juaristi, J. I.}, title = {Femtosecond-laser induced dynamics of CO on Ru(0001): Deep insights from a hot-electron friction model including surface motion}, series = {Physical review : B, Condensed matter and materials physics}, volume = {94}, journal = {Physical review : B, Condensed matter and materials physics}, publisher = {American Physical Society}, address = {College Park}, issn = {2469-9950}, doi = {10.1103/PhysRevB.94.165447}, pages = {17}, year = {2016}, abstract = {A Langevin model accounting for all six molecular degrees of freedom is applied to femtosecond-laser induced, hot-electron driven dynamics of Ru(0001)(2 x 2): CO. In our molecular dynamics with electronic friction approach, a recently developed potential energy surface based on gradient-corrected density functional theory accounting for van der Waals interactions is adopted. Electronic friction due to the coupling of molecular degrees of freedom to electron-hole pairs in the metal are included via a local density friction approximation, and surface phonons by a generalized Langevin oscillator model. The action of ultrashort laser pulses enters through a substrate-mediated, hot-electron mechanism via a time-dependent electronic temperature (derived from a two-temperature model), causing random forces acting on the molecule. The model is applied to laser induced lateral diffusion of CO on the surface, "hot adsorbate" formation, and laser induced desorption. Reaction probabilities are strongly enhanced compared to purely thermal processes, both for diffusion and desorption. Reaction yields depend in a characteristic (nonlinear) fashion on the applied laser fluence, as well as branching ratios for various reaction channels. Computed two-pulse correlation traces for desorption and other indicators suggest that aside from electron-hole pairs, phonons play a non-negligible role for laser induced dynamics in this system, acting on a surprisingly short time scale. Our simulations on precomputed potentials allow for good statistics and the treatment of long-time dynamics (300 ps), giving insight into this system which hitherto has not been reached. We find generally good agreement with experimental data where available and make predictions in addition. A recently proposed laser induced population of physisorbed precursor states could not be observed with the present low-coverage model.}, language = {en} } @article{KuleszaTitovDalyetal.2016, author = {Kulesza, Alexander Jan and Titov, Evgenii and Daly, Steven and Wlodarczyk, Radoslaw and Megow, J{\"o}rg and Saalfrank, Peter and Choi, Chang Min and MacAleese, Luke and Antoine, Rodolphe and Dugourd, Philippe}, title = {Excited States of Xanthene Analogues: Photofragmentation and Calculations by CC2 and Time-Dependent Density Functional Theory}, series = {ChemPhysChem : a European journal of chemical physics and physical chemistry}, volume = {17}, journal = {ChemPhysChem : a European journal of chemical physics and physical chemistry}, publisher = {Wiley-VCH}, address = {Weinheim}, issn = {1439-4235}, doi = {10.1002/cphc.201600650}, pages = {3129 -- 3138}, year = {2016}, abstract = {Action spectroscopy has emerged as an analytical tool to probe excited states in the gas phase. Although comparison of gas-phase absorption properties with quantum-chemical calculations is, in principle, straightforward, popular methods often fail to describe many molecules of interest-such as xanthene analogues. We, therefore, face their nano-and picosecond laser-induced photofragmentation with excited-state computations by using the CC2 method and time-dependent density functional theory (TDDFT). Whereas the extracted absorption maxima agree with CC2 predictions, the TDDFT excitation energies are blueshifted. Lowering the amount of Hartree-Fock exchange in the DFT functional can reduce this shift but at the cost of changing the nature of the excited state. Additional bandwidth observed in the photofragmentation spectra is rationalized in terms of multiphoton processes. Observed fragmentation from higher-lying excited states conforms to intense excited-to-excited state transitions calculated with CC2. The CC2 method is thus suitable for the comparison with photofragmentation in xanthene analogues.}, language = {en} } @article{BoeseSaalfrank2016, author = {Boese, Adrian Daniel and Saalfrank, Peter}, title = {CO Molecules on a NaCl(100) Surface: Structures, Energetics, and Vibrational Davydov Splittings at Various Coverages}, series = {The journal of physical chemistry : C, Nanomaterials and interfaces}, volume = {120}, journal = {The journal of physical chemistry : C, Nanomaterials and interfaces}, publisher = {American Chemical Society}, address = {Washington}, issn = {1932-7447}, doi = {10.1021/acs.jpcc.6b03726}, pages = {12637 -- 12653}, year = {2016}, abstract = {In this work, we study the adsorption of CO from low to high coverage at a defect-free NaCl(100) surface by means of duster and periodic models, using highly accurate wave function-based QM:QM embedding as well as density functional theory. At low coverages, the most accurate methods predict a zero-point-corrected adsorption energy of around 13 kJ/mol, and the CO molecules are found to be oriented perpendicular to the surface. At higher coverages, lower-energy phases with nonparallel/upright, tilted orientations emerge. Besides the well-known p(2 x 1)/antiparallel phase (T/A), we find other tilted phases (tilted/irregular, T/I; tilted/spiral, T/S) as local minima. Vibrational frequencies for CO adsorbed on NaCl(100) and Davydov splittings of the C-O stretch vibration are also determined. The IR spectra are characteristic fingerprints for the relative orientation of CO molecules and may therefore be used as sensitive probes to distinguish parallel/upright from various tilted adsorption phases.}, language = {en} } @article{TitovGranucciGoetzeetal.2016, author = {Titov, Evgenii and Granucci, Giovanni and Goetze, Jan Philipp and Persico, Maurizio and Saalfrank, Peter}, title = {Dynamics of Azobenzene Dimer Photoisomerization: Electronic and Steric Effects}, series = {The journal of physical chemistry letters}, volume = {7}, journal = {The journal of physical chemistry letters}, publisher = {American Chemical Society}, address = {Washington}, issn = {1948-7185}, doi = {10.1021/acs.jpciett.6b01401}, pages = {3591 -- 3596}, year = {2016}, abstract = {While azobenzenes readily photoswitch in solution, their photoisomerization in densely packed self-assembled monolayers (SAMs) can be suppressed. Reasons for this can be steric hindrance and/or electronic quenching, e.g., by exciton coupling. We address these possibilities by means of nonadiabatic molecular dynamics with trajectory surface hopping calculations, investigating the trans -> cis isomerization of azobenzene after excitation into the pi pi* absorption band. We consider a free monomer, an isolated dimer and a dimer embedded in a SAM-like environment of additional azobenzene molecules, imitating in this way the gradual transition from an unconstrained over an electronically coupled to an electronically coupled and sterically hindered, molecular switch. Our simulations reveal that in comparison to the single molecule the quantum yield of the trans -> cis photoisomerization is similar for the isolated dimer, but greatly reduced in the sterically constrained situation. Other implications of dimerization and steric constraints are also discussed.}, language = {en} } @article{LoncaricAlducinSaalfranketal.2016, author = {Loncaric, Ivor and Alducin, Maite and Saalfrank, Peter and Inaki Juaristi, J.}, title = {Femtosecond laser pulse induced desorption: A molecular dynamics simulation}, series = {Nature climate change}, volume = {382}, journal = {Nature climate change}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0168-583X}, doi = {10.1016/j.nimb.2016.02.051}, pages = {114 -- 118}, year = {2016}, abstract = {In recent simulations of femtosecond laser induced desorption of molecular oxygen from the Ag(110) surface, it has been shown that depending on the properties (depth and electronic environment) of the well in which 02 is adsorbed, the desorption can be either induced dominantly by hot electrons or via excitations of phonons. In this work we explore whether the ratios between the desorption yields from different adsorption wells can be tuned by changing initial surface temperature and laser pulse properties. We show that the initial surface temperature is an important parameter, and that by using low initial surface temperatures the electronically mediated process can be favored. In contrast, laser properties seem to have only a modest influence on the results. (C) 2016 Elsevier B.V. All rights reserved.}, language = {en} } @article{ScholzLindnerLoncaricetal.2019, author = {Scholz, Robert and Lindner, Steven and Loncaric, Ivor and Tremblay, Jean Christophe and Juaristi, J. and Alducin, Maite and Saalfrank, Peter}, title = {Vibrational response and motion of carbon monoxide on Cu(100) driven by femtosecond laser pulses: Molecular dynamics with electronic friction}, series = {Physical review : B, Condensed matter and materials physics}, volume = {100}, journal = {Physical review : B, Condensed matter and materials physics}, number = {24}, publisher = {American Physical Society}, address = {College Park}, issn = {2469-9950}, doi = {10.1103/PhysRevB.100.245431}, pages = {20}, year = {2019}, abstract = {Carbon monoxide on copper surfaces continues to be a fascinating, rich microlab for many questions evolving in surface science. Recently, hot-electron mediated, femtosecond-laser pulse induced dynamics of CO molecules on Cu(100) were the focus of experiments [Inoue et al., Phys. Rev. Lett. 117, 186101 (2016)] and theory [Novko et al., Phys. Rev. Lett. 122, 016806 (2019)], unraveling details of the vibrational nonequilibrium dynamics on ultrashort (subpicoseconds) timescales. In the present work, full-dimensional time-resolved hot-electron driven dynamics are studied by molecular dynamics with electronic friction (MDEF). Dissipation is included by a friction term in a Langevin equation which describes the coupling of molecular degrees of freedom to electron-hole pairs in the copper surface, calculated from gradient-corrected density functional theory (DFT) via a local density friction approximation (LDFA). Relaxation due to surface phonons is included by a generalized Langevin oscillator model. The hot-electron induced excitation is described via a time-dependent electronic temperature, the latter derived from an improved two-temperature model. Our parameter-free simulations on a precomputed potential energy surface allow for excellent statistics, and the observed trends are confirmed by on-the-fly ab initio molecular dynamics with electronic friction (AIMDEF) calculations. By computing time-resolved frequency maps for selected molecular vibrations, instantaneous frequencies, probability distributions, and correlation functions, we gain microscopic insight into hot-electron driven dynamics and we can relate the time evolution of vibrational internal CO stretch-mode frequencies to measured data, notably an observed redshift. Quantitatively, the latter is found to be larger in MDEF than in experiment and possible reasons are discussed for this observation. In our model, in addition we observe the excitation and time evolution of large-amplitude low-frequency modes, lateral CO surface diffusion, and molecular desorption. Effects of surface atom motion and of the laser fluence are also discussed.}, language = {en} } @article{MartinazzoNestSaalfranketal.2006, author = {Martinazzo, Rocco and Nest, Mathias and Saalfrank, Peter and Tantardini, Gian Franco}, title = {A local coherent-state approximation to system-bath quantum dynamics}, series = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, volume = {125}, journal = {The journal of chemical physics : bridges a gap between journals of physics and journals of chemistr}, number = {19}, publisher = {American Institute of Physics}, address = {Melville}, issn = {0021-9606}, doi = {10.1063/1.2362821}, pages = {16}, year = {2006}, abstract = {A novel quantum method to deal with typical system-bath dynamical problems is introduced. Subsystem discrete variable representation and bath coherent-state sets are used to write down a multiconfigurational expansion of the wave function of the whole system. With the help of the Dirac-Frenkel variational principle, simple equations of motion-a kind of Schrodinger-Langevin equation for the subsystem coupled to (pseudo) classical equations for the bath-are derived. True dissipative dynamics at all times is obtained by coupling the bath to a secondary, classical Ohmic bath, which is modeled by adding a friction coefficient in the derived pseudoclassical bath equations. The resulting equations are then solved for a number of model problems, ranging from tunneling to vibrational relaxation dynamics. Comparison of the results with those of exact, multiconfiguration time-dependent Hartree calculations in systems with up to 80 bath oscillators shows that the proposed method can be very accurate and might be of help in studying realistic problems with very large baths. To this end, its linear scaling behavior with respect to the number of bath degrees of freedom is shown in practice with model calculations using tens of thousands of bath oscillators.}, language = {en} } @article{LoncaricFuchselJuaristietal.2017, author = {Loncaric, Ivor and Fuchsel, Gernot and Juaristi, J. I. and Saalfrank, Peter}, title = {Strong Anisotropic Interaction Controls Unusual Sticking and Scattering of CO at Ru(0001)}, series = {Physical review letters}, volume = {119}, journal = {Physical review letters}, publisher = {American Physical Society}, address = {College Park}, issn = {0031-9007}, doi = {10.1103/PhysRevLett.119.146101}, pages = {5}, year = {2017}, abstract = {Complete sticking at low incidence energies and broad angular scattering distributions at higher energies are often observed in molecular beam experiments on gas-surface systems which feature a deep chemisorption well and lack early reaction barriers. Although CO binds strongly on Ru(0001), scattering is characterized by rather narrow angular distributions and sticking is incomplete even at low incidence energies. We perform molecular dynamics simulations, accounting for phononic (and electronic) energy loss channels, on a potential energy surface based on first-principles electronic structure calculations that reproduce the molecular beam experiments. We demonstrate that the mentioned unusual behavior is a consequence of a very strong rotational anisotropy in the molecule-surface interaction potential. Beyond the interpretation of scattering phenomena, we also discuss implications of our results for the recently proposed role of a precursor state for the desorption and scattering of CO from ruthenium.}, language = {en} } @article{AndrianovSaalfrank2006, author = {Andrianov, Ivan and Saalfrank, Peter}, title = {Free vibrational relaxation of H adsorbed on a Si(100) surface investigated with the multi-configurational time-dependent Hartree method}, series = {Chemical physics letters}, volume = {433}, journal = {Chemical physics letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0009-2614}, doi = {10.1016/j.cplett.2006.11.067}, pages = {91 -- 96}, year = {2006}, abstract = {The results of a quantum-mechanical study of vibrational relaxation of hydrogen adsorbed on a Si(100) surface with the multi-configurational time-dependent Hartree (MCTDH) method are presented. A two-dimensional subsystem is coupled non-linearly to a bath of harmonic oscillators (phonons of the Si bulk), and the relaxation of subsystem vibrations proceeds primarily via a two-phonon process. Characteristic times of the system evolution agree well with our previous perturbation theory study. The vibrational population decay is non-exponential, exhibiting pronounced recurrences due to finite bath size. The dependence of the lifetimes of the vibrational levels on the bath size and on the coupling details is investigated.}, language = {en} }