Dear Reader

After years of misusing this blog for a website, I finally made a real website for myself. You can find it here (www.janmewes.de) since I moved the redirect today.

Have some nice and healthy holidays and a happy new year!

Jan

once a blog about workaday quantum chemistry

Dear Reader

After years of misusing this blog for a website, I finally made a real website for myself. You can find it here (www.janmewes.de) since I moved the redirect today.

Have some nice and healthy holidays and a happy new year!

Jan

Schmitz, Schwerdtfeger, Grimme, **Mewes**, *JACS, ***2021****/22** (just accepted)

“It's Complicated: On Relativistic Effects and Periodic Trends in the Melting and Boiling Points of the Group 11 Coinage Metals”

“It's Complicated: On Relativistic Effects and Periodic Trends in the Melting and Boiling Points of the Group 11 Coinage Metals”

Kunze, Hansen, Grimme, **Mewes**, *JPCL, ***2021*** *

“PCM-ROKS for the Description of Charge-Transfer States in Solution: Singlet-Triplet Gaps with Chemical Accuracy from Open-Shell Self-Consistent Reaction-Field Calculations”

“PCM-ROKS for the Description of Charge-Transfer States in Solution: Singlet-Triplet Gaps with Chemical Accuracy from Open-Shell Self-Consistent Reaction-Field Calculations”

“Comment on ‘The Nature of Chalcogen‐Bonding‐Type Tellurium–Nitrogen Interactions’: Fixing the Description of Finite-Temperature Effects Restores the Agreement Between Experiment and Theory”

Exclusively Relativistic: Periodic Trends in the Melting and Boiling Points of Group 12

Grimme, Hansen, Ehlert, **Mewes**, *JCP*, **2021**

r²SCAN-3c: A "Swiss Army Knife" Composite Electronic-Structure Method

r²SCAN-3c: A "Swiss Army Knife" Composite Electronic-Structure Method

Accurate Elemental Boiling Points From First Principles

Scholz, Massoth, Bursch,

Smits**, Mewes**, Jerabek, Schwerdtfeger, *Angewandte, ***2020**:

Caldeweyher, **Mewes**, Ehlert, Grimme, *PCCP, ***2020**:

Extension and evaluation of the D4 London-dispersion correction for periodic solids

*(HOT Article)*

**Mewes**, Smits, Kresse, Schwerdtfeger, *Angewandte, ***2019**:

Copernicium: A Relativistic Noble Liquid

*(Open Access)*

**Mewes**, Smits, Jerabek, Schwerdtfeger, *Angewandte, ***2019:**

Oganesson is a Semiconductor: On the Relativistic Band-Gap Narrowing in the Heaviest Noble-Gas Solids

*(Hot Paper, Open Access, auf Deutsch)*

Trombach, Ehlert, Grimme, Schwerdtfeger,**Mewes**, *PCCP**, ***2019:**

Exploring the chemical nature of super-heavy main-group elements by means of efficient plane-wave density-functional theory

*(HOT Article, Frontcover)*

**Mewes**, *PCCP**, ***2018: **

Modeling TADF in organic emitters requires careful consideration of the environment and going beyond the Franck–Condon approximation*(HOT Article) *

__Professional:__* *

*Jan 2021 – current***Scientific Consultant, **Lanthanoid-based OLEDs

beeOLED GmbH,**Dresden**

*Jun 2019 – current*
**Postdoc/Senior Researcher** (one-year returning Scholarship, AvH)

with Prof.**Stefan** **Grimme**,

Mulliken-Center for Theoretical Chemistry,**Bonn**

* Sep 2016 – May 2019 *
**Feodor-Lynen Postdoc **(Scholarship of the Alexander-von-Humboldt Foundation)

with Prof.**Peter** **Schwerdtfeger**,

New Zealand Institute for Advanced Study (NZIAS),**Auckland**

* Mar 2015 – Aug 2016 *
**Postdoc** (Scholarship of HGS Mathcomp)

