Hydrogen‐Bond Cooperative Effects in Small Cyclic Water Clusters as Revealed by the Interacting Quantum Atoms Approach

The cooperative effects of hydrogen bonding in small water clusters (H₂O)_n (n=3–6) have been studied by using the partition of the electronic energy in accordance with the interacting quantum atoms (IQA) approach. The IQA energy splitting is complemented by a topological analysis of the electron density (ρ( r )) compliant with the quantum theory of atoms‐in‐molecules (QTAIM) and the calculation of electrostatic interactions by using one‐ and two‐electron integrals, thereby avoiding convergence issues inherent to a multipolar expansion. The results show that the cooperative effects of hydrogen bonding in small water clusters arise from a compromise between: 1) the deformation energy (i.e., the energy necessary to modify the electron density and the configuration of the nuclei of the isolated water molecules to those within the water clusters), and 2) the interaction energy (E_int) of these contorted molecules in (H₂O)_n. Whereas the magnitude of both deformation and interaction energies is enhanced as water molecules are added to the system, the augmentation of the latter becomes dominant when the size of the cluster is increased. In addition, the electrostatic, classic, and exchange components of E_int for a pair of water molecules in the cluster (H₂O)_n?1 become more attractive when a new H₂O unit is incorporated to generate the system (H₂O)_n with the last‐mentioned contribution being consistently the most important part of E_int throughout the hydrogen bonds under consideration. This is opposed to the traditional view, which regards hydrogen bonding in water as an electrostatically driven interaction. Overall, the trends of the delocalization indices, δ(Ω,Ω′), the QTAIM atomic charges, the topology of ρ( r ), and the IQA results altogether show how polarization, charge transfer, electrostatics, and covalency contribute to the cooperative effects of hydrogen bonding in small water clusters. It is our hope that the analysis presented in this paper could offer insight into the different intra‐ and intermolecular interactions present in hydrogen‐bonded systems.     

Continuity of strong solutions of the Reaction–Diffusion equation in initial data

This paper is concerned with the continuity of strong solutions of the Reaction–Diffusion equation with respect to initial conditions. It is proved that the solutions of the equation are continuous in H¹$H^1$ with respect to initial data provided the nonlinearity satisfies a dissipative condition with polynomial growth order p≥2$p\ge 2$.     

Mechanical properties of anodized coatings over molten aluminum alloy

A method to measure interfacial mechanical properties at high temperatures and in a controlled atmosphere has been developed to study anodized aluminum surface coatings at temperatures where the interior aluminum alloy is molten. This is the first time that the coating strength has been studied under these conditions. We have investigated the effects of ambient atmosphere, temperature, and surface finish on coating strength for samples of aluminum alloy 7075. Surprisingly, the effective Young's modulus or strength of the coating when tested in air was twice as high as when samples were tested in an inert nitrogen or argon atmosphere. Additionally, the effective Young's modulus of the anodized coating increased with temperature in an air atmosphere but was independent of temperature in an inert atmosphere. The effect of surface finish was also examined. Sandblasting the surface prior to anodization was found to increase the strength of the anodized coating with the greatest enhancement noted for a nitrogen atmosphere. Machining marks were not found to significantly affect the strength.     

New ferrocenyl-chalcones and bichalcones: Synthesis and characterization

Ferrocenyl-chalcones and their bichalcone analogues were characterized by IR and NMR spectroscopy, as well as electrochemically. Their UV–visible spectra were recorded, and the electronic transitions were assigned by time-dependent DFT calculations. The single-crystal X-ray structures were determined for two ferrocenyl bichalcones.     

A theoretical study of hydration effects on the prototropic tautomerism of selenouracils

The prototropic tautomerism of 2-, 4-selenouracil and 2,4-diselenouracil has been studied using density functional theory (DFT) methods, at the B3LYP/6-311 + G(3df,2p)//B3LYP/6-31G(d,p) level. The relative stability order of selenouracil tautomers does not resemble that of uracil tautomers, but it is similar to that of thiouracils, even though the energy gaps between the different tautomers of selenouracils are smaller than for thiouracils. The tautomerism activation barriers are high enough as to conclude that only the oxo-selenone or the diselenone structures should be found in the gas phase. The specific interaction with one water molecule reduces these barriers by a half, but still the oxo-selenone form is always the most stable tautomer. The addition of a second water molecule has a relatively small effect, as well as bulk effects, evaluated by means of a continuum-polarized model. For isolated 2- and 4-selenouracils, the more favorable tautomerization process corresponds to a hydrogen transfer towards the selenium atom, the activation barriers for transfer towards the oxygen atom being much higher. This situation changes when specific and bulk effects are included, and the latter process becomes the more favorable one. For 2,4-diselenouracil the more favorable tautomerization, in the gas phase, corresponds to the H shift from N1 to the Se atom at C2, while solvation effects favor the transfer from N3 to the Se atom at C4.     

On the Critical Effect of the Metal (Mo vs. W) on the [3+2] Cycloaddition Reaction of M3S4 Clusters with Alkynes: Insights from Experiment and Theory

