Comprehensive Benchmarking of a Density-Dependent Dispersion Correction

Active catalyst construction for CO2 recycling via catalytic synthesis of N-doped carbon on supported Cu

Bridging homogeneous and heterogeneous catalysis is a long-term pursuit in the field of catalysis. Herein, we report our results in integration of nano- and molecular catalysis via catalytic synthesis of nitrogen doped carbon layers on AlOx supported nano-Cu which can finely tune the catalytic performance of the supported copper catalyst. This synthetic catalytic material, which can be generated in situ by the reaction of CuAlOx and 1,10-Phen in the presence of hydrogen, could be used for controllable synthesis of N,N-dimethylformamide (DMF) from dimethylamine and CO₂/H₂ via blocking reaction pathways of further catalytic hydrogenation of DMF to N(CH₃)₃. Detailed characterizations and DFT calculations reveal that the presence of N-doped layered carbon on the surface of the nano-Cu particles results in higher activation energy barriers during the conversion of DMF to N(CH₃)₃. Our primary results could promote merging of homogeneous catalysis and heterogeneous catalysis and CO₂ recycling.

Bridging homogeneous and heterogeneous catalysis is a long-term pursuit in the field of catalysis. Here, the authors present results on integration of nano- and molecular catalysis via catalytic synthesis of nitrogen doped carbon layers on AlO_x supported nano-Cu which can finely tune the catalytic performance of the supported copper catalyst.

Can Density Functional Theory Be Trusted for High-Order Electric Properties? The Case of Hydrogen-Bonded Complexes

This work reports on an extensive assessment of the performance of a wide palette of density functional approximations in predicting the (high-order) electric properties of hydrogen-bonded complexes. To this end, we compute the electronic and vibrational contributions to the electric polarizability and the first and second hyperpolarizabilities, using the CCSD(T)/aug-cc-pVTZ level of theory as reference. For all the studied properties, the average absolute errors below 20% can only be obtained using the CAM-B3LYP functional, while LC-BLYP and MN15 are shown to be only slightly less accurate (average absolute errors not exceeding 30%). Among Minnesota density functionals, i.e., M06, M06-2X, and MN15, we only recommend the latter one, which quite accurately predicts the electronic and vibrational (hyper)polarizabilities. We also analyze the optimal tuning of the range-separation parameter μ for the LC-BLYP functional, finding that this approach does not bring any systematic improvement in the predictions of electronic and vibrational (hyper)polarizabilities and the accuracy of computed properties is largely system-dependent. Finally, we report huge errors in predicting the vibrational second hyperpolarizability by ωB97X, M06, and M06-2X functionals. Based on the explicit evaluation of anharmonic terms contributing to the second hyperpolarizability, this failure is traced down to a poor determination of third- and fourth-order energy derivatives with respect to normal modes. These results reveal serious flaws of some density functional approximations and suggest caution in selecting the appropriate functional to calculate not only electronic and vibrational (hyper)polarizabilities but also other molecular properties that contain vibrational anharmonic contributions.

Comparing the Self-Assembly of Sexiphenyl-Dicarbonitrile on Graphite and Graphene on Cu(111)

A comparative study on the self-assembly of sexiphenyl-dicarbonitrile on highly oriented pyrolytic graphite and single-layer graphene on Cu(111) is presented. Despite an overall low molecule-substrate interaction, the close-packed structures exhibit a peculiar shift repeating every four to five molecules. This shift has hitherto not been reported for similar systems and is hence a unique feature induced by the graphitic substrates.

Water Splitting Reaction at Polar Lithium Niobate Surfaces

Water splitting is a highly promising, environmentally friendly approach for hydrogen production. It is often discussed in the context of carbon dioxide free combustion and storage of electrical energy after conversion to chemical energy. Since the oxidation and reduction reactions are related to significant overpotentials, the search for suitable catalysts is of particular importance. Ferroelectric materials, for example, lithium niobate, attracted considerable interest in this respect. Indeed, the presence of surfaces with different polarizations and chemistries leads to spatial separation of reduction and oxidation reactions, which are expected to be boosted by the electrons and holes available at the positive and negative surfaces, respectively. Employing the density functional theory and a simplified thermodynamic approach, we estimate the overpotentials related to the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER) on both polar LiNbO3 (0001) surfaces. Our calculations performed for ideal surfaces in vacuum predict the lowest overpotential for the hydrogen evolution reaction (0.4 V) and for the oxygen evolution reaction (1.2 V) at the positive and at the negative surfaces, respectively, which are lower than (or comparable with) commonly employed catalysts. However, calculations performed to model the aqueous solution in which the reactions occur reveal that the presence of water substantially increases the required overpotential for the HER, even inverting the favorable polarization direction for oxidation and reduction reactions. In aqueous solution, we predict an overpotential of 1.2 V for the HER at the negative surface and 1.1 V for the OER at the positive surface.