with Prof.**Andreas** **Dreuw**

Interdisciplinary Center for Scientific Computing, **Heidelberg**

__Education:__

Extension and evaluation of the D4 London-dispersion correction for periodic solids

Copernicium: A Relativistic Noble Liquid

Oganesson is a Semiconductor: On the Relativistic Band-Gap Narrowing in the Heaviest Noble-Gas Solids

Trombach, Ehlert, Grimme, Schwerdtfeger,

Exploring the chemical nature of super-heavy main-group elements by means of efficient plane-wave density-functional theory

beeOLED GmbH,

with Prof.

Mulliken-Center for Theoretical Chemistry,

with Prof.

New Zealand Institute for Advanced Study (NZIAS),

with Prof.

Oct 2010 – Feb 2015

with Prof.

Interdisciplinary Center for Scientific Computing (IWR), Heidelberg,

Fellow of the HGS MathComp

Thesis: „Development and Application of Methods for the Description of Photochemical Processes in Condensed Phase“.

with Prof.

Goethe-University, Frankfurt

Thesis: „The Mechanism of Photodecarboxylation of Nitrophenylacetate and the Implications for Ortho-Nitrobenzyl Caged Compounds“.

with Prof.

Goethe-University, Frankfurt

Thesis (German): „Energietransfer in Bakteriellen Lichtsammelkomplexen“.

After almost three years without a new post, I have to admit that as an early-career scientist and father of two (1 and 4), the time for running a research blog is very limited. Thus, I will turn this website into a somewhat more static source of information.

so long,

Jan

so long,

Jan

This post is a follow up on my last post, in which I discussed the accuracy of the ADC/SS-PCM approach for a number of push-pull systems. So if you haven't read part one yet, please do so before you continue.

Part one ended with the finding that solvent effects were significantly over-estimated for the type of molecule under investigation. More specifically, the shift of the fluorescence energy between non-polar cyclohexane (chx) and polar acetonitrile (acn) solution calculated with ADC(2)/SS-PCM was exactly twice as large as the experimental numbers - for almost all of the systems. Even more surprisingly, the theoretically ill-defined one-shot approach (quick reminder: one-shot means that I use the result from the first step of the solvent-field iterations) seemingly provides a much better agreement with the experiment than the converged solvent-field:

I ended that first part with the promise to explain why I'm sure that this is a coincidence, what/how we can learn from it, and what underlying cause is responsible for the problem. So let's go:

The fact that the error is so suspiciously constant just means that it is very systematic. Systematic errors are much nicer than statistic ones, because its very straightforward to correct them. In this particular case, the correction is simply dividing the shifts by two. Inspecting the plot of the results above, you will find that the results from the converged solvent field divided by two (green line) reproduce the experiment even better than the one-shot results. Or in other words: Although the one-shot approach is much closer to the experiment in absolute numbers, the converged results follow the experimental trend more closely, i.e., their statistical error is smaller. The good news is that after all, this shows that the solvent-field iterations indeed improve the results compared to the one-shot approach. This was just difficult to see because of a large systematic error. Eventually, we are left with the remaining question:

As I already mentioned on a number of occasions, ADC(2) seems to over-stabilize states with large density shifts, such as CT states. I think this due to too large orbital-relaxation effects. Orbital relaxation is the response of the molecular density to the primary electron transfer during an excitation. While this is quite difficult to explain in words, it is quite apparent from a comparison of the the electron and hole densities of an excitation (initial electron transfer) to the total attachment and detachment densities (including orbital relaxation). Here I have calculated and visualized all of these for ZMSO2-14:

