ABSTRACT
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
ABSTRACT
As a typical antibiotic pollutant, tetracycline (TC) is producing increasing threats to the ecosystem and human health, and exploring convenient means for monitoring of TC is needed. Here, we proposed alkali-etched imprinted Mn-based Prussian blue analogues featuring superior oxidase-mimetic activity and precise recognition for the colorimetric sensing of TC. Simply etching Mn-based Prussian blue analogues (Mn-PBAs) with NaOH could expose the sites and surfaces to significantly improve their catalytic activity. Density functional theory calculations were employed to screen the molecularly imprinted polymer (MIP) layer for target identification. Consequently, the designed Mn-PBANaOH@MIP possessed the rich channels for substrates to get in touch with the active Mn-PBANaOH core, showing an excellent catalytic capacity to trigger the chromogenic oxidation of 3,3',5,5'-tetramethylbenzidine (TMB) without the use of H2O2. If TC was introduced, it would be recognized selectively by the MIP shell and masked the channels for TMB access, resulting in the obstruction of the chromogenic reaction. According to this mechanism, selective optical detection of TC was achieved, and performance stability, reusability, and reliability as well as practicability were also verified, promising potential for TC monitoring in complex matrices. Our work not only presents an effective way to enhance the enzyme-like activity of Prussian blue analogues but also provides a facile approach for TC sensing. Additionally, the work will inspire the exploration of molecularly imprinted nanozymes for various applications.
Subject(s)
Colorimetry , Oxidoreductases , Humans , Colorimetry/methods , Hydrogen Peroxide , Ecosystem , Reproducibility of Results , Sodium Hydroxide , Tetracycline , Anti-Bacterial AgentsABSTRACT
Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. The present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be nonzero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T2-S1 state mixing in pyrazine, spin-forbidden versus spin-allowed 4f-5d transitions in the Tb3+ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy)3].
ABSTRACT
A comprehensive ab initio study of periodic actinide-ligand bonding trends for trivalent actinides is performed. Relativistic density functional theory (DFT) and complete active-space (CAS) self-consistent field wavefunction calculations are used to dissect the chemical bonding in the [AnCl6]3-, [An(CN)6]3-, [An(NCS)6]3-, [An(S2PMe2)3], [An(DPA)3]3-, and [An(HOPO)]- series of actinide (An = U-Es) complexes. Except for some differences for the early actinide complexes with DPA, bond orders and excess 5f-shell populations from donation bonding show qualitatively similar trends in 5fâ¯n active-space CAS vs DFT calculations. The influence of spin-orbit coupling on donation bonding is small for the tested systems. Along the actinide series, chemically soft vs chemically harder ligands exhibit clear differences in bonding trends. There are pronounced changes in the 5f populations when moving from Pu to Am or Cm, which correlate with previously noted "breaks" in chemical trends. Bonding involving 5f becomes very weak beyond Cm/Bk. We propose that Cm(III) is a borderline case among the trivalent actinides that can be meaningfully considered to be involved in ground-state 5f covalent bonding.
ABSTRACT
Accurate Pt-An bond dissociation enthalpies (BDEs) for PtAnOn+ (An = U, Np, Pu and n = 0-2) and the corresponding enthalpies for the Pt + OAnOn+ substitution reactions have been studied for the first time using an accurate composite coupled cluster approach. Analogous O-AnOn+ bond dissociation enthalpies are also presented. To make the study possible, new correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are developed and reported for Pt and Au, with accompanying benchmark calculations of their atomic ionization potentials to demonstrate the effectiveness of the new basis sets. For the charged PtAnOn+ species (n = 1, 2), a low-spin state (LSS) for which the Pt-An σ bond is doubly occupied is studied together with a high-spin state (HSS) obtained by unpairing the σ bond orbital and placing one electron into the An 5f shell. The relative energies of the two spin states have been compared and qualitatively assessed via natural population and natural bond analyses. The enthalpies for the Pt substitution reactions, i.e., Pt + OAnOn+ â PtAnOn+ + O, are calculated to range from about 14-62 kcal/mol, and the Pt-AnOn+ bond dissociation enthalpies range from about 78-149 kcal/mol for the ground electronic states. For the PtAnO+ species, the LSSs were all predicted to be the ground state, whereas the PtAnO2+ molecules all favored the HSSs. The prediction for PtUO2+ is consistent with previous theoretical findings. The natural bond orbital analyses indicate a triple bond between An and O, with a double to quadruple bond between the An and Pt.
ABSTRACT
Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab initio wavefunction calculations at the multireference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron-nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on nonrelativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538-549] and therefore had limited applicability. Results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.
ABSTRACT
Predicting accurate nuclear magnetic resonance chemical shieldings relies upon cancellation of different types of errors between the theoretically calculated shielding constant of the analyte of interest and the reference. Often, the intrinsic error in computed shieldings due to basis sets, approximations in the Hamiltonian, description of the wave function, and dynamic effects is nearly identical between the analyte and reference, yet if the electronic structure or sensitivity to local environment differs dramatically, this cannot be taken for granted. Detailed prior work has examined the octahedral trivalent cation Al(H2O)6 3+, accounting for ab initio intrinsic errors. However, the use of this species as a reference for the chemically distinct tetrahedral anion Al(OH)4 - requires an understanding of how these errors cancel in order to define the limits of accurately predicting Al27 chemical shielding in Al(OH)4 -. In this work, we estimate the absolute shielding of the Al27 nucleus in Al(OH)4 - at the coupled cluster level (515.1 ± 5.3 ppm). Shielding sensitivity to the choice of method approximation and atomic basis sets used has been evaluated. Solvent and thermal effects are assessed through ensemble averaging techniques using ab initio molecular dynamics. The contribution of each type of intrinsic error is assessed for the Al(H2O)6 3+ and Al(OH)4 - ions, revealing significant differences that fundamentally hamper the ability to accurately calculate the Al27 chemical shift of Al(OH)4 - from first principles.
ABSTRACT
Thermochemical and spectroscopic properties for actinyl complexes involving UO22+/1+/0 and NpO22+/1+/0 with N2 and CO, together with the UO2-O2, UO2+-O2, and UO2+-NO complexes, have been studied for the first time using an accurate composite coupled cluster approach. Two general bonding motifs were investigated, end-on (η1) and side-on (η2) relative to the metal center of the actinyls. For end-on CO complexes, both C-coordinated (An-C) and O-coordinated (An-O) structures were examined, with the former always being lower in energy. All of the η1 complexes were calculated to be stable, with dissociation energies ranging from 2 to 36 kcal/mol, except for that of UO2+-O2 (the η1 orientation for UO2+-NO was not amenable to single reference coupled cluster). In agreement with a previous study, the η2 structure for UO2+-O2 was calculated to be relatively strongly bound, by 22.3 kcal/mol in this work. The closely related NO complex, however, had a calculated dissociation energy of just 4.0 kcal/mol. The binding energy of O2 to neutral UO2 in a η2 orientation was calculated to be very strong, 75.4 kcal/mol, and strongly resembled a UO2+(O2-) complex at equilibrium. The N-N and C-O bonds were found to be somewhat activated for all the side-on (η2) neutral An(IV) complexes, with stretching frequencies of N2 or CO being red-shifted by as much as 480 cm-1 with a 0.06 Šbond length elongation. Dissociation energies for the η1 complexes are strongly correlated with the extent of electron transfer from ligand to actinyl. The nature of bonding in the actinyl complexes is examined using natural resonance theory (NRT). The correlation between bonding motif and small molecule activation is in agreement with experiments in condensed phases.
ABSTRACT
A series of cerium(iv) mixed-ligand guanidinate-amide complexes, {[(Me3Si)2NC(N i Pr)2] x CeIV[N(SiMe3)2]3-x }+ (x = 0-3), was prepared by chemical oxidation of the corresponding cerium(iii) complexes, where x = 1 and 2 represent novel complexes. The Ce(iv) complexes exhibited a range of intense colors, including red, black, cyan, and green. Notably, increasing the number of the guanidinate ligands from zero to three resulted in significant redshift of the absorption bands from 503 nm (2.48 eV) to 785 nm (1.58 eV) in THF. X-ray absorption near edge structure (XANES) spectra indicated increasing f occupancy (n f) with more guanidinate ligands, and revealed the multiconfigurational ground states for all Ce(iv) complexes. Cyclic voltammetry experiments demonstrated less stabilization of the Ce(iv) oxidation state with more guanidinate ligands. Moreover, the Ce(iv) tris(guanidinate) complex exhibited temperature independent paramagnetism (TIP) arising from the small energy gap between the ground- and excited states with considerable magnetic moments. Computational analysis suggested that the origin of the low energy absorption bands was a charge transfer between guanidinate π orbitals that were close in energy to the unoccupied Ce 4f orbitals. However, the incorporation of sterically hindered guanidinate ligands inhibited optimal overlaps between Ce 5d and ligand N 2p orbitals. As a result, there was an overall decrease of ligand-to-metal donation and a less stabilized Ce(iv) oxidation state, while at the same time, more of the donated electron density ended up in the 4f shell. The results indicate that incorporating guanidinate ligands into Ce(iv) complexes gives rise to intense charge transfer bands and noteworthy electronic structures, providing insights into the stabilization of tetravalent lanthanide oxidation states.
ABSTRACT
Gas phase actinyl cation-cation interactions (CCIs) were studied by an accurate composite coupled cluster thermochemical approach for the first time. A number of CCI dimers were constructed from the monomers UO22+, UO2+, NpO22+, NpO2+, PuO2+, and AmO2+. All CCI dimers studied were calculated to be thermodynamically unstable, with dissociation energies ranging from -60 to -90 kcal mol-1, but in many cases kinetic stability was indicated by calculated local minima with well depths as large as â¼15 kcal mol-1. Most of the dimers studied involved a T-shaped geometry, although one side-on dimer, (UO2+)2, was included since it was amenable to coupled cluster methods. In the T-shaped isomers the most stable dimers were calculated to arise when the oxo-group of an An(v) actinyl cation was oriented towards the metal center of an An(vi) actinyl cation. For both mixed-valent An(vi)/An(v) and mono-valent An(v) dimers, the stability as estimated from the depth of the calculated local minimum decreased in the donor series U(v) > Np(v) > Pu(v) > Am(v). These trends correlate well with experimental trends in condensed phase CCIs. A rationale for the bonding in CCIs was investigated by carrying out charge transfer analyses using the natural bond orbital (NBO) method. Augmenting the usual Lewis acid-base explanation, CCIs are the direct result of a competition between charge transfer stabilization, which can be as much as 0.11e or 30.7 kcal mol-1 at equilibrium, and Coulombic repulsive destabilization.
ABSTRACT
New correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are reported for the actinide elements Ac and Np through Lr. These complete the series of sets reported previously for Th-U [K. A. Peterson, J. Chem. Phys. 142, 074105 (2015); M. Vasiliu et al., J. Phys. Chem. A 119, 11422 (2015)]. The new sets range in size from double- to quadruple-zeta and encompass both those optimized for valence (6s6p5f7s6d) and outer-core electron correlations (valence + 5s5p5d). The final sets have been contracted for both the DKH3 and eXact 2-component (X2C) Hamiltonians, yielding cc-pVnZ-DK3/cc-pVnZ-X2C sets for valence correlation and cc-pwCVnZ-DK3/cc-pwCVnZ-X2C sets for outer-core correlation (n = D, T, Q in each case). In order to test the effectiveness of the new basis sets, both atomic and molecular benchmark calculations have been carried out. In the first case, the first three atomic ionization potentials (IPs) of all the actinide elements Ac-Lr have been calculated using the Feller-Peterson-Dixon (FPD) composite approach, primarily with the multireference configuration interaction (MRCI) method. Excellent convergence towards the respective complete basis set (CBS) limits is achieved with the new sets, leading to good agreement with experiment, where these exist, after accurately accounting for spin-orbit effects using the 4-component Dirac-Hartree-Fock method. For a molecular test, the IP and atomization energy (AE) of PuO2 have been calculated also using the FPD method but using a coupled cluster approach with spin-orbit coupling accounted for using the 4-component MRCI. The present calculations yield an IP0 for PuO2 of 159.8 kcal/mol, which is in excellent agreement with the experimental electron transfer bracketing value of 162 ± 3 kcal/mol. Likewise, the calculated 0 K AE of 305.6 kcal/mol is in very good agreement with the currently accepted experimental value of 303.1 ± 5 kcal/mol. The ground state of PuO2 is predicted to be the Σ0g+5 state.
ABSTRACT
Gas-phase acidities and aqueous solution pKa's are predicted for MO2(OH)2, where the center atom M is a main Group 6, 16, and U atom using the Feller-Peterson-Dixon approach based on coupled cluster CCSD(T) calculations with additional corrections. The gas-phase acidities of the MO2(OH)2 compounds are essentially the same for elements (M) of the same group, 304-310 kcal/mol at 298 K. All of the Group 6 compounds are 5-6 kcal/mol less acidic in the gas phase than H2SO4. The gas-phase acidity of UO2(OH)2 is calculated to be up to 338.0 kcal/mol, â¼10% less acidic in the gas phase than the other MO2(OH)2 acids. The most acidic molecule in aqueous solution is predicted to be H2SO4. Overall, for the Group 16 compounds, the pKa's increase going down the group, with H2PoO4 predicted to be slightly more acidic than nitric acid. H2CrO4 is the most acidic of the Group 6 transition metal compounds. The aqueous acidities of H2MoO4 and H2WO4 are comparable and about 3 pKa units less acidic than H2CrO4 and comparable in acidity to HNO3. H2UO4 is not acidic at all in aqueous solution with a pKa near 20 pKa units and is also not predicted to readily undergo hydrolysis reactions.