Q-Chem 5.0 with ADC/SS-PCM has been released!

While I'm still in the process of writing the second post about boron-subtituted aromatics, let me quickly mention that Q-Chem 5.0 has recently been released. 

Why is this special, you may ask? It is, because it contains the ADC/SS-PCM approach that I have developed during my PhD and that has been discussed numerous times in this blog. With this approach, you can accurately model your favourite excited states and transitions from (emission) and to (absorption) them in solution at the ADC level of theory at up to third order in perturbation theory, i.e., ADC(3), a very accurate benchmark method.

You can also employ e.g. the very efficient riSOS-ADC(2) approach and investigate quite large molecules (say: materials) with up to 500 basis functions. Since ADC is, in contrast to TD-DFT, an ab-initio method, it does provide complete and physically sound description of the electrostatics of the systems as well as its interaction with the solvent/dielectric environment. We compared the ADC2/SS-PCM and ADC3/SS-PCM approaches already in the publication presenting the method.

I'm currently in the final steps of publishing another article about a project with my friend and colleague Felix (find his blog in the menu), which turned out to be a nice showcase for the model. Turns out I just made this table-of-contents graphic, which I want to share with you:

It shows the excitation-hole (in blue) and excited electron (in red) of the lowest excited singlet state (S1) of an already quite large push-push system that we studied in the article. In the top left corner, it shows the excited electron computed with the SS-PCM simulating the non-polar solvent cyclohexane, and in the bottom right corner with parameters for the polar solvent acetonitrile. From this quite intuitive visualization it becomes clear how the electron-hole or in other word charge separation (the distance between the blue and red blops) increases as it is stabilized by the polar solvent and, more importantly, how this affects the properties of the system. I'll post a link to the article as soon as it appears online. (*Link)

A brief description of the theory behind the ADC/SS-PCM model and its capabilities can be found in the Q-Chem 5.0 online manual, which you can find here.

Please note that very unfortunately, due to a last minute change of some defaults of the PCM solvent model that was not communicated very well, the description of ADC/SS-PCM in the in the manual is incomplete concerning one detail: For all calculations with the model the line "ChargeSeparation Marcus" has to be included in the $pcm block of the input file. This bug will be fixed with the next release (5.0.1) in July.

So long,
Jan

A story about boron-substituted aromatics - part one

Another project I have been working on last year was concerned with the quantum-chemical investigation of a bunch of boron-containing aromatics that have been synthesized and studied in the group of Prof. Wagner in Frankfurt. This project is interesting for a number of reasons, the most story-worthy of which is that Prof. Wagners lecture "Allgemeine und Anorganische Chemie" marked the very beginning of my academic career in October 2005. This cooperation was a nice opportunity to "turn the spit" and explain aspects of photochemistry to the Professor who introduced me to chemistry as a scientific discipline some 10 years ago. But also from a scientific point of view this project is exciting. These boron compounds do in general have a number of quite unique photophysical properties (some e.g. exhibit temperature-dependent delayed fluorescence with a curios solvent-dependency), and the project turned out to provide the opportunities to try out a couple of methods I had wanted to try out for a while. The paper itself can be found here. In the following I want to shed some light onto aspects that didn't get too much attention in the brevity of a scientific article.

Two of the four investigated molecules can be derived from anthracene and pentacene by replacing the middle carbons with boron atoms and are named accordingly: Dibora-anthracene (DBA) and dibora-pentacene (DBP). Diborinine (DBI) consists of two antiaromatic biphenylene subunits bridged via the dibora-benzylic motif, while iso-DBP is an asymetric variant of DBP, in which the boron-containing ring is shifted to one side.