While the electron (top left) and hole (top right) densities appear to be rather independent (positions and sizes of the visible blobs are uncorrelated), they show that in this excited state an electron from the nitrogen lone-pair is excited into a pi* orbital. The attachment (bottom left) and detachment (bottom right) densities are, in contrast to that, clearly not independent, but correlated. Wherever there is a positive contribution in the hole/attachment density (e.g. the extra electron in the pi-system in the left part of the molecule), one can find negative blue (orbital-relaxation) contributions at the same place in the sigma system in the detachment density. Vice versa, wherever density vanishes, e.g. from the nitrogen lone-pair, there are positive relaxation contributions in the attachment density, which are clearly not there in the primary electron density. This response of the density to the initially excited electron/excitation-hole is exactly what I understand to be orbital relaxation. To have an even more intuitive picture, I should probably go ahead and plot just the difference between the electron/hole and attachment/detachment densities, but for now the pictures above must suffice.

However, while these pictures are nice to illustrate the effect, one needs a more quantitative measure to see how strong orbital relaxation is in these molecules, and to compare it to other molecules. For this purpose, its advisable to look at the so-called promotion numbers, which can be seen as integrals over the attachment/detachment densities, or in other words the number of electrons that is involved in/shifted around during the excitation. A typical locally excited state would have a promotion number of around 1.5, meaning that in addition to the initially excited electron, a total of 0.5 electrons make up the orbital relaxation. For charge-transfer states which do in general have stronger orbital relaxation (because more charge is shifted around), these numbers are closer to or maybe even larger than two. Let us now take a look how the promotion numbers of ZMS-14 and ZMSO2-14 evolve during the solvent-field iterations in chx, diethylether (eth), dichloromethane (dcm) and acn (please focus only on the reddish and blueish lines for now):

Apparently, there is a significant increase of the promotion numbers and thus the amount of orbital relaxation during the solvent-field iterations in polar solvents, even for the weaker push-pull systems. The promotion numbers increase to unphysically large values of 2.5 and larger for ZMSO2-14 in polar solvents and reach 2.2 even in chx. In general, I was surprised that the promotion number change at all during the iterations, which is something I have never observed before. I have investigated a number of large push-pull systems designed for TADF with the very same methodology and small systems such as DMABN with ADC2 and ADC3/SS-PCM, for which I obtained much better agreement for the calculated chx to acn shifts. Also, in all of these examples, the promotion numbers hardly ever change during the iterations, as evident from the the representative example termed F1 shown in the figure above as greenish lines (the molecule is taken form this article).

A similar problem with unphysically large promotion numbers accommodated by faulty excitation energies at the ADC(2) level of theory was reported by my friend and colleague Felix in an investigation of iridium complexes. In his case, the problem appeared in the gas-phase calculation and was corrected at the ADC(3) level of theory.

After all, I tend to think this is an ADC(2) problem that is amplified/unveiled during the solvent-field iterations. Consequently, one should carefully check the evolution of the promotion numbers during the solvent field iterations (they can be obtained by activating "state_analysis = true" and "nto_pairs = 2" in the $rem block of the input file). If there are large changes, one should be extra careful with the results of the method.

Dear colleagues,

in the framework of a recent project published **in this article**,
I for the first time had the chance to systematically study the
performance of ADC/SS-PCM for a group of click-chemistry generated
molecules with varying charge-transfer character in their emitting
excited state.

The common motif of these molecule is that they consist of an electron-pushing triphenylamin group (right side of the molecule in the figure below) and an electron-pulling group (left side), whose strength can be modulated by oxidizing the contained sulphur atom. In the lowest excited singlet state of all of these systems, an electron is excited from the electron-pushing triphenylamine to the electron pulling part of the molecule. While the excitation-hole turned out to be essentially identical in all of these systems and all solvents, the location and structure of the excited electron differed significantly depending on the molecule and environment.

The figure below shows the structures of the molecules as well as the location of the excited electron computed at the ADC2/SS-PCM/SV level of theory for the non-polar solvent cyclohexane on the left and for the polar solvent acetonitrile in the right.

Since the discussion of the performance of the methodology was (as usual) quite brief in the article, I will share some more methodological insights here and in the following posts.

