Realistic modeling of ruthenium-catalyzed transfer hydrogenation

We report the first computational study of a fully atomistic model of the ruthenium-catalyzed transfer hydrogenation of formaldehyde and the reverse reaction in an explicit methanol solution. Using ab initio molecular dynamics techniques, we determined the thermodynamics, mechanism, and electronic structure along the reaction path. To assess the effect of the solvent quantitatively, we make a direct comparison with the gas-phase reaction. We find that the energy profile in solution bears little resemblance to the profile in the gas phase and a distinct solvation barrier is found: the activation barriers in both directions are lowered and the concerted hydride and proton transfer in the gas phase are converted into a sequential mechanism in solution with the substrate appearing as methoxide-like intermediate. Our results indicate that besides the metal-ligand bifunctional mechanism, as proposed by Noyori, also a concerted solvent-mediated mechanism is feasible. Our study gives a new perspective of the active role a solvent can have in transition-metal-catalyzed reactions.     

Polyoxometalates with internal cavities: redox activity, basicity, and cation encapsulation in [Xn+P5W30O110](15-n)- Preyssler complexes, with X = Na+, Ca2+, Y3+, La3+, Ce3+, and Th4+

The Preyssler anion, of general formula [Xn+P5W30O110](15-n)-, is the smallest polyoxometalate (POM) with an internal cavity allowing cation exchange with the solution. The Preyssler anion features a rich chemistry evidenced by its ability to accept electrons at low potentials, to selectively capture various metal cations, and to undergo acid-base reactions. A deep understanding of these topics is herein provided by means of DFT calculations on the title series of compounds. We evaluate the energetics of the release/encapsulation process for several Xn+ cations and identify the effect of the encapsulated ion on the properties of the Preyssler anion. We revisit the relationship between the internal cation charge and the electrochemical behavior of the POM. A linear dependence between the first one-electron reduction energies and the encapsulated Xn+ charge is found, with a slope of 48 mV per unit charge. The protonation also shifts the reduction potential to more positive values, but the effect is much larger. In connection to this, the last proton's pKa = 2 for the Na+ derivative was estimated to be in reasonable agreement with experiment. The electronic structure of lanthanide derivatives is more complex than conventional POM structures. The reduction energy for the CeIV-Preyssler + 1e- --> CeIII-Preyssler process was computed to be more exothermic than that of very oxidant species such as S2Mo18O624-.     

Oxidative properties of FeO2+: electronic structure and solvation effects

An electronic structure analysis is provided of the action of solvated FeO(2+), [FeO(H(2)O)(5)](2+), as a hydroxylation catalyst. It is emphasized that the oxo end of FeO(2+) does not form hydrogen bonds (as electron donor and H-bond acceptor) with H-bond donors nor with aliphatic C-H bonds, but it activates C-H bonds as an electron acceptor. It is extremely electrophilic, to the extent that it can activate even such poor electron donors as aliphatic C-H bonds, the C-H bond orbital acting as electron donor in a charge transfer type of interaction. Lower lying O-H bonding orbitals are less easily activated. The primary electron accepting orbital in a water environment is the 3sigma*alpha orbital, an antibonding combination of Fe-3d(z(2)) and O-2p(z), which is very low-lying relative to the pi*alpha compared with, for example, the sigma* orbital in O(2) relative to its pi*. This is ascribed to relatively small Fe-3d(z(2)) with O-2p(z) overlap, due to the nodal structure of the 3d(z(2)).The H-abstraction barrier is very low in the gas phase, but it is considerably enhanced in water solvent. This is shown to be due to strong screening effects of the dielectric medium, leading to relative destabilization of the levels of the charged [FeO(H(2)O)(5)](2+) species compared to those of the neutral substrate molecules, making it a less effective electron acceptor. The solvent directly affects the orbital interactions responsible for the catalytic reaction.     

Tensile deformation behaviour of the polymer phase of flexible polyurethane foams and polyurethane elastomers

The change in micromorphology of the polymer phase (single strut) of a flexible polyurethane foam during deformation has been investigated by attenuated total reflection infrared spectroscopy - linear dichroism and by atomic force microscopy. Deformation and, therefore, orientation take place mainly in the soft rubbery phase. This two-phase elastic deformation process has been translated into a mathematical model, which correctly predicts the shape of a single-strut stress-strain curve. The theory also predicts the ultimate shape of stress-strain curves of polyurethane elastomers at various hard phase contents and of low-density polyethylene at various temperatures. Deviations from the elastic behaviour could be ascribed to yielding in combination with the rubbery behaviour.     

Gas phase rate coefficients and activation energies for the reaction of butanal and 2‐methyl‐propanal with nitrate radicals

Rate coefficients have been determined for the reaction of butanal and 2‐methyl‐propanal with NO₃ using relative and absolute methods. The relative measurements were accomplished by using a static reactor with long‐path FTIR spectroscopy as the analytical tool. The absolute measurements were made using fast‐flow–discharge technique with detection of NO₃ by optical absorption. The resulting average coefficients from the relative rate experiments were k = (1.0 ± 0.1) × 10⁻¹⁴ and k = (1.2 ± 0.2) × 10⁻¹⁴ (cm³ molecule⁻¹ s⁻¹) for butanal and 2‐methyl‐propanal, respectively. The results from the absolute measurements indicated secondary reactions involving NO₃ radicals and the primary formed acyl radicals. The prospect of secondary reactions was investigated by means of mathematical modeling. Calculations indicated that the unwanted NO₃ radical reactions could be suppressed by introducing molecular oxygen into the flow tube. The rate coefficients from the absolute rate experiments with oxygen added were and k = (1.2 ± 0.1) × 10⁻¹⁴ and = (0.9 ± 0.1) × 10⁻¹⁴ (cm³ molecule⁻¹ s⁻¹) for butanal and 2‐methyl‐propanal. The temperature dependence of the reactions was studied in the range between 263 and 364 K. Activation energies for the reactions were determined to 12 ± 2 kJ mole⁻¹ and 14 ± 1 kJ mole⁻¹ for butanal and 2‐methyl‐propanal, respectively. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 294–303, 2000     

