On-Surface Synthesis and Determination of the Open-Shell Singlet Ground State of Tridecacene.

The character of the electronic structure of acenes has been the subject of longstanding discussion. However, convincing experimental evidence of their open-shell character has so far been missing. Here, we present the on-surface synthesis of tridecacene molecules by thermal annealing of octahydrotridecacene on a Au(111) surface. We characterized the electronic structure of the tridecacene by scanning probe microscopy, which reveals the presence of an inelastic signal at 126 meV. We attribute the inelastic signal to spin excitation from the singlet diradical ground state to the triplet excited state. To rationalize the experimental findings, we carried out many-body ab initio calculations as well as model Hamiltonians to take into account the effect of the metallic substrate. Moreover, we provide a detailed analysis of how the dynamic electron correlation and virtual charge fluctuation between the molecule and metallic surface reduces the singlet-triplet band gap. Thus, this work provides the first experimental confirmation of the magnetic character of tridecacene.

Hilbert space multireference coupled cluster tailored by matrix product states.

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants-the state universal one and the Brillouin-Wigner and Mukherjee's state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

TAMM: Tensor algebra for many-body methods.

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

Interplay between π-Conjugation and Exchange Magnetism in One-Dimensional Porphyrinoid Polymers.

The synthesis of novel polymeric materials with porphyrinoid compounds as key components of the repeating units attracts widespread interest from several scientific fields in view of their extraordinary variety of functional properties with potential applications in a wide range of highly significant technologies. The vast majority of such polymers present a closed-shell ground state, and, only recently, as the result of improved synthetic strategies, the engineering of open-shell porphyrinoid polymers with spin delocalization along the conjugation length has been achieved. Here, we present a combined strategy toward the fabrication of one-dimensional porphyrinoid-based polymers homocoupled via surface-catalyzed [3 + 3] cycloaromatization of isopropyl substituents on Au(111). Scanning tunneling microscopy and noncontact atomic force microscopy describe the thermal-activated intra- and intermolecular oxidative ring closure reactions as well as the controlled tip-induced hydrogen dissociation from the porphyrinoid units. In addition, scanning tunneling spectroscopy measurements, complemented by computational investigations, reveal the open-shell character, that is, the antiferromagnetic singlet ground state (S = 0) of the formed polymers, characterized by singlet-triplet inelastic excitations observed between spins of adjacent porphyrinoid units. Our approach sheds light on the crucial relevance of the π-conjugation in the correlations between spins, while expanding the on-surface synthesis toolbox and opening avenues toward the synthesis of innovative functional nanomaterials with prospects in carbon-based spintronics.

On-Surface Synthesis and Characterization of [7]Triangulene Quantum Ring.

The ability to engineer geometrically well-defined antidots in large triangulene homologues allows for creating an entire family of triangulene quantum rings (TQRs) with tunable high-spin ground state, crucial for next-generation molecular spintronic devices. Herein, we report the synthesis of an open-shell [7]triangulene quantum ring ([7]TQR) molecule on Au(111) through the surface-assisted cyclodehydrogenation of a rationally designed kekulene derivative. Bond-resolved scanning tunneling microscopy (BR-STM) unambiguously imaged the molecular backbone of a single [7]TQR with a triangular zigzag edge topology, which can be viewed as [7]triangulene decorated with a coronene-like antidot in the center. Additionally, dI/dV mapping reveals that both inner and outer zigzag edges contribute to the edge-localized and spin-polarized electronic states of [7]TQR. Both experimental results and spin-polarized density functional theory calculations indicate that [7]TQR retains its open-shell septuple ground state (S = 3) on Au(111). This work demonstrates a new route for the design of high-spin graphene quantum rings for future quantum devices.

On-Surface Strain-Driven Synthesis of Nonalternant Non-Benzenoid Aromatic Compounds Containing Four- to Eight-Membered Rings.

