Effect of absolute laser phase on reaction paths in laser-induced chemical reactions

Potential surfaces, dipole moments, and polarizabilities are calculated by ab initio methods [unrestricted MP2(full)/6-311++G(2d,2p)] along the reaction paths of the F+CH4 and Cl+CH4 reaction systems. It is found that in general dipole moments and polarizabilities exhibit peaks near the transition state. In the case of X=F these peaks are on the products side and in the case of X=Cl they are on the reactants side indicating an early transition state in the case of fluorine and a late transition state in the case of chlorine. An analysis of the geometric changes along the reaction paths reveals a one-to-one correspondence between the peaks in the electric properties and peaks in the rate of change of certain internal geometric coordinates along the reaction path. Interaction with short infrared intense laser fields pulses leads to the possibility of interferences between the dipole and polarizability laser-molecule interactions as a function of laser phase. The larger dipole moment in the Cl+CH4 reaction can lead to the creation of deep wells (instead of energy barriers) and new strongly bound states in the transition state region. This suggests possible coherent control of the reaction path as a function of the absolute phase of the incident field, by significant modification of the potential surfaces along the reaction path and, in particular, in the transition state region.     

Hydrogen–hydrogen bonding in biphenyl revisited

The evidence for the stabilizing nature of the H–H bonding in planar biphenyl is succinctly reviewed. The stabilizing nature of the H–H bonding is revealed through a comparison of the atomic energy of every atom in planar biphenyl with the same atom in the twisted equilibrium structure. It is shown that the barrier to rotation via the planar transition state is the net resultant of a stabilisation of the four ortho-hydrogen atoms (by 8 kcal/mol each), a stabilisation of the two para-carbon atoms (by 3 kcal/mol each) and by the dominant destabilisation of the two carbon atoms joining the two rings—the two junction carbon atoms—(by 22 kcal/mol each). The energetic stabilisation of the four ortho-hydrogen atoms is further shown to be in large proportion due to the formation of the hydrogen–hydrogen interatomic surface. Furthermore, neither the “bond order” between the two junction carbon atoms nor the total electron delocalisation between the two rings exhibit a significant change in going from the planar to the twisted equilibrium geometry. These findings are in contrast with the classical view of a balance between “steric non-bonded repulsion” and better electron delocalisation as a function of the twist dihedral angle. Similar conclusions have been recently reached by Pacios and Gómez through a study of the electrostatic potential at the position of the hydrogen nuclei. We dedicate this article to Professor TM Krygowski on the occasion of his 70th birthday wishing him a long and productive life.     

Twisted amides: X-ray crystallographic and theoretical study of two acylated glycolurils with aromatic substituents

X-ray crystallography and theoretical analysis were applied to explore the molecular basis for the efficient and selective Claisen-like condensations of diacylglycolurils. The crystal structures of 1-acetyl-6-benzoyl-3,4,7,8-tetramethylglycoluril (4b), and of 1-(3′-oxo-3′-phenylpropionyl)-3,4,7,8-tetramethylglycoluril (5b), the product of base-promoted intramolecular condensation of 4b, were obtained by X-ray diffraction. The acetyl (Ac) group in 4b is essentially coplanar with the attached tetrahydroimidazolone ring of the glycoluril core (τ=7°), while the benzoyl (Bz) group is twisted by τ=45° relative to a plane through the ring to which it is bonded. Product 5b contains a flat amide (τ=7°). Ab initio energy optimizations of the experimental structures for 4b and 5b give optimized geometries which are not dramatically altered, suggesting that crystal packing effects are small. An atoms-in-molecules study of the delocalization of the Fermi hole reveals that electrons in the Bz C=O group of 4b are delocalized into the phenyl ring as well as into the urea moiety of the glycoluril core. This effect stabilizes the Bz over the Ac carbonyl group, and accounts for selective twisting of the Bz group. The Laplacian of the electron density reveals a non-bonded valence shell charge concentration at O of the Ac group, corresponding to a lone-pair region, aligned with a charge depletion in the valence shell of the Bz C=O carbon [(C15–O16C18)=113°]. The angle of approach [(O16C18=O19)] is 100°, equal to the angle for ideal nucleophilic attack on a carbonyl group. Oxygen atom O16 is thus poised to attack C18; only the O16C18 distance (3.248 Å) seems to prevent reaction. These results suggest that the same distance restraint may prevent O-acylation in the enolate intermediate 6b derived from4b. By contrast, the transition state for C-acylation, leading from 6b towards product 5b requires a different geometry, which may explain the observed selectivity for C-acylation in this enolate. The results show that, as 4b is converted to 5b, amide torsional strain is relieved, which may account for the high reactivity of 4b and the efficiency and irreversibility of this condensation process. This study provides a starting point for quantitative correlation of substrate structure in diacylglycolurils with kinetic data for the rearrangement reaction.     

Structures,g-tensors,and hyperfine coupling constants of L-α-alanine radicals in radiation dosimetry: An ab initio molecular dynamics simulation study