Interactions of multiple spark-generated bubbles with phase differences

This paper aims to study the complex interaction between multiple bubbles, and to provide a summary and physical explanation of the phenomena observed during the interaction of two bubbles. High-speed photography is utilized to observe the experiments involving multiple spark-generated bubbles. Numerical simulations corresponding to the experiments are performed using the Boundary Element Method (BEM). The bubbles are typically between 3 and 5 mm in radius and are generated either in-phase (at the same time) or with phase differences. Complex phenomena are observed such as bubble splitting, and high-speed jetting inside a bubble caused by another collapsing bubble nearby (termed the ‘catapult’ effect). The two-bubble interactions are broadly classified in a graph according to two parameters: the relative inter-bubble distance and the phase difference (a new parameter introduced). The BEM simulations provide insight into the physics, such as bubble shape changes in detail, and jet velocities. Also presented in this paper are the experimental results of three bubble interactions. The interesting and complex observations of multiple bubble interaction are important for a better understanding of real life applications in medical ultrasonic treatment and ultrasonic cleaning. Many of the three bubble interactions can be explained by isolating bubble pairs and classifying their interaction according to the graph for the two bubble case. This graph can be a useful tool to predict the behavior of multiple bubble interactions.     

Dynamic Interplay between Defective UiO-66 and Protic Solvents in Activated Processes

UiO-66, composed by Zr-oxide inorganic bricks [Zr₆(μ₃-O)₄(μ₃-OH)₄] and organic terephthalate linkers, is one of the most studied metal–organic frameworks (MOFs) due to its exceptional thermal, chemical, and mechanical stability. Thanks to its high connectivity, the material can withstand structural deformations during activation processes such as linker exchange, dehydration, and defect formation. These processes do alter the zirconium coordination number in a dynamic way, creating open metal sites for catalysis and thus are able to tune the catalytic properties. In this work, it is shown, by means of first-principle molecular-dynamics simulations at operating conditions, how protic solvents may facilitate such changes in the metal coordination. Solvent can induce structural rearrangements in the material that can lead to undercoordinated but also overcoordinated metal sites. This is demonstrated by simulating activation processes along well-chosen collective variables. Such enhanced MD simulations are able to track the intrinsic dynamics of the framework at realistic conditions.     

Assessment of a simple correction for the long-range charge-transfer problem in time-dependent density-functional theory

The failure of the time-dependent density-functional theory to describe long-range charge-transfer (CT) excitations correctly is a serious problem for calculations of electronic transitions in large systems, especially if they are composed of several weakly interacting units. The problem is particularly severe for molecules in solution, either modeled by periodic boundary calculations with large box sizes or by cluster calculations employing extended solvent shells. In the present study we describe the implementation and assessment of a simple physically motivated correction to the exchange-correlation kernel suggested in a previous study [O. Gritsenko and E. J. Baerends J. Chem. Phys. 121, 655 (2004)]. It introduces the required divergence in the kernel when the transition density goes to zero due to a large spatial distance between the "electron" (in the virtual orbital) and the "hole" (in the occupied orbital). A major benefit arises for solvated molecules, for which many CT excitations occur from solvent to solute or vice versa. In these cases, the correction of the exchange-correlation kernel can be used to automatically "clean up" the spectrum and significantly reduce the computational effort to determine low-lying transitions of the solute. This correction uses a phenomenological parameter, which is needed to identify a CT excitation in terms of the orbital density overlap of the occupied and virtual orbitals involved. Another quantity needed in this approach is the magnitude of the correction in the asymptotic limit. Although this can, in principle, be calculated rigorously for a given CT transition, we assess a simple approximation to it that can automatically be applied to a number of low-energy CT excitations without additional computational effort. We show that the method is robust and correctly shifts long-range CT excitations, while other excitations remain unaffected. We discuss problems arising from a strong delocalization of orbitals, which leads to a breakdown of the correction criterion.     

Spatial sign preprocessing: a simple way to impart moderate robustness to multivariate estimators

The spatial sign is a multivariate extension of the concept of sign. Recently multivariate estimators of covariance structures based on spatial signs have been examined by various authors. These new estimators are found to be robust to outlying observations. From a computational point of view, estimators based on spatial sign are very easy to implement as they boil down to a transformation of the data to their spatial signs, from which the classical estimator is then computed. Hence, one can also consider the transformation to spatial signs to be a preprocessing technique, which ensures that the calibration procedure as a whole is robust. In this paper, we examine the special case of spatial sign preprocessing in combination with partial least squares regression as the latter technique is frequently applied in the context of chemical data analysis. In a simulation study, we compare the performance of the spatial sign transformation to nontransformed data as well as to two robust counterparts of partial least squares regression. It turns out that the spatial sign transform is fairly efficient but has some undesirable bias properties. The method is applied to a recently published data set in the field of quantitative structure-activity relationships, where it is seen to perform equally well as the previously described best linear model for these data.