TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted

Time Resolved Photoelectron Spectroscopy of Thioflavin T Photoisomerization: A Simulation Study

The excited state isomerization of thioflavin T (ThT) is responsible for the quenching of its fluorescence in a non-restricted environment. The fluorescence quantum yield increases substantially upon binding to amyloid fibers. Simulations reveal that the variation of the twisting angle between benzothiazole and benzene groups (ϕ(1)) is responsible for the sub-picosecond fluorescence quenching. The evolution of the twisting process can be directly probed by photoelectron emission with energies ε ≥ 1.0 eV before the molecule reaches the ϕ(1)-twisted configuration (~300 fs)

Non-adiabatic Dynamics Using Time-Dependent Density Functional Theory: assessing the Coupling Strengths

Light-driven reactions constitute an important class of processes in physics, chemistry, and biology. The development of accurate and efficient computational tools for the study of excited states dynamics has thus become of primary importance. Recently, we have developed a non-adiabatic ab initio molecular dynamics (AIMD) method, which combines Tully's surface hopping with TDDFT. Here, we investigate some fundamental aspects of this AIMD scheme that deal with the accuracy of TDDFT potential energy surfaces in regions of strong coupling and with the associated intensity of the non-adiabatic couplings (NACs). Of particular interest is the coupling between the ground and the first excited singlet state, which constitute a potential pitfall for all non-adiabatic AND based on TDDFT. To this end, we have investigated the excited state dynamics of protonated formaldimine (CH2NH2+) with particular emphasis on the analysis of the NAC strengths in regions of the configurational space close to surface hopping points and conical intersections. A good agreement of the structural and dynamical properties is found with respect to state averaged multiconfiguration self consistent field results. (C) 2009 Elsevier B.V. All rights reserved

Direct photolysis of carbonyl compounds dissolved in cloud and fog droplets

Gas phase photolysis is an important tropospheric sink for many carbonyl compounds, however the significance of direct photolysis of carbonyl compounds dissolved in cloud and fog droplets is uncertain. We develop a theoretical approach to assess the importance of aqueous photolysis for a series of carbonyls that possess carboxyl and hydroxyl functional groups by comparison with rates of other atmospheric processes. We use computationally and experimentally derived Henry's law parameters, hydration equilibrium parameters, aqueous hydroxyl radical (OH) rate constants, and optical extinction coefficients to identify types of compounds that will not have competitive aqueous photolysis rates. We also present molecular dynamics simulations of atmospherically relevant carbonyl compounds designed to estimate gas and aqueous phase extinction coefficients. In addition, experiments designed to measure the photolysis rate of glyceraldehyde, an atmospherically relevant water soluble organic compound, reveal that aqueous quantum yields are highly molecule-specific and cannot be extrapolated from measurements of structurally similar compounds. We find that only three out of the 92 carbonyl compounds investigated, pyruvic acid, 3-oxobutanoic acid, and 3-oxopropanoic acid, may have aqueous photolysis rates that exceed the rate of oxidation by dissolved OH. For almost all carbonyl compounds lacking α, β conjugation, atmospheric removal by direct photolysis in cloud and fog droplets can be neglected

Direct photolysis of carbonyl compounds dissolved in cloud and fog~droplets

Gas-phase photolysis is an important tropospheric sink for many carbonyl compounds; however the significance of direct photolysis of these compounds dissolved in cloud and fog droplets is uncertain. We develop a theoretical approach to assess the importance of aqueous photolysis for a series of carbonyls that possess carboxyl and hydroxyl functional groups by comparison with rates of other atmospheric processes. We use computationally and experimentally derived effective Henry's law constants, hydration equilibrium parameters, aqueous hydroxyl radical (OH) rate constants, and optical extinction coefficients to identify types of compounds that will (or will not) have competitive aqueous photolysis rates. We also present molecular dynamics simulations designed to estimate gas- and aqueous-phase extinction coefficients of unstudied atmospherically relevant compounds found in d-limonene and isoprene secondary organic aerosol. In addition, experiments designed to measure the photolysis rate of glyceraldehyde, an atmospherically relevant water-soluble organic compound, reveal that aqueous quantum yields are highly molecule-specific and cannot be extrapolated from measurements of structurally similar compounds. We find that only two out of the 92 carbonyl compounds investigated, pyruvic acid and acetoacetic acid, may have aqueous photolysis rates that exceed the rate of oxidation by dissolved OH. For almost all carbonyl compounds lacking α,β-conjugation that were investigated, atmospheric removal by direct photolysis in cloud and fog droplets can be neglected under typical atmospheric conditions

Mixed time-dependent density-functional theory/classical photodynamics study of oxirane photochemistry

We present a mixed time-dependent density-functional theory (TDDFT)/classical trajectory surface hopping (SH) study of the photochemical ring opening in oxirane. Previous preparatory work limited to the symmetric CC ring-opening pathways of oxirane concluded that the Tamm-Dancoff approximation (TDA) is important for improving the performance of TDDFT away from the equilibrium geometry. This observation is supported by the present TDDFT TDA/SH calculations which successfully confirm the main experimentally derived Gomer-Noyes mechanism for the photochemical CO ring opening of oxirane and, in addition, provide important state-specific information not easily accessible from experiments. In particular, we find that, while one of the lowest two excited states is photochemically relatively inert, excitation into the other excited state leads predominantly to rapid ring opening, cyclic-C2H4O→H2CH2O•. This is followed by hopping to the electronic ground state where hot (4000 K) dynamics leads to further reactions, namely, H2CH2O•→CH3CHO→H3+HO and CH4+CO. We note that, in the dynamics, we are not limited to following minimum energy pathways and several surface hops may actually be needed before products are finally reached. The performance of different functionals is then assessed by comparison of TDDFT and diffusion Monte Carlo potential energy curves along a typical TDDFT TDA/SH reaction path. Finally, although true (S0,S1) conical intersections are expected to be absent in adiabatic TDDFT, we show that the TDDFT TDA is able to approximate a conical intersection in this syste