Prediction of new inorganic molecules with quantum chemical methods

Strong influence of weak hydrogen bonding on actinide-phosphonate complexation: accurate predictions from DFT followed by experimental validation

Among the varied classes of weak hydrogen bond, the CHO type is one of immense interest as it governs the finer structures of biological and chemical molecules, hence determining their functionalities. In the present work, this weak hydrogen bond has been shown to strongly influence the complexation behaviour of uranyl nitrate [UO2(NO3)2] with diamyl-H-phosphonate (DAHP) and its branched isomer disecamyl-H-phosphonate (DsAHP). The structures of the bare ligands and complexes have been optimized by density functional theory (DFT) calculations. Surprisingly, despite having the same chemical composition the branched UO2(NO3)2·2DsAHP complex shows a remarkably higher stability (by ∼14 kcal mol-1) compared to the UO2(NO3)2·2DAHP complex. Careful inspection of the optimized structures reveals the existence of multiple CHO hydrogen-bonding interactions between the nitrate oxygens or U[double bond, length as m-dash]O oxygens and the α-hydrogens in the alkyl chains of the ligands. Comparatively stronger such bonds are found in the UO2(NO3)2·2DsAHP complex. The binding free energies associated with the complexes are computed and favoured superior binding energetics for the more stable UO2(NO3)2·2DsAHP complex. Calculations involving diisoamyl-H-phosphonate (DiAHP) and its complexes have also been performed. Theoretical predictions are experimentally tested by carrying out the extraction of U(vi) from nitric acid media using these ligands. DAHP, DsAHP and DiAHP are synthesised, characterised by NMR and evaluated for their physicochemical properties viz. viscosity, density and aqueous solubility. It was experimentally discovered that indeed DsAHP complexation with uranyl nitrate is more favoured. H-phosphonates are generically classified as acidic extractants owing to the formation of an enol tautomer at lower acidities, hence complexing the metal ion by proton exchange. Our experiments interestingly reveal a neutral ligand characteristic for DsAHP alone which is generically an acidic extractant. Furthermore, the enol tautomer of H-phosphonates that governs their extraction profiles at low acidities is also explored by DFT and the anomalous pH dependent complexation trend of DsAHP could be successfully explained. The extractions of Pu(iv) and Th(iv) have also been carried out in addition to U(vi). Solvent extraction behaviour of Am(iii) was also studied with all three ligands; the positive binding energies computed for the Am(iii) complexation corroborate with our experimental results on the poor extraction of Am(iii).

Theoretical study of the electronic spectra of neutral and cationic NpO and NpO2

The electronic spectra of neutral NpO and NpO2 as well as of their mono- (NpO(+), NpO2(+)) and dications (NpO(2+), NpO2(2+)) were studied using multiconfigurational relativistic quantum chemical calculations at the complete active space self-consistent field/CASPT2 level of theory taking into account spin-orbit coupling. The active space included 16 orbitals: all the 7s, 6d, and 5f orbitals of neptunium together with selected orbitals of oxygen. The vertical excitation energies on the ground state geometries have been computed up to ca. 35,000 cm(-1). The gas-phase electronic spectra were evaluated on the basis of the computed Einstein coefficients at 298 K and 3000 K. The computed vertical transition energies show good agreement with previous condensed-phase results on NpO2(+) and NpO2(2+).

Structure and other molecular properties of actinide trichlorides AnCl3 (An = Th-Cm)

The ground-state molecular properties of the trichlorides of light actinides (An = Th-Cm) have been predicted by state-of-the-art quantum chemical calculations. The ground electronic states have been determined by multireference calculations at the CASPT2 level including both scalar and spin-orbit relativistic effects. These studies supported the expected single-configuration character of ThCl3 and CmCl3 with their well-defined 6dσ/7s hybrid and 5f(7) configurations, respectively. In contrast, the intermediate actinides (PaCl3-AmCl3) with partly filled 5f shells have numerous very low-lying excited states and consequently a mixed character of the spin-orbit ground states. Apart from the planar ThCl3 the ground-state molecular geometries proved to be pyramidal with C(3v) symmetry. The gradually decreasing An-Cl bond distances reveal the actinide contraction known for the atomic and ionic radii of these actinide atoms. Other ground-state molecular properties as vibrational frequencies and natural charges have been obtained by density functional theory calculations using the B3LYP exchange-correlation functional in conjunction with small-core relativistic energy-consistent pseudopotentials for the actinides.

The generalized active space concept in multiconfigurational self-consistent field methods

A multiconfigurational self-consistent field method based on the concept of generalized active space (GAS) is presented. GAS wave functions are obtained by defining an arbitrary number of active spaces with arbitrary occupation constraints. By a suitable choice of the GAS spaces, numerous ineffective configurations present in a large complete active space (CAS) can be removed, while keeping the important ones in the CI space. As a consequence, the GAS self-consistent field approach retains the accuracy of the CAS self-consistent field (CASSCF) ansatz and, at the same time, can deal with larger active spaces, which would be unaffordable at the CASSCF level. Test calculations on the Gd atom, Gd(2) molecule, and oxoMn(salen) complex are presented. They show that GAS wave functions achieve the same accuracy as CAS wave functions on systems that would be prohibitive at the CAS level.

