Refine
Year of publication
Document Type
- Article (61)
- Postprint (2)
- Part of a Book (1)
- Doctoral Thesis (1)
- Other (1)
Language
- English (66)
Is part of the Bibliography
- yes (66) (remove)
Keywords
- DFT (2)
- STM (2)
- charge localisation (2)
- cluster model (2)
- surface chemistry (2)
- Basis sets (1)
- Born-Oppenheimer MD (1)
- Chelates (1)
- Chemical calculations (1)
- Delta-Kohn-Sham (1)
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.
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.
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.
Hot localised charge carriers on the Si(111)-7×7 surface are modelled by small charged clusters. Such resonances induce non-local desorption, i.e. more than 10 nm away from the injection site, of chlorobenzene in scanning tunnelling microscope experiments. We used such a cluster model to characterise resonance localisation and vibrational activation for positive and negative resonances recently. In this work, we investigate to which extent the model depends on details of the used cluster or quantum chemistry methods and try to identify the smallest possible cluster suitable for a description of the neutral surface and the ion resonances. Furthermore, a detailed analysis for different chemisorption orientations is performed. While some properties, as estimates of the resonance energy or absolute values for atomic changes, show such a dependency, the main findings are very robust with respect to changes in the model and/or the chemisorption geometry.
Utilizing suitable precursor molecules, a thermally activated and surface-assisted synthesis results in the formation of defect-free graphene nanoribbons (GNRs), which exhibit electronic properties that are not present in extended graphene. Most importantly, they have a band gap in the order of a few electron volts, depending on the nanoribbon width. In this study, we investigate the electronic structure changes during the formation of GNRs, nitrogen-doped (singly and doubly N-doped) as well as non-N-doped chevron-shaped CGNRs on Au(111). Thus we determine the optical gaps of the precursor molecules, the intermediate nonaromatic polymers, and finally the aromatic GNRs, using high-resolution electron energy loss spectroscopy and density functional theory calculations. As expected, we find no influence of N-doping on the size of the optical gaps. The gap of the precursor molecules is around 4.5 eV. Polymerization leads to a reduction of the gap to a value of 3.2 eV due to elongation and thus enhanced delocalization. The CGNRs exhibit a band gap of 2.8 eV, thus the gap is further reduced in the nanoribbons, since they exhibit an extended delocalized pi-electron system.
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.
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.
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.
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.
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.
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∥ and g⊥, 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 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.
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.
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.
Fluoroionophores of fluorophore-spacer-receptor format were prepared for detection of PdCl2 by fluorescence enhancement. The fluorophore probes 1-13 consist of a fluorophore group, in alkyl spacer and a dithiomaleonitrile PdCl2 receptor. First, varying the length of the alkylene spacer (compounds 1-3) revealed, dominant through-space pathway for oxidative photoinduced electron transfer (PET) in CH2-bridged dithiomaleonitrile fluoroionophores. Second. fluorescent probes 4-9 containing two anthracene or pyrene fragments connected through CH2 bridges to the dithiomaleonitrile unit were synthesized. Modulation of the oxidation potential (E-Ox) through electron-withdrawing or -donating groups on the anthracene moiety regulates file thermodynamic driving force for oxidative PET (Delta G(PET)) in bis(anthrylmethylthio)maleonitriles and therefore the fluorescence quantum yields (Phi(f)), too. The new concept was confirmed and transferred to pyrenyl ligands, and fluorescence enhancements (FE) greater than 3.2 in the presence of PdCl2 were achieved by 7 and 8 (FE=5.4 and 5.2). Finally, for comparison, monofluorophore ligands 10-13 were synthesized.
Fluoroionophores of fluorophore-spacer-receptor format were prepared for detection of PdCl2 by fluorescence enhancement. The fluorescent probes 1-13 consist of a fluorophore group, an alkyl spacer and a dithiomaleonitrile PdCl2 receptor. First, varying the length of the alkylene spacer (compounds 1-3) revealed a dominant through-space pathway for oxidative photoinduced electron transfer (PET) in CH2-bridged dithiomaleonitrile fluoroionophores. Second, fluorescent probes 4-9 containing two anthracene or pyrene fragments connected through CH2 bridges to the dithiomaleonitrile unit were synthesized. Modulation of the oxidation potential (EOx) through electron-withdrawing or -donating groups on the anthracene moiety regulates the thermodynamic driving force for oxidative PET (GPET) in bis(anthrylmethylthio)maleonitriles and therefore the fluorescence quantum yields (f), too. The new concept was confirmed and transferred to pyrenyl ligands, and fluorescence enhancements (FE) greater than 3.2 in the presence of PdCl2 were achieved by 7 and 8 (FE=5.4 and 5.2). Finally, for comparison, monofluorophore ligands 10-13 were synthesized.
In this paper, we present quantum dynamical calculations on electron correlation dynamics in atoms and molecules using explicitly time-dependent ab initio configuration interaction theory. The goals are (i) to show that in which cases it is possible to switch off the electronic correlation by ultrashort laser pulses, and (ii) to understand the temporal evolution and the time scale on which it reappears. We characterize the appearance of correlation through electron-electron scattering when starting from an uncorrelated state, and we identify pathways for the preparation of a Hartree-Fock state from the correlated, true ground state. Exemplary results for noble gases, alkaline earth elements, and selected molecules are provided. For Mg we show that the uncorrelated state can be prepared using a shaped ultrashort laser pulse.
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.
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.
In this paper, we perform many-electron dynamics using the time-dependent configuration-interaction method in its reduced density matrix formulation (rho-TDCI). Dissipation is treated implicitly using the Lindblad formalism. To include the effect of ionization on the state-resolved dynamics, we extend a recently introduced heuristic model for ionizing states to the rho-TDCI method, which leads to a reduced density matrix evolution that is not norm-preserving. We apply the new method to the laser-driven excitation of H(2) in a strongly dissipative environment, for which the state-resolve lifetimes are tuned to a few femtoseconds, typical for dynamics of adsorbate at metallic surfaces. Further testing is made on the laser-induced intramolecular charge transfer in a quinone derivative as a model for a molecular switch. A modified scheme to treat ionizing states is proposed to reduce the computational burden associated with the density matrix propagation, and it is thoroughly tested and compared to the results obtained with the former model. The new approach scales favorably (similar to N(2)) with the number of configurations N used to represent the reduced density matrix in the rho-TDCI method, as compared to a N(3) scaling for the model in its original form.