Alanine is used as a transfer standard dosimeter for gamma ray and electron beam calibration. An important factor affecting its dosimetric response is humidity which can lead to errors in absorbed dose calculations. Ab initio molecular dynamics calculations were performed to determine the environmental effects on the electron paramagnetic resonance (EPR) parameters of L-α-alanine radicals in acidic and alkaline solutions. A new result, not dissimilar to the closed-shell amino acid molecule alanine, is that the non-zwitterionic form of the alanine radical is the stable form in the gas phase while the zwitterionic neutral alanine radical is not a stable structure in the gas phase. Geometric and EPR parameters of radicals in both gas and solution phases are found to be dependent on hydrogen bonding of water molecules with the polar groups and on dynamic solvation. Calculations on the optimized free radicals in the gas phase revealed that for the neutral radical, hydrogen bonding to water molecules drives a decrease in the magnitudes of g-tensor components g _xx and g _yy without affecting neither g _zz component nor the hyperfine coupling constants (HFCCs). The transfer from the gas to solution phase of the alanine radical anion is accompanied with an increase in the spin density on the carboxylic group's oxygen atoms. However, for the neutral radical, this transfer from gas to solution phase is accompanied with the decrease in the spin density on oxygen atoms. Calculated isotropic HFCCs and g-tensor of all radicals are in good agreement with experiment in both acidic and alkaline solutions.     

Comparison of multicomponent fuel droplet vaporization experiments in forced convection with the Sirignano model

The vaporization of multicomponent fuel droplets was studied experimentally in a heated flow and the results were compared to the model proposed by Abramzon and Sirignano. The droplet was suspended on a permanent holder which was set up in a thermal wind-tunnel. This wind-tunnel was fitted with a video recording system and an infra-red camera. The period during which the droplet was suspended on the holder before the opening of the hot air flow damper was recorded. This first sequence corresponds to the droplet vaporization in natural convection, whose initial experiment conditions, especially diameter, temperature, composition of the droplet, are well known. Then the damper was turn on, and the sequence of forced convection begun. The initial diameter of the droplet was recorded by the video system. The other initial conditions of this second sequence cannot be determined experimentally. The distribution of temperature in the droplet and the surface temperature, the mass fraction distribution in the droplet and the surface mass fraction were unknown. These unknown parameters were determined by coupling our experiment with a model using “the film concept” in natural convection. Experimental results were compared with the calculations and found satisfactory, in natural convection as well as in forced convection initiated by this method. The method was tested in the case of a fuel mixture droplets (heptane–decane) for different initial concentrations and variable durations of the sequence in natural convection.     

The localized electrons detector as an ab initio representation of molecular structures

The local single particle momentum is proposed as a localized‐electrons detector (LED) that provides a direct three‐dimensional representation of bonding interactions in molecules. It is given exclusively in terms of the electron density and its gradient. We show that the graphical representation of bonding interactions given by LED is consistent with the local curvatures of the electron density as given by the eigenvalues of the Hessian matrix, according to a local symmetry classification of the critical points here introduced. LED consistently complements the topological analysis of the electron density given by the quantum theory of atoms in molecules, by providing a graphical representation of the symmetry of the bonding interactions in molecular systems. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2418–2425, 2010     

Remarks on a lemma by Jacques-Louis Lions

Let Ω   be a bounded and connected open subset of R^N${\mathbb{R}}^N$ with a Lipschitz-continuous boundary ∂Ω, the set Ω being locally on one side of ∂Ω  . It is shown in this Note that a fundamental characterization of the space L²(Ω)$L^2(\Omega )$ due to Jacques-Louis Lions is in effect equivalent to a variety of other properties. One of the keys for establishing these equivalences is a specific “approximation lemma”, itself one of these equivalent properties.     

The stationary Oseen equations in an exterior domain: An approach in weighted Sobolev spaces

In this work, we study the linearized Navier–Stokes equations in an exterior domain of R³${\mathbb{R}}^3$ at the steady state, that is, the Oseen equations. We are interested in the existence and the uniqueness of weak, strong and very weak solutions in L^p$L^p$-theory which makes our work more difficult. Our analysis is based on the principle that linear exterior problems can be solved by combining their properties in the whole space R³${\mathbb{R}}^3$ and the properties in bounded domains. Our approach rests on the use of weighted Sobolev spaces.     

Stability of a Rayleigh-Bénard Poiseuille flow for yield stress fluids—Comparison between Bingham and regularized models

A linear stability analysis of a Rayleigh-Bénard Poiseuille flow is performed for yield stress fluids whether we use the Bingham or regularized models. A fundamental difference between those models is that the effective viscosity is not defined in the plug zone for the Bingham model, while it is defined in the whole domain for the regularized models. For these models, the viscosity depends highly on a parameter ? near the axis and increases drastically in an intermediate region. The convergence of the critical conditions between the simple and the Bingham models is not obtained. However, we show that the Bercovier and Papanastasiou models can tend to the exact Bingham results.