The initial motivation for the project was to find an explanation for the surprisingly different fluorescence quantum yields of these closely related molecules. While isoDBP and DBP exhibit high quantum yields of 90% and 50%, respectively, DBA is hardly fluorescent (2%) and DBI completely dark. But let us postpone this question to one of the upcoming posts about this project and start  off with another question that popped up while I was inspecting the absorption spectra of these compounds in various solvents: DBI and to some extent also DBA exhibit a huge influence of some solvents onto the shape of the absorption spectra. Specifically, the main peaks of the absorption spectrum of DBI in non-polar cyclohexane (CHX) are blue-shifted and almost vanish in acetonitrile (ACN), while this is not the case in other polar solvents. We concluded that there has to exist an specific interaction between the lewis-basic solvent ACN and the lewis-acidic boron-atoms in DBI, or in other words that ACN behaves as a non-innocent solvent.

As a first step to confirm this hypothesis, I computed the structures and electronic energies of one mono- and two di-adducts (cis and trans) with a cheap and fast (exploratory) methodology (B3LYP-D3BJ/SVP, today I would use PBEh-3c). The resulting structures are shown together with the structure of bare DBI in the following picture:

Calculated ground-state minimum structures of bare DBI (a), the ACN-monoadduct (b), the trans-diadduct (c) and the cis-diadduct (d). Already from these cheap, approximate calculations it was evident that only the mono- and trans-diadduct are energetically feasible.

For the energetically feasible adducts of DBI and DBA, I subsequently calculated vibrationally resolved absorption spectra, which can be directly compared to the experimental spectra. For this purpose, I used the same methodology, namely TD-B3LYP with a small SVP basis, since it is fast and due to a surprisingly stable inherent error compensation typically also quite accurate (but also not reliable).

The calculation of vibrationally resolved spectra is one of the things I wanted to try out for a while now and hence, I'll briefly explain the hows and whys in the following. To make this type of calculation possible for such a large molecule, one has to employ a drastic-sounding approximation with a very long name: Independent-mode-shifted-harmonic-oscillator (IMDHO)-approximation. This is necessary to reduce the untraceable problem of knowing the shape of all potential-energy surfaces of all excited-states of interest around the initial (ground-state) geometry to something actually doable. Using IMDHO, the problem at hand is reduced to knowing the vibrational frequencies (normal modes) of the ground-state and the gradients along them in each of the excited states of interest. The assumptions required to get there are:

1. All modes can be approximated as a harmonic oscillator. (Common and not too bad for the ground state, crude for excited states)
2. The vibrational frequencies in the excited states are so similar to those of the ground state than we can treat them as being identical (Very crude)

From these assumptions, it follows that the only difference between the ground- and excited-state PES is the position of the respective minima, which is why its called "shifted"HO. Since we are in the harmonic approximation, all we need to know this shift w.r.t. to ground-state minimum are the gradients along these normal modes in the excited states of interest. Hence, to obtain the vibrationally resolved absorption spectrum within the IMDHO approximation, one first computes the ground state vibrational frequencies (normal-modes), and subsequently the gradient along these normal-modes for each of the excited-states. Having these quantities, we (actually the handy tool ORCA_ASA does that for us) can build the approximate PES and compute the coupling between the electronic and vibrational excitations, yielding the spectra, which are shown together with the absorption spectra recorded in various mixtures of CHX and ACN in the following picture:

Comparison of measured (top) and calculated (bottom) absorption spectra of DBI and DBA and their respective ACN mono- and diadducts. Calculations of the vibrationally resolved spectra were carried out at the (TD)B3LYP/SVP level of theory employing the so-called independent-mode-shifted-harmonic-oscillator (IMDHO)-approximation. Considering the crude IMDHO approximation and the poor level of electronic structure theory, the agreement between experiment and calculation is surprisingly good.

Due to the very reasonable agreement between experiment and calculation, in particular for DBI, we can doubtlessly identify the ACN-monoadduct in the spetra recorded in ACN. However, it is difficult to say this with certainty since the trans-diadduct is essentially dark (no absorption) and would not show up in the spectra. Thus, to further investigate the issue I conducted a thorough thermochemical analysis of the formation of the ACN adducts for DBA, DBI and DBP, which I will present and discuss in the upcoming post.

So long!
Jan