Counterion effects on the mesomorphic and electrochemical properties of guanidinium salts

Ionic liquid crystals (ILCs) combine the ion mobility of ionic liquids with the order and self-assembly of thermotropic mesophases. To understand the role of the anion in ILCs, wedge-shaped arylguanidinium salts with tetradecyloxy side chains were chosen as benchmark systems and their liquid crystalline self-assembly in the bulk phase as well as their electrochemical behavior in solution were studied depending on the anion. Differential scanning calorimetry (DSC), polarizing optical microscopy (POM) and X-ray diffraction (WAXS, SAXS) experiments revealed that for spherical anions, the phase width of the hexagonal columnar mesophase increased with the anion size, while for non-spherical anions, the trends were less clear cut. Depending on the anion, the ILCs showed different stability towards electrochemical oxidation and reduction with the most stable being the PF6 based compound. Cyclic voltammetry (CV) and density functional theory (DFT) calculations suggest a possible contribution of the guanidinium cation to the oxidation processes.

Optical response and multi-exciton effects in 2D PTCDA aggregates with local excitation

3,4,9,10-Perylenetetracarboxylic-dianhydride (PTCDA) aggregates have unique optical properties and are model materials for studying exciton energy transfer (EET) in planar stacked molecular aggregates. In the framework of density matrix theory, a hierarchy of molecular transition operator expectation values could be constructed to derive the equations of motion of multi-exciton states. Realistic parameters for PTCDA molecules are used to study EET and the optical response of two-dimensional aggregates upon local excitation. Our simulations show that information about the dark state can be obtained with local field excitation and the inter-chain coupling results in a red-shift of the lowest excitonic energy level. Configuration effects, inter-chain detuning and multi-exciton states are discussed. The calculated lowest excitonic energy level of a 2D PTCDA aggregate is qualitatively consistent with the lowest experimental absorption peak of a PTCDA film. The obtained results are valuable for the study of aggregates in optical nanocavities and for the design of photoelectric devices.

Approximation Schemes to Include Nuclear Motion in Laser-Driven Ab Initio Electron Dynamics: Application to High Harmonic Generation

The quantum mechanical description of many-electron dynamics in molecules driven by short laser pulses is at the heart of theoretical attochemistry. In addition to the formidable time-dependent electronic structure problem, the field faces the challenge that nuclear motion, ideally also treated quantum mechanically, may not be negligible, but scales enormously in effort. As a consequence, most first-principles calculations on ultrafast electron dynamics in molecules are done within the fixed-nuclear approximation. For laser-pulse excitation in H₂⁺, where an "exact" treatment of the coupled nuclear-electron dynamics is possible, it has been shown that nuclear motion can have a nonnegligible impact on high harmonic generation (HHG) spectra (Witzorky et al., J. Chem. Theor. Comput. 2021, 17, 7353-7365). It is not so clear, however, how to include (quantum) nuclear motion also for more complicated molecules, with more electrons and/or nuclei, in particular when the electronic structure is described by correlated, multistate wavefunction methods such as the time-dependent configuration interaction (TD-CI). In this work, we suggest a scheme in which the Born-Oppenheimer potential energy surfaces of a molecule are approximated by model potentials (harmonic and asymptotic, as an expansion in 1/R), obtained from only a few ab initio calculations, with the prospect to treat complex molecular systems. The method is tested successfully for HHG by few-cycle laser pulses for the "exact" H₂⁺ reference. It is then applied for diatomic molecules with more electrons and for a two-dimensional model of the water molecule using TD-CIS (S = single) for the electronic structure part.

Tetrahalidometallate(II) Ionic Liquids with More than One Metal: The Effect of Bromide versus Chloride