Optoelectronic Properties of TiS2: A Never Ended Story Tackled by Density Functional Theory and Many-Body Methods

Herein is reported a thorough computational investigation on the bulk TiS₂ material with the CdI₂ structure type and the ideal 1:2 Ti:S stoichiometry. Computations were performed using some of the most refined models, e.g., a hybrid functional together with dispersion effects (Grimme's), the GW ansatz, and the Bethe-Salpether equation for the optical properties. We showed that switching from Perdew-Berke-Enzerhof (PBE) to PBE0 leads to a gap opening. Moreover, our results demonstrate unambiguously that van der Waals interactions must be properly treated with dispersion effects in order to retrieve the experimental crystal structure and the appropriate c/ a ratio. Indeed, the calculations prove that when one uses a highly accurate computational protocol, the bulk hexagonal TiS₂ is a semiconductor with a small gap, whereas using the generalized gradient approximation (GGA) PBE functional leads to a semimetal. Furthermore, the band structure is significantly modified when dispersion parameters are taken into account. Pressure effects were also investigated, and they fully describe the previously simulated electronic transition behavior of the material, e.g., TiS₂ becomes semimetallic under strain.

Assessment of van der Waals inclusive density functional theory methods for adsorption and selective dehydrogenation of formic acid on Pt(111) surface

In this work, we studied the adsorption and catalytic dehydrogenation of formic acid (HCOOH) on Pt(111) surface using different van der Waals inclusive density functional theory (DFT) methods.

Carbon chain growth by formyl coupling over the Cu/γ-AlOOH(001) surface in syngas conversion

Carbon chain growth in syngas-to-ethanol conversion over the Cu/γ-AlOOH(001) surface is mainly ascribed to formyl coupling.

Theory and practice of modeling van der Waals interactions in electronic-structure calculations

The accurate description of long-range electron correlation, most prominently including van der Waals (vdW) dispersion interactions, represents a particularly challenging task in the modeling of molecules and materials. vdW forces arise from the interaction of quantum-mechanical fluctuations in the electronic charge density. Within (semi-)local density functional approximations or Hartree-Fock theory such interactions are neglected altogether. Non-covalent vdW interactions, however, are ubiquitous in nature and play a key role for the understanding and accurate description of the stability, dynamics, structure, and response properties in a plethora of systems. During the last decade, many promising methods have been developed for modeling vdW interactions in electronic-structure calculations. These methods include vdW-inclusive Density Functional Theory and correlated post-Hartree-Fock approaches. Here, we focus on the methods within the framework of Density Functional Theory, including non-local van der Waals density functionals, interatomic dispersion models within many-body and pairwise formulation, and random phase approximation-based approaches. This review aims to guide the reader through the theoretical foundations of these methods in a tutorial-style manner and, in particular, highlight practical aspects such as the applicability and the advantages and shortcomings of current vdW-inclusive approaches. In addition, we give an overview of complementary experimental approaches, and discuss tools for the qualitative understanding of non-covalent interactions as well as energy decomposition techniques. Besides representing a reference for the current state-of-the-art, this work is thus also designed as a concise and detailed introduction to vdW-inclusive electronic structure calculations for a general and broad audience.

Theoretical insight into the origin of the electrochemical promotion of ethylene oxidation on ruthenium oxide

In EPOC the adsorption energies of the reactants and products are modified by applying an electrochemical potential to the catalyst. DFT computations unveil that the C–C and O–O bond dissociation are accelerated under positive and negative potential, respectively.

Can microsolvation effects be estimated from vacuum computations? A case-study of alcohol decomposition at the H2O/Pt(111) interface

Activation and reaction energies of alcohol decomposition at Pt(111) are barely modified by a PCM, in contrast to adding a single water molecule, whose effect can be predicted based on vacuum computations.

Can single-walled carbon nanotube diameter be defined by catalyst particle diameter?

The need of designing and controlling single-walled carbon nanotube (SWCNT) properties is a challenge in a growing nanomaterials-related industry. Recently, great progress has been made experimentally to selectively control SWCNT diameter and chirality. However, there is not yet a complete understanding of the synthesis process and there is a lack of mathematical models that explain nucleation and diameter selectivity of stable carbon allotropes. Here, in-situ analysis of chemical vapor deposition SWCNT synthesis confirms that the nanoparticle to nanotube diameter ratio varies with the catalyst particle size. It is found that the tube diameter is larger than that of the particle below a specific size (d_c ≈ 2nm) and above this value is smaller than particle diameters. To explain these observations, we develop a statistical mechanics based model that correlates possible energy states of a nascent tube with the catalyst particle size. This model incorporates the equilibrium distance between the nucleating SWCNT layer and the metal catalyst (e.g. Fe, Co, Ni) evaluated with density functional theory (DFT) calculations. The theoretical analysis explains and predicts the observed correlation between tube and solid particle diameters during growth of supported SWCNTs. This work also brings together previous observations related to the stability condition for SWCNT nucleation. Tests of the model against various published data sets and our own experimental results show good agreement, making it a promising tool for evaluating SWCNT synthesis processes.

Energy Decomposition Analysis for Metal Surface-Adsorbate Interactions by Block Localized Wave Functions

The energy decomposition analysis based on block localized wave functions (BLW-EDA) allows one to gain physical insight into the nature of chemical bonding, decomposing the interaction energy in (1) a "frozen" term, accounting for the attraction due to electrostatic and dispersion interactions, modulated by Pauli repulsion, (2) the variationally assessed polarization energy, and (3) the charge transfer. This method has so far been applied to gas- and condensed-phase molecular systems. However, its standard version is not compatible with fractionally occupied orbitals (i.e., electronic smearing) and, as a consequence, cannot be applied to metallic surfaces. In this work, we propose a simple and practical extension of BLW-EDA to fractionally occupied orbitals, termed Ensemble BLW-EDA. As illustrative examples, we have applied the developed method to analyze the nature of the interaction of various adsorbates on Pt(111), ranging from physisorbed water to strongly chemisorbed ethylene. Our results show that polarization and charge transfer both contribute significantly at the adsorption minimum for all studied systems. The energy decomposition analysis provides details with respect to competing adsorption sites (e.g., CO on atop vs hollow sites) and elucidates the respective importance of polarization and charge transfer for the increased adsorption energy of H₂S compared to H₂O. Our development will enable a deeper understanding of the impact of charge transfer on catalytic processes in general.

Entrapment of THF-Stabilized Iridacyclic IrIII Silylenes from Double H-Si Bond Activation and H2 Elimination

The reaction of H₃ SiR (R=Ph, nBu) with cationic η₅ -C₅ Me₅ - (Cp*) and benzo[h]quinolinyl-based iridacycle [1 b]⁺ gives rise to new [(IrH)→SiRH₂ ]⁺ adducts. In the presence of THF these adducts readily undergo elimination of H₂ gas at subambient temperature to form THF-stabilized metallacyclic Ir^III silylene complexes, which were characterized in situ by NMR spectroscopy, trapped in minute amounts by reactive crystallization, and structurally characterized by XRD. Theoretical investigations (static DFT-D reaction-energy profiling, ETS-NOCV) support the promoting role of THF in the H₂ elimination step and the consolidation of the Ir-to-Si interaction in the spontaneous (ΔG<0) formation of Ir silylenes in the presence of THF. Mechanistic insights indicate that the Ir silylene species arising from the [1 b]⁺ /phenylsilane system are relevant catalytic species in the hydrodefluorination of fluoroalkanes.

Range-Separated Double-Hybrid Functional from Nonempirical Constraints

On the basis of our previous developments in the field of nonempirical double hybrids, we present here a new exchange-correlation functional based on a range-separated model for the exchange part and integrating a nonlocal perturbative correction to the electron correlation contribution. Named RSX-QIDH, the functional is free from any kind of empirical parametrization. Its range-separation parameter is set to recover the total energy of the hydrogen atom, thus eliminating the self-interaction error for this one-electron system. Subsequent tests on some relevant benchmark data sets confirm that the self-interaction error is particularly low for RSX-QIDH. This new functional provides also correct dissociation profiles for charged rare-gas dimers and very accurate ionization potentials directly from Kohn-Sham orbital energies. Above all, these good results are not obtained at the expense of other properties. Indeed, further tests on standard benchmarks show that RSX-QIDH is competitive with the more empirical ωB97X-2 double hybrid and outperforms the parent LC-PBE long-range corrected hybrid, thus underlining the important role of the nonlocal perturbative correlation.

Few-Layer PdSe2 Sheets: Promising Thermoelectric Materials Driven by High Valley Convergence

Herein, we report a comprehensive study on the structural and electronic properties of bulk, monolayer, and multilayer PdSe2 sheets. First, we present a benchmark study on the structural properties of bulk PdSe2 by using 13 commonly used density functional theory (DFT) functionals. Unexpectedly, the most commonly used van der Waals (vdW)-correction methods, including DFT-D2, optB88, and vdW-DF2, fail to provide accurate predictions of lattice parameters compared to experimental data (relative error > 15%). On the other hand, the PBE-TS series functionals provide significantly improved prediction with a relative error of <2%. Unlike hexagonal two-dimensional materials like graphene, transition metal dichalcogenides, and h-BN, the conduction band minimum of monolayer PdSe2 is not located along the high symmetry lines in the first Brillouin zone; this highlights the importance of the structure-property relationship in the pentagonal lattice. Interestingly, high valley convergence is found in the conduction and valence bands in monolayer, bilayer, and trilayer PdSe2 sheets, suggesting promising application in thermoelectric cooling.

Extension of the D3 dispersion coefficient model

A new model, termed D4, for the efficient computation of molecular dipole-dipole dispersion coefficients is presented. As in the related, well established D3 scheme, these are obtained as a sum of atom-in-molecule dispersion coefficients over atom pairs. Both models make use of dynamic polarizabilities obtained from first-principles time-dependent density functional theory calculations for atoms in different chemical environments employing fractional atomic coordination numbers for interpolation. Different from the D3 model, the coefficients are obtained on-the-fly by numerical Casimir-Polder integration of the dynamic, atomic polarizabilities α(iω). Most importantly, electronic density information is now incorporated via atomic partial charges computed at a semi-empirical quantum mechanical tight-binding level, which is used to scale the polarizabilities. Extended statistical measures show that errors for dispersion coefficients with the proposed D4 method are significantly lower than with D3 and other, computationally more involved schemes. Alongside, accurate isotropic charge and hybridization dependent, atom-in-molecule static polarizabilities are obtained with an unprecedented efficiency. Damping function parameters are provided for three standard density functionals, i.e., TPSS, PBE0, and B3LYP, allowing evaluation of the new DFT-D4 model for common non-covalent interaction energy benchmark sets.

Force Field for Water over Pt(111): Development, Assessment, and Comparison

Metal/water interfaces are key in many natural and industrial processes, such as corrosion, atmospheric, or environmental chemistry. Even today, the only practical approach to simulate large interfaces between a metal and water is to perform force-field simulations. In this work, we propose a novel force field, GAL17, to describe the interaction of water and a Pt(111) surface. GAL17 builds on three terms: (i) a standard Lennard-Jones potential for the bonding interaction between the surface and water, (ii) a Gaussian term to improve the surface corrugation, and (iii) two terms describing the angular dependence of the interaction energy. The 12 parameters of this force field are fitted against a set of 210 adsorption geometries of water on Pt(111). The performance of GAL17 is compared to several other approaches that have not been validated against extensive first-principles computations yet. Their respective accuracy is evaluated on an extended set of 802 adsorption geometries of H₂O on Pt(111), 52 geometries derived from icelike layers, and an MD simulation of an interface between a c(4 × 6) Pt(111) surface and a water layer of 14 Å thickness. The newly developed GAL17 force field provides a significant improvement over previously existing force fields for Pt(111)/H₂O interactions. Its well-balanced performance suggests that it is an ideal candidate to generate relevant geometries for the metal/water interface, paving the way to a representative sampling of the equilibrium distribution at the interface and to predict solvation free energies at the solid/liquid interface.

Noncovalent Molecular Electronics

Molecular electronics covers several distinctly different conducting architectures, including organic semiconductors and single-molecule junctions. The noncovalent interactions, abundant in the former, are also often found in the latter, i.e., the dimer junctions. In the present work, we draw the parallel between the two types of noncovalent molecular electronics for a range of π-conjugated heteroaromatic molecules. In silico modeling allows us to distill the factors that arise from the chemical nature of their building blocks and from their mutual arrangement. We find that the same compounds are consistently the worst and the best performers in the two types of electronic assemblies, emphasizing the universal imprint of the underlying chemistry of the molecular cores on their diverse charge transport characteristics. The interplay between molecular and intermolecular factors creates a spectrum of noncovalent conductive architectures, which can be manipulated using the design strategies based upon the established relationships between chemistry and transport.

Ordering of Air-Oxidized Decanethiols on Au(111)

Self-assembled monolayers (SAMs) of alkanethiols on gold are a commonly used platform for nanotechnology owing to their ease of preparation and high surface coverage. Unfortunately, the gold-sulfur bond is oxidized at ambient conditions which alters the stability and structure of the monolayer. We show using scanning tunneling microscopy and X-ray photoelectron spectroscopy that decanethiolate molecules oxidize into decanesulfonates that organize into a hitherto unknown striped phase. Air-exposed SAMs oxidize, as can be determined by a shift of the S 2p peak and the appearance of O 1s photoelectrons as part of the decanethiol monolayer transforms into a lamellae-like decanesulfonate structure when exposed to air. The herringbone structure of the Au(111) surface is preserved, indicating that the interaction between the molecules and the surface is rather weak as these findings are substantiated by density functional theory calculations.

GAtor: A First-Principles Genetic Algorithm for Molecular Crystal Structure Prediction

We present the implementation of GAtor, a massively parallel, first-principles genetic algorithm (GA) for molecular crystal structure prediction. GAtor is written in Python and currently interfaces with the FHI-aims code to perform local optimizations and energy evaluations using dispersion-inclusive density functional theory (DFT). GAtor offers a variety of fitness evaluation, selection, crossover, and mutation schemes. Breeding operators designed specifically for molecular crystals provide a balance between exploration and exploitation. Evolutionary niching is implemented in GAtor by using machine learning to cluster the dynamically updated population by structural similarity and then employing a cluster-based fitness function. Evolutionary niching promotes uniform sampling of the potential energy surface by evolving several subpopulations, which helps overcome initial pool biases and selection biases (genetic drift). The various settings offered by GAtor increase the likelihood of locating numerous low-energy minima, including those located in disconnected, hard to reach regions of the potential energy landscape. The best structures generated are re-relaxed and re-ranked using a hierarchy of increasingly accurate DFT functionals and dispersion methods. GAtor is applied to a chemically diverse set of four past blind test targets, characterized by different types of intermolecular interactions. The experimentally observed structures and other low-energy structures are found for all four targets. In particular, for Target II, 5-cyano-3-hydroxythiophene, the top ranked putative crystal structure is a Z' = 2 structure with P1̅ symmetry and a scaffold packing motif, which has not been reported previously.

Origin of the size-dependence of the equilibrium van der Waals binding between nanostructures

Nanostructures can be bound together at equilibrium by the van der Waals (vdW) effect, a small but ubiquitous many-body attraction that presents challenges to density functional theory. How does the binding energy depend upon the size or number of atoms in one of a pair of identical nanostructures? To answer this question, we treat each nanostructure as a whole object, not as a collection of atoms. Our calculations start from an accurate static dipole polarizability for each considered nanostructure, and an accurate equilibrium center-to-center distance for the pair (the latter from experiment or from the vdW-DF-cx functional). We consider the competition in each term -C_2k/d^2k (k = 3, 4, 5) of the long-range vdW series for the interaction energy, between the size dependence of the vdW coefficient C_2k and that of the 2kth power of the center-to-center distance d. The damping of these vdW terms can be negligible, but in any case, it does not affect the size dependence for a given term in the absence of non-vdW binding. To our surprise, the vdW energy can be size-independent for quasi-spherical nanoclusters bound to one another by vdW interaction, even with strong nonadditivity of the vdW coefficient, as demonstrated for fullerenes. We also show that, for low-dimensional systems, the vdW interaction yields the strongest size-dependence, in stark contrast to that of fullerenes. We illustrate this with parallel planar polycyclic aromatic hydrocarbons. The size dependences of other morphologies or bonding types lie between, as shown by sodium clusters.

Metastable Structures in Cluster Catalysis from First-Principles: Structural Ensemble in Reaction Conditions and Metastability Triggered Reactivity

Reactivity studies on catalytic transition metal clusters are usually performed on a single global minimum structure. With the example of a Pt₁₃ cluster under a pressure of hydrogen, we show from first-principle calculations that low energy metastable structures of the cluster can play a major role for catalytic reactivity and that hence consideration of the global minimum structure alone can severely underestimate the activity. The catalyst is fluxional with an ensemble of metastable structures energetically accessible at reaction conditions. A modified genetic algorithm is proposed to comprehensively search for the low energy metastable ensemble (LEME) structures instead of merely the global minimum structure. In order to reduce the computational cost of density functional calculations, a high dimensional neural network potential is employed to accelerate the exploration. The presence and influence of LEME structures during catalysis is discussed by the example of H covered Pt₁₃ clusters for two reactions of major importance: hydrogen evolution reaction and methane activation. The results demonstrate that although the number of accessible metastable structures is reduced under reaction condition for Pt₁₃ clusters, these metastable structures can exhibit high activity and dominate the observed activity due to their unique electronic or structural properties. This underlines the necessity of thoroughly exploring the LEME structures in catalysis simulations. The approach enables one to systematically address the impact of isomers in catalysis studies, taking into account the high adsorbate coverage induced by reaction conditions.

How do London Dispersion Interactions Impact the Photochemical Processes of Molecular Switches?

In the last two decades, linear-response time-dependent density functional theory (LR-TDDFT) has become one of the most widely used approaches for the computation of the excited-state properties of atoms and molecules. Despite its success in describing the photochemistry and the photophysics of a vast majority of molecular systems, its domain of applicability has been limited by several substantial drawbacks. Commonly identified problems of LR-TDDFT include the correct description of Rydberg states, charge-transfer excited states, doubly excited states, and nearly degenerate states. In addition to these widely recognized shortcomings, the approximate functionals used in LR-TDDFT are unable to fully describe London dispersion interactions. In this work, we aim at understanding the impact of van der Waals interactions on the properties of chemical systems beyond their electronic ground state. For this purpose, we compare the results of excited-state energy profiles and dynamic trajectories for the prototypical cis-stilbene molecule with its 3-3',5-5'-tetra-tert-butyl derivative. While the explicit treatment of London dispersion interactions results in negligible changes for the cis-stilbene, we show that these attractive forces have a substantial influence on the energetics and structural evolution of the substituted derivative. In the latter case, intramolecular dispersion interactions impact the outcome of the simulation qualitatively, leading to an increased preference for the photocyclization pathway. The methodological consequences of this work are not uniquely applicable to the illustrative stilbene case. In fact, this molecule is representative of a whole class of chemical situations, where dispersion forces dominate the interactions between the unexcited substituents of a photoexcited chromophore. This is, for instance, a common situation in organic photovoltaics where donor molecules are usually functionalized with long alkyl side chains to improve solubility and assembly.

Supported gold–nickel nano-alloy as a highly efficient catalyst in levulinic acid hydrogenation with formic acid as an internal hydrogen source

The surface Au–Ni nano-alloy was very efficiently used for the first time for the hydrogenation of levulinic acid with formic acid as an internal hydrogen source.

Unsupported shaped cobalt nanoparticles as efficient and recyclable catalysts for the solvent-free acceptorless dehydrogenation of alcohols

Unsupported cobalt nanoparticles are active and chemoselective catalysts for the solvent-free acceptorless dehydrogenation of alcohols.

Direct n-octanol amination by ammonia on supported Ni and Pd catalysts: activity is enhanced by “spectator” ammonia adsorbates

DFT calculations highlight the role of co-adsorbed ammonia in catalytic activity in the amination of alcohols by ammonia.

A β-Carbon elimination strategy for convenient in situ access to cyclopentadienyl metal complexes

The electronic and steric properties of tailored cyclopentadienyl (Cp) ligands are powerful handles to modulate the catalytic properties of their metal complexes. This requires the individual preparation, purification and storage of each ligand/metal combination. Alternative, ideally in situ, complexation protocols would be of high utility. We disclose a new approach to access Cp metal complexes. Common metal precursors rapidly react with cyclopentadienyl carbinols via β-carbon eliminations to directly give the Cp-metal complexes. An advantage of this is the direct and flexible use of storable pre-ligands. No auxiliary base is required and the Cp complexes can be prepared in situ in the reaction vessel for subsequent catalytic transformations.