Product Integration Rules at Clenshaw-Curtis and Related Points: A Robust Implementation

Product integration rules generalizing the Fej?r, Clenshaw-Curtis,and Filippi quadrature rules respectively are derived for integralswith trigonometric and hyperbolic weight factors. The Chebyshevmoments of the weight functions are found to be given by well-conditionedexpressions, in terms of hypergeometric functions ₀F₁. An a priori error estimator is discussed which is shown bothto avoid wasteful invocation of the integration rule and toincrease significantly the robustness of the automatic quadratureprocedure. Then, specializing to extended Clenshaw-Curtis (ECC) rules,three types of a posteriori error estimates are considered andthe existence of a great risk of their failure is demonstratedby large scale validation tests. An empirical error estimator,superseding them for slowly varying integrands, is found toresult in a spectacular increase in the output reliability. Finally, enhancements in the control of the interval subdivisionstrategy aiming at increasing code robustness is discussed.Comparison with the code DQAWO of QUADPACK, with about a hundredthousand solved integrals, is illustrative of the increasedrobustness and error estimate reliability of our computer codeimplementation of the ECC rules.     

Thermodynamic Complexation of Dopamine with Zinc(II) in Media with Different Dielectric Constants

The complexation of zinc(II) with dopamine has been investigated by spectrophotometric measurements in mixed solvent system at an ionic strength of 0.2 mol•dm－3 sodium chloride, employed at at (15±0.1), (25±0.1), (35±0.1) ℃ at inin a pH range of ca. 6 to ca. 7 with a high ratio of ligand to metal. The effect of solvent systems on protonation and complexation are was was discussed. Linear relationships are werewere observed by plotting lg K versus 1/D, where K and D show stability and dielectric constants, respectively.     

Conformational analysis of 18-azacrown-6 and its bonding with late first transition series divalent metals: insight from DFT combined with NPA and QTAIM analyses

Density functional theory calculations, together with quantum theory of atoms in molecules (QTAIM) analyses, have been performed to investigate 18-azacrown-6 complexes of the high-spin late first transition series divalent metal ions in the gas phase and, in some cases, in aqueous solution simulated by a polarizable continuum model. Six intramolecular H-H bonding interactions in the meso-complexes are found to arise from folding of the ligand upon its electrostatic interaction with the metal ions, which are largely absent in the lowest-energy C(2h) conformer of the free ligand. The ligand-to-metal charge transfer obtained from QTAIM analysis, among other things, is found to be an important factor that controls the stability of these complexes. The inter-relationship between the ligand preorganization energy, the zero-point corrected formation energy of the metal complexes, and the H-H bonding pair distances, as well as the dependence of the electron density and the total energy density at the H-H bond critical points on the H-H bonding pair distances, provides a physical basis for understanding and explaining the stabilizing nature of these closed-shell interactions, which are often viewed as steric clashes that lead to complex destabilization.     

Can a dipole-bound electron form a pseudo-atom? An atoms-in-molecules study of the hydrated electron

Non-nuclear local maxima, or attractors, of electron density are a rare but very interesting feature of the electron density distribution in molecules and solids. Recently, non-nuclear attractors (NNAs) and the corresponding pseudoatoms of electron density have been identified with the quantum theory of atoms in molecules for some anionic clusters formed by several polar solvent molecules and an excess electron bound in either a solvated-electron or dipole-bound fashion. This contribution reports a detailed study of the topology of the electron density for a series of dipole-bound water cluster anions, as calculated with Hartree-Fock, M?ller-Plesset perturbation theory, and coupled-cluster methods together with basis sets augmented with extra diffuse basis functions to accommodate the excess electron. For dipole-bound clusters, electron densities obtained with insufficient inclusion of electron correlation effects and tight basis sets feature a well-pronounced pseudoatom due to the excess electron, which ultimately disappears when a higher level of electronic structure theory and a more diffuse basis set are used. On the other hand, for solvated-electron clusters, where the excess electron is surrounded by solvent molecules, the existence of NNAs does not seem to be an artifact of the method employed, but rather a genuine feature of the electron density distribution. Pseudoatoms of electron density thus appear to be an exclusive feature of confined environments and are unlikely to be found on the tip of a cluster dipole or on solid surfaces.     