The central observation, which also made it into the paper, was that the solvent effects are systematically overestimated for this group of molecules, leading to too low emission energies in polar solvents (here acetonitrile, ACN). In the non-polar solvent cyclohexane (CHX) the error seems to cancel out with the systematic overestimation of the emission energies by the quantum-chemical method (ADC2 with a small, non-polarized split-valence basis), leading to a very reasonable agreement.

The common motif of these molecule is that they consist of an electron-pushing triphenylamin group (right side of the molecule in the figure below) and an electron-pulling group (left side), whose strength can be modulated by oxidizing the contained sulphur atom. In the lowest excited singlet state of all of these systems, an electron is excited from the electron-pushing triphenylamine to the electron pulling part of the molecule. While the excitation-hole turned out to be essentially identical in all of these systems and all solvents, the location and structure of the excited electron differed significantly depending on the molecule and environment.

The figure below shows the structures of the molecules as well as the location of the excited electron computed at the ADC2/SS-PCM/SV level of theory for the non-polar solvent cyclohexane on the left and for the polar solvent acetonitrile in the right.

Since the discussion of the performance of the methodology was (as usual) quite brief in the article, I will share some more methodological insights here and in the following posts.

The central observation, which also made it into the paper, was that the solvent effects are systematically overestimated for this group of molecules, leading to too low emission energies in polar solvents (here acetonitrile, ACN). In the non-polar solvent cyclohexane (CHX) the error seems to cancel out with the systematic overestimation of the emission energies by the quantum-chemical method (ADC2 with a small, non-polarized split-valence basis), leading to a very reasonable agreement.

Something
that we only mentioned but did not show in the article was that the
spin-opposite scaled variant of ADC(2), SOS-ADC(2) [short: SOS(2)]
improves on this systematic error. After the article was published I
took a closer look at this and recomputed all of the numbers also with
SOS(2). Turns out this method pretty much nails the emission energies in
ACN, with those in CHX getting a little worse:

While
these trends are most certainly interesting to know for future
investigations (use SOS-ADC/SV for polar solvents and ADC/SV for
non-polar ones), the question that remains is if this systematic
overestimation of the solvent-stabilization is
due to A) the quantum-chemical methodology, B) the
solvent-model, C) the nature of these molecules or D) a combination of
all three.

Speaking for A) is that **I have observed before how ADC(2) yields too low energies for charge-transfer states**
compared to locally excited states. Since in this case the picture was
very consistent, including the fact that the problem is corrected by
ADC(3), I'm pretty confident that A) contributes at least to some
extent. The question that really bothers me, however, is if B) is also
the case, because so far, I am pretty sure that the self-consistent
treatment of solvent polarization is the one and only physically correct
way to do it. You can find my arguments for that **in the respective article (same as above).**

To
further look into the issue, let us eliminate the inherent error of the
quantum-chemical methodology as far as possible by looking at the CHX
to ACN shifts in the emission energy. Let us furthermore include the
emission energies from the first solvent-field iteration in the
comparison, i.e. the ones calculated with the solvent field obtained for
the excited-state computed in the relaxed ground-state solvent field
(ptSS-calculation).

Seeing
this plot for the first time literally shocked me. The agreement of the
one-shot approach with the experimental data is so apparent and
convincing that I immediately starting looking for the fundamental flaw
in the self-consistent approach. The latter apparently overestimated the
shifts by about a factor of two. Did I say about? For ADC(2), its
pretty much exactly a factor of two:

Can this be a coincidence? Exactly a factor of two? The lines are pretty much on top of each other. Yet, I'm pretty sure (as sure as it gets for a scientist) that this is indeed a coincidence, for a number of convincing reasons, which I will elaborate in the next post. And the best thing will be, that we can still learn from this coincidence!

So long,

Jan

In my last post I described how I calculated vibrationally resolved UV/vis absorption spectra for a couple of boronaromatics and their acetonitrile (ACN) adducts and how this helped to resolve the mystery of their odd solvatochromism (solvent dependent absorption spectra). The one question that remained open was: Do DBI and DBA exclusively form 1:1 monoadducts with ACN or also the 1:2 diadducts? Since the latter are predicted to have essentially no absorption in the relevant spectral range, i.e., between 300 and 700 nm, I could neither confirm nor exclude their formation.

