Florian Schiffmann

Ab initio molecular dynamics using hybrid density functionals

Ab initio molecular dynamics simulations with hybrid density functionals have so far found little application due to their computational cost. In this work, an implementation of the Hartree-Fock exchange is presented that is specifically targeted at ab initio molecular dynamics simulations of medium sized systems. We demonstrate that our implementation, which is available as part of the CP2K/Quickstep program, is robust and efficient. Several prescreening techniques lead to a linear scaling cost for integral evaluation and storage. Integral compression techniques allow for in-core calculations on systems containing several thousand basis functions. The massively parallel implementation respects integral symmetry and scales up to hundreds of CPUs using a dynamic load balancing scheme. A time-reversible multiple time step scheme, exploiting the difference in computational efficiency between hybrid and local functionals, brings further time savings. With extensive simulations of liquid water, we demonstrate the ability to perform, for several tens of picoseconds, ab initio molecular dynamics based on hybrid functionals of systems in the condensed phase containing a few thousand Gaussian basis functions.

An atomistic picture of the regeneration process in dye sensitized solar cells

A highly efficient mechanism for the regeneration of the cis-bis(isothiocyanato)bis(2,2′-bipyridyl-4,4′-dicarboxylato)-ruthenium(II) sensitizing dye (N3) by I- in acetonitrile has been identified by using molecular dynamics simulation based on density functional theory. Barrier–free complex formation of the oxidized dye with both I- and , and facile dissociation of and from the reduced dye are key steps in this process. In situ vibrational spectroscopy confirms the reversible binding of I2 to the thiocyanate group. Additionally, simulations of the electrolyte near the interface suggest that acetonitrile is able to cover the (101) surface of anatase with a passivating layer that inhibits direct contact of the redox mediator with the oxide, and that the solvent structure specifically enhances the concentration of I- at a distance which further favors rapid dye regeneration.

Synthesis of two-dimensional analogues of copolymers by site-to-site transmetalation of organometallic monolayer sheets

Monolayer sheets have gained attention due to the unique properties derived from their two-dimensional structure. One of the key challenges in sheet modification/synthesis is to exchange integral parts while keeping them intact. We describe site-to-site transmetalation of Zn(2+) in the netpoints of cm(2)-sized, metal-organic sheets by Fe(2+), Co(2+), and Pb(2+). This novel transformation was done both randomly and at predetermined patterns defined by photolithography to create monolayer sheets composed of different netpoints. All transmetalated sheets are mechanically strong enough to be spanned over 20 × 20 μm(2) sized holes. Density functional theory calculations provide both a model for the molecular structure of an Fe(2+)-based sheet and first insights into how transmetalation proceeds. Such transmetalated sheets with random and patterned netpoints can be considered as two-dimensional analogues of linear copolymers. Their nanoscale synthesis presents an advance in monolayer/polymer chemistry with applications in fields such as surface coating, molecular electronics, device fabrication, imaging, and sensing.

A Polar Corundum Oxide Displaying Weak Ferromagnetism at Room Temperature

Combining long-range magnetic order with polarity in the same structure is a prerequisite for the design of (magnetoelectric) multiferroic materials. There are now several demonstrated strategies to achieve this goal, but retaining magnetic order above room temperature remains a difficult target. Iron oxides in the +3 oxidation state have high magnetic ordering temperatures due to the size of the coupled moments. Here we prepare and characterize ScFeO3 (SFO), which under pressure and in strain-stabilized thin films adopts a polar variant of the corundum structure, one of the archetypal binary oxide structures. Polar corundum ScFeO3 has a weak ferromagnetic ground state below 356 K—this is in contrast to the purely antiferromagnetic ground state adopted by the well-studied ferroelectric BiFeO3.

Combining linear-scaling DFT with subsystem DFT in born-oppenheimer and ehrenfest molecular dynamics simulations: From molecules to a virus in solution

In this work, methods for the efficient simulation of large systems embedded in a molecular environment are presented. These methods combine linear-scaling (LS) Kohn-Sham (KS) density functional theory (DFT) with subsystem (SS) DFT. LS DFT is efficient for large subsystems, while SS DFT is linear scaling with a smaller prefactor for large sets of small molecules. The combination of SS and LS, which is an embedding approach, can result in a 10-fold speedup over a pure LS simulation for large systems in aqueous solution. In addition to a ground-state Born-Oppenheimer SS+LS implementation, a time-dependent density functional theory-based Ehrenfest molecular dynamics (EMD) using density matrix propagation is presented that allows for performing nonadiabatic dynamics. Density matrix-based EMD in the SS framework is naturally linear scaling and appears suitable to study the electronic dynamics of molecules in solution. In the LS framework, linear scaling results as long as the density matrix remains sparse during time propagation. However, we generally find a less than exponential decay of the density matrix after a sufficiently long EMD run, preventing LS EMD simulations with arbitrary accuracy. The methods are tested on various systems, including spectroscopy on dyes, the electronic structure of TiO2 nanoparticles, electronic transport in carbon nanotubes, and the satellite tobacco mosaic virus in explicit solution.