Towards quantum-chemical method development for arbitrary basis functions

We present the design of a flexible quantum-chemical method development framework, which supports employing any type of basis function. This design has been implemented in the light-weight program package molsturm, yielding a basis-function-independent self-consistent field scheme. Versatile interfaces, making use of open standards like python, mediate the integration of molsturm with existing third-party packages. In this way both rapid extension of the present set of methods for electronic structure calculations as well as adding new basis function types can be readily achieved. This makes molsturm well-suitable for testing novel approaches for discretising the electronic wave function and allows comparing them to existing methods using the same software stack. This is illustrated by two examples, an implementation of coupled-cluster doubles as well as a gradient-free geometry optimisation, where in both cases, an arbitrary basis functions could be used. molsturm is open-source and can be obtained from https://molsturm.org.Comment: 15 pages and 7 figure

Quantum chemistry with Coulomb Sturmians:Construction and convergence of Coulomb Sturmian basis sets at the Hartree-Fock level

The first discussion of basis sets consisting of exponentially decaying Coulomb Sturmian functions for modelling electronic structures is presented. The proposed basis set construction selects Coulomb Sturmian functions using separate upper limits to their principle, angular momentum and magnetic quantum numbers. Their common Coulomb Sturmian exponent is taken as a fourth parameter. The convergence properties of such basis sets are investigated for second and third row atoms at the Hartree-Fock level. Thereby important relations between the values of the basis set parameters and the physical properties of the electronic structure are recognised. For example, an unusually large limit for the angular momentum quantum number in unrestricted Hartree-Fock calculations can be linked to the breaking of spherical symmetry in such cases. Furthermore, a connection between the optimal, i.e. minimum-energy, Coulomb Sturmian exponent and the average Slater exponents values obtained by Clementi and Raimondi (E. Clementi and D. L. Raimondi, J. Chem. Phys. 38, 2686 (1963)) is made. These features of Coulomb Sturmian basis sets emphasise their ability to correctly reproduce the physical features of Hartree-Fock wave functions.Comment: 16 pages, 14 figures, supporting inf

Study of internal stresses in a TWIP steel analyzing transient and permanent softening during reverse shear tests.

Recent Bauschinger-type tests conducted on a twinning-induced plasticity (TWIP) steel highlights the important contribution of internal stresses to work hardening [1]. Along this line we present Bauschinger experiments in a Fe-22Mn wt.%-0.6C wt.% TWIP steel. The mechanical behaviour upon load reversal shows transient and permanent softening effects. Determination of the internal stress from the magnitude of the permanent softening yields a contribution to work hardening of the order of 20%. Analysis of the transient softening, during strain reversal, indicates that internal stress are consistent with reported data on high carbon spheroidized steels.Acknowledgements The authors would like to acknowledge the financial support by the German Research Foundation within the framework of the SFB 761 ‘‘steel ab initio’’ and the CICYT grant MAT2009-14452 awarded by the Spanish Ministry of Science and Innovation.Peer reviewe

Advances in Molecular Quantum Chemistry Contained in the Q-Chem 4 Program Package

A summary of the technical advances that are incorporated in the fourth major release of the Q-Chem quantum chemistry program is provided, covering approximately the last seven years. These include developments in density functional theory methods and algorithms, nuclear magnetic resonance (NMR) property evaluation, coupled cluster and perturbation theories, methods for electronically excited and open-shell species, tools for treating extended environments, algorithms for walking on potential surfaces, analysis tools, energy and electron transfer modelling, parallel computing capabilities, and graphical user interfaces. In addition, a selection of example case studies that illustrate these capabilities is given. These include extensive benchmarks of the comparative accuracy of modern density functionals for bonded and non-bonded interactions, tests of attenuated second order Møller–Plesset (MP2) methods for intermolecular interactions, a variety of parallel performance benchmarks, and tests of the accuracy of implicit solvation models. Some specific chemical examples include calculations on the strongly correlated Cr2 dimer, exploring zeolite-catalysed ethane dehydrogenation, energy decomposition analysis of a charged ter-molecular complex arising from glycerol photoionisation, and natural transition orbitals for a Frenkel exciton state in a nine-unit model of a self-assembling nanotube

Human cytomegalovirus immediate-early 1 protein rewires upstream STAT3 to downstream STAT1 signaling switching an IL6-type to an IFNγ-like response

MN and CP were supported by the Wellcome Trust (www.wellcome.ac.uk) Institutional Strategic Support Fund and CP was supported by the Deutsche Forschungsgemeinschaft (PA 815/2-1; www.dfg.de).The human cytomegalovirus (hCMV) major immediate-early 1 protein (IE1) is best known for activating transcription to facilitate viral replication. Here we present transcriptome data indicating that IE1 is as significant a repressor as it is an activator of host gene expression. Human cells induced to express IE1 exhibit global repression of IL6- and oncostatin M-responsive STAT3 target genes. This repression is followed by STAT1 phosphorylation and activation of STAT1 target genes normally induced by IFNγ. The observed repression and subsequent activation are both mediated through the same region (amino acids 410 to 445) in the C-terminal domain of IE1, and this region serves as a binding site for STAT3. Depletion of STAT3 phenocopies the STAT1-dependent IFNγ-like response to IE1. In contrast, depletion of the IL6 receptor (IL6ST) or the STAT kinase JAK1 prevents this response. Accordingly, treatment with IL6 leads to prolonged STAT1 instead of STAT3 activation in wild-type IE1 expressing cells, but not in cells expressing a mutant protein (IE1dl410-420) deficient for STAT3 binding. A very similar STAT1-directed response to IL6 is also present in cells infected with a wild-type or revertant hCMV, but not an IE1dl410-420 mutant virus, and this response results in restricted viral replication. We conclude that IE1 is sufficient and necessary to rewire upstream IL6-type to downstream IFNγ-like signaling, two pathways linked to opposing actions, resulting in repressed STAT3- and activated STAT1-responsive genes. These findings relate transcriptional repressor and activator functions of IE1 and suggest unexpected outcomes relevant to viral pathogenesis in response to cytokines or growth factors that signal through the IL6ST-JAK1-STAT3 axis in hCMV-infected cells. Our results also reveal that IE1, a protein considered to be a key activator of the hCMV productive cycle, has an unanticipated role in tempering viral replication.Publisher PDFPeer reviewe

Software for the frontiers of quantum chemistry:An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange–correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear–electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an “open teamware” model and an increasingly modular design