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.
Modern laser technology and ultrafast spectroscopies have pushed the timescales for detecting and manipulating dynamical processes in molecules from the picosecond over femtosecond domains, to the attosecond regime (1 as = 10(-18) s). This way, real-time dynamics of electrons after their photoexcitation can be probed and manipulated. In particular, experiments are moving more and more from atomic and solid state systems to molecules, opening the fields of "molecular electron dynamics" and "attosecond chemistry." Also on the theory side, powerful quantum dynamical tools have been developed to rationalize experiments on ultrafast electron dynamics in molecular species. <br /> In this contribution, we concentrate on theoretical aspects of ultrafast electron dynamics in molecules, mostly driven by lasers. The dynamics will be described with the help of wavefunction-based ab initio methods such as time-dependent configuration interaction (TD-CI) or the multiconfigurational time-dependent Hartree-Fock (MCTDHF) methods. Besides a survey of the methods and their extensions toward, e.g., treatment of ionization, laser pulse optimization, and open quantum systems, two specific examples of applications will be considered: The creation and/or dynamical fate of electronic wavepackets, and the nonlinear optical response to laser pulse excitation in molecules by high harmonic generation (HHG).
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.
The response of the hydrogen molecular ion, H-2(+), to few-cycle laser pulses of different intensities is simulated. To treat the coupled electron-nuclear motion, we use adiabatic potentials computed with Gaussian-type basis sets together with a heuristic ionization model for the electron and a grid representation for the nuclei. Using this mixed-basis approach, the time-dependent Schrodinger equation is solved, either within the Born-Oppenheimer approximation or with nonadiabatic couplings included. The dipole response spectra are compared to all-grid-based solutions for the three-body problem, which we take as a reference to benchmark the Gaussian-type basis set approaches. Also, calculations employing the fixed-nuclei approximation are performed, to quantify effects due to nuclear motion. For low intensities and small ionization probabilities, we get excellent agreement of the dynamics using Gaussian-type basis sets with the all-grid solutions. Our investigations suggest that high harmonic generation (HHG) and high-frequency response, in general, can be reliably modeled using Gaussian-type basis sets for the electrons for not too high harmonics. Further, nuclear motion destroys electronic coherences in the response spectra even on the time scale of about 30 fs and affects HHG intensities, which reflect the electron dynamics occurring on the attosecond time scale. For the present system, non-Born-Oppenheimer effects are small. The Gaussian-based, nonadiabatically coupled, time-dependent multisurface approach to treat quantum electron-nuclear motion beyond the non-Born-Oppenheimer approximation can be easily extended to approximate wavefunction methods, such as time-dependent configuration interaction singles (TD-CIS), for systems where no benchmarks are available.
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.
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.
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.
Quantum chemical approach to atomic manipulation of chlorobenzene on the Si(111)-7 x 7 surface
(2017)
We present a cluster model to describe the localization of hot charge carriers on the Si(111)-7 x 7 surface, which leads to (nonlocal) desorption of chlorobenzene molecules in scanning tunneling microscope (STM) manipulation experiments. The localized charge carriers are modeled by a small cluster. By means of quantum chemical calculations, this cluster model explains many experimental findings from STM manipulation. We show that the negative charge is mainly localized in the surface, while the positive one also resides on the molecule. Both resonances boost desorption: In the negative resonance the adatom is elevated; in the positive one the chemisorption bond between the silicon surface adatom and chlorobenzene is broken. We find normal modes promoting desorption matching experimental low-temperature activation energies for electron-and hole-induced desorption.
We use clusters for the modeling of local ion resonances caused by low energy charge carriers in STM-induced desorption of benzene derivates from Si(111)-7 x 7. We perform Born-Oppenheimer molecular dynamics for the charged systems assuming vertical transitions to the charged states at zero temperature, to rationalize the low temperature activation energies, which are found in experiment for chlorobenzene. Our calculations suggest very similar low temperature activation energies for toluene and benzene. For the cationic resonance transitions to physisorption are found even at 0 K, while the anion remains chemisorbed during the propagations. Further, we also extend our previous static quantum chemical investigations to toluene and benzene. In addition, an in depth analysis of the ionization potentials and electron affinities, which are used to estimate resonance energies, is given.
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.
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.
Non-Born-Oppenheimer quantum dynamics of H-2(+) and HD+ excited by single one-cycle laser pulses linearly polarized along the molecular (z) axis have been studied within a three-dimensional model, including the internuclear distance R and electron coordinates z and rho, by means of the numerical solution of the time-dependent Schrodinger equation on the timescale of about 200 fs. Laser carrier frequencies corresponding to the wavelengths of lambda(l) = 400 and 50 nm have been used and the amplitudes of the pulses have been chosen such that the energies of H-2(+) and HD+ are above the dissociation threshold after the ends of the laser pulses. It is shown that excitation of H-2(+) and HD+ above the dissociation threshold is accompanied by formation of vibrationally "hot" and "cold" ensembles of molecules. Dissociation of vibrationally "hot" molecules does not prevent the appearance of post-laser-pulse electronic oscillations, parallel z oscillations, and transversal rho oscillations. Moreover, dissociation of "hot" molecules does not influence characteristic frequencies of electronic z and rho oscillations. The main difference between the laser-induced quantum dynamics of homonuclear H-2(+) and its heteronuclear isotope HD+ is that fast post-laser-pulse electronic z oscillations in H-2(+) are regularly shaped with the period of tau(shp) approximate to 30 fs corresponding to nuclear oscillations in H-2(+), while electronic z oscillations in HD+ arise as "echo pulses" of its initial excitation and appear with the period of tau(echo) approximate to 80 fs corresponding to nuclear motion in HD+. Accordingly, corresponding power spectra of nuclear motion contain strong low-frequency harmonics at omega(shp) = 2 pi/tau(shp) in H2(+) and omega(echo) = 2 pi/tau(echo) in HD+. Power spectra related to both electronic and nuclear motion have been calculated in the acceleration form. Both higher- and lower-order harmonics are generated at the laser wavelength lambda(l) = 400 nm, while only lower-order harmonics are well pronounced at lambda(l) = 50 nm. It is also shown that a rationalized harmonic order, defined in terms of the frequency of the laser-induced electronic z oscillations, agrees with the concept of inversion symmetry for electronic motion in diatomic molecules.
We combine the stochastic pulse optimization (SPO) scheme with the time-dependent configuration interaction singles method in order to control the high frequency response of a simple molecular model system to a tailored femtosecond laser pulse. For this purpose, we use H-2 treated in the fixed nuclei approximation. The SPO scheme, as similar genetic algorithms, is especially suited to control highly non-linear processes, which we consider here in the context of high harmonic generation. Here, we will demonstrate that SPO can be used to realize a "non-harmonic" response of H2 to a laser pulse. Specifically, we will show how adding low intensity side frequencies to the dominant carrier frequency of the laser pulse and stochastically optimizing their contribution can create a high-frequency spectral signal of significant intensity, not harmonic to the carrier frequency. At the same time, it is possible to suppress the harmonic signals in the same spectral region, although the carrier frequency is kept dominant during the optimization. (C) 2016 AIP Publishing LLC.
We present and discuss the results of crystallographic and electron paramagnetic resonance (EPR) spectroscopic analyses of five tetrachloridocuprate(II) complexes to supply a useful tool for the structural characterisation of the [CuCl4](2-) moiety in the liquid state, for example in ionic liquids, or in solution. Bis(benzyltriethylammonium)-, bis(trimethylphenylammonium)-, bis(ethyltriphenylphosphonium)-, bis(benzyltriphenylphosphonium)-, and bis(tetraphenylarsonium) tetrachloridocuprate(II) were synthesised and characterised by elemental, IR, EPR and X-ray analyses. The results of the crystallographic analyses show distorted tetrahedral coordination geometry of all [CuCl4](2-) anions in the five complexes and prove that all investigated complexes are stabilised by hydrogen bonds of different intensities. Despite the use of sterically demanding ammonium, phosphonium and arsonium cations to obtain the separation of the paramagnetic Cu(II) centres for EPR spectroscopy no hyperfine structure was observed in the EPR spectra but the principal values of the electron Zeeman tensor, g(parallel to) and g(perpendicular to), could be determined. With these EPR data and the crystallographic parameters we were able to carry out a correlation study to anticipate the structural situation of tetrachloridocuprates in different physical states. This correlation is in good agreement with DFT calculations.
We present a novel laser pulse control for the chiroptical switch 1-(2-cis-fluoroethenyl)-2-fluoro-3,5-dibromobenzene mounted on adamantane, where the latter imitates a linker group or part of a solid surface. This molecular device offers three switching states: a true achiral "off"-state and two chiral "on"-states of opposite handedness. Due to the alignment of its chiral axis along the surface normal several defined orientations of the switch have to be considered for an efficient stereocontrol strategy. In addition to these different initial conditions, coupled torsional degrees of freedom around the chiral axis make the quest for highly stereoselective laser pulses a challenge. The necessary flexibility in pulse accomplished by employing the iterative stochastic pulse optimization method we presented recently. Still, the complexity of the system dictates a combined treatment by fast molecular dynamics and computationally intensive quantum dynamics. Although quantum effects are found to be of importance, the pulses optimized within the classical treatment allow us to turn on the chirality of the switch, achieving high enantioselectivity in the quantum treatment for all orientations at the same time.
Motivated by recent atomic manipulation experiments, we report quantum chemical calculations for chemi- and physisorption minima of chlorobenzene on the Si(111)-7x7 surface. A density functional theory cluster approach is applied, using the B3LYP hybrid functional alongside Grimme's empirical dispersion corrections (D3). We were able to identify chemisorption sites of binding energies of 1.6 eV and physisorption energies of 0.6 eV, both in encouraging agreement with the trend of experimental data. The cluster approach opens up the possibility of a first-principles based dynamical description of STM manipulation experiments on this system, the interpretation of which involves both the chemi- and physisorbed states. However, we found that special care has to be taken regarding the choice of clusters, basis sets, and the evaluation of the dispersion corrections.
Based on the analysis of optical absorption spectra, it has recently been speculated that the excitonic coupling between individual azobenzene-functionalized alkanethiols arranged in a self-assembled monolayer (SAM) on a gold surface could be strong enough to hinder collective trans-cis isomerization-on top of steric hindrance [Gahl et al., J. Am. Chem. Soc., 2010, 132, 1831]. Using models of SAMs of increasing complexity (dimer, linear N-mers, and two-dimensionally arranged N-mers) and density functional theory on the (TD-) B3LYP/6-31G* level, we determine optical absorption spectra, the nature and magnitude of excitonic couplings, and the corresponding spectral shifts. It is found that at inter-monomer distances of about 20 angstrom and above, TD-B3LYP excitation frequencies (and signal intensities) can be well described by the frequently used point-dipole approximation. Further, calculated blue shifts in optical absorption spectra account for the experimental observations made for azobenzene/gold SAMs, and hint to the fact that they can indeed be responsible for reduced switching probability in densely packed self-assembled structures.
The formation of CuCl nanoplatelets from the ionic liquid precursor (ILP) butylpyridinium tetrachlorocuprate [C4Py](2)[CuCl4] using ascorbic acid as a reducing agent was investigated. In particular, electron paramagnetic resonance (EPR) spectroscopy was used to evaluate the interaction between ascorbic acid and the Cu(II) ion before reduction to Cu(I). EPR spectroscopy suggests that the [CuCl4](2-) ion in the neat IL is a distorted tetrahedron, consistent with DFT calculations. Addition of ascorbic acid leads to the removal of one chloride from the [CuCl4](2-) anion, as shown by DFT and the loss of symmetry by EPR. DFT furthermore suggests that the most stable adduct is formed when only one hydroxyl group of the ascorbic acid coordinates to the Cu(II) ion.
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.