The synthesis of polycyclic aromatic hydrocarbons containing various non-benzenoid rings remains a big challenge facing contemporary organic chemistry despite a considerable effort made over the last decades. Herein, we present a novel route, employing on-surface chemistry, to synthesize nonalternant polycyclic aromatic hydrocarbons containing up to four distinct kinds of non-benzenoid rings. We show that the surface-induced mechanical constraints imposed on strained helical reactants play a decisive role leading to the formation of products, energetically unfavorable in solution, with a peculiar ring current stabilizing the aromatic character of the π-conjugated system. Determination of the chemical and electronic structures of the most frequent product reveals its closed-shell character and low band gap. The present study renders a new route for the synthesis of novel nonalternant polycyclic aromatic hydrocarbons or other hydrocarbons driven by internal stress imposed by the surface not available by traditional approaches of organic chemistry in solution.

Massively parallel quantum chemical density matrix renormalization group method.

We present, to the best of our knowledge, the first attempt to exploit the super-computer platform for quantum chemical density matrix renormalization group (QC-DMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candidates for QC-DMRG calculations. In case of the largest calculation, which is the nitrogenase FeMo cofactor cluster with the active space comprising 113 electrons in 76 orbitals and bond dimension equal to 6000, our parallel approach scales up to approximately 2000 CPU cores.

Perturbative triples correction to domain-based local pair natural orbital variants of Mukherjee's state specific coupled cluster method.

In this article we report an implementation of the perturbative triples correction to Mukherjee's state-specific multireference coupled cluster method based on the domain-based pair natural orbital approach (DLPNO-MkCC). We tested the performance of DLPNO-MkCCSD(T) in calculations involving tetramethyleneethane and isomers of naphthynes. These tests show that more than 97% of triples energy was recovered with respect to the canonical MkCCSD(T) method, which together with the DLPNO-MkCCSD part accounts for about 99.70-99.85% of the total correlation energy. The applicability of the method was demonstrated on calculations of singlet-triplet gaps for several large systems: triangulene, dynemicin A, and a beryllium complex.

Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen).

We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.

Strong fluctuations in aboveground population size do not limit genetic diversity in populations of an endangered biennial species.

Assessing genetic diversity within populations of rare species and understanding its determinants are crucial for effective species protection. While a lot is known about the relationships between genetic diversity, fitness, and current population size, very few studies explored the effects of past population size. Knowledge of past population size may, however, improve our ability to predict future population fates. We studied Gentianella praecox subsp. bohemica, a biennial species with extensive seed bank. We tested the effect of current, past minimal and maximal population size, and harmonic mean of population sizes within the last 15 years on genetic diversity and fitness. Maximum population size over the last 15 years was the best predictor of expected heterozygosity of the populations and was significantly related to current population size and management. Plant fitness was significantly related to current as well as maximum population size and expected heterozygosity. The results suggested that information on past population size may improve our understanding of contemporary genetic diversity across populations. They demonstrated that despite the strong fluctuations in population size, large reductions in population size do not result in immediate loss of genetic diversity and reduction of fitness within the populations. This is likely due to the seed bank of the species serving as reservoir of the genetic diversity of the populations. From a conservation point of view, this suggests that the restoration of small populations of short-lived species with permanent seed bank is possible as these populations may still be genetically diverse.

Reduced-cost sparsity-exploiting algorithm for solving coupled-cluster equations.

We present an algorithm for reducing the computational work involved in coupled-cluster (CC) calculations by sparsifying the amplitude correction within a CC amplitude update procedure. We provide a theoretical justification for this approach, which is based on the convergence theory of inexact Newton iterations. We demonstrate by numerical examples that, in the simplest case of the CCD equations, we can sparsify the amplitude correction by setting, on average, roughly 90% nonzero elements to zeros without a major effect on the convergence of the inexact Newton iterations.

Perturbative universal state-selective correction for state-specific multi-reference coupled cluster methods.