To answer this question, I started to conduct what was planned as a "brief" thermochemical analysis of formation of the respective mono- and diadducts. However, since the first results were inconclusive the I ended up employing a hierarchy of methods of increasing sophistication (and computational demand), which only eventually (at the very highest level of theory) cumulated in a good enough agreement with the experimental observations. I think the convergence of the results with respect to level of theory and employed basis-sets is quite instructive (without having much experience in quantum-thermochemistry), but at the same time too long and theoretical to be included in the article, and hence provides a nice topic for this blog. Lets get to it!

The initial plan was to just collect and investigate the results from the calculations I had already conducted to model the vibrationally resolved spectra. These include all optimizations and normal-mode analyses for naked DBI, DBA and DBP as well as their monoadducts and trans-diadducts, such that solely the respective calculations for ACN were missing. Having completed the latter and putting together all the numbers (computed at the same level of theory as the vibronic spectra, i.e., B3LYP-D3BJ/SVP) afforded the following picture:

To answer this question, I started to conduct what was planned as a "brief" thermochemical analysis of formation of the respective mono- and diadducts. However, since the first results were inconclusive the I ended up employing a hierarchy of methods of increasing sophistication (and computational demand), which only eventually (at the very highest level of theory) cumulated in a good enough agreement with the experimental observations. I think the convergence of the results with respect to level of theory and employed basis-sets is quite instructive (without having much experience in quantum-thermochemistry), but at the same time too long and theoretical to be included in the article, and hence provides a nice topic for this blog. Lets get to it!

The initial plan was to just collect and investigate the results from the calculations I had already conducted to model the vibrationally resolved spectra. These include all optimizations and normal-mode analyses for naked DBI, DBA and DBP as well as their monoadducts and trans-diadducts, such that solely the respective calculations for ACN were missing. Having completed the latter and putting together all the numbers (computed at the same level of theory as the vibronic spectra, i.e., B3LYP-D3BJ/SVP) afforded the following picture:

The good news is that at least the ordering suits the experimental obervsations with the DBI monoadduct attaining the lowest energy, followed by DBA and eventually DBP. The bad news is that all values are about 20 kJ/mol too low to explain the experiment. The negative values of the free energies of formation ("+RT" values) for
all mono-adducts suggest a quantitative formation for all of
them, even if only traces of ACN are present. Remembering the experimental observations (no changes in the spectrum of DBP even in pure ACN, weak influence for DBA at high ACN concentrations, major changes only for DBI), the values are apparently systematically too low.

Considering that thermochemical calculations are (in contrast to the excited-state calculations) quite sensitive to basis-set size and that I have employed the small SVP basis, this should not really be a surprise. The magnitude of the effect, however, did surprise me. The underlying problem is that small incomplete basis sets are susceptible to the so-called basis-set superposition error (BSSE). In a nutshell, BSSE is the result of an artificially constrained electronic wavefunction. Consequently, any calculation for a larger system (e.g. adducts, dimers) attain energies that are systematically lowered compared to respective calculation for the respective subs systems (their isolated constituents, mononers), since in the supermolecular calculation the wavefunction of a fragment A can also use the basis-functions on fragment B and*vice versa*.

Although there are corrections for the BSSE, the straightforward approach is typically to use a larger (augmented) basis set, which is what I did. To save time in the calculations with the larger (def2-TZVP) basis set, I skipped the optimization of the the geometries and used the ones from before, which is a quite common thing to do since geometries are not as sensitive as energies to basis-set size. To save even more time I used another common trick: I did not redo the normal-mode analyses for the ZPE and RT corrections but used the ones obtained with SVP, since they are also much less basis-set dependent. Such a combined approach is abbreviated B3LYP-D3BJ/def2-TZVP//SVP (behind the double-slash follows the level used for the geometries/corrections).

