Search

Multiply-bonded main group compounds of groups 13-15 are attracting significant interest not only because they provide fundamental insight into the nature of metal-metal bonding, but also for their potential in small molecule bond activation and catalysis. This includes dialumenes, neutral Al(I) compounds that contain AlAl double bonds, which display high reactivity owing to their intrinsic diradicaloid character. The electronic structure of the simplest dialumene, Al 2 H 2 , is here analyzed up to a practical Full-CI limit using DMRG and selected CI methods for the bond dissociation energy (BDE), geometry and properties of the electron density (difference density, ELF). Acquiring Full-CI reference values for the simplest dialumene (but possessing the highest diradical character) allows for a rigorous benchmarking of simpler correlated wavefunction theory (WFT) methods and density functional methods in treating the electronic structure of such systems. Single-reference coupled cluster theory using a RHF reference is found to reliably converge to the Full-CI limit and CCSD(T) is fully capable of capturing the diradical character, while multi-reference methods offer no clear advantages. Density functional methods struggle to fully describe the electronic structure complexity although non-hybrid functionals such as TPSS come close. Solving the inverse Kohn-Sham problem for a Full-CI-quality density revealed minimal density-driven errors in the TPSS-calculated BDE unlike high-percentage hybrids such as M06-2X. No density functional, however, predicts accurate relative energies. Fractional-occupation density plots at the TPSS level correlate well with WFT-based diradical character metrics, a useful result for determining diradical character in larger systems.

Hydrides are present in the reduced states of the iron-molybdenum cofactor (FeMoco) of Mo nitrogenase and are believed to play a key mechanistic role in the dinitrogen reduction reaction catalyzed by the enzyme. Two hydrides are present in the E4 state according to 1 H ENDOR and there is likely a single hydride in the E2 redox state. The 2-hydride E4 state has been experimentally observed to bind N 2 and it has been speculated that E3 may bind N 2 as well. However, the E3 state has not been directly observed and very little is known about its molecular and electronic structure or reactivity. In recent computational studies, we have explored the energy surfaces of the E2 and E4 by QM/MM modelling, and found that the most stable hydride isomers contain bridging or partially bridging hydrides with an open protonated belt sulfide-bridge. In this work we systematically explore the energy surface of the E3 redox state, comparing single hydride and two-hydride isomers with varying coordination and bridging vs. terminal sulfhydryl groups. We also include a model featuring a triply protonated carbide. The results are only mildly dependent on the QM-region size. The three most stable E3 isomers at the r 2 SCAN level of theory have in common: an open belt sulfide-bridge (terminal sulfhydryl group on Fe6) and either 2 bridging hydrides (between Fe2 and Fe6), 1 bridging-1-terminal hydride (around Fe2 and Fe6) or a dihydrogen ligand bound at the Fe2 site. Analyzing the functional dependency of the results, we find that functionals previously found to predict accurate structures of spin-coupled Fe/Mo dimers and FeMoco (TPSSh, B97-D3, r 2 SCAN, and B3LYP*) are in generally good agreement about the stability of these 3 E3 isomers. However, B3LYP*, similar to its parent B3LYP method, predicts a triply protonated carbide isomer as the most stable isomer, an unlikely scenario in view of the lack of experimental evidence for carbide protonation occurring in reduced FeMoco states. Distinguishing further between the 3 hydride isomers is difficult and this flexible coordination nature of hydrides suggests that multiple hydride isomers could be present during experimental conditions. N 2 binding was explored and resulted in geometries with 2 bridging hydrides and N 2 bound to either Fe 2 or Fe 6 with a local low-spin state on the Fe. N 2 binding is predicted to be mildly endothermic, similar to the E2 state, and it seems unlikely that the E3 state is capable of binding N 2 .