Multiconfigurational Second-Order Perturbation Theory Restricted Active Space (RASPT2) Method for Electronic Excited States: A Benchmark Study

The recently developed second-order perturbation theory restricted active space (RASPT2) method has been benchmarked versus the well-established complete active space (CASPT2) approach. Vertical excitation energies for valence and Rydberg excited states of different groups of organic (polyenes, acenes, heterocycles, azabenzenes, nucleobases, and free base porphin) and inorganic (nickel atom and copper tetrachloride dianion) molecules have been computed at the RASPT2 and multistate (MS) RASPT2 levels using different reference spaces and compared with CASPT2, CCSD, and experimental data in order to set the accuracy of the approach, which extends the applicability of multiconfigurational perturbation theory to much larger and complex systems than previously. Relevant aspects in multiconfigurational excited state quantum chemistry such as the valence-Rydberg mixing problem in organic molecules or the double d-shell effect for first-row transition metals have also been addressed.

Hydration of lanthanide chloride salts: a quantum chemical and classical molecular dynamics simulation study

We present the results of a quantum chemical and classical molecular dynamics simulation study of some solutions containing chloride salts of La(3+), Gd(3+), and Er(3+) at various concentrations (from 0.05 to 5 M), with the purpose of understanding their structure and dynamics and analyzing how the coordination varies along the lanthanide series. In the La-Cl case, nine water molecules surround the central La(3+) cation in the first solvation shell, and chloride is present only in the second shell for all solutions but the most concentrated one (5 M). In the Gd(3+) case, the coordination number is ∼8.6 for the two lowest concentrations (0.05 and 0.1 M), and then it decreases rapidly. In the Er(3+) case, the coordination number is 7.4 for the two lowest concentrations (0.05 and 0.1 M), and then it decreases. The counterion Cl(-) is not present in the first solvation shell in the La(3+) case for most of the solutions, but it becomes progressively closer to the central cation in the Gd(3+) and Er(3+) cases, even at low concentrations.

Generating Cu(II)-oxyl/Cu(III)-oxo species from Cu(I)-alpha-ketocarboxylate complexes and O2: in silico studies on ligand effects and C-H-activation reactivity

A mechanism for the oxygenation of Cu(I) complexes with alpha-ketocarboxylate ligands that is based on a combination of density functional theory and multireference second-order perturbation theory (CASSCF/CASPT2) calculations is elaborated. The reaction proceeds in a manner largely analogous to those of similar Fe(II)-alpha-ketocarboxylate systems, that is, by initial attack of a coordinated oxygen molecule on a ketocarboxylate ligand with concomitant decarboxylation. Subsequently, two reactive intermediates may be generated, a Cu-peracid structure and a [CuO](+) species, both of which are capable of oxidizing a phenyl ring component of the supporting ligand. Hydroxylation by the [CuO](+) species is predicted to proceed with a smaller activation free energy. The effects of electronic and steric variations on the oxygenation mechanisms were studied by introducing substituents at several positions of the ligand backbone and by investigating various N-donor ligands. In general, more electron donation by the N-donor ligand leads to increased stabilization of the more Cu(II)/Cu(III)-like intermediates (oxygen adducts and [CuO](+) species) relative to the more Cu(I)-like peracid intermediate. For all ligands investigated, the [CuO](+) intermediates are best described as Cu(II)-O(*-) species with triplet ground states. The reactivity of these compounds in C-H abstraction reactions decreases with more electron-donating N-donor ligands, which also increase the Cu-O bond strength, although the Cu-O bond is generally predicted to be rather weak (with a bond order of about 0.5). A comparison of several methods to obtain singlet energies for the reaction intermediates indicates that multireference second-order perturbation theory is likely more accurate for the initial oxygen adducts, but not necessarily for subsequent reaction intermediates.

The restricted active space followed by second-order perturbation theory method: theory and application to the study of CuO2 and Cu2O2 systems

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.

Experimental and theoretical evidence for U(C6H6) and Th(C6H6) complexes

The vibrational spectra of UBz and ThBz have been measured in solid argon. Complementary quantum chemical calculations have allowed the assignments of the vibrational spectra. According to the calculations, AcBz are stable molecules, as well as other species like BzAcBz and BzAc2Bz. Experimentally, there is no evidence for the sandwich compounds BzAcBz and BzAc2Bz due to the limitations in the reagent concentrations.

Quantum Chemical Characterization of the Bonding of N-Heterocyclic Carbenes to Cp2MI Compounds [M = Ce(III), U(III)]

The binding of N-heterocyclic carbenes to Ce(III) and U(III) compounds is characterized by quantum chemical methods. Density functional methods are in qualitative agreement with experiment that binding to U(III) is more favorable than to Ce(III); after correcting for basis-set superposition error, quantitative agreement with experiment is achieved with a multireference second-order perturbation theory approach accounting for relativistic effects. The small computed (and observed) preference derives from a combination of several small effects, including differences in electronic binding energies, rovibrational partition functions, and solvation free energies. Prospects for ligand modification to improve the differentiation between lanthanides and actinides are discussed on the basis of computational predictions.