Whereas the cluster [Mo₃S₄(acac)₃(py)₃]⁺ ([ 1 ]⁺, acac=acetylacetonate, py=pyridine) reacts with a variety of alkynes, the cluster [W₃S₄(acac)₃(py)₃]⁺ ([ 2 ]⁺) remains unaffected under the same conditions. The reactions of cluster [ 1 ]⁺ show polyphasic kinetics, and in all cases clusters bearing a bridging dithiolene moiety are formed in the first step through the concerted [3+2] cycloaddition between the C?C atoms of the alkyne and a Mo(μ‐S)₂ moiety of the cluster. A computational study has been conducted to analyze the effect of the metal on these concerted [3+2] cycloaddition reactions. The calculations suggest that the reactions of cluster [ 2 ]⁺ with alkynes feature ΔG^≠ values only slightly larger than its molybdenum analogue, however, the differences in the reaction free energies between both metal clusters and the same alkyne reach up to approximately 10 kcal mol^?1, therefore indicating that the differences in the reactivity are essentially thermodynamic. The activation strain model (ASM) has been used to get more insights into the critical effect of the metal center in these cycloadditions, and the results reveal that the change in reactivity is entirely explained on the basis of the differences in the interaction energies E_int between the cluster and the alkyne. Further decomposition of the E_int values through the localized molecular orbital‐energy decomposition analysis (LMO‐EDA) indicates that substitution of the Mo atoms in cluster [ 1 ]⁺ by W induces changes in the electronic structure of the cluster that result in weaker intra‐ and inter‐fragment orbital interactions.     

Multiphysics modelling of the separation of suspended particles via frequency ramping of ultrasonic standing waves

A model was developed to determine the local changes of concentration of particles and the formations of bands induced by a standing acoustic wave field subjected to a sawtooth frequency ramping pattern. The mass transport equation was modified to incorporate the effect of acoustic forces on the concentration of particles. This was achieved by balancing the forces acting on particles. The frequency ramping was implemented as a parametric sweep for the time harmonic frequency response in time steps of 0.1 s. The physics phenomena of piezoelectricity, acoustic fields and diffusion of particles were coupled and solved in COMSOL Multiphysics? (COMSOL AB, Stockholm, Sweden) following a three step approach. The first step solves the governing partial differential equations describing the acoustic field by assuming that the pressure field achieves a pseudo steady state. In the second step, the acoustic radiation force is calculated from the pressure field. The final step allows calculating the locally changing concentration of particles as a function of time by solving the modified equation of particle transport. The diffusivity was calculated as function of concentration following the Garg and Ruthven [1] equation which describes the steep increase of diffusivity when the concentration approaches saturation. However, it was found that this steep increase creates numerical instabilities at high voltages (in the piezoelectricity equations) and high initial particle concentration. The model was simplified to a pseudo one-dimensional case due to computation power limitations. The predicted particle distribution calculated with the model is in good agreement with the experimental data as it follows accurately the movement of the bands in the centre of the chamber.     

A local lagrangian analysis of passive particle advection in a gas flow field

A local analysis is performed to study the departure from passive advection by small inertial particles based on a Lagrangian framework. The analysis considers heavy particles immersed in a gaseous flow and is restricted to short times, making it relevant to the PIV technique. A necessary (but not sufficient condition) for passive particle advection of inertial particles is that the quantity Λ_maxτ_p be much smaller than one, where Λ_max is the largest modulus of the eigenvalues corresponding to the velocity gradient tensor. This allows for the inertial and passive time scales to match beyond the initial transient, and consequently for the respective trajectories to remain relatively close. Due to this important role regarding advection behavior, Λ_maxτ_p is offered as a definition of a local Stokes number, St_Λ. Since this quantity is a field quantity, it directly provides indication of when and where passive advection by particles can be expected. A departure equation is obtained in one-dimension, where the influence of initial velocity and gravity are explicitly shown. If the flow is irrotational, the higher dimensional analysis reduces to a set of decoupled one-dimensional equations acting along each respective eigenvector of the velocity gradient tensor. A similar expression is found for the case of a purely temporal flow field.     

Intermolecular Weak Interactions in HTeXH Dimers (X=O,S, Se,Te): Hydrogen Bonds,Chalcogen–Chalcogen Contacts and Chiral Discrimination

A theoretical study of the HTeXH (X=O, S, Se and Te) monomers and homodimers was carried out by means of second‐order Møller‐Plesset perturbation theory (MP2) computational methods. In the case of monomers, the isomerization energy from HTeXH to H₂Te=X and H₂X=Te (X=O, S, Se, and Te) and the rotational transition‐state barriers were obtained. Due to the chiral nature of these compounds, homo and heterochiral dimers were found. The electron density of the complexes was characterized with the atoms‐in‐molecules (AIM) methodology, finding a large variety of interactions. The charge transfer within the dimers was analyzed by means of natural bond orbitals (NBO). The density functional theory‐symmetry adapted perturbation theory (DFT‐SAPT) method was used to compute the components of the interaction energies. Hydrogen bonds and chalcogen–chalcogen interactions were characterized and their influence analyzed concerning the stability and chiral discrimination of the dimers.