In this work, we report an extension of our previous development of the universal state-selective (USS) multireference coupled-cluster (MRCC) formalism. It was shown [Brabec et al., J. Chem. Phys. 136, 124102 (2012)] and [Banik et al., J. Chem. Phys. 142, 114106 (2015)] that the USS(2) approach significantly improves the accuracy of Brillouin-Wigner and Mukherjee MRCC formulations, however, the numerical and storage costs associated with calculating highly excited intermediates pose a significant challenge, which can restrict the applicability of the USS(2) method. Therefore, we introduce a perturbative variant of the USS(2) approach (USS(pt)), which substantially reduces numerical overhead of the full USS(2) correction while preserving its accuracy. Since the new USS(pt) implementation calculates the triple and quadruple projections in on-the-fly manner, the memory bottleneck associated with the need of storing expensive recursive intermediates is entirely eliminated. On the example of several benchmark systems, we demonstrate accuracies of USS(pt) and USS(2) approaches and their efficiency in describing quasidegenerate electronic states. It is also shown that the USS(pt) method significantly alleviates problems associated with the lack of invariance of MRCC theories upon the rotation of active orbitals.

Iterative universal state selective correction for the Brillouin-Wigner multireference coupled-cluster theory.

As a further development of the previously introduced a posteriori Universal State-Selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011); J. Brabec et al., ibid. 136, 124102 (2012)], we suggest an iterative form of the USS correction by means of correcting effective Hamiltonian matrix elements. We also formulate USS corrections via the left Bloch equations. The convergence of the USS corrections with excitation level towards the full configuration interaction (FCI) limit is also investigated. Various forms of the USS and simplified diagonal USS corrections at the singles and doubles and perturbative triple levels are numerically assessed on several model systems and on the ozone and tetramethyleneethane molecules. It is shown that the iterative USS correction can successfully replace the previously developed a posteriori Brillouin-Wigner coupled cluster size-extensivity correction, while it is not sensitive to intruder states and performs well also in other cases when the a posteriori one fails, like, e.g., for the asymmetric vibration mode of ozone.

Highly entangled polyradical nanographene with coexisting strong correlation and topological frustration.

Open-shell nanographenes exhibit unconventional π-magnetism arising from topological frustration or strong electron-electron interaction. However, conventional design approaches are typically limited to a single magnetic origin, which can restrict the number of correlated spins or the type of magnetic ordering in open-shell nanographenes. Here we present a design strategy that combines topological frustration and electron-electron interactions to fabricate a large fully fused 'butterfly'-shaped tetraradical nanographene on Au(111). We employ bond-resolved scanning tunnelling microscopy and spin-excitation spectroscopy to resolve the molecular backbone and reveal the strongly correlated open-shell character, respectively. This nanographene contains four unpaired electrons with both ferromagnetic and anti-ferromagnetic interactions, harbouring a many-body singlet ground state and strong multi-spin entanglement, which is well described by many-body calculations. Furthermore, we study the magnetic properties and spin states in the nanographene using a nickelocene magnetic probe. The ability to imprint and characterize many-body strongly correlated spins in polyradical nanographenes paves the way for future advancements in quantum information technologies.

Role of Spin Polarization and Dynamic Correlation in Singlet-Triplet Gap Inversion of Heptazine Derivatives.

The new generation of proposed light-emitting molecules for organic light-emitting diodes (OLEDs) has raised considerable research interest due to its exceptional feature─a negative singlet-triplet (ST) gap violating Hund's multiplicity rule in the excited S1 and T1 states. We investigate the role of spin polarization in the mechanism of ST gap inversion. Spin polarization is associated with doubly excited determinants of certain types, whose presence in the wave function expansion favors the energy of the singlet state more than that of the triplet. Using a perturbation theory-based model for spin polarization, we propose a simple descriptor for prescreening of candidate molecules with negative ST gaps and prove its usefulness for heptazine-type molecules. Numerical results show that the quantitative effect of spin polarization decreases linearly with the increasing highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) exchange integral. Comparison of single- and multireference coupled-cluster predictions of ST gaps shows that the former methods provide good accuracy by correctly balancing the effects of doubly excited determinants and dynamic correlation. We also show that accurate ST gaps may be obtained using a complete active space model supplemented with dynamic correlation from multireference adiabatic connection theory.