Fifteen N-butylpyridinium salts - five monometallic [C₄ Py]₂ [MBr₄ ] and ten bimetallic [C₄ Py]₂ [M_0.5 ^a M_0.5 ^b Br₄ ] (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 °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 °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.

Gaussian-Type Orbital Calculations for High Harmonic Generation in Vibrating Molecules: Benchmarks for H2. J Chem Theory Comput 2021;17:7353-7365. [PMID: 34747605 DOI: 10.1021/acs.jctc.1c00837]

The response of the hydrogen molecular ion, H₂⁺, 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 Schrödinger 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.

Many-electron dynamics in laser-driven molecules: wavefunction theory vs. density functional theory

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 (π-pulses), i.e., bound electron motion, and (ii) large-amplitude electron motion leading to High Harmonic Generation (HHG). Test systems are a H2 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 π-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.

PSIXAS: A Psi4 plugin for efficient simulations of X-ray absorption spectra based on the transition-potential and Δ-Kohn-Sham method

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 Δ-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.

Desorption induced by low energy charge carriers on Si(111)-7 × 7: First principles molecular dynamics for benzene derivates

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 × 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. © 2018 Wiley Periodicals, Inc.

Electron correlation dynamics in atoms and molecules

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.

STM-switching of organic molecules on semiconductor surfaces: an above threshold density matrix model for 1,5 cyclooctadiene on Si(100)

The scanning tunnelling microscope (STM)-induced switching of a single cyclooctadiene molecule between two stable conformations chemisorbed on a Si(100) surface is investigated using an above threshold model including a neutral ground state and an ionic excited state potential. Switching was recently achieved experimentally with an STM operated at cryogenic temperatures (Nacci et al 2008 Phys. Rev. B 77 121405(R)) and rationalized by a below threshold model using just a single potential energy surface (Nacci et al 2009 Nano Lett. 9 2997). In the present paper, we show that experimental key findings on the inelastic electron tunnelling (IET) switching can also be rationalized using an above threshold density matrix model, which includes, in addition to the neutral ground state potential, an anionic or cationic excited potential. We use one and two-dimensional potential energy surfaces. Furthermore, the influence of two key parameters of the density matrix description, namely the electronic lifetime of the ionic resonance and the vibrational lifetimes, on the ground state potential are discussed.

Concept of a single temperature for highly nonequilibrium laser-induced hydrogen desorption from a ruthenium surface

Laser-induced condensed phase reactions are often interpreted as nonequilibrium phenomena that go beyond conventional thermodynamics. Here, we show by Langevin dynamics and for the example of femtosecond-laser desorption of hydrogen from a ruthenium surface that light adsorbates thermalize rapidly due to ultrafast energy redistribution after laser excitation. Despite the complex reaction mechanism involving hot electrons in the surface region, all desorption product properties are characterized by equilibrium distributions associated with a single, unique temperature. This represents an example of ultrahot chemistry on the subpicosecond time scale.

Optical absorption and excitonic coupling in azobenzenes forming self-assembled monolayers: a study based on density functional theory

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 Å 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.

Microscopic model of the optical absorption of carbon nanotubes functionalized with molecular spiropyran photoswitches

The adsorption of molecules to the surface of carbon nanostructures opens a new field of hybrid systems with distinct and controllable properties. We present a microscopic study of the optical absorption in carbon nanotubes functionalized with molecular spiropyran photoswitches. The switching process induces a change in the dipole moment leading to a significant coupling to the charge carriers in the nanotube. As a result, the absorption spectra of functionalized tubes reveal a considerable redshift of transition energies depending on the switching state of the spiropyran molecule. Our results suggest that carbon nanotubes are excellent substrates for the optical readout of spiropyran-based molecular switches. The gained insights can be applied to other noncovalently functionalized one-dimensional nanostructures in an externally induced dipole field.

Stochastic approach to laser-induced ultrafast dynamics: the desorption of H(2)/D(2) from Ru(0001)

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.

Ab initio Electron Dynamics with the Multi-Configuration Time-Dependent Hartree-Fock Method

This paper gives an overview of the theory and some selected applications of the Multi-Configuration Time-Dependent Hartree-Fock (MCTDHF) method. This method can be seen as an explicitly time-dependent version of the Complete Active Space Self-Consistent Field (CASSCF) method of standard quantum chemistry. Strengths and shortcomings are discussed, and some comparisons with other methods of correlated electron dynamics are given.

Current versus temperature-induced switching in a single-molecule tunnel junction: 1,5 cyclooctadiene on Si(001)

The biconformational switching of single cyclooctadiene molecules chemisorbed on a Si(001) surface was explored by quantum chemical and quantum dynamical calculations and low-temperature scanning tunneling microscopy experiments. The calculations rationalize the experimentally observed switching driven by inelastic electron tunneling (IET) at 5 K. At higher temperatures, they predict a controllable crossover behavior between IET-driven and thermally activated switching, which is fully confirmed by experiment.

Stereoselective isomerization of an ensemble of adsorbed molecules with multiple orientations: stochastic laser pulse optimization for selective switching between achiral and chiral atropisomers

We present quantum dynamical simulations for the laser driven isomerization of an ensemble of surface mounted stereoisomers with multiple orientations. The model system 1-(2-cis-fluoroethenyl)-2-fluorobenzene supports two chiral and one achiral atropisomers upon torsion around the C-C single bond connecting phenyl ring and ethylene group. An infrared picosecond pulse is used to excite the internal rotation around the chiral axis, thereby controlling the chirality of the molecule. In order to selectively switch the molecules--independent of their orientation on a surface--from their achiral to either their left- or right-handed form, a stochastic pulse optimization algorithm is applied. The stochastic pulse optimization is performed for different sets of defined orientations of adsorbates corresponding to the rotational symmetry of the surface. The obtained nonlinearly polarized laser pulses are highly enantioselective for each orientation.

From stochastic pulse optimization to a stereoselective laser pulse sequence: simulation of a chiroptical molecular switch mounted on adamantane

Quantum dynamical simulations for the laser-controlled isomerization of 1-(2-cis-fluoroethenyl)-2-fluorobenzene mounted on adamantane are reported based on a one-dimensional electronic ground-state potential and dipole moment calculated by density functional theory. The model system 1-(2-cis-fluoroethenyl)-2-fluorobenzene supports two chiral and one achiral atropisomers upon torsion around the C-C single bond connecting the phenyl ring and ethylene group. The molecule itself is bound to an adamantyl frame which serves as a model for a linker or a surface. Due to the C3 symmetry of the adamantane molecule, the molecular switch can have three equivalent orientations. An infrared picosecond pulse is used to excite the internal rotation around the chiral axis, thereby controlling the chirality of the molecule. In order to selectively switch the molecules--independent of their orientations-- from their achiral to either their left- or right-handed form, a stochastic pulse optimization algorithm is applied. A subsequent detailed analysis of the optimal pulse allows for the design of a stereoselective laser pulse sequence of analytical form. The developed control scheme of elliptically polarized laser pulses is enantioselective and orientation-selective.

Molecular response properties from explicitly time-dependent configuration interaction methods

In this paper we report the calculation of molecular electric response properties with the help of explicitly time-dependent configuration interaction (TD-CI) methods. These methods have the advantage of being applicable (within the limitations of the time-dependent Schrodinger equation) to time-dependent perturbations of arbitrary shape and strength. Three variants are used to solve the time-dependent electronic Schrodinger equation, namely, the TD-CIS (inclusion of single excitations only), TD-CISD (inclusion of single and double excitations), and TD-CIS(D) (single excitations and perturbative treatment of double excitations) methods and applied for illustration to small molecules, H(2) and H(2)O. In the calculation, slowly varying off-resonant electric fields are applied to the molecules and linear (polarizabilities) and nonlinear (hyperpolarizabilities, harmonic generation) response properties are determined from the time-dependent dipole moments.