Considering that thermochemical calculations are (in contrast to the excited-state calculations) quite sensitive to basis-set size and that I have employed the small SVP basis, this should not really be a surprise. The magnitude of the effect, however, did surprise me. The underlying problem is that small incomplete basis sets are susceptible to the so-called basis-set superposition error (BSSE). In a nutshell, BSSE is the result of an artificially constrained electronic wavefunction. Consequently, any calculation for a larger system (e.g. adducts, dimers) attain energies that are systematically lowered compared to respective calculation for the respective subs systems (their isolated constituents, mononers), since in the supermolecular calculation the wavefunction of a fragment A can also use the basis-functions on fragment B and

Although there are corrections for the BSSE, the straightforward approach is typically to use a larger (augmented) basis set, which is what I did. To save time in the calculations with the larger (def2-TZVP) basis set, I skipped the optimization of the the geometries and used the ones from before, which is a quite common thing to do since geometries are not as sensitive as energies to basis-set size. To save even more time I used another common trick: I did not redo the normal-mode analyses for the ZPE and RT corrections but used the ones obtained with SVP, since they are also much less basis-set dependent. Such a combined approach is abbreviated B3LYP-D3BJ/def2-TZVP//SVP (behind the double-slash follows the level used for the geometries/corrections).

Electronic energies obtained at the B3LYP-D3BJ/def2-TZVP//SVP level of theory. |

However, I wasn't quite satisfied with the agreement and more importantly just curious how ab-initio methods would compare to the DFT results. So I eventually conducted additional SCS-MP2 and CEPA/1 (a coupled-cluster variant) calculations. To also improve on the structures, I re-optimized them at the SCS-MP2/SVP level of theory and later refined them at the SCS-MP2/def2-TZVP level of theory (this was JUST possible for the DBI diadduct). For the final energies at the SCS-MP2 and CEPA/1 levels of theory, I even conducted complete basis-set (CBS) extrapolations using the def2-SVP and def2-TZVP sets. This technique estimates the energy that would be obtained with a hypothetical, complete set of basis functions from the differences between the energies of two limited sets. The free energies of formation for the monoadducts of DBI DBA and DBP are summarized in this final figure:

Summary of the free energies of formation of the monoadducts at (from left to right) increasing levels of theory. |

The systematic overbinding is improved but not eliminated at the mixed approach with the larger basis (def2-TZVP//SVP), and becomes slightly worse at the fully consistent SCS-MP2/def2-TZVP level of theory. Note that the difference between the mixed SVP//def2-TZVP and fully consistent def2-TZVP approaches is quite small, whereas the mixed approach is MUCH cheaper.

The CBS extrapolation, which I would have expected to correct the systematic errors, actually worsens the agreement with the experiment. All energies of formation are again systematically reduced, such that even the DBP monoadduct is suggested to be stable at this level. I had heard that MP2 tends to overbind organic molecules, but this is a larger error than I expected. In particular since I have used the spin-component scaled variant, which should be less prone to these problems.

Only at the CEPA/1/CBS level of theory do the calculated free energies of formation ultimately agree with the observed behaviour: DBI showing quantitative monoadduct formation is weakly bound, DBA showing weak but significant monoadduct formation at high MeCN concentrations is energy-neutral, and DBP, which does not show any signs of adduct formation, is strongly endothermic. To make the coupled-cluster calculations possible for these already quite large systems, I employed the very handy domain-localized pair natural orbital (DLPNO) approximation implemented in Orca, and still burned a lot of computer time.

At the time, I found it quite instructive (and ultimately satisfying) to see how these numbers eventually converge to agree with the experimentally observed behaviour at the highest level of theory, and how all of the methods show the issues I had heard about, but never observed before. I hope you has a similar experience and can take something home. If you have any questions please leave a comment.

The next post will be about the investigation of the mechanism underlying the different fluorescence quantum yields of DBI, DBA and DBP.

So long!

Jan

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:

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

Abonnieren
Posts (Atom)