| 1. |
- Aquilante, Francesco, et al.
(author)
-
Modern quantum chemistry with [Open]Molcas
- 2020
-
In: Journal of Chemical Physics. - : AIP Publishing. - 0021-9606 .- 1089-7690. ; 152:21
-
Journal article (peer-reviewed)abstract
- MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.
|
|
| 2. |
- De Kerpel, Jan O A, et al.
(author)
-
Geometric and Electronic Structure of Co(II)-Substituted Azurin
- 1999
-
In: The Journal of Physical Chemistry Part B. - : American Chemical Society (ACS). - 1520-5207 .- 1520-6106. ; 103:39, s. 8375-8382
-
Journal article (peer-reviewed)abstract
- The molecular and electronic structures of Co(II)-substituted azurin have been investigated using several realistic models of the metal coordination sphere. The geometry of the models was optimized using the hybrid density functional B3LYP method and compared to the structures obtained for similar Cu(II) models. It is found that Co(II) prefers a distorted tetrahedral structure with four strong bonds to two histidine nitrogens, the cysteine sulphur, and the backbone carbonyl group. This is in contrast to Cu(II), where two weak axial bonds to methionine and the backbone oxygen are found, combined with three strong bonds to the histidines and cysteine in the equatorial plane of a trigonal bipyramidal structure. The optimal structure of the models conforms with experimental crystal data, indicating that the active-site structure in these proteins is determined by the preferences of the metal ion and its ligand and not by protein strain. The electronic structure and spectrum of the Co(imidazole) 2(SH)(SH) 2(HCONH 2) + model have been investigated in detail using multiconfigurational second-order perturbation theory based on a complete active-space wavefunction (CASPT2). Nine ligand-field transitions and six S cys → Co charge-transfer transitions have been calculated, and all experimentally observed absorption bands in the absorption spectrum of Co(II) azurin have been assigned. It is shown that the Co-S cys bond is more ionic than the Cu-S cys bond and that this causes the blue shift and weakening of the charge-transfer states in the spectrum of Co(II)-substituted azurin compared to native copper protein.
|
|
| 3. |
- De Kerpel, Jan O A, et al.
(author)
-
Theoretical study of the structural and spectroscopic properties of stellacyanin
- 1998
-
In: The Journal of Physical Chemistry Part B. - : American Chemical Society (ACS). - 1520-5207 .- 1520-6106. ; 102:23, s. 4638-4647
-
Journal article (peer-reviewed)abstract
- The electronic spectrum of the azurin Met121Gln mutant, a model of the blue copper protein stellacyanin, has been studied by ab initio multiconfigurational second-order perturbation theory (the CASPT2 method), including the effect of the protein and solvent by point charges. The six lowest electronic transitions have been calculated and assigned with an error of less than 2400 cm-1. The ground-state singly occupied orbital is found to be a predominantly π antibonding orbital involving Cu3d and Scys3pπ. However, it also contains a significant amount (18%) of Cu-Scys σ antibonding character. This σ interaction is responsible for the appearance in the absorption spectrum of a band at 460 nm, with a significantly higher intensity than observed for other, axial, type 1 copper proteins (i.e., plastocyanin or azurin). The π-σ mixing is caused by the axial glutamine ligand binding at a much shorter distance to copper than the corresponding methionine ligand in the normal blue copper proteins. This explains why, based on its spectral properties, stellacyanin is classified among the rhombic type 1 copper proteins, although its structure is clearly trigonal, as it is for the axial proteins. Calculations have also been performed on structures where the glutamine model coordinates to the copper ion via the deprotonated N∈ atom instead of the O∈ atom. However, the resulting transition energies do not resemble the experimental spectrum obtained at normal or elevated pH. Thus, the results do not confirm the suggestion that the coordinating atom of glutamine changes at high pH.
|
|
| 4. |
- Delcey, Mickaël G., 1988-, et al.
(author)
-
Accurate calculations of geometries and singlet-triplet energy differences for active-site models of [NiFe] hydrogenase
- 2014
-
In: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 16:17, s. 7927-7938
-
Journal article (peer-reviewed)abstract
- We have studied the geometry and singlet-triplet energy difference of two mono-nuclear Ni2+ models related to the active site in [NiFe] hydrogenase. Multiconfigurational second-order perturbation theory based on a complete active-space wavefunction with an active space of 12 electrons in 12 orbitals, CASPT2(12,12), reproduces experimental bond lengths to within 1 pm. Calculated singlet-triplet energy differences agree with those obtained from coupled-cluster calculations with single, double and (perturbatively treated) triple excitations (CCSD(T)) to within 12 kJ mol(-1). For a bimetallic model of the active site of [NiFe] hydrogenase, the CASPT2(12,12) results were compared with the results obtained with an extended active space of 22 electrons in 22 orbitals. This is so large that we need to use restricted active-space theory (RASPT2). The calculations predict that the singlet state is 48-57 kJ mol(-1) more stable than the triplet state for this model of the Ni-Sl(a) state. However, in the [NiFe] hydrogenase protein, the structure around the Ni ion is far from the square-planar structure preferred by the singlet state. This destabilises the singlet state so that it is only similar to 24 kJ mol(-1) more stable than the triplet state. Finally, we have studied how various density functional theory methods compare to the experimental, CCSD(T), CASPT2, and RASPT2 results. Semi-local functionals predict the best singlet-triplet energy differences, with BP86, TPSS, and PBE giving mean unsigned errors of 12-13 kJ mol(-1) (maximum errors of 25-31 kJ mol(-1)) compared to CCSD(T). For bond lengths, several methods give good results, e. g. TPSS, BP86, and M06, with mean unsigned errors of 2 pm for the bond lengths if relativistic effects are considered.
|
|
| 5. |
- Dong, Geng, et al.
(author)
-
H2 binding to the active site of [NiFe] hydrogenase studied by multiconfigurational and coupled-cluster methods
- 2017
-
In: Physical Chemistry Chemical Physics. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 19:16, s. 10590-10601
-
Journal article (peer-reviewed)abstract
- [NiFe] hydrogenases catalyse the reversible conversion of molecular hydrogen to protons and electrons. This seemingly simple reaction has attracted much attention because of the prospective use of H2 as a clean fuel. In this paper, we have studied how H2 binds to the active site of this enzyme. Combined quantum mechanical and molecular mechanics (QM/MM) optimisation was performed to obtain the geometries, using both the TPSS and B3LYP density-functional theory (DFT) methods and considering both the singlet and triplet states of the Ni(ii) ion. To get more accurate energies and obtain a detailed account of the surroundings, we performed calculations with 819 atoms in the QM region. Moreover, coupled-cluster calculations with singles, doubles, and perturbatively treated triples (CCSD(T)) and cumulant-approximated second-order perturbation theory based on the density-matrix renormalisation group (DMRG-CASPT2) were carried out using three models to decide which DFT methods give the most accurate structures and energies. Our calculations show that H2 binding to Ni in the singlet state is the most favourable by at least 47 kJ mol-1. In addition, the TPSS functional gives more accurate energies than B3LYP for this system.
|
|
| 6. |
- Dong, Geng, et al.
(author)
-
Reaction Mechanism of [NiFe] Hydrogenase Studied by Computational Methods
- 2018
-
In: Inorganic Chemistry. - : American Chemical Society (ACS). - 0020-1669 .- 1520-510X. ; 57:24, s. 15289-15298
-
Journal article (peer-reviewed)abstract
- [NiFe] hydrogenases catalyze the reversible conversion of molecular hydrogen to protons and electrons. This seemingly simple reaction has attracted much attention because of the prospective use of H2 as a clean fuel. In this paper, we have studied the full reaction mechanism of this enzyme with various computational methods. Geometries were obtained with combined quantum mechanical and molecular mechanics (QM/MM) calculations. To get more accurate energies and obtain a detailed account of the surroundings, we performed big-QM calculations with 819 atoms in the QM region. Moreover, QM/MM thermodynamic cycle perturbation calculations were performed to obtain free energies. Finally, density matrix renormalisation group complete active space self-consistent field calculations were carried out to study the electronic structures of the various states in the reaction mechanism. Our calculations indicate that the Ni-L state is not involved in the reaction mechanism. Instead, the Ni-C state is reduced by one electron and then the bridging hydride ion is transferred to the sulfur atom of Cys546 as a proton and the two electrons transfer to the Ni ion. This step turned out to be rate-determining with an energy barrier of 58 kJ/mol, which is consistent with the experimental rate of 750 ± 90 s-1 (corresponding to ∼52 kJ/mol). The cleavage of the H-H bond is facile with an energy barrier of 33 kJ/mol, according to our calculations. We also find that the reaction energies are sensitive to the size of the QM system, the basis set, and the density functional theory method, in agreement with previous studies.
|
|
| 7. |
- Galván, Ignacio Fdez., et al.
(author)
-
OpenMolcas : From Source Code to Insight
- 2019
-
In: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 15:11, s. 5925-5964
-
Journal article (peer-reviewed)abstract
- In this Article we describe the OpenMolcas environment and invite the computational chemistry community to collaborate. The open-source project already includes a large number of new developments realized during the transition from the commercial MOLCAS product to the open-source platform. The paper initially describes the technical details of the new software development platform. This is followed by brief presentations of many new methods, implementations, and features of the OpenMolcas program suite. These developments include novel wave function methods such as stochastic complete active space self-consistent field, density matrix renormalization group (DMRG) methods, and hybrid multiconfigurational wave function and density functional theory models. Some of these implementations include an array of additional options and functionalities. The paper proceeds and describes developments related to explorations of potential energy surfaces. Here we present methods for the optimization of conical intersections, the simulation of adiabatic and nonadiabatic molecular dynamics, and interfaces to tools for semiclassical and quantum mechanical nuclear dynamics. Furthermore, the Article describes features unique to simulations of spectroscopic and magnetic phenomena such as the exact semiclassical description of the interaction between light and matter, various X-ray processes, magnetic circular dichroism, and properties. Finally, the paper describes a number of built-in and add-on features to support the OpenMolcas platform with postcalculation analysis and visualization, a multiscale simulation option using frozen-density embedding theory, and new electronic and muonic basis sets.
|
|
| 8. |
- Guo, Meiyuan, et al.
(author)
-
Fingerprinting Electronic Structure of Heme Iron by Ab Initio Modeling of Metal L-Edge X-ray Absorption Spectra
- 2019
-
In: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 15:1, s. 477-489
-
Journal article (peer-reviewed)abstract
- The capability of the multiconfigurational restricted active space approach to identify electronic structure from spectral fingerprints is explored by applying it to iron L-edge X-ray absorption spectroscopy (XAS) of three heme systems that represent the limiting descriptions of iron in the Fe-O-2 bond, ferrous and ferric [Fe(P)(ImH)(2)](0/1+) (P = porphine, ImH = imidazole), and Fe-II(P). The level of agreement between experimental and simulated spectral shapes is calculated using the cosine similarity, which gives a quantitative and unbiased assignment. Further dimensions in fingerprinting are obtained from the L-edge branching ratio, the integrated absorption intensity, and the edge position. The results show how accurate ab initio simulations of metal L-edge XAS can complement calculations of relative energies to identify unknown species in chemical reactions.
|
|
| 9. |
- Malmqvist, Per-Åke, et al.
(author)
-
The restricted active space followed by second-order perturbation theory method: Theory and application to the study of CuO(2) and Cu(2)O(2) systems.
- 2008
-
In: Journal of Chemical Physics. - : AIP Publishing. - 0021-9606 .- 1089-7690. ; 128:20
-
Journal article (peer-reviewed)abstract
- A multireference second-order perturbation theory using a restricted active space self-consistent field wave function as reference (RASPT2/RASSCF) is described. This model is particularly effective for cases where a chemical system requires a balanced orbital active space that is too large to be addressed by the complete active space self-consistent field model with or without second-order perturbation theory (CASPT2 or CASSCF, respectively). Rather than permitting all possible electronic configurations of the electrons in the active space to appear in the reference wave function, certain orbitals are sequestered into two subspaces that permit a maximum number of occupations or holes, respectively, in any given configuration, thereby reducing the total number of possible configurations. Subsequent second-order perturbation theory captures additional dynamical correlation effects. Applications of the theory to the electronic structure of complexes involved in the activation of molecular oxygen by mono- and binuclear copper complexes are presented. In the mononuclear case, RASPT2 and CASPT2 provide very similar results. In the binuclear cases, however, only RASPT2 proves quantitatively useful, owing to the very large size of the necessary active space.
|
|
| 10. |
- Olsson, Mats H M, et al.
(author)
-
On the relative stability of tetragonal and trigonal Cu(II) complexes with relevance to the blue copper proteins
- 1998
-
In: Journal of Biological Inorganic Chemistry. - : Springer Science and Business Media LLC. - 0949-8257 .- 1432-1327. ; 3:2, s. 109-125
-
Journal article (peer-reviewed)abstract
- The role of the cysteine thiolate ligand for the unusual copper coordination geometry in the blue copper proteins has been studied by comparing the electronic structure, geometry, and energetics of a number of small Cu(II) complexes. The geometries have been optimised with the density functional B3LYP method, and energies have been calculated by multi- configurational second-order perturbation theory (the CASPT2 method). Most small inorganic Cu(II) complexes assume a tetragonal geometry, where four ligands make σ bonds to a Cu 3d orbital. If a ligand lone-pair orbital instead forms a π bond to the copper ion, it formally occupies two ligand positions in a square coordination, and the structure becomes trigonal. Large, soft, and polarisable ligands, such as SH- and SeH-, give rise to covalent copper-ligand bonds and structures close to a tetrahedron, which might be trigonal or tetragonal with approximately the same stability. On the other hand, small and hard ligands, such as NH3, OH2, and OH-, give ionic bonds and flattened tetragonal structures. It is shown that axial type 1 (blue) copper proteins have a trigonal structure with a π bond to the cysteine sulphur atom, whereas rhombic type 1 and type 2 proteins have a tetragonal structure with σ bonds to all strong ligands. The soft cysteine ligand is essential for the stabilisation of a structure that is close to a tetrahedron (either trigonal or tetragonal), which ensures a low reorganisation energy during electron transfer.
|
|