Resonance description of S-nitrosothiols: insights into reactivity

A resonance representation of the electronic structure of S-nitrosothiols as a combination of the conventional R-S-N=O structure, a zwitterionic structure R-S+=N-O-, and a RS-/NO+ ion pair is proposed. The resonance forms are employed to predict and rationalize the structural and conformational properties of RSNOs, their interaction with Lewis acids, and their overall chemical reactivity.     

Accurate ab initio potential for the Na+...I* complex

High-level ab initio calculations employing the multireference configuration interaction and coupled clusters methods with a correlation-consistent sequence of basis sets have been used to obtain accurate potential energy curves for the complex of the sodium cation with the iodine atom. Potential curves for the first two electronic Lambda-S states have very different characters: the potential for the 2pi state has a well depth of approximately 10 kcal/mol, while the 2sigma state is essentially unbound. This difference is rationalized in terms of the anisotropic interaction of the quadrupole moment of the iodine atom with the sodium cation, which is stabilizing in the case of the 2pi state and destabilizing in the case of the 2sigma state. The effects of spin-orbit coupling have been accounted for with both ab initio and semiempirical approaches, which have been found to give practically the same results. Inclusion of spin-orbit interactions does not affect the X(omega = 32) ground state, which retains its 2pi character, but it results in two omega = 12 spin-orbit states, with mixed 2sigma and 2pi characters and binding energies roughly half of that of the ground spin-orbit state. Complete basis set (CBS) extrapolations of potential curves, binding energies, and equilibrium geometries were also performed, and used to calculate a number of rovibronic parameters for the Na+...I* complex and to parameterize model potentials. The final CBS-extrapolated and zero-point vibrational energy-corrected binding energy is 10.2 kcal/mol. Applications of the present results for simulations of NaI photodissociation femtosecond spectroscopy are discussed.     

Lattice Boltzmann simulation of natural convection in nanofluid-filled 2D long enclosures at presence of magnetic field

In this paper, the effects of a magnetic field on natural convection flow in filled long enclosures with Cu/water nanofluid have been analyzed by lattice Boltzmann method. This study has been carried out for the pertinent parameters in the following ranges: the Rayleigh number of base fluid, Ra = 10³–10⁵, the volumetric fraction of nanoparticles between 0 and 6 %, the aspect ratio of the enclosure between A = 0.5 and 2. The Hartmann number has been varied from Ha = 0 to 90 with interval 30 while the magnetic field is considered at inclination angles of θ = 0°, 30°, 60° and 90°. Results show that the heat transfer decreases by the increment of Hartmann number for various Rayleigh numbers and the aspect ratios. Heat transfer decreases with the growth of the aspect ratio but this growth causes the effect of the nanoparticles to increase. The magnetic field augments the effect of the nanoparticles at high Rayleigh numbers (Ra = 10⁵). The effect of the nanoparticles rises for high Hartmann numbers when the aspect ratio increases. The rise in the magnetic field inclination improves heat transfer at aspect ratio of A = 0.5.     

Substituent effects and the role of negative hyperconjugation in siloxycarbene rearrangements

Substituent effects and the role of negative hyperconjugation in 1,2-silyl migration and decarbonylation of methoxy(substituted-siloxy)carbenes have been investigated using quantum chemical calculations and natural bond orbital analysis. It has been found that sigma-electron-withdrawing substituents generally lower the barriers for 1,2-silyl migration and decarbonylation, consistent with symmetry-forbidden concerted rearrangements involving intramolecular front-side nucleophilic attack by the carbene lone pair at silicon and by the methoxy oxygen at silicon, respectively. However, while good linear Hammett correlations are obtained for 1,2-silyl migration, those obtained for decarbonylation are poor. In addition, there appears to be a relationship between the extent of pertinent hyperconjugative interactions in the siloxycarbene conformers and the ease of intramolecular reactivity. As a matter of fact, the finding that 1,2-silyl migration is more favorable than decarbonylation seems to be primarily related to stronger negative hyperconjugation between the carbene lone pair and the O-Si antibonding orbital, compared to that between the methoxy oxygen n(sigma) lone pair and the O-Si antibonding orbital. Moreover, the activation enthalpies for 1,2-silyl migration decrease linearly with stronger negative hyperconjugation, although no such correlation could be established for decarbonylation.