A Nonorthogonal Configuration Interaction Approach to Singlet Fission in Perylenediimide Compounds

Perylenediimide molecules constitute a family of chromophores that undergo singlet fission, a process in which an excited singlet state converts into lower energy triplets on two neighboring molecules, potentially increasing the efficiency of organic solar cells. Here, the nonorthogonal configuration interaction method is applied to study the effect of the different crystal packing of various perylenediimide derivatives on the relative energies of the singlet and triplet states, the intermolecular electronic couplings, and the relative rates for singlet fission. The analysis of the wave functions and electronic couplings reveals that charge transfer states play an important role in the singlet fission mechanism. Dimer conformations where the PDI molecules are at large displacements along the long axis and short on the short axis are posed as the most favorable for singlet fission. The role of the substituent at the imide group has been inspected concluding that, although it has no effect in the energies, for some conformations it significantly influences the electronic couplings, and therefore, replacing this substituent with hydrogen may introduce artifacts in the computational modeling of the PDI molecules.

GronOR: Scalable and Accelerated Nonorthogonal Configuration Interaction for Molecular Fragment Wave Functions

GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.

On the role of dynamic electron correlation in non-orthogonal configuration interaction with fragments

Two different approaches have been implemented to include the effect of dynamic electron correlation in the Non-Orthogonal Configuration Interaction for Fragments (NOCI-F) method. The first is based on shifting the diagonal matrix elements of the NOCI matrix, while the second incorporates the dynamic correlation explicitly in the fragment wave functions used to construct the many-electron basis functions of the NOCI. The two approaches are illustrated for the calculation of the electronic coupling relevant in singlet fission and the coupling of spin moments in organic radicals. Comparison of the calculated diabatic couplings, the NOCI energies and wave functions shows that dynamic electron correlation is not only efficiently but also effectively incorporated by the shifting approach and can largely affect the coupling between electronic states. Also, it brings the NOCI coupling of the spin moments in close agreement with benchmark calculations.

GronOR: Massively parallel and GPU-accelerated non-orthogonal configuration interaction for large molecular systems

GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented.

Ab initio modelling of magnetic anisotropy in Sr3NiPtO6

First principles calculations in the framework of Density Functional Theory (DFT) and wavefunction-based correlated methods have been performed to investigate in detail the magnetic anisotropy in Sr3NiPtO6. This material is known for the easy-plane anisotropy with a large anisotropy constant of about 7.5-9.3 meV. We find that by properly choosing the onsite Coulomb repulsion and exchange parameters, DFT can correctly explain the easy-plane magnetocrystalline anisotropy of the material, but the magnitude of the anisotropy constant is underestimated. On the other hand, a quantitative agreement with respect to experiments, both in the magnitude and direction of the magnetic anisotropy, can be recovered by using the wavefunction-based approach which is able to fully describe the multiplet physics. We also show that the presence of structural distortions of the local NiO6 coordination is crucial for stabilizing the magnetic anisotropy in this compound.

Magnetic interactions in calcium and sodium ladder vanadates

Magnetic interactions in ladder vanadates are determined with quantum chemical computational schemes using the embedded cluster model approach to represent the material. The available experimental data for calcium vanadate is accurately reproduced and the nature of the interladder interaction is established to be ferromagnetic. An analysis of the main contributions to the magnetic couplings is presented and the role of the covalently bonded apex oxygen is elucidated. In the sodium vanadate, the ground state configuration of the rungs is V-3d1-O-2p5-V-3d1. We show that with this configuration good intrachain coupling constants are obtained for the high-temperature phase. The interchain coupling in NaV2O5 is predicted to be approximately 34 K, ferromagnetic in nature.

Theoretical characterization of the ground and optically excited states of alpha'-NaV2O5