Open-shell transition metal chemistry presents challenges to contemporary electronic structure methods, based on either density functional or wavefunction theory. While CCSD(T) is the well-trusted gold standard for maingroup thermochemistry, the accuracy and robustness of the method is less clear for open-shell transition metal chemistry, requiring benchmarking of CCSD(T)-based protocols against either higher-level theory or experiment. Ionization energies (IEs) of metallocenes provide an interesting test case with metallocenes being common redox reagents as well as playing roles as redox mediators and cocatalysts in redox catalysis. Using highly accurate ZEKE-MATI experimental measurements of gas phase adiabatic (5.3275 ± 0.0006 eV) and vertical (5.4424 ± 0.0006 eV) ionization energies of cobaltocene, we systematically assessed the accuracy of the local coupled-cluster method DLPNO-CCSD(T) with respect to geometry, reference determinant, basis set size and extrapolation schemes, PNO cut-off and extrapolation, local triples approximation, relativistic effects and core-valence correlation. We show that PNO errors are controllable via the recently introduced PNO extrapolation schemes and that the expensive iterative triples (T 1 ) contribution can be made more manageable by calculating it as a smaller-basis/smaller PNO-cutoff correction. The reference determinant turns out to be a critical aspect in these calculations with the HF determinant resulting in large DLPNO-CCSD(T) errors, likely due to the qualitatively flawed molecular orbital spectrum. The BP86 functional on the other hand was found to provide reference orbitals giving small DLPNO-CCSD(T) errors, likely due to more realistic orbitals as suggested by the more consistent MO spectrum compared to HF. A protocol including complete basis set extrapolations with correlation-consistent basis sets, complete PNO space extrapolations, iterative triples- and core-valence correlation corrections was found to give errors of -0.07 eV and -0.03 eV for adiabatic- and vertical-IE of cobaltocene, respectively, giving close to chemical accuracy for both properties. A computationally efficient DLPNO-CCSD(T) protocol was devised and tested against adiabatic ionization energies of 6 different metallocenes (V, Cr, Mn, Fe, Co, Ni). For the other metallocenes, the iterative triples (T 1 ) and PNO extrapolation contributions turn out to be even more important. The results give errors close to the experimental uncertainty, similar to recent auxiliary-field quantum Monte Carlo results. The quality of the reference determinant orbitals is identified as the main source of uncertainty in CCSD(T) calculations of metallocenes.

In this study, we present experimental and theoretical results on dissociative electron attachment and dissociative ionisation for the potential FEBID precursor cis-Pt(CO) 2 Cl 2 . UHV surface studies have shown that high purity platinum deposits can be obtained from cis-Pt(CO) 2 Cl 2 . The efficiency and energetics of ligand removal through these processes are discussed and experimental appearance energies are compared to calculated thermochemical thresholds. The present results demonstrate the potential effectiveness of electron-induced reactions in the deposition of this FEBID precursor, and these are discussed in conjunction with surface science studies on this precursor and the design of new FEBID precursors.

We present a combined theoretical and experimental study on the ionization and primary fragmentation channels of the mono-halogenated biphenyls; 2-chlorobiphenyl, 2-bromobiphenyl and 2-iodobiphenyl. The ionization energies (IEs) of the 2-halobiphenyls and the appearance energies (AEs) of the principal fragments are determined through electron impact ionization, while quantum mechanical calculations at the coupled cluster level of theory are used to elucidate the observed processes and the associated dynamics. The primary fragmentation channels are the direct loss of the halogen upon ionization, the loss of the respective hydrogen halides (HX) as well as loss of the hydrogen halide and an additional hydrogen. We find that the dihedral angle strongly influences the relative potential energy of the neutral and the cation on their respective ground state surfaces, an effect caused by the strong influence of the nuclear motion on the conjugation between the phenyl rings. For the principal dissociative ionization channels from the mono-halogenated biphenyls we reason that these can not be described as statistical decay from the ground state cation, but must rather be understood as direct, state-selective processes from specific excited cationic states characterized through local ionization of either the halogenated or the non-substituted phenyl ring.

