Probing the Isolobal Relation between Cp'''NiP3 and White Phosphorus by Experimental Charge Density Analysis.

An in-depth analysis of the description of bonding within Cp'''Ni-cyclo-P3 (Cp'''=1,2,4-tri-tert-butylcyclopentadienyl, [Ni]P3) employing X-ray diffraction based multipolar modeling, density functional theory (DFT) as well as an "experimental wavefunction" obtained from X-ray restrained wavefunction (XRW) fitting is presented. The results are compared to DFT calculations on white phosphorus - an isolobal analogue to [Ni]P3. A complementary bonding analysis shows insights into the reactivity of [Ni]P3. The isolobal principle is reflected in every aspect of our analysis and the employed methods seamlessly predict the differences in reactivity of [Ni]P3 and P4. Crystallographic modeling, solid-state NMR, and DFT calculations describe the dynamic behavior of the cyclo-P3 unit in the title molecule.

Reversing the Bond Length Alternation Order in Conjugated Polyenes by Substituent Effects.

We have synthesised several push-pull substituted conjugated polyenes and determined their accurate C-C bond lengths and charge-density distributions by utilising quantum crystallographic techniques. In a series of alkene, dienes, and triene bearing two (trifluoromethyl)sulfonyl (triflyl) groups on the terminal carbon atom, unique reversal of the bond-length alternation (BLA) order has been observed. This is a pronounced aberration from the molecular structure predicted by the Lewis-structure-based neutral resonance structure. Such reversal of BLA order has not been observed in push-pull compounds bearing conventional electron-withdrawing groups such as carbonyl and cyano groups instead of triflyl groups. Bonding behaviour of both normal and reversed bond length alternating systems has been revealed by complementary bonding analysis using several bond descriptors based on the experimentally fitted wavefunctions.

The effects of experimentally obtained electron correlation and polarization on electron densities and exchange-correlation potentials.

In X-ray constrained wavefunction (XCW) fitting, external information, such as electron correlation and polarization, is included into a single-determinantal isolated-molecule wavefunction. In a first step, we show that the extraction of these two physical effects by XCW fitting is complete and accurate by comparing to theoretical reference calculations. In a second step, we show that fitting to data from single-crystal x-ray diffraction measurements provides the same results qualitatively and how the physical effects can be separated, although always inherently convolved in the experiment. We further demonstrate that exchange-correlation potentials are systematically affected by XCW fitting in a physically meaningful way, which could be exploited for method development in quantum chemistry, subject to some remaining challenges that we also outline.

Aurophilic Interactions Studied by Quantum Crystallography.

This is the first use of a wave-function-based crystallographic method to characterize aurophilic interactions from X-ray diffraction data. Theoretical calculations previously suggested the importance of electron correlation and dispersion forces, but no influence of relativistic corrections to the Au...Au interaction energy was found. In this study, we confirm the importance of relativistic corrections in the characterization of aurophilic interactions in addition to electron correlation and dispersion.

Similarities and Differences between Crystal and Enzyme Environmental Effects on the Electron Density of Drug Molecules.

The crystal interaction density is generally assumed to be a suitable measure of the polarization of a low-molecular weight ligand inside an enzyme, but this approximation has seldomly been tested and has never been quantified before. In this study, we compare the crystal interaction density and the interaction electrostatic potential for a model compound of loxistatin acid (E64c) with those inside cathepsin B, in solution, and in vacuum. We apply QM/MM calculations and experimental quantum crystallography to show that the crystal interaction density is indeed very similar to the enzyme interaction density. Less than 0.1 e are shifted between these two environments in total. However, this difference has non-negligible consequences for derived properties.

On the Role of Hydrogen Bonding in Gas-Phase SN2 Reactions at Silicon.

The shape of the potential energy surface (PES) of gas-phase SN2 reactions at silicon is determined by the type of nucleophile, the leaving group, and substituents which remain bonded to silicon. In this study, we present PES scans along the reaction coordinate of six symmetrical SN2 reactions: X- + SiR3X → XSiR3 + X-, where X = Cl or F and R = H, Me, or OMe. While the fluorine systems and the ClSiH3Cl system only give single-well PESs, ClSiMe3Cl and ClSi(OMe)3Cl give triple- and double-well PESs with stable pre- and post-reaction complexes. A complementary bonding analysis (energy decomposition analysis, quantum theory of atoms in molecules, and natural bond orbitals) reveals that the leaving group (X-) is stabilized by hydrogen bonding in the XSiMe3X and XSi(OMe)3X systems. It is shown that this so far neglected stabilizing contribution, along with σ-hole bonding, is responsible for the shapes of the PESs of ClSiMe3Cl and ClSi(OMe)3Cl in the gas phase.

A Small Cationic Organo-Copper Cluster as Thermally Robust Highly Photo- and Electroluminescent Material.

Organic light-emitting diodes (OLEDs) are revolutionizing display applications. In this aspect, luminescent complexes of precious metals such as iridium, platinum, or ruthenium still playing a significant role. Emissive compounds of earth-abundant copper with equivalent performance are desired for practical, large-scale applications such as solid-state lighting and displays. Copper(I)-based emitters are well-known to suffer from weak spin-orbit coupling and a high reorganization energy upon photoexcitation. Here we report a cationic organo-copper cluster [Cu4(PCP)3]+ (PCP = 2,6-(PPh2)2C6H3) that features suppressed nonradiative decays, giving rise to a robust narrow-band green luminophore with a photoluminescent (PL) efficiency up to 93%. PL decay kinetics corroborated by DFT calculations reveal a complex emission mechanism involving contributions of both thermally activated delayed fluorescence and phosphorescence. This robust compound was solution-processed into a thin film in prototype OLEDs with external quantum efficiency up to 11% and a narrow emission bandwidth (65 nm fwhm).

Non-Oxido-Vanadium(IV) Metalloradical Complexes with Bidentate 1,2-Dithienylethene Ligands: Observation of Reversible Cyclization of the Ligand Scaffold in Solution.

Derivatives of 1,2-dithienylethene (DTE) have superb photochromic properties due to an efficient reversible photocyclization reaction of their hexatriene structure and, thus, have application potential in materials for optoelectronics and (multi-responsive) molecular switches. Transition-metal complexes bearing switchable DTE motifs commonly incorporate their coordination site rather distant from the hexatriene system. In this work the redox active ligand 1,2-bis(2,5-dimethylthiophen-3-yl)ethane-1,2-dione is described, which reacts with [V(TMEDA)2 Cl2 ] to give a rare non-oxido vanadium(IV) species 3(M,M/P,P). This blue complex has two bidentate en-diolato ligands which chelate the VIV center and give rise to two five-membered metallacycles with the adjacent hexatriene DTE backbone bearing axial chirality. Upon irradiation with UVA light or prolonged heating in solution, the blue compound 3(M,M/P,P) converts into the purple atropisomer 4(para,M/para,P). Both complexes were isolated and structurally characterized by single-crystal X-ray diffraction analysis (using lab source and synchrotron radiation). The antiparallel configuration (M or P helicity) present in both 3(M,M/P,P) and 4(para,M/para,P) is a prerequisite for (reversible) 6π cyclization reactions. A CW EPR spectroscopic study reveals the metalloradical character for 3(M,M/P,P) and 4(para,M/para,P) and indicates dynamic reversible cyclization of the DTE backbone in complex 3(M,M/P,P) at ambient temperature in solution.

Proximity Enforced Agostic Interactions Involving Closed-Shell Coinage Metal Ions.

A proximity enforcing diarylsilane ligand is reported, which gives rise to unusual Si-H···M interactions with the d10 metal ions Cu+ and Ag+ upon complexation. These interactions are studied in detail both experimentally and computationally and can be classified to be weakly agostic in nature for the Si-H···Cu interaction. The Si-H···Ag interaction has more signatures of an electrostatic contact.

Chemical Bonding in Polarised Push-Pull Ethylenes.

1,1-Diamino-2,2-bis(triflyl)ethylenes with both twisted and planar structures around the partial "C=C" bond were synthesised. Bonding properties in these compounds were analysed by an experimental approach using high-resolution X-ray diffraction data treated with X-ray wavefunction refinement (XWR). In the twisted compound, a dominant contribution of the charge-separated resonance structure was revealed. On the contrary, the nearly planar compound still showed π-bonding character, however, with a considerable contribution of the charge-separated resonance structure.

Nature of Bonding in Donor-Acceptor Interactions Exemplified by Complexes of N-Heterocyclic Carbenes with 1,2,5-Telluradiazoles.

Comprehensive structural, spectroscopic, and quantum chemical analyses of new donor-acceptor complexes between N-heterocyclic carbenes and 1,2,5-telluradiazoles and a comparison with previously known complexes involving tellurenyl cations showed that the dative C-Te bonds cannot be solitarily described with only one Lewis formula. Canonical Lewis formulas that denote covalency and arrows emphasizing ionicity complement each other in varying extents. Evaluation of the relative weights of these resonance forms requires proper bonding description with a well-balanced toolbox of analytical methods. If for conciseness only, one resonance form is used, it must be the most significant one according to the analytical evaluation. If unclear, all significant resonance forms should be displayed.

Superelectrophilic Behavior of an Anion Demonstrated by the Spontaneous Binding of Noble Gases to [B12 Cl11 ].

It is common and chemically intuitive to assign cations electrophilic and anions nucleophilic reactivity, respectively. Herein, we demonstrate a striking violation of this concept: The anion [B12 Cl11 ]- spontaneously binds to the noble gases (Ngs) xenon and krypton at room temperature in a reaction that is typical of "superelectrophilic" dications. [B12 Cl11 Ng]- adducts, with Ng binding energies of 80 to 100 kJ mol-1 , contain B-Ng bonds with a substantial degree of covalent interaction. The electrophilic nature of the [B12 Cl11 ]- anion is confirmed spectroscopically by the observation of a blue shift of the CO stretching mode in the IR spectrum of [B12 Cl11 CO]- and theoretically by investigation of its electronic structure. The orientation of the electric field at the reactive site of [B12 Cl11 ]- results in an energy barrier for the approach of polar molecules and facilitates the formation of Ng adducts that are not detected with reactive cations such as [C6 H5 ]+ . This introduces the new chemical concept of "dipole-discriminating electrophilic anions."

Refinement of X-ray and electron diffraction crystal structures using analytical Fourier transforms of Slater-type atomic wavefunctions in Olex2.

An implementation of Slater-type spherical scattering factors for X-ray and electron diffraction for elements in the range Z = 1-103 is presented within the software Olex2. Both high- and low-angle Fourier behaviour of atomic electron density and electrostatic potential can thus be addressed, in contrast to the limited flexibility of the four Gaussian plus constant descriptions which are currently the most widely used method for calculating atomic scattering factors during refinement. The implementation presented here accommodates the increasing complexity of the electronic structure of heavier elements by using complete atomic wavefunctions without any interpolation between precalculated tables or intermediate fitting functions. Atomic wavefunctions for singly charged ions are implemented and made accessible, and these show drastic changes in electron diffraction scattering factors compared with the neutral atom. A comparison between the two different spherical models of neutral atoms is presented as an example for four different kinds of X-ray and two electron diffraction structures, and comparisons of refinement results using the existing diffraction data are discussed. A systematic but slight improvement in R values and residual densities can be observed when using the new scattering factors, and this is discussed relative to effects on the atomic displacement parameters and atomic positions, which are prominent near the heavier elements in a structure.

A tetrahedral neptunium(V) complex.

Neptunium is an actinide element sourced from anthropogenic production, and, unlike naturally abundant uranium, its coordination chemistry is not well developed in all accessible oxidation states. High-valent neptunium generally requires stabilization from at least one metal-ligand multiple bond, and departing from this structural motif poses a considerable challenge. Here we report a tetrahedral molecular neptunium(V) complex ([Np5+(NPC)4][B(ArF5)4], 1-Np) (NPC = [NPtBu(pyrr)2]-; tBu = C(CH3)3; pyrr = pyrrolidinyl (N(C2H4)2); B(ArF5)4 = tetrakis(2,3,4,5,6-pentafluourophenyl)borate). Single-crystal X-ray diffraction, solution-state spectroscopy and density functional theory studies of 1-Np and the product of its proton-coupled electron transfer (PCET) reaction, 2-Np, demonstrate the unique bonding that stabilizes this reactive ion and establishes the thermochemical and kinetic parameters of PCET in a condensed-phase transuranic complex. The isolation of this four-coordinate, neptunium(V) complex reveals a fundamental reaction pathway in transuranic chemistry.

Aspherical atom refinements on X-ray data of diverse structures including disordered and covalent organic framework systems: a time-accuracy trade-off.

Aspherical atom refinement is the key to achieving accurate structure models, displacement parameters, hydrogen-bond lengths and analysis of weak interactions, amongst other examples. There are various quantum crystallographic methods to perform aspherical atom refinement, including Hirshfeld atom refinement (HAR) and transferable aspherical atom model (TAAM) refinement. Both HAR and TAAM have their limitations and advantages, the former being more accurate and the latter being faster. With the advent of non-spherical atoms in Olex2 (NoSpherA2), it is now possible to overcome some limitations, like treating disorder, twinning and network structures, in aspherical refinements using HAR, TAAM or both together. TAAM refinement in NoSpherA2 showed significant improvement in refinement statistics compared with independent atom model (IAM) refinements on a diverse set of X-ray diffraction data. The sensitivity of TAAM towards poor data quality and disorder was observed in terms of higher refinement statistics for such structures. A comparison of IAM with TAAM and HAR in NoSpherA2 indicated that the time taken by TAAM refinements was of the same order of magnitude as that taken by IAM, while in HAR the time taken using a minimal basis set was 50 times higher than for IAM and rapidly increased with increasing size of the basis sets used. The displacement parameters for hydrogen and non-hydrogen atoms were very similar in both HAR and TAAM refinements. The hydrogen-bond lengths were slightly closer to neutron reference values in the case of HAR with higher basis sets than in TAAM. To benefit from the advantages of each method, a new hybrid refinement approach has been introduced, allowing a combination of IAM, HAR and TAAM in one structure refinement. Refinement of coordination complexes involving metal-organic compounds and network structures such as covalent organic frameworks and metal-organic frameworks is now possible in a hybrid mode such as IAM-TAAM or HAR-TAAM, where the metal atoms are treated via either the IAM or HAR method and the organic part via TAAM, thus reducing the computational costs without compromising the accuracy. Formal charges on the metal and ligand can also be introduced in hybrid-mode refinement.

An investigation into the Brønsted acidity of the perfluorinated alkoxy silanes {(F3C)3CO}3SiH and {(F6C5)3CO}2Si(Cl)H.

The perfluorinated alkoxy silanes {(F3C)3CO}3SiH (1) and {(F5C6)3CO}2Si(Cl)H (2) were prepared and fully characterized. Despite the high calculated Brønsted acidities, all attempts to deprotonate 1 and 2 to give the conjugate silanide ions failed due to the exceptionally short and strong Si-H bonds. In the solid state, the Si-H units are not involved in any intermolecular interactions, but instead the crystal packing consists of exceptionally short and strong F⋯F interactions. The cohesive energies are entirely comprised of London dispersion interactions, similarly as in the crystal structures of noble gases.

Influence of modelling disorder on Hirshfeld atom refinement results of an organo-gold(I) compound.

Details of the validation of disorder modelling with Hirshfeld atom refinement (HAR) for a previously investigated organo-gold(I) compound are presented here. The impact of refining disorder on HAR results is discussed using an analysis of the differences of dynamic structure factors. These dynamic structure factor differences are calculated from thermally smeared quantum mechanical electron densities based on wavefunctions that include or exclude electron correlation and relativistic effects. When disorder is modelled, the electron densities stem from a weighted superposition of two (or more) different conformers. Here this is shown to impact the relative importance of electron correlation and relativistic effect estimates expressed by the structure factor magnitudes. The role of disorder modelling is also compared with the effect of the treatment of hydrogen anisotropic displacement parameter (ADP) values and atomic anharmonicity of the gold atom. The analysis of ADP values of gold and disordered carbon atoms showed that the effect of disorder significantly altered carbon ADP values and did not influence those of the gold atom.

Refinement of anomalous dispersion correction parameters in single-crystal structure determinations.

Correcting for anomalous dispersion is part of any refinement of an X-ray dif-fraction crystal structure determination. The procedure takes the inelastic scattering in the diffraction experiment into account. This X-ray absorption effect is specific to each chemical compound and is particularly sensitive to radiation energies in the region of the absorption edges of the elements in the compound. Therefore, the widely used tabulated values for these corrections can only be approximations as they are based on calculations for isolated atoms. Features of the unique spatial and electronic environment that are directly related to the anomalous dispersion are ignored, although these can be observed spectroscopically. This significantly affects the fit between the crystallographic model and the measured intensities when the excitation wavelength in an X-ray diffraction experiment is close to an element's absorption edge. Herein, we report on synchrotron multi-wavelength single-crystal X-ray diffraction, as well as X-ray absorption spectroscopy experiments which we performed on the mol-ecular compound Mo(CO)6 at energies around the molybdenum K edge. The dispersive (f') and absorptive (f'') terms of the anomalous dispersion can be refined as independent parameters in the full-matrix least-squares refinement. This procedure has been implemented as a new feature in the well-established OLEX2 software suite. These refined parameters are in good agreement with the independently recorded X-ray absorption spectrum. The resulting crystallographic models show significant improvement compared to those employing tabulated values.

Relativistic Hirshfeld atom refinement of an organo-gold(I) compound.

The main goal of this study is the validation of relativistic Hirshfeld atom refinement (HAR) as implemented in Tonto for high-resolution X-ray diffraction datasets of an organo-gold(I) compound. The influence of the relativistic effects on statistical parameters, geometries and electron density properties was analyzed and compared with the influence of electron correlation and anharmonic atomic motions. Recent work in this field has indicated the importance of relativistic effects in the static electron density distribution of organo-mercury compounds. This study confirms that differences in electron density due to relativistic effects are also of significant magnitude for organo-gold compounds. Relativistic effects dominate not only the core region of the gold atom, but also influence the electron density in the valence and bonding region, which has measurable consequences for the HAR refinement model parameters. To study the effects of anharmonic motion on the electron density distribution, dynamic electron density difference maps were constructed. Unlike relativistic and electron correlation effects, the effects of anharmonic nuclear motion are mostly observed in the core area of the gold atom.

Vanishing of the atomic form factor derivatives in non-spherical structural refinement - a key approximation scrutinized in the case of Hirshfeld atom refinement.

When calculating derivatives of structure factors, there is one particular term (the derivatives of the atomic form factors) that will always be zero in the case of tabulated spherical atomic form factors. What happens if the form factors are non-spherical? The assumption that this particular term is very close to zero is generally made in non-spherical refinements (for example, implementations of Hirshfeld atom refinement or transferable aspherical atom models), unless the form factors are refinable parameters (for example multipole modelling). To evaluate this general approximation for one specific method, a numerical differentiation was implemented within the NoSpherA2 framework to calculate the derivatives of the structure factors in a Hirshfeld atom refinement directly as accurately as possible, thus bypassing the approximation altogether. Comparing wR2 factors and atomic parameters, along with their uncertainties from the approximate and numerically differentiating refinements, it turns out that the impact of this approximation on the final crystallographic model is indeed negligible.