The character of the ground and optically excited states was investigated by quantum chemical calculations. We propose a rung ground state with V 3d(1)(xy)-O 2p(1)(y)-V 3d(1)(xy) character, instead of the conventional picture of one unpaired electron shared by 2 V ions. The unpaired electron on O is low-spin coupled to the V d electrons and spin density is predicted to be localized on vanadium. The absorption peak at 0.9 eV is assigned to a state with similar orbital occupations but a different spin coupling scheme, resulting in spin density localized on the bridging oxygen.

Atomic many-body effects for the p-shell photoelectron spectra of transition metals

Ab initio theoretical results for the 2p- and 3p-hole states of an Mn(2+) ion are reported in order to determine the importance of atomic contributions to the photoelectron spectra of bulk MnO. A combined treatment of relativity and electron correlation reveals important physical effects that have been neglected in virtually all previous work. The many-body and relativistic effects included in the atomic model are able, without any ad hoc empirical parameters, to explain most of the features of the MnO photoelectron spectra. In particular, it is not necessary to invoke charge transfer to explain the complex p-level spectra.

Specificity of multiple homologous genomic regions in Spodoptera exigua nucleopolyhedrovirus DNA replication

The region upstream of the Spodoptera exigua multicapsid nucleopolyhedrovirus (SeMNPV) ubiquitin gene contains four near-identical 68-bp-long palindromic repeats. This region, named Sehr6 and located at map unit (m.u.) 88 of the SeMNPV genome on pSeEcoRI-2.2, showed structural homology to previously identified homologous regions (hrs) in a number of other baculoviruses. Hrs function as enhancers of transcription and as putative origins (oris) of baculovirus DNA replication. Five additional hrs (Sehr1-Sehr5) were identified on the SeMNPV genome by Southern blot hybridization with an 18-bp-long oligonucleotide complementary to a sequence conserved within the arms of the four palindromic repeats of Sehr6. Sehr1-Sehr6 were dispersed on the SeMNPV genome at m.u. 8.0, 30.0, 38.5, 51.0, 77.0 and 88.0, respectively. Sequence analysis of these hrs confirmed the presence of palindromic repeats, highly similar to those found in pSeEcoRI-2.2. The number of palindromes varied from one (Sehr4) to nine (Sehr1) per hr. The Sehrs are all present in non-coding regions of the SeMNPV genome and also contain multiple putative transcription recognition sequences. Plasmids containing either of the Sehrs replicated in an SeMNPV-dependent DNA replication assay. The Sehrs were unable to replicate in an AcMNPV-dependent DNA replication assay. This was in contrast to the previously observed SeMNPV non-hr type ori, which replicated in the presence of both AcMNPV and SeMNPV. These data suggest that the replication of SeMNPV and the role of hrs in this process is highly specific.

Identification and functional analysis of a non-hr origin of DNA replication in the genome of Spodoptera exigua multicapsid nucleopolyhedrovirus

The genome of Spodoptera exigua multicapsid nucleopolyhedrovirus (SeMNPV) was screened for the presence of putative origins of DNA replication (oris). Using a transient DNA replication assay, several fragments were identified that underwent SeMNPV-dependent DNA replication in Spodoptera frugiperda cells (Sf-AE-21). Preliminary sequence data revealed the presence of multiple copies of homologous repeats (hrs). Restriction fragment Xbal-F2 showed a distinct sequence reminiscent of Autographa californica and Orgyla pseudotsugata MNPV (AcMNPV and OpMNPV) non-hr oris. Deletion analysis of this fragment indicated that the essential sequences of this putative non-hr ori mapped within a region of 800 bp. Sequence analysis of this region showed a unique distribution of six different (im)-perfect palindromes, several polyadenylation motifs and the occurrence of multiple direct repeats. No sequence homology or similarities to other reported baculovirus oris were detected. The spatial and modular distribution of these motifs are similar to those of the non-hr oris of AcMNPV and OpMNPV. Comparison of baculovirus non-hr and consensus eukaryotic oris revealed no consensus ori but indicated that each of the non-hrs studied so far is unique. From the structural similarity, however, it was concluded that the SeMNPV Xbal-F2 ori represents a baculovirus non-hr type ori. In addition, evidence is provided that SeMNPV renders more specificity to baculovirus DNA replication than AcMNPV.

