Prediction of the reduction potential in transition‐metal containing complexes: How expensive? For what accuracy?

Accurate computationally derived reduction potentials are important for catalyst design. In this contribution, relatively inexpensive density functional theory methods are evaluated for computing reduction potentials of a wide variety of organic, inorganic, and organometallic complexes. Astonishingly, SCRF single points on B3LYP optimized geometries with a reasonably small basis set/ECP combination works quite well‐‐B3LYP with the BS1 [modified‐LANL2DZ basis set/ECP (effective core potential) for metals, LANL2DZ(d,p) basis set/LANL2DZ ECP for heavy nonmetals (Si, P, S, Cl, and Br), and 6‐31G(d') for other elements (H, C, N, O, and F)] and implicit PCM solvation models, SMD (solvation model based on density) or IEFPCM (integral equation formalism polarizable continuum model with Bondi atomic radii and α = 1.1 reaction field correction factor). The IEFPCM‐Bondi‐B3LYP/BS1 methodology was found to be one of the least expensive and most accurate protocols, among six different density functionals tested (BP86, PBEPBE, B3LYP, B3P86, PBE0, and M06) with thirteen different basis sets (Pople split‐valence basis sets, correlation consistent basis sets, or Los Alamos National Laboratory ECP/basis sets) and four solvation models (SMD, IEFPCM, IPCM, and CPCM). The MAD (mean absolute deviation) values of SCRF‐B3LYP/BS1 of 49 studied species were 0.263 V for SMD and 0.233 V for IEFPCM‐Bondi; and the linear correlations had respectable R ² values (R ² = 0.94 for SMD and R ² = 0.93 for IEFPCM‐Bondi). These methodologies demonstrate relatively reliable, convenient, and time‐saving functional/basis set/solvation model combinations in computing the reduction potentials of transition metal complexes with moderate accuracy. © 2017 Wiley Periodicals, Inc.     

Can configuration entropy losses be predicted from the binding affinities of hydrogen-bonded complexes with varying numbers of single bonds?

The design of new supramolecular complexes often depends on reducing entropic contributions to improve binding. However, few studies provide reliable values for the cost of entropic contributions. Here, the binding affinities of a series of six alpha,omega-diamides to alpha,omega-dicarboxylates are calculated using a predominant states method and an implicit solvent model based upon finite difference solutions of the Poisson-Boltzmann equation. The calculations are able to reproduce the observed increase in binding free energy as the number of single bonds increases. However, calculations show that the increase in binding free energy is not due to an increase in entropy. Instead, the increase is due to the changing ability of the alpha,omega-diamides to form internal hydrogen bonds that must be disrupted to bind to the dicarboxylate receptors. This suggests that interpreting experimental free-energy trends to give rotational entropy contributions may be problematic.     

Can metal ions be used as gas‐phase disulfide bond cleavage reagents? A survey of coinage metal complexes of model peptides containing an intermolecular disulfide bond

The role that a metal ion can have in promoting disulfide bond cleavage has been assessed by surveying the tandem mass spectra of the following metal complexes of model peptides containing an intermolecular disulfide bond: [M--H+Cu(II)](+); [M--H+Cu(II)(bipy)](+); [M+Ag(I)](+); and [M+Au(I)(PMe(3))](+). In comparison to previously studied protonated peptides, these binary and ternary metal complexes generally yield more abundant S--S and/or C--S bond cleavage. In general, [M--H+Cu(II)](+) ions cleave the adjacent C--S bond more readily, while the [M+Au(I)(PMe(3))](+) ion cleaves the S--S bond more readily. The ternary metal complex [M--H+Cu(II)(bipy)](+), on the other hand, fragments by exclusive loss of the bipyridyl ligand for the larger model peptides studied. Of all coinage metal systems studied, Me(3)PAu(+) is superior in promoting disulfide bond cleavage.     

Can the pi-facial selectivity of solvation be predicted by atomistic simulation?

This work is concerned with the rationalization and prediction of solvent and temperature effects in nucleophilic addition to alpha-chiral carbonyl compounds leading to facial diastereoselectivity. We study, using molecular dynamics simulations, the facial solvation of (R)-2-phenyl-propionaldehyde in n-pentane and n-octane at a number of temperatures and compare it with experimental selectivity data for the nBuLi addition leading to syn- and anti-(2R)-2-phenyl-3-heptanol, which give nonlinear Eyring plots with the presence of inversion temperatures. We have found from simulations that the facial solvation changes with temperature and alkane. Moreover, by introducing a suitable molecular chirality index we have been able to predict break temperatures (T(CI)) for the two solvents within less than 20 degrees of the inversion temperatures experimentally observed in the diastereoselective nBuLi addition. We believe this could lead to a viable approach for predicting inversion temperatures and other subtle solvent effects in a number of stereoselective reactions.     

Can nanofiltration be fully predicted by a model?

There is an increasing need for model-based tools to design membrane processes for new industrial applications or to optimise existing membrane installations. The advantage of such tools is that costs can be saved by reducing the number of expermiments. In this study, the requirements for a membrane filtration model, suitable for practical use, are summarised. It is investigated to what extent it is possible to set-up such a model with the current available literature and knowledge. A membrane filtration model has been set-up based on the Maxwell–Stefan transport equations. A Freundlich equation is used to describe the membrane charge by means of adsorption of ions. With the model the permeate flux and rejections of multi-component liquid feeds can be calculated as a function of membrane properties (mean pore size, porosity, thickness, surface charge characteristic) and feed pressure. With two NF-membranes (Desal 5DK and a prototype capillary type 2 membrane) rejection experiments have been carried out with glucose, single salt solutions (NaCl, CaCl₂, Na₂SO₄) and ternary ion mixtures of these salts. With the model the experimental flux-rejection curves can be fitted reasonably well. However, each salt mixture needs its own set of fitted parameters for the membrane charge isotherms. Furthermore, the fitted membrane charges are in contradiction with values from the literature obtained by electrokinetic measurements. Obviously, the membrane charge parameters have lost their physical meaning and are used to compensate for physical phenomena not included in the model. Extending the model with an electrostatic free energy term will be a step forward in development. Further research is needed to fulfil all requirements for the wide scope of industrial applications.     

Theoretical calculations: Can Gibbs free energy for intermolecular complexes be predicted efficiently and accurately?

The theoretical study has been performed to refine the procedure for calculations of Gibbs free energy with a relative accuracy of less than 1 kcal/mol. Three benchmark intermolecular complexes are examined via several quantum-chemical methods, including the second-order Moller-Plesset perturbation (MP2), coupled cluster (CCSD(T)), and density functional (BLYP, B3LYP) theories augmented by Dunnings correlation-consistent basis sets. The effects of electron correlation, basis set size, and anharmonicity are systematically analyzed, and the results are compared with available experimental data. The results of the calculations suggest that experimental accuracy can be reached only by extrapolation of MP2 and CCSD(T) total energies to the complete basis set. The contribution of anharmonicity to the zero point energy and TDeltaSint values is fairly small. The new, economic way to reach chemical accuracy in the calculations of the thermodynamic parameters of intermolecular interactions is proposed. In addition, interaction energy (De) and free energy change (DeltaA) for considered species have been evaluated by Carr-Parrinello molecular dynamics (CPMD) simulations and static BLYP-plane wave calculations. The free energy change along the reaction paths were determined by the thermodynamic integration/"Blue Moon Ensemble" technique. Comparison between obtained values, and available experimental and conventional ab initio results has been made. We found that the accuracy of CPMD simulations is affected by several factors, including statistical uncertainty and convergence of constrained forces (TD integration), and the nature of DFT (density functional theory) functional. The results show that CPMD technique is capable of reproducing interaction and free energy with an accuracy of 1 kcal/mol and 2-3 kcal/mol respectively.     

Mechanism of GTP: Can all of the available data be accommodated?

Group transfer polymerization (GTP)^a of acrylic monomers is a living system mediated by a trialkylsilyl capped growing chain end. The fact that it operates at temperatures as high as 100° C differentiates GTP from living anionic polymerization, which at best operates at 25° C for hindered methacrylates. To accommodate all of the mechanistic data available it appears that at least two mechanistic pathways are required: an associative process for mild nucleophilic catalysts and a dissociative process for strong nucleophilic catalysts.     

Can we understand the different coordinations and structures of closed‐shell metal cation‐water clusters?

We present a model potential for studying M(q+)(H(2)O)(n=1,9) clusters where M stands for either Na(+), Cs(+), Ca(2+), Ba(2+), or La(3+). The potential energy surfaces (PES) are explored by the Monte Carlo growth method. The results for the most significant equilibrium structures of the PES as well as for energetics are favorably compared to the best ab initio calculations found in the literature and to experimental results. Most of these complexes have a different coordination number in cluster compared to experimental results in solution or solid phase. An interpretation of the coordination number in clusters is given. In order to well describe the transition between the first hydration sphere and the second one we show that an autocoherent treatment of the electric field is necessary to correctly deal with polarization effects. We also explore the influence of the cation properties (charge, size, and polarizability) on both structures and coordination number in clusters, as well as the meaning of the second hydration sphere. Such an approach shows that the leading term in the interaction energy for a molecule in the second hydration sphere is an electrostatic attraction to the cation and not a hydrogen bond with the water molecules in the first hydration sphere.     

Can radical cations of the constituents of nucleic acids be formed in the gas phase using ternary transition metal complexes?

Electrospray ionization (ESI) tandem mass spectrometry (MS/MS) of ternary transition metal complexes of [M(L(3))(N)](2+) (where M = copper(II) or platinum(II); L(3) = diethylenetriamine (dien) or 2,2':6',2'-terpyridine (tpy); N = the nucleobases: adenine, guanine, thymine and cytosine; the nucleosides: 2'deoxyadenosine, 2'deoxyguanosine, 2'deoxythymine, 2'deoxycytidine; the nucleotides: 2'deoxyadenosine 5'-monophosphate, 2'deoxyguanosine 5'-monophosphate, 2'deoxythymine 5'-monophosphate, 2'deoxycytidine 5'-monophosphate) was examined as a means of forming radical cations of the constituents of nucleic acids in the gas phase. In general, sufficient quantities of the ternary complexes [M(L(3))(N)](2+) could be formed for MS/MS studies by subjecting methanolic solutions of mixtures of a metal salt [M(L(3))X(2)] (where M = Cu(II) or Pt(II); L(3) = dien or tpy; X = Cl or NO(3)) and N to ESI. The only exceptions were thymine and its derivatives, which failed to form sufficient abundances of [M(L(3))(N)](2+) ions when: (a) M = Pt(II) and L(3) = dien or tpy; (b) M = Cu(II) and L(3) = dien. In some instances higher oligomeric complexes were formed; e.g., [Pt(tpy)(dG)(n)](2+) (n = 1-13). Each of the ternary complexes [M(L(3))(N)](2+) was mass-selected and then subjected to collision-induced dissociation (CID) in a quadrupole ion trap. The types of fragmentation reactions observed for these complexes depend on the nature of all three components (metal, auxiliary ligand and nucleic acid constituent) and can be classified into: (i) a redox reaction which results in the formation of the radical cation of the nucleic acid constituent, N(+.); (ii) loss of the nucleic acid constituent in its protonated form; and (iii) fragmentation of the nucleic acid constituent. Only the copper complexes yielded radical cations of the nucleic acid constituent, with [Cu(tpy)(N)](2+) being the preferred complex due to suppression, in this case, of the loss of the nucleobase in its protonated form. The yields of the radical cations of the nucleobases from the copper complexes follow the order of their ionization potentials (IPs): G (lowest IP) > A > C > T (highest IP). Sufficient yields of the radical cations of each of the nucleobases allowed their CID reactions (in MS(3) experiments) to be compared to their even-electron counterparts.     

Can the transition from tunneling to hopping in molecular junctions be predicted by theoretical calculation?

The electron transport mechanism changes from tunneling to hopping as molecular length increases. To validate the theoretical simulation after the transition point and clarify influence of electronic structures on the transition, we calculated the conductance of a series of conjugated molecules by density functional theory together with the nonequilibrium Green's function. We found that the highest occupied molecular orbital energy level, transmission spectrum, and the reorganization energy are good indicators for the transition of the electron transport mechanism. The calculated resistances of short junctions (<50 Å, before the transition point) are consistent with the experimental result, following the tunneling mechanism. However, the theoretical predication failed for long molecules, indicating the limitation of the theoretical framework of elastic scattering when the electron transport mechanism changes to hopping. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Do b‐ions occur from vibrational excitation upon H‐desorption in electron capture dissociation?

Electron capture dissociation (ECD) has recently been shown in some cases to produce abundant N-terminal b-ion peptide fragments. These product ions are usually only observed when activation occurs via vibrational excitation as in collision-induced dissociation (CID). Here, we show that occurrence of b-ions in the ECD spectra of synthetic peptides are correlated with low gas-phase basicity and that the observed b-ion fragments are N-terminal products. Furthermore, all ECD spectra containing b-ions also had abundant losses of hydrogen and ammonia from the charge-reduced species.     

Fluorescence of New o‐Carborane Compounds with Different Fluorophores: Can it be Tuned?

Two sets of o‐carborane derivatives incorporating fluorene and anthracene fragments as fluorophore groups have been successfully synthesized and characterized, and their photophysical properties studied. The first set, comprising fluorene‐containing carboranes 6 – 9 , was prepared by catalyzed hydrosilylation reactions of ethynylfluorene with appropriate carboranylsilanes. The compound 1‐[(9,9‐dioctyl‐fluorene‐2‐yl)ethynyl]carborane ( 11 ) was synthesized by the reaction of 9,9‐dioctyl‐2‐ethynylfluorene and decaborane (B₁₀H₁₄). Furthermore, reactions of the lithium salt of 11 with 1 equivalent of 4‐(chloromethyl)styrene or 9‐(chloromethyl)anthracene yielded compounds 12 and 13 . Members of the second set of derivatives, comprising anthracene‐containing carboranes, were synthesized by reactions of monolithium or dilithium salts of 1‐Me‐1,2‐C₂B₁₀H₁₁, 1‐Ph‐1,2‐C₂B₁₀H₁₁, and 1,2‐C₂B₁₀H₁₂ with 1 or 2 equivalents of 9‐(chloromethyl)anthracene, respectively, to produce compounds 14 – 16 . In addition, 2 equivalents of the monolithium salts of 1‐Me‐1,2‐C₂B₁₀H₁₁ (Me‐o‐carborane) and 1‐Ph‐1,2‐C₂B₁₀H₁₁ (Ph‐o‐carborane) were reacted with 9,10‐bis(chloromethyl)anthracene to produce compounds 17 and 18 , respectively. Fluorene derivatives 6 – 9 exhibit moderate fluorescence quantum yields (32–44 %), whereas 11 – 13 , in which the fluorophore is bonded to the C_cluster (C_c), show very low emission intensity (6 %) or complete fluorescence quenching. The anthracenyl derivatives containing the Me‐o‐carborane moiety exhibit notably high fluorescence emissions, with ?_F=82 and 94 %, whereas their Ph‐o‐carborane analogues are not fluorescent at all. For these compounds, we have observed a correlation between the C_c?C_c bond length and the fluorescence intensity in CH₂Cl₂ solution, comparable to that observed for previously reported styrene‐containing carboranes. Thus, our hypothesis is that for systems of this type the fluorescence may be tuned and even predicted by changing the substituent on the adjacent C_c.     

Can Enantioselectivity be Computed in Enthalpic Barrierless Reactions? The Case of CuI‐Catalyzed Cyclopropanation of Alkenes

An extensive computational study has been carried out on different catalytic systems for cyclopropanation reactions based on copper. Most DFT schemes used present drawbacks that preclude the calculation of accurate absolute kinetic properties (energy barriers) of such systems, excepting the M05 and M06 suites of density functionals. On the other hand, there is a wide range of DFT methods capable of reproducing relative energy values, which can be easily translated into selectivities. Most of the theoretical levels used tend to overestimate activation barriers, allowing the location of the transition state (TS) on the potential-energy surface (PES) of the most reactive systems, which are probably artifacts of the method. However, after a thorough analysis of the calculated PES, and the origin of the energy differences obtained for the different alkene approaches in chiral systems, it is found that energy differences are almost constant over a wide range of geometries covering the reaction channel zone in which the true TS on the Gibbs free-energy surface (GFES) lies. Therefore, many computational schemes can still be used confidently to explain and predict enantioselectivities in these systems.     

Can be Azo Dyes Obtained by Grinding under Solvent‐Free Conditions?

The solid‐solid reactions of some electron‐donors with sulfanilic acid in the presence of solid sodium nitrite afford azo dyes by self‐catalyzed diazotization of sulfanilic acid (2) under solvent‐free conditions with moderate yields. Also the reactions of some electron‐donors with diazotization of o‐nitroaniline (5), m‐nitroaniline (6) and p‐nitroaniline (7) in the presence of solid sodium nitrite catalyzed by p‐toluenesulfonic acid (PTSA) afford azo dyes under solvent‐free conditions in good yields. This new method totally avoids the use of acids, alkalies, and toxic and/or expensive solvents in diazotization and diazo coupling reactions.     

Can Coordination‐Driven Supramolecular Self‐Assembly Reactions Be Conducted from Fully Aliphatic Linkers?

The reaction between a preassembled Cu^I bimetallic molecular clip with a short intermetallic distance and a series of fully aliphatic cyano‐capped ditopic linkers with increasing lengths was investigated. It is shown that, depending on the length of the ditopic linkers, the rational design of unprecedented supramolecular compact metallacycles containing fully aliphatic walls is possible. The specific preorganized molecular arrangement of the molecular clip used favors stabilizing interlinker London dispersion interactions, which allow, as the length of the linkers increases, the selective formation of discrete compact metallacycles at the expense of 1D coordination polymers. The generalizability of this approach was demonstrated by the reaction of fully aliphatic cyano‐capped linkers with two other types of preassembled Cu^I bimetallic molecular clips that also had short intermetallic distances.     

How Can a Hydrophobic MOF be Water‐Unstable? Insight into the Hydration Mechanism of IRMOFs

We report an ab initio molecular dynamics study of the hydration process in a model IRMOF material. At low water content (one molecule per unit cell), water physisorption is observed on the zinc cation but the free?bound equilibrium strongly favors the free state. This is consistent with the hydrophobic nature of the host matrix and its type‐V isotherm observed in a classical Monte Carlo simulation. At higher loading, a water cluster can be formed at the Zn₄O site and this is shown to stabilize the water‐bound state. This structure rapidly transforms into a linker‐displaced state, where water has fully displaced one arm of a linker and which corresponds to the loss of the material’s fully ordered structure. Thus an overall hydrophobic MOF material can also become water unstable, a feature that has not been fully understood until now.     

Synthesis,structures and electrochemical properties of four organic–inorganic hybrid polyoxometalates constructed from polyoxotungstate clusters and transition metal complexes

Four novel organic–inorganic hybrid compounds [Cu₅ ^I(4,4′-bpy)₃(2,2′-bpy)₄][BW₁₂O₄₀] · H₂O (1), [Ni_0.5(2,2′-bpy)_1.25][Ni(2,2′-bpy)₃][Ni(2,2′-bpy)₂(H₂O)(SiW₁₁^VIW^VO₄₀)] · 0.5H₂O (2), [H₂bpy]₂[Zn(2,2′-bpy)₃]₂[Si₂W₁₈O₆₂] · 1.5H₂O (3) and [Cu^II(2,2′-bpy)₂]₂[SiW₁₂O₄₀] · 2H₂O (4) (2,2′-bpy = 2,2′-bipyridine, 4,4′-bpy = 4,4′-bipyridine) have been synthesized under hydrothermal conditions and characterized by elemental analysis, IR spectroscopy, thermal gravimetric analysis, electrochemical measurements and single-crystal X-ray diffraction. Compound (1) is a novel [BW₁₂O₄₀]⁵⁻ polyoxoanion bisupported by copper(I) coordination cations with mixed 2,2′-bpy and 4,4′-bpy ligands. Compound (2) is constructed from the [SiW₁₁^VIW^VO₄₀]⁵⁻ polyoxoanions supported by [Ni(2,2′–bpy)₂]²⁺. Compound (3) is composed of a novel [Si₂W₁₈O₆₂]⁸⁻ cluster and [Zn(2,2′–bpy)₃]²⁺ complexes, which held together into a three-dimensional (3D) supramolecular network through hydrogen-bonding interactions. Compound (4) shows a 2D layer framework constructed from a bisupporting Keggin polyoxoanion cluster and [Cu(2,2′–bpy)₂]²⁺ coordination polymer fragments, resulting in 3D networks via supramolecular interactions. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Weighing synthetic polymers of ultra‐high molar mass and polymeric nanomaterials: What can we learn from charge detection mass spectrometry?

Advances in soft ionization techniques for mass spectrometry (MS) of polymeric materials make it possible to determine the masses of intact molecular ions exceeding megadaltons. Interfacing MS with separation and fragmentation methods has additionally led to impressive advances in the ability to structurally characterize polymers. Even if the gap to the megadalton range has been bridged by MS for polymers standards, the MS‐based analysis for more complex polymeric materials is still challenging. Charge detection mass spectrometry (CDMS) is a single‐molecule method where the mass and the charge of each ion are directly determined from individual measurements. The entire molecular mass distribution of a polymer sample can be thus accurately measured. Described in this perspective paper is how molecular weight distribution as well as charge distribution can provide new insights into the structural and compositional studies of synthetic polymers and polymeric nanomaterials in the megadalton to gigadalton range of molecular weight. The recent multidimensional CDMS studies involving couplings with separation and dissociation techniques will be presented. And, finally, an outlook for the future avenues of the CDMS technique in the field of synthetic polymers of ultra‐high molar mass and polymeric nanomaterials will be provided.     

How Can Theory Predict the Selectivity of Palladium‐Catalyzed Cross‐Coupling of Pristine Aromatic Molecules?

The new approach for palladium‐catalyzed cross‐coupling of two non‐activated aromatic compounds (D. R. Stuart, K. Fagnou, Science 2007 , 316, 1172) was studied theoretically. The energetic span model (S. Kozuch, S. Shaik, Acc. Chem. Res. 2011 , 44, 101, and references therein) was employed to analyze the kinetic behavior of the catalytic cycle. The computed energy profile, combined with the energetic span model, accounts for the experimental selectivity, which favors the hetero‐coupling of benzene with indole. This selectivity is driven by a fine balance of the entropic contributions and the high ratio of concentrations used for benzene over indole. This analysis may allow future theoretical predictions of how different aromatic compounds can be effectively coupled.