The production of alloyed nanostructures presents a unique problem in focused electron beam induced deposition (FEBID). Deposition of such structures has historically involved the mixing of two or more precursor gases in situ or via multiple channel gas injection systems, thereby making the production of precise, reproducible alloy compositions difficult. Promising recent efforts to address this problem have involved the use of multi-centred, heterometallic FEBID precursor species. In this vein, we present here a study of low-energy electron interactions with cyclopentadienyl iron dicarbonyl manganese pentacarbonyl ((η 5 -Cp)Fe(CO) 2 Mn(CO) 5 ), a bimetallic species with a polyhapto ligand (Cp) and seven terminal carbonyl ligands. Gas phase studies and coupled cluster calculations of observed low-energy electron-induced reactions were conducted in order to predict the performance of this precursor in FEBID. In dissociative electron attachment, we find single CO loss and cleavage of the Fe-Mn bond, leading to the formation of [Mn(CO) 5 ] - , to be the two dominant channels. Contributions through further CO loss from the intact core and the formation of [Mn(CO) 4 ] - are minor channels. In dissociative ionization (DI), the fragmentation is significantly more extensive and the DI spectra are dominated by fragments formed through the loss of 5 and 6 CO ligands, and fragments formed through cleavage of the Fe-Mn bond accompanied by substantial CO loss. The gas phase fragmentation channels observed are discussed in relation to the underlying processes and their energetics, and in context to related surface studies and the likely performance of this precursor in FEBID.

X-ray absorption spectroscopy (XAS) is a widely used experimental technique capable of selectively probing the local structure around an absorbing atomic species in molecules and materials. When applied to heavy elements, however, the quantitative interpretation can be challenging due to the intrinsic spectral broadening arising from the decrease in the core-hole lifetime. In this work we have used high-energy resolution fluorescence detected XAS (HERFD-XAS) to investigate a series of molybdenum complexes. The sharper spectral features obtained by HERFD-XAS measurements enable a clear assignment of the features present in the pre-edge region. Time-dependent density functional theory (TDDFT) has been previously shown to predict K-pre-edge XAS spectra of first row transition metal compounds with a reasonable degree of accuracy. Here we extend this approach to molybdenum K-edge HERFD-XAS and present the necessary calibration. Modern pure and hybrid functionals are utilized and relativistic effects are accounted for using either the Zeroth Order Regular Approximation (ZORA) or the second order Douglas-Kroll-Hess (DKH2) scalar relativistic approximations. We have found that both the predicted energies and intensities are in excellent agreement with experiment, independent of the functional used. The model chosen to account for relativistic effects also has little impact on the calculated spectra. This study provides an important calibration set for future applications of molybdenum HERFD-XAS to complex catalytic systems.

(51)V NMR parameters have been calculated for VOCl(3), the reference compound in (51)V NMR spectroscopy, in order to capture environmental effects in both the neat liquid and the solid state. Using a combination of periodic geometry optimizations and Car-Parrinello molecular dynamics simulations with embedded cluster NMR calculations, we are able to test the ability of current computational approaches to reproduce (51)V NMR properties (isotropic shifts, anisotropic shifts and quadrupole coupling constants) in the gas, liquid and solid states, for direct comparison with liquid and solid-state experimental data. The results suggest that environmental effects in the condensed phases can be well captured by an embedded cluster approach and that the remaining discrepancy with experiment may be due to the approximate density functionals in current use. The predicted gas-to-liquid shift on the isotropic shielding constant is small, validating the common practice to use a single VOCl(3), molecule as reference in (51)V NMR computations.

Density functional methods were evaluated in their ability to predict relative conformational energies of a test set of monosubstituted cyclohexanes and six-membered heterocycles. It is shown that while popular density functionals like B3LYP are unreliable for predicting accurate conformational energies for the axial/equatorial equilibrium of monosubstituted cyclohexanes, 1-silacyclohexanes and tetrahydropyrans, density functionals that take into account dispersion interactions like M06-2X and B2PLYP-D result in energy differences close to CCSD(T)/CBS results. Using the M06-2X density functional, we have then investigated the conformational properties of a large number of monosubstituted silacyclohexanes, with the number of silicon atoms ranging from 1 to 6. Our calculations suggest remarkably different conformational properties when compared to cyclohexane. The carbon/silicon exchange in a cyclohexane ring often has systematic, yet counterintuitive effects on the conformational properties. Dispersion interactions are shown to be especially important for accurate relative energy calculations of polysilacyclohexanes.