Characterization of the ecdysteroid UDP-glucosyltransferase gene of a single nucleocapsid nucleopolyhedrovirus of Buzura suppressaria

A putative ecdysteroid UDP-glucosyltransferase (egt) gene was identified in the single nucleocapsid nucleopolyhedrovirus of Buzura suppressaria (BusuNPV). This is the first egt gene to be characterized in a SNPV, suggesting that egt genes are prevalent in nucleopolyhedroviruses and possibly in all baculoviruses. The open reading frame (ORF) of the gene is 1539 nucleotides (nt) long, encoding a putative protein (EGT) of 513 amino acids (aa) with a M of 58922. The 5' noncoding region contains three possible TATA boxes. A polyadenylation signal, AATAAA, was found downstream of the translation stop codon. A putative signal peptide of 16 residues was present at the N-terminus of the EGT. The BusuNPV egt gene has a high degree of nucleotide and amino acid sequence homology to multiple nucleocapsid (M) NPV egt genes, the highest being to the Spodoptera exigua MNPV egt. A phylogenetic tree of eleven known EGTs was constructed using maximum parsimony analysis.

A hybrid baculovirus-bacteriophage T7 transient expression system

A hybrid recombinant baculovirus-bacteriophage T7 expression system was developed for transient expression in insect cells of plasmids with foreign genes provided with a T7 promoter. The coding sequence for T7 RNA polymerase, with or without a nuclear localization signal, was inserted into the genome of Autographa californica nuclear polyhedrosis virus. Recombinant viruses stably expressed T7 RNA polymerase in insect cells. Upon transfection of infected insect cells with plasmids containing the genes for chloramphenicol acetyltransferase (CAT), the hepatitis B virus precore-, core- or e- antigens under control of the T7 promoter, transient expression of these genes was detected by ELISA. The results obtained indicate that this baculovirus/T7 system provides a simple and widely applicable tool for transient gene expression studies.

Nucleotide sequence and genetic organization of a 7.3 kb region (map unit 47 to 52.5) of Autographa californica nuclear polyhedrosis virus fragment EcoRI-C

The nucleotide sequence and genetic organization of a 7297 bp region within the EcoRI-C fragment of Autographa californica multiple nucleocapsid nuclear polyhedrosis virus (AcMNPV) are presented. Eight putative open reading frames were found and their respective amino acid sequences compared with a number of data libraries. ORF 1227 corresponded with gp41 and its predicted protein sequence was found to be 55 amino acids longer at its C terminus than reported previously. Moreover the main part of the ORF 1227 product, including the additional 55 amino acids, showed a high degree of homology with protein p40 of Helicoverpa zea single nucleocapsid nuclear polyhedrosis virus (HzSNPV). Three other ORFs in the analysed AcMNPV region showed homology with ORFs in the HzSNPV sequence, indicating that the general organization of this region is similar in both viruses. However one ORF found in the AcMNPV sequence was absent from the corresponding HzSNPV sequence.

Glycoprotein H of pseudorabies virus is essential for entry and cell-to-cell spread of the virus

To study the function of the envelope glycoprotein gH of pseudorabies virus, a gH null mutant was constructed. A premature translation termination codon was introduced in the gH gene by linker insertion mutagenesis, and a mutant virus was rescued by using a cell line that expresses the wild-type protein. Mutant virus isolated from complementing cells was unable to form plaques on noncomplementing cells, indicating that gH is essential in the life cycle of the virus. Immunological staining and electron microscopy showed that the mutant virus produced noninfectious progeny and was unable to spread from infected to uninfected cells by cell-cell fusion. Thus, similar to gH of herpes simplex virus, gH of pseudorabies virus is required for entry and cell-to-cell spread.