Refine
Year of publication
Document Type
- Article (49)
- Part of a Book (1)
- Doctoral Thesis (1)
- Other (1)
- Postprint (1)
Language
- English (53)
Is part of the Bibliography
- yes (53)
Keywords
- surface chemistry (2)
- Basis sets (1)
- Born-Oppenheimer MD (1)
- Chelates (1)
- Chemical calculations (1)
- Delta-Kohn-Sham (1)
- Electron dynamics (1)
- Ionization (1)
- Lasers (1)
- Ligand design (1)
Institute
- Institut für Chemie (53) (remove)
Fifteen N-butylpyridinium salts - five monometallic [C4Py](2)[MBr4] and ten bimetallic [C4Py](2)[(M0.5M0.5Br4)-M-a-Br-b] (M=Co, Cu, Mn, Ni, Zn) - were synthesized, and their structures and thermal and electrochemical properties were studied. All the compounds are ionic liquids (ILs) with melting points between 64 and 101 degrees C. Powder and single-crystal X-ray diffraction show that all ILs are isostructural. The electrochemical stability windows of the ILs are between 2 and 3 V. The conductivities at room temperature are between 10(-5) and 10(-6) S cm(-1). At elevated temperatures, the conductivities reach up to 10(-4) S cm(-1) at 70 degrees C. The structures and properties of the current bromide-based ILs were also compared with those of previous examples using chloride ligands, which illustrated differences and similarities between the two groups of ILs.
High Harmonic Generation (HHG) is a nonlinear optical process that provides a tunable source for high-energy photons and ultrashort laser pulses. Recent experiments demonstrated that HHG spectroscopy may also be used as an analytical tool to discriminate between randomly oriented configurational isomers of polyatomic organic molecules, namely, between the cis- and trans-forms of 1,2-dichloroethene (DCE) [M. C. H. Wong et al., Phys. Rev. A 84, 051403 (2011)]. Here, we suggest as an economic and at the same time a reasonably accurate method to compute HHG spectra for polyatomic species, Time-Dependent Configuration Interaction Singles (TD-CIS) theory in combination with extended atomic orbital bases and different models to account for ionization losses. The HHG spectra are computed for aligned and unaligned cis- and trans-DCE. For the unaligned case, a coherent averaging over possible rotational orientations is introduced. Furthermore, using TD-CIS, possible differences between the HHG spectra of cis- and trans-DCE are studied. For aligned molecules, spectral differences between cis and trans emerge, which can be related to their different point group symmetries. For unaligned, randomly oriented molecules, we also find distinct HHG spectra in partial agreement with experiment. In addition to HHG response in the frequency space, we compute time-frequency HHG spectra to gain insight into which harmonics are emitted at which time. Further differences between the two isomers emerge, suggesting time-frequency HHG as another tool to discriminate configurational isomers.
With recent experimental advances in laser-driven electron dynamics in polyatomic molecules, the need arises for their reliable theoretical modelling. Among efficient, yet fairly accurate methods for many-electron dynamics are Time-Dependent Configuration Interaction Singles (TD-CIS) (a Wave Function Theory (WFT) method), and Real-Time Time-Dependent Density Functional Theory (RT-TD-DFT), respectively. Here we compare TD-CIS combined with extended Atomic Orbital (AO) bases, TD-CIS/AO, with RT-TD-DFT in a grid representation of the Kohn-Sham orbitals, RT-TD-DFT/Grid. Possible ionization losses are treated by complex absorbing potentials in energy space (for TD-CIS/AO) or real space (for RT-TD-DFT), respectively. The comparison is made for two test cases: (i) state-to-state transitions using resonant lasers (pi-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H-2 molecule and cis- and trans-1,2-dichlorethene, C2H2Cl2, (DCE). From time-dependent electronic energies, dipole moments and from HHG spectra, the following observations are made: first, for bound state-to-state transitions enforced by pi-pulses, TD-CIS nicely accounts for the expected population inversion in contrast to RT-TD-DFT, in agreement with earlier findings. Secondly, when using laser pulses under non-resonant conditions, dipole moments and lower harmonics in HHG spectra are obtained by TD-CIS/AO which are in good agreement with those obtained with RT-TD-DFT/Grid. Deviations become larger for higher harmonics and at low laser intensities, i.e., for low-intensity HHG signals. We also carefully test effects of basis sets for TD-CIS/AO and grid size for RT-TD-DFT/Grid, different exchange-correlation functionals in RT-TD-DFT, and absorbing boundaries. Finally, for the present examples, TD-CIS/AO is observed to be at least an order of magnitude more computationally efficient than RT-TD-DFT/Grid.
Angle-resolved two-photon photoemission and high-resolution electron energy loss spectroscopy are employed to derive the electronic structure of a subnanometer atomically precise quasi-one-dimensional graphene nanoribbon (GNR) on Au(111). We resolved occupied and unoccupied electronic bands including their dispersion and determined the band gap, which possesses an unexpectedly large value of 5.1 eV. Supported by density functional theory calculations for the idealized infinite polymer and finite size oligomers, an unoccupied nondispersive electronic state with an energetic position in the middle of the band gap of the GNR could be identified. This state resides at both ends of the ribbon (end state) and is only found in the finite sized systems, i.e., the oligomers.
Near edge X-ray absorption fine structure (NEXAFS) simulations based on the conventional configuration interaction singles (CIS) lead to excitation energies, which are systematically blue shifted. Using a (restricted) open shell core hole reference instead of the Hartree Fock (HF) ground state orbitals improves (Decleva et al., Chem. Phys., 1992, 168, 51) excitation energies and the shape of the spectra significantly. In this work, we systematically vary the underlying SCF approaches, that is, based on HF or density functional theory, to identify best suited reference orbitals using a series of small test molecules. We compare the energies of the K edges and NEXAFS spectra to experimental data. The main improvement compared to conventional CIS, that is, using HF ground state orbitals, is due to the electrostatic influence of the core hole. Different SCF approaches, density functionals, or the use of fractional occupations lead only to comparably small changes. Furthermore, to account for bigger systems, we adapt the core-valence separation for our approach. We demonstrate that the good quality of the spectrum is not influenced by this approximation when used together with the non-separated ground state wave function. Simultaneously, the computational demands are reduced remarkably. (C) 2016 Wiley Periodicals, Inc.
Near edge X-ray absorption fine structure (NEXAFS) spectra and their pump-probe extension (PP-NEXAFS) offer insights into valence- and core-excited states. We present PSIXAS, a recent implementation for simulating NEXAFS and PP-NEXAFS spectra by means of the transition-potential and the Delta-Kohn-Sham method. The approach is implemented in form of a software plugin for the Psi4 code, which provides access to a wide selection of basis sets as well as density functionals. We briefly outline the theoretical foundation and the key aspects of the plugin. Then, we use the plugin to simulate PP-NEXAFS spectra of thymine, a system already investigated by others and us. It is found that larger, extended basis sets are needed to obtain more accurate absolute resonance positions. We further demonstrate that, in contrast to ordinary NEXAFS simulations, where the choice of the density functional plays a minor role for the shape of the spectrum, for PP-NEXAFS simulations the choice of the density functional is important. Especially hybrid functionals (which could not be used straightforwardly before to simulate PP-NEXAFS spectra) and their amount of "Hartree-Fock like" exact exchange affects relative resonance positions in the spectrum.
Tetrahalidocuprates(II) show a high degree of structural flexibility. We present the results of crystallographic and electron paramagnetic resonance (EPR) spectroscopic analyses of four new tetrabromidocuprate(II) compounds and compare the results with previously reported data. The cations in the new compounds are the sterically demanding benzyltriphenylphosphonium, methyltriphenylphosphonium, tetraphenylphosphonium, and hexadecyltrimethylammonium ions; they were used to achieve a reasonable separation of the paramagnetic Cu(II) ions for EPR spectroscopy. X-Ray crystallography shows that in all four complexes the [CuBr4](2-) units have a distorted tetrahedral coordination geometry which is in agreement with DFT calculations. The EPR hyperfine structure was not resolved. This is due to the exchange broadening resulting from still incomplete separation of the paramagnetic Cu(II) centres. Nevertheless, the principal values of the electron Zeemann tensor (g(parallel to) and g(perpendicular to)) of the complexes could be determined. A correlation of structural (X-ray) parameters with the spin density at the copper centres (DFT) is well reflected in the EPR spectra of the bromidocuprates. This enables the correlation of X-ray and EPR parameters to predict the structure of tetrabromidocuprates in physical states other than the crystalline state. As a result, we provide a method to structurally characterize [CuBr4](2-) in, for example, ionic liquids or in solution, which has important implications for e.g. catalysis or materials science.
We report quantum chemical calculations, mostly based on density functional theory, on azobenzene and substituted azobenzenes as neutral molecules or ions, in ground and excited states. Both the cis and trans configurations are computed as well as the activation energies to transform one isomer into the other and the possible reaction paths and reaction surfaces along the torsion and inversion modes. All calculations are done for the isolated species, but results are discussed in light of recent experiments aiming at the switching of surface mounted azobenzenes by scanning tunneling microscopes.
An electronic friction approach based on Langevin dynamics is used to describe the multidimensional (six-dimensional) dynamics of femtosecond laser induced desorption of H-2 and D-2 from a H(D)-covered Ru(0001) surface. The paper extends previous reduced-dimensional models, using a similar approach. In the present treatment forces and frictional coefficients are calculated from periodic density functional theory (DFT) and essentially parameter-free, while the action of femtosecond laser pulses on the metal surface is treated by using the two-temperature model. Our calculations shed light on the performance and validity of various adiabatic, non-adiabatic, and Arrhenius/Kramers type kinetic models to describe hot-electron mediated photoreactions at metal surfaces. The multidimensional frictional dynamics are able to reproduce and explain known experimental facts, such as strong isotope effects, scaling of properties with laser fluence, and non-equipartitioning of vibrational, rotational, and translational energies of desorbing species. Further, detailed predictions regarding translations are made, and the question for the controllability of photoreactions at surfaces with the help of vibrational preexcitation is addressed.
Stochastic approach to laser-induced ultrafast dynamics : the desorption of H-2/D-2 from Ru(0001)
(2010)
The desorption of molecular hydrogen and deuterium induced by femtosecond-laser pulses is studied theoretically for the so-called DIMET (Desorption Induced by Multiple Electronic Transitions) process. These investigations are based on nonadiabatic classical Monte Carlo trajectory (CMCT) simulations on a ground and an excited state potential energy surface, including up to all six adsorbate degrees of freedom. The focus is on the hot-electron mediated energy transfer from the surface to the molecule and back, and the energy partitioning between the different degrees of freedom of the desorbing molecules. We first validate for a two-mode model comprising the desorption mode and the internal vibrational coordinate, the classical Monte Carlo trajectory method by comparing with Monte Carlo wavepacket (MCWP) calculations arising from a fully quantum mechanical open-system density matrix treatment. We then proceed by extending the CMCT calculations to include all six nuclear degrees of freedom of the desorbing molecule. This allows for a detailed comparison between theory and experiment concerning isotope effects, energy partitioning (translational, vibrational, and rotational energies and their distributions), and the dependence of these properties on the laser fluence. The most important findings are as follows. (i) CMCT agrees qualitative with the MCWP scheme. (ii) The basic experimental features such as the large isotope effect, the non-linear increase of yield with laser fluence, translationally hot products (in the order of several 1000 K) and non-equipartitioning of translational and internal energies (E-trans > E- vib > E-rot) are well reproduced. (iii) Predictions concerning a strong angular dependence of translational energies at large observation angles are also made.
Selective excitation of molecule-surface vibrations in H2 and D2 dissociatively adsorbed on Ru(0001)
(2012)
In this contribution we report about the selective vibrational excitation of H2 and D2 on Ru(0001) as an example for nonadiabatic coupling of an open quantum system to a dissipative environment. We investigate the possibility of achieving state-selective vibrational excitations of H2 and D2 adsorbed on a Ru(0001) surface using picosecond infrared laser pulses. The systems behavior is explored using pulses that are rationally designed and others that are optimized using a time-local variant of Optimal Control Theory. The effects of dissipation on the laser-driven dynamics are studied using the reduced-density matrix formalism. The non-adiabatic couplings between adsorbate and surface are computed perturbatively, for which our recently introduced state-resolved anharmonic rate model is used. It is shown that mode- and state-selective excitation can be achieved in the absence of dissipation when using optimized laser pulses. The inclusion of dissipation in the model reduces the state selectivity and the population transfer yield to highly excited states. In this case, mode activation is most effectively realized by a rational pulse of carefully chosen duration rather than by a locally optimized pulse.
We investigate the recombinative desorption of hydrogen and deuterium from a Ru(0001) surface initiated by femtosecond laser pulses. We adopt a quantum mechanical two-state model including three molecular degrees of freedom to describe the dynamics within the desorption induced by electronic transition (DIET) limit. The energy distributions as well as the state-resolved and ensemble properties of the desorbed molecules are analyzed in detail by using the time-energy method. Our results shed light on the experimentally observed 1) large isotopic effects regarding desorption yields and translational energies and 2) the nonequal energy partitioning into internal and translational modes. In particular, it is shown that a single temperature is sufficient to characterize the energy distributions for all degrees of freedom. Further, we confirm that quantization effects play an important role in the determination of the energy partitioning.
We use quantum chemical cluster models together with constrained density STM Ph CI functional theory (DFT) and ab initio molecular dynamics (AIMD) for open system to simulate tip and rationalize nonlocal scanning tunneling microscope (STM) manipulation experiments for Philh ci chlorobenzene (PhCl) on a Si(111)-7 X 7 surface. We consider three different processes, namely, the electron-induced dissociation of the carbon-chlorine bond for physisorbed PhCl molecules at low temperatures and the electron- or hole-induced desorption of chemisorbed PhCl at 300 K. All processes can be induced nonlocally, i.e., up to several nanometers (nm) away from the injection site, in STM experiments. We rationalize and explain the experimental findings regarding the STM-induced dissociation using constrained DFT. The coupling of STM-induced ion resonances to nuclear degrees of freedom is simulated with AIMD using the Gadzuk averaging approach for open systems. From this data, we predict a 4 fs lifetime for the cationic resonance. For the anion model, desorption could not be observed. In addition, the same cluster models are used for transition-state theory calculations, which are compared to and validated against time-lapse STM experiments.
We report explicitly time-dependent coupled cluster singles doubles (TD-CCSD) calculations, which simulate the laser-driven correlated many-electron dynamics in molecular systems. Small molecules, i.e., HF, H(2)O, NH(3), and CH(4), are treated mostly with polarized valence double zeta basis sets. We determine the coupled cluster ground states by imaginary time propagation for these molecules. Excited state energies are obtained from the Fourier transform of the time-dependent dipole moment after an ultrashort, broadband laser excitation. The time-dependent expectation values are calculated from the complex cluster amplitudes using the corresponding configuration interaction singles doubles wave functions. Also resonant laser excitations of these excited states are simulated, in order to explore the limits for the numerical stability of our current TD-CCSD implementation, which uses time-independent molecular orbitals to form excited configurations.
In this paper we report time-dependent configuration interaction singles calculations modeling the laser- induced current through a metal-insulator-metal (MIM) contact. We compare our results to recent experiments [D. Diesing, M. Merschdorf, A. Thon, and W. Pfeiffer, Appl. Phys. B (to be published)]. We use two jellium slabs separated by a vacuum region in a one-dimensional model to describe the MIM contact. The contact is coupled to ultrashort (fs) laser pulses by the semiclassical dipole approximation. We discuss simulated two-pulse correlation spectra in comparison to experimental results
In this paper we report dynamical simulations of laser-driven, coupled nuclear-electron dynamics for a molecule- surface system. Specifically, the laser desorption of a small molecule (NO) from a metal slab (Pt) in the so-called DIET limit (Desorption Induced by Electronic Transitions), is studied. The excitation of the metal electrons by a laser pulse followed by the formation of a negative ion resonance, its subsequent decay, and the simultaneous desorption of the molecule are all treated within a single quantum mechanical model. This model is based on an earlier theory of Harris and others [S. M. Harris, S. Holloway, and G. R. Darling, J. Chem. Phys. 102, 8235 (1995)], according to which a nuclear degree of freedom is coupled to an electronic one, both propagated on a single non-Born-Oppenheimer potential energy surface. The goals of the present contribution are (i) to make a conceptual connection of this model to the frequently adopted nonadiabatic "multi-state" models of photodesorption, (ii) to understand details of the desorption mechanism, (iii) to explicitly account for the laser pulse, and (iv) to study the photodesorption as a function of the thickness of the metal film, and the laser parameters. As an important methodological aspect we also present a highly efficient numerical scheme to propagate the wave packet in a problem-adapted diabatic basis
Ultrafast electronic excitations of small sodium clusters and the onset of electron thermalization
(2009)
In this paper we report simulations of the ultrafast laser excitation and relaxation of the correlated valence electrons of a Na-8 cluster. The aim is twofold: first, while the total energy stays constant when the exciting laser pulse is over, we observe that the entropy computed from the reduced one electron density matrix rises on a much longer time scale. We discuss whether this can be understood as the onset of the thermalization of a finite system. Second, we describe this process with eight different methods of wavefunction-based electronic structure theory, which have been adapted for an explicitly time-dependent description. Their respective advantages and limitations for the simulation of the excitation and subsequent relaxation are explained.
In this paper, we report simulations of laser-driven many-electron dynamics by means of the time-dependent configuration interaction singles (TD-CIS) approach. Photoionization is included by a heuristic model within calculations employing standard Gaussian basis sets. Benzo[g]-N-methyl-quinolinium-7-hydroxylate (BMQ7H) serves as a test system to generate predefined wave packets, i.e. a superposition between the ground and fifth excited state, in a large molecule. For this molecule, these two states have a very similar geometry, which enables us to use the fixed nuclei approximation. Furthermore, this geometric stability would also prevent a dephasing of the electron wave packet due to nuclear dynamics in an experimental realization of our simulations. We also simulate the possible detection of such a wave packet by ultra short probe laser pulses, i.e. pump-probe spectra.