Communication: Application of state-specific multireference coupled cluster methods to core-level excitations.

The concept of the model space underlying multireference coupled-cluster (MRCC) formulations is a powerful tool to deal with complex correlation effects for various electronic states. Here, we demonstrate that iterative state-specific MRCC methods (SS-MRCC) based on properly defined model spaces can be used to describe core-level excited states even when Hartree-Fock orbitals are utilized. We show that the SS-MRCC models with single and double excitations are comparable in accuracy to high-level single reference equation-of-motion coupled cluster (EOMCC) formalism.

Universal state-selective corrections to multi-reference coupled-cluster theories with single and double excitations.

The recently proposed universal state-selective (USS) corrections [K. Kowalski, J. Chem. Phys. 134, 194107 (2011)] to approximate multi-reference coupled-cluster (MRCC) energies can be commonly applied to any type of MRCC theory based on the Jeziorski-Monkhorst [B. Jeziorski and H. J. Monkhorst, Phys. Rev. A 24, 1668 (1981)] exponential ansatz. In this paper we report on the performance of a simple USS correction to the Brillouin-Wigner and Mukherjee's MRCC approaches employing single and double excitations (USS-BW-MRCCSD and USS-Mk-MRCCSD). It is shown that the USS-BW-MRCCSD correction, which employs the manifold of single and double excitations, can be related to a posteriori corrections utilized in routine BW-MRCCSD calculations. In several benchmark calculations we compare the USS-BW-MRCCSD and USS-Mk-MRCCSD results with the results obtained with the full configuration interaction method.

Implementation of the multireference Brillouin-Wigner and Mukherjee's coupled cluster methods with non-iterative triple excitations utilizing reference-level parallelism.

In this paper we discuss the performance of the non-iterative state-specific multireference coupled cluster (SS-MRCC) methods accounting for the effect of triply excited cluster amplitudes. The corrections to the Brillouin-Wigner and Mukherjee's MRCC models based on the manifold of singly and doubly excited cluster amplitudes (BW-MRCCSD and Mk-MRCCSD, respectively) are tested and compared with exact full configuration interaction results for small systems (H(2)O, N(2), and Be(3)). For the larger systems (naphthyne isomers) the BW-MRCC and Mk-MRCC methods with iterative singles, doubles, and non-iterative triples (BW-MRCCSD(T) and Mk-MRCCSD(T)) are compared against the results obtained with single reference coupled cluster methods. We also report on the parallel performance of the non-iterative implementations based on the use of processor groups.

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron-Sulfur Cubanes.

Iron-sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. To provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III)4S4(SMe)4]-2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

Volume of the amygdala is reduced in patients with narcolepsy - a structural MRI study.

OBJECTIVE: Based on the clinical observation that patients suffering from narcolepsy with cataplexy (NC) have cataplectic attacks when they experience positive emotions, it is therefore hypothesised that the abnormal processing of external emotional input through the limbic system, or motor dysregulation induced by emotions, takes place during these episodes. To date, imaging studies have failed to reveal consistent brain abnormalities in NC patients. METHODS: Considering the discrepancies in reported structural or functional abnormalities of the hypothalamus, amygdala, and nucleus accumbens, we used the MRI volumetry to determine the volumes of the amygdala and nucleus accumbens in a group of eleven patients with NC (5 males and 6 females, mean age 41.7 years ± 17.7). This data was compared to an equal number of examinations in healthy volunteers matched for age and gender. RESULTS: We found a decrease in the amygdalar volume of NC patients in both raw (p<0.001) and relative (p<0.01) data sets. The difference in amygdalar volume between healthy volunteers and NC patients was about 17%. In contrast to the amygdala, we did not find any differences in the volumes of nucleus accumbens. CONCLUSION: In the present MRI volumetric study, we found bilateral gray matter loss in the amygdala only.