Dataset Information

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Huckel-Mobius Interconversions in Expanded Porphyrins.

ABSTRACT: Localized orbital coupled cluster theory has recently emerged as a nonempirical alternative to DFT for large systems. Intuitively, one might expect such methods to perform less well for highly delocalized systems. In the present work, we apply both canonical CCSD(T) approximations and a variety of localized approximations to a set of flexible expanded porphyrins-macrocycles that can switch between Hückel, figure-eight, and Möbius topologies under external stimuli. Both minima and isomerization transition states are considered. We find that Möbius(-like) structures have much stronger static correlation character than the remaining structures, and that this causes significant errors in DLPNO-CCSD(T) and even DLPNO-CCSD(T1) approaches, unless TightPNO cutoffs are employed. If sub-kcal mol-1 accuracy with respect to canonical relative energies is required even for Möbius-type systems (or other systems plagued by strong static correlation), then Nagy and Kallay's LNO-CCSD(T) method with "tight" settings is the suitable localized approach. We propose the present POLYPYR21 data set as a benchmark for localized orbital methods or, more broadly, for the ability of lower-level methods to handle energetics with strongly varying degrees of static correlation.

SUBMITTER: Sylvetsky N

PROVIDER: S-EPMC7304861 | biostudies-literature | 2020 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Hückel-Möbius Interconversions in Expanded Porphyrins.

Sylvetsky Nitai N Banerjee Ambar A Alonso Mercedes M Martin Jan M L JML

Journal of chemical theory and computation 20200511 6

Localized orbital coupled cluster theory has recently emerged as a nonempirical alternative to DFT for large systems. Intuitively, one might expect such methods to perform less well for highly delocalized systems. In the present work, we apply both canonical CCSD(T) approximations and a variety of localized approximations to a set of flexible expanded porphyrins-macrocycles that can switch between Hückel, figure-eight, and Möbius topologies under external stimuli. Both minima and isomerization t ...[more]

PMID: 32338891

Similar Datasets

Project description:Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. DFT remains the workhorse for modeling such extended macrocycles, when taking into account their size and huge conformational flexibility. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, π···π stacking, steric effects, ring strain, and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wave function methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a stronger degree of static correlation than the Hückel and twisted-Hückel structures, and as a result the relative energies of singly twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between the extrapolated MP2/CBS energies estimated from cc-pV{T,Q}Z, diffuse-augmented heavy-aug-cc-pV{T,Q}Z and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n-particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol-1 with the CCSD(T) relative energies. Regarding DFT methods, the range-separated double hybrids, such as ωB97M(2) and B2GP-PLYP, outperform other functionals with RMSDs of 0.6 and 0.8 kcal mol-1, respectively. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended π-systems, the errors drop down to 1.9 kcal mol-1 for the revised revDOD-PBEP86-NL, which eliminates the same-spin correlation energy. Minnesota meta-GGA functionals with high fractions of exact exchange (M06-2X and M08-HX) also perform reasonably well, outperforming more robust and significantly less empirically parametrized functionals like SCAN0-D3.

Project description:The syntheses, interconversions, and spectroscopic properties of a set of iron carbonyl clusters containing an interstitial carbide are reported. This includes the low temperature X-ray structures of the six-iron clusters (Y)2[Fe6(μ6-C)(μ2-CO)4(CO)12] (1a-c; where Y = NMe4, NEt4, PPh4); the five-iron cluster [Fe5(μ5-C)(CO)15] (3); and the novel formulation of the five-iron cluster (NMe4)2[Fe5(μ5-C)(μ2-CO)(CO)13] (4). Also included in this set is the novel charge-neutral cluster, [Fe6(μ6-C)(CO)18] (2), for which we were unable to obtain a crystallographic structure. As synthetic proof for the identity of 2, we performed a closed loop of interconversions within a family of crystallographically defined species (1, 3, and 4): [Fe6]2- → [Fe6]0 → [Fe5]0 → [Fe5]2- → [Fe6]2-. The structural, spectroscopic, and electronic properties of this "missing link" cluster 2 were investigated by IR, Raman, XPS, and Mössbauer spectroscopies-as well as by DFT calculations. A single νCO feature (1965 cm-1) in the IR spectrum of 2, as well as a prominent Raman feature (νsymm = 1550 cm-1), are consistent with the presence of terminal carbonyls and a {(μ6-C)Fe6} arrangement of iron centers around the central carbide. The XPS of 2 exhibits a higher energy Fe 2p3/2 feature (707.4 eV) as compared to that of 1 (705.5 eV), consistent with the two-electron oxidation induced by treatment of 1 with two equivalents of [Fc](PF6) under CO atmosphere (for the two added CO ligands). DFT calculations indicate two axial and four equatorial Fe sites in 1, all of which have the same or similar oxidation states, for example, two Fe(0) and four Fe(+0.5). These assignments are supported by Mössbauer spectra for 1, which exhibit two closely spaced quadrupole doublets with δ = 0.076 and 0.064 mm s-1. The high-field Mössbauer spectrum of 2 (4.2 K) exhibits three prominent quadrupole doublets with δ = -0.18, -0.11, and +0.41 mm s-1. This indicates three pairs of chemically equivalent Fe sites. The first two pairs arise from irons of a similar oxidation state, while the last pair arises from irons in a different oxidation state, indicating a mixed-valent cluster. Variable field Mössbauer spectra for 2 were simulated assuming these two groups and a diamagnetic ground state. Taken together, the Mössbauer results and DFT calculations for 2 indicate two axial Fe(II) sites and four equatorial sites of lower valence, probably Fe(0). In the DFT optimized pentagonal bipyramidal structure for 2, the Fe(II)-Ccarbide distances are compressed (∼1.84 Å), while the Fe(0)-Ccarbide distances are elongated (∼2.05 Å). Analysis of the formulations for 1 (closo-square bipyramid) and 2 (nido-pentagonal bipyramid) is considered in the context of the textbook electron-counting rules of 14n+2 and 14n+4 for closo and nido clusters, respectively. This redox-dependent intracluster disproportionation of Fe oxidation states is concluded to arise from changes in bonding to the central carbide. A similar phenomenon may be promoted by the central carbide of the FeMoco cluster of nitrogenase, which may in turn stimulate N2 reduction.

Dataset Information

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Huckel-Mobius Interconversions in Expanded Porphyrins.

Publications

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Hückel-Möbius Interconversions in Expanded Porphyrins.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets