Reaction mechanism and tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution: a combined Monte Carlo and quantum mechanics study

A combined Monte Carlo and quantum mechanical study was carried out to analyze the tautomeric equilibrium of 2-mercaptopyrimidine in the gas phase and in aqueous solution. Second- and fourth-order M?ller-Plesset perturbation theory calculations indicate that in the gas phase thiol (Pym-SH) is more stable than the thione (Pym-NH) by ca. 8 kcal/mol. In aqueous solution, thermodynamic perturbation theory implemented on a Monte Carlo NpT simulation indicates that both the differential enthalpy and Gibbs free energy favor the thione form. The calculated differential enthalpy is DeltaH(SH)(-->)(NH)(solv) = -1.7 kcal/mol and the differential Gibbs free energy is DeltaG(SH)(-->)(NH)(solv) = -1.9 kcal/mol. Analysis is made of the contribution of the solute-solvent hydrogen bonds and it is noted that the SH group in the thiol and NH group in the thione tautomers act exclusively as a hydrogen bond donor in aqueous solution. The proton transfer reaction between the tautomeric forms was also investigated in the gas phase and in aqueous solution. Two distinct mechanisms were considered: a direct intramolecular transfer and a water-assisted mechanism. In the gas phase, the intramolecular transfer leads to a large energy barrier of 34.4 kcal/mol, passing through a three-center transition state. The proton transfer with the assistance of one water molecule decreases the energy barrier to 17.2 kcal/mol. In solution, these calculated activation barriers are, respectively, 32.0 and 14.8 kcal/mol. The solvent effect is found to be sizable but it is considerably more important as a participant in the water-assisted mechanism than the solvent field of the solute-solvent interaction. Finally, the calculated total Gibbs free energy is used to estimate the equilibrium constant.     

Modeling protic to dipolar aprotic solvent rate acceleration and leaving group effects in S(N)2 reactions: A theoretical study of the reaction of acetate ion with ethyl halides in aqueous and dimethyl sulfoxide solutions

The S(N)2 reactions between acetate ions and ethyl chloride, ethyl bromide, and ethyl iodide in aqueous and dimethyl sulfoxide (DMSO) solutions were theoretically investigated at an ab initio second-order M?ller-Plesset perturbation level of theory for geometry optimizations and at a fourth-order M?ller-Plesset perturbation level for energy calculations. The solvent effect was included by the polarizable continuum model using the Pliego and Riveros parametrization for DMSO and the Luque et al. scale factor for the water solution. The calculated DeltaG() values of 24.9, 20.0, and 18.5 kcal mol(-1) in a DMSO solution for ethyl chloride, ethyl bromide, and ethyl iodide are in good agreement with the estimated experimental values of 22.3, 20.0, and 16.6 kcal mol(-1), respectively. In an aqueous solution, the theoretical Delta G++ barriers of 26.9, 23.1, and 22.1 kcal mol(-1) are also in good agreement with the estimated experimental values of 26.1, 25.2, and 24.7 kcal mol(-1), respectively. The present ab initio calculations are reliable to predict the absolute and relative reactivities of ethyl halides in a DMSO solution, but in the aqueous phase, the results are less accurate. The protic to dipolar aprotic solvent rate acceleration is theoretically predicted, although this effect is underestimated. We suggest that further improvement of the present results could be obtained by including liquid-phase optimization in both solvents and treating specific solvation by water molecules for the reaction in the aqueous phase.     

Model of uptake of OH radicals on nonreactive solids

A model of adsorption and recombination of OH radicals was developed for nonreactive solid surfaces of atmospheric interest. A parametrization of this heterogeneous mechanism was carried out to determine the role of the catalytic properties of these solid surfaces, taking into account the adsorption energy, defects, surface diffusion, and chemical reactions in the gas-solid interface. The uptake process was simulated for diffusion-controlled chemical reactions on the surface on the basis of Langmuir-Hinshelwood and Eley-Rideal mechanisms. Using an analytical approach and the Monte Carlo technique, we show the dependencies of the uptake probability of the heterogeneous reactions on the OH concentration and adsorption energy. The model is employed in the analysis of the empirically derived uptake coefficient for water ice, Al(2)O(3), NaCl, NH(4)NO(3), NH(4)HSO(4), and (NH(4))(2)SO(4). We found the following values for the free energy of adsorption of OH radicals: E(ice) = 7.3-7.6 kcal/mol, E(Al)(2)(O)(3) = 11-11.7 kcal/mol, E(NH)(4)(NO)(3) = 10.2 kcal/mol, E(NaCl) = 10.2 kcal/mol, E(NH)(4)(HSO)(4) = 9.8 kcal/mol, and E((NH)(4))(2)(SO)(4) = 9.8 kcal/mol. The atmospheric implications of the catalytic reactions of OH with adsorbed reactive molecules are discussed. The results of the modeling of the uptake process showed that the heterogeneous decay rate can exceed the corresponding gas-phase reaction rate under atmospheric conditions.     

A microscopic model of chemical reactions with reorganization. High-energy approximation

A dynamical model of a chemical reaction, accompanied by reorganization of the immediate environment of the isolated chemical subsystem, is proposed. The model enables studying the emergence of nonequilibrium distribution functions as a combined result of the interaction within the dynamical subsystem and the energy exchange with a subsystem of inactive degrees of freedom (thermal bath). The study is based on the quasiclassical high-energy approximation for nonadiabatic effects in the energy exchange within the dynamical subsystem, for strong and weak coupling of the oscillator mode with the thermal bath. Such an approximation allows for the important statement that nonequilibrium effects in thermal reactions are absent if the initial translational distribution along the reaction coordinate and the initial vibrational distribution in transversal degrees of freedom are Boltzmann-like with the same temperature. The results obtained in the absence of the initial equilibrium distribution have been used for interpreting the kinetics of endothermic plasmochemical reactions proceeding under nonequilibrium conditions.     

Role of solvation in the reduction of the uranyl dication by water: a density functional study

We have studied the solvation of uranyl, UO(2)(2+), and the reduced species UO(OH)(2+) and U(OH)(2)(2+) systematically using three levels of approximation: direct application of a continuum model (M1); explicit quantum-chemical treatment of the first hydration sphere (M2); a combined quantum-chemical/continuum model approach (M3). We have optimized complexes with varying numbers of aquo ligands (n = 4-6) and compared their free energies of solvation. Models M1 and M2 have been found to recover the solvation energy only partially, underestimating it by approximately 100 kcal/mol or more. With our best model M3, the calculated hydration free energy Delta(h)G degrees of UO(2)(2+) is about -420 kcal/mol, which shifts to about -370 kcal/mol when corrected for the expected error of the model. This value agrees well with the experimentally determined interval, -437 kcal/mol < Delta(h)G degrees < -318 kcal/mol. Complexes with 5 and 6 aquo ligands have been found to be about equally favored with models M2 and M3. The same solvation models have been applied to a two-step reduction of UO(2)(2+) by water, previously theoretically studied in the gas phase. Our results show that the solvation contribution to the reaction free energy, about 60 kcal/mol, dominates the endoergicity of the reduction.     

Ab initio study of the SN1Ar and SN2Ar reactions of benzenediazonium ion with water. On the conception of "unimolecular dediazoniation" in solvolysis reactions

The nucleophilic substitution of N2 in benzenediazonium ion 1 by one H2O molecule to form protonated phenol 2 has been studied with ab initio (RHF, MP2, QCISD(T)//MP2) and hybrid density functional (B3LYP) methods. Three mechanisms were considered: (a) the unimolecular process SN1Ar with steps 1 --> Ph+ + N2 and Ph+ + H2O --> 2, (b) the bimolecular process SN2Ar with precoordination 1 + H2O --> 1 x H2O, SN reaction 1 x H2O --> [TS]++ --> 2 x N2 and dissociation of the postcoordination complex 2 x N2 --> 2 + N2, and (c) the direct bimolecular process SN2Ar that bypasses precoordination and involves just the SN reaction 1 + H2O --> [TS]++ --> 2 + N2. The SN2Ar reactions proceed by way of a Cs symmetric SN2Ar transition state structure that is rather loose, contains essentially a phenyl cation weakly bound to N2 and OH2, and is analogous to the transition state structures of front-side nucleophilic replacement at saturated centers. In solvolysis reactions, all of these processes follow first-order kinetics, and the electronic relaxation is essentially the same. It is argued that "unimolecular dediazoniations" have to proceed by way of SN2Ar transition state structures because strict SN1Ar reactions cannot be realized in solvolyses, despite the fact that the Gibbs free energy profile favors the strict SN1Ar process over the SN2Ar reaction by 6.7 kcal/mol. It is further argued that the direct SN2Ar process is the best model for the solvolysis reaction for dynamic reasons, and its Gibbs free energy of activation is 19.3 kcal/mol and remains higher than the SN1Ar value. Even though the SN1Ar and SN2Ar models provide activation enthalpies and SKIE values that closely match the experimental data, the analysis leads us to the unavoidable conclusion that this agreement is fortuitous. While the experiments do show that the solvent effect on the activation energy is about the same for all solvents, they do not show the absence of a solvent effect. The ab initio results presented here suggest that the solvent effect on the direct SN2Ar dediazoniation is approximately 12 kcal/mol, and computation of solvent effects with the isodensity polarized continuum model (IPCM) support this conclusion.     

Dissociation dynamics of thiolactic acid at 193 nm: detection of the nascent OH product by laser-induced fluorescence

Electronically excited thiolactic acid (2-mercaptopropionic acid), H(3)C-CH(SH)-COOH, undergoes the C-OH bond cleavage on excitation to the S(2) state at 193 nm, generating the primary product OH (v,J), which is detected by laser-induced fluorescence technique in a collisionless condition of flow system. The partitioning of the available energy between vibrational, rotational, and translational degrees of freedom of nascent photofragments is obtained from relative intensities of ro-vibronic lines in laser-induced fluorescence spectrum of OH, and their Doppler profiles. The rotational population of OH (v(")=0) is characterized by rotational temperature of 408+/-25 K. OH is produced in a vibrationally cold state, i.e., mostly in v(")=0. The average translational energy of OH (v(")=0,J(")) is found to be 21.5+/-2.0 kcal/mol, which implies 25.6 kcal/mol of energy in relative translation of photoproducts corresponding to the f(t) value of approximately 0.6. The observed high translational energy is due to the presence of a barrier in the exit channel, implying that the C-OH bond scission takes place on an electronically excited potential energy surface. The observed partitioning of the available energy between various degrees of the photofragments is theoretically modeled, and the hybrid model, with 26.0 kcal/mol of barrier in the exit channel, is found to explain the measured data quite well. The experimental results are also supported with ab initio molecular orbital calculations for both the ground and the excited electronic states. Time-dependent density functional theory is used to understand the nature of various electronic transitions connecting the lower excited states. Potential energy curves as a function of the C-OH bond length of thiolactic acid suggest distinct exit barriers in the S(1), T(1), and T(2) states. But, we could locate the transition state structure for OH formation in the S(1) state alone. Thus, although thiolactic acid is excited to the S(2) state at 193 nm, it undergoes internal conversion to S(1) where it dissociates to yield OH. In addition to the OH channel from excited electronic states, we studied theoretically all probable dissociation channels occurring on the ground electronic state of thiolactic acid.     

Dichlorocarbene addition to cyclopropenes: a computational study

Reaction paths for addition of dichlorocarbene to 1,2-disubstituted cyclopropenes were calculated using hybrid density functional theory (B3LYP/6-31G) in the gas phase and in the presence of a continuum solvation model corresponding to acetonitrile. In both the gas phase and acetonitrile, :CCl2-cyclopropene addition follows an asymmetric, non-least-motion approach. Barriers to addition range from 0 to 2 kcal/mol. The reactions proceed in concerted fashion in both the gas phase and solution to yield 1,3-dienes or bicyclobutanes. The reaction pathway on this complex potential energy surface of this reaction appears to bifurcate, and the product distribution is believed to be controlled by reaction dynamics. At the present level of theory, there appears to be no minimum on the potential energy surface corresponding to a dipolar intermediate.     

Dynamics of H and D abstraction in the reaction of Cl atom with butane-1,1,1,4,4,4-d6

We report the primary (D-atom) and secondary (H-atom) abstraction dynamics of chlorine atom reaction with butane-1,1,1,4,4,4-d(6). The H- and D-atom abstraction channels were studied over a range of collision energies: 10.4 kcal mol(-1) and 12.9 kcal mol(-1); 5.2 kcal mol(-1) to 12.8 kcal mol(-1), respectively, using crossed molecular beam dc slice ion imaging techniques. Single photon ionization at 157 nm was used to probe the butyl radical products resulting from the H- and D-atom abstraction reactions. These two channels manifest distinct dynamics principally in the translational energy distributions, while the angular distributions are remarkably similar. The reduced translational energy distribution for the primary abstraction showed marked variation with collision energy in the backward direction, while the secondary abstraction showed this variation in the forward direction.     

A Theoretical Study of Ene Reactions in Solution: A Solution‐Phase Translational Entropy Model

Several density functional theory (DFT) methods, such as CAM‐B3LYP, M06, ωB97x, and ωB97xD, are used to characterize a range of ene reactions. The Gibbs free energy, activation enthalpy, and entropy are calculated with both the gas‐ and solution‐phase translational entropy; the results obtained from the solution‐phase translational entropies are quite close to the experimental measurements, whereas the gas‐phase translational entropies do not perform well. For ene reactions between the enophile propanedioic acid (2‐oxo‐1,3‐dimethyl ester) and π donors, the two‐solvent‐involved explicit+implicit model can be employed to obtain accurate activation entropies and free‐energy barriers, because the interaction between the carbonyl oxygen atom and the solvent in the transition state is strengthened with the formation of C?C and O?H bonds. In contrast, an implicit solvent model is adequate to calculate activation entropies and free‐energy barriers for the corresponding reactions of the enophile 4‐phenyl‐1,2,4‐triazoline‐3,5‐dione.     

Theoretical modeling of hydroxyl-radical-induced lipid peroxidation reactions

The OH-radical-induced mechanism of lipid peroxidation, involving hydrogen abstraction followed by O2 addition, is explored using the kinetically corrected hybrid density functional MPWB1K in conjunction with the MG3S basis set and a polarized continuum model to mimic the membrane interior. Using a small nonadiene model of linoleic acid, it is found that hydrogen abstraction preferentially occurs at the mono-allylic methylene groups at the ends of the conjugated segment rather than at the central bis-allylic carbon, in disagreement with experimental data. Using a full linoleic acid, however, abstraction is correctly predicted to occur at the central carbon, giving a pentadienyl radical. The Gibbs free energy for abstraction at the central C11 is approximately 8 kcal/mol, compared to 9 kcal/mol at the end points (giving an allyl radical). Subsequent oxygen addition will occur at one of the terminal atoms of the pentadienyl radical fragment, giving a localized peroxy radical and a conjugated butadiene fragment, but is associated with rather high free energy barriers and low exergonicity at the CPCM-MPWB1K/MG3S level. The ZPE-corrected potential energy surfaces obtained without solvent effects, on the other hand, display considerably lower barriers and more exergonic reactions.     

Hydration process as an activation of trans- and cisplatin complexes in anticancer treatment. DFT and ab initio computational study of thermodynamic and kinetic parameters

The thermodynamic and kinetic aspects of hydration reactions of cis-/transplatin were explored. The polarizable continuum model was used for estimation of solvent effects. Using the B3LYP/6-31+G(d) method, the structures were optimized and vibrational frequencies estimated. Interaction energies and activation barriers were determined at the CCSD(T)/6-31++G(d,p) level within the COSMO approach. An associative mechanism was assumed with a trigonal-bipyramidal structure of the transition state. Within the applied model, all the hydration reactions are slightly endothermic. The Gibbs energies of cisplatin hydration amount to 7.0 and 14.2 kcal/mol for the chloride and ammonium replacement, respectively. Analogous values for the transplatin reactions are 6.8 and 11.9 kcal/mol. The determined rate constants are by several (three to four) orders of magnitude larger for the dechlorination process than for deammination. The cisplatin dechlorination rate constant was established as 1.3 x 10(-4) s(-1) in excellent accord with the experiment.     

Separable cooperative and localized translational motions of water confined by a chemically heterogeneous environment

We report quasi-elastic neutron scattering experiments at two resolutions that probe timescales of picoseconds to nanoseconds for the hydration dynamics of water, confined in a concentrated solution of N-acetyl-leucine-methylamide (NALMA) peptides in water over a temperature range of 248 K to 288 K. The two QENS resolutions used allow for a clean separation of two observable translational components, and ultimately two very different relaxation processes, that become evident when analyzed under a combination of the jump diffusion model and the relaxation cage model. The first translational motion is a localized beta-relaxation process of the bound surface water, and exhibits an Arrhenius temperature dependence and a large activation energy of approximately 8 kcal mol(-1). The second non-Arrhenius translational component is a dynamical signature of the alpha-relaxation of more fluid water, exhibiting a glass transition temperature of approximately 116 K when fit to the Volger Fulcher Tamman functional form. These peptide solutions provide a novel experimental system for examining confinement in order to understand the dynamical transition in bulk supercooled water by removing the unwanted interface of the confining material on water dynamics.     

Structure and conformations of N-(fluoroformyl)imidosulfurous dichloride,FC(O)N=SCl2

The IR (gas) and Raman (liquid) spectra of FC(O)NSCl(2) demonstrate the presence of a conformational mixture in both phases. According to a gas electron diffraction study, the main conformer (94(8)%) possesses a syn-syn structure (C(O)F group synperiplanar with respect to the SCl(2) bisector and the C=O bond synperiplanar to the N=S bond). Quantum chemical calculations (HF, B3LYP and MP2 with 6-31G basis set, and MP2/6-311(2df)) predict a syn-anti structure for the second conformer. Analysis of the IR (gas) spectrum results in a contribution of 5(1)% of the minor form, corresponding to a Gibbs free energy difference DeltaG degrees = G degrees (syn-anti) - G degrees (syn-syn) = 1.75(15) kcal/mol. This value is reproduced very well by quantum chemical calculations, which include electron correlation effects (DeltaG degrees = 1.28-1.56 kcal/mol). The HF approximation overestimates this energy difference (DeltaG degrees = 3.24 kcal/mol).     

One- and two-dimensional translational energy distributions in the iodine-loss dissociation of 1,2-C2H4I2+ and 1,3-C3H6I2+: what does this mean?

Threshold photoelectron photoion coincidence (TPEPICO) has been used to study the sequential photodissociation reaction of internal energy selected 1,2-diiodoethane cations: C(2)H(4)I(2)(+) → C(2)H(4)I(+) + I → C(2)H(3)(+) + I + HI. In the first I-loss reaction, the excess energy is partitioned between the internal energy of the fragment ion C(2)H(4)I(+) and the translational energy. The breakdown diagram of C(2)H(4)I(+) to C(2)H(3)(+), i.e., the fractional ion abundances below and above the second dissociation barrier as a function of the photon energy, yields the internal energy distribution of the first daughter, whereas the time-of-flight peak widths yield the released translational energy in the laboratory frame directly. Both methods indicate that the kinetic energy release in the I-loss step is inconsistent with the phase space theory (PST) predicted two translational degrees of freedom, but is well-described assuming only one translational degree of freedom. Reaction path calculations partly confirm this and show that the reaction coordinate changes character in the dissociation, and it is, thus, highly anisotropic. For comparison, data for the dissociative photoionization of 1,3-diiodopropane are also presented and discussed. Here, the reaction coordinate is expected to be more isotropic, and indeed the two degrees of freedom assumption holds. Characterizing kinetic energy release distributions beyond PST is crucial in deriving accurate dissociative photoionization onset energies in sequential reactions. On the basis of both experimental and theoretical grounds, we also suggest a significant revision of the 298 K heat of formation of 1,2-C(2)H(4)I(2)(g) to 64.5 ± 2.5 kJ mol(-1) and that of CH(2)I(2)(g) to 113.5 ± 2 kJ mol(-1) at 298 K.     

Gibbs energy of aggregation: An integral thermodynamic characteristic of the self-organization of equilibrium liquid systems

The Gibbs energy of aggregation, a new integral thermodynamic characteristic of the self-organization of equilibrium liquid systems due to specific interactions such as hydrogen bonds, is introduced. This function is defined as the difference between the Gibbs energy of the associated liquid system and the Gibbs energy of the hypothetical nonideal liquid system consisting of monomers at the same temperature and pressure. General expressions for molar Gibbs energies of aggregation of a pure liquid and a binary solution, taking into account the contributions from specific and universal (dipole, dispersion) interactions, are derived in the context of the quasi-chemical approach. An expression for the dependence of the Gibbs energy of aggregation of a pure liquid on the degree of chain aggregation is derived using the athermal associated solution model. In the same approximation, an expression for the dependence of the Gibbs energy of aggregation for a binary solution of an associating component in an inert solvent on the solution’s composition and the degree of the chain association of a pure liquid is derived. Results from numerical calculations of the Gibbs energies of aggregation of a pure liquid and a binary solution are reported. A correlation is made between the Gibbs energy of aggregation of a pure liquid and the standard Gibbs energy of formation of the dimer. It is shown that the Gibbs energy of aggregation versus the standard Gibbs energy of dimerization dependence is nonlinear in the range of low degrees of aggregation.     

Importance of correlated motions on the low barrier rotational potentials of crystalline molecular gyroscopes

The energetic and structural changes taking place upon rotation of the central phenylene of 1,4-bis(3,3,3-triphenylpropynyl)benzene in the solid state were computed using molecular mechanics calculations. Pseudopolymorphic crystals of a benzene clathrate (1A) and a desolvated form (1B) were analyzed with models that account for varying degrees of freedom within the corresponding lattices. The calculated rotational barriers in a rigid lattice approximation, 78 kcal/mol for 1A and 72 kcal/mol for 1B, are about 5 times greater than those previously measured by variable-temperature 13C CPMAS NMR and quadrupolar echo 2H NMR line-shape analysis: 12.8 kcal/mol for 1A and 14.6 kcal/mol for 1B. The potential energy barriers calculated with a model that restricts whole body rotation and translational motions but allows for internal rotations give results that are near the experimental free-energy barriers. The calculated barriers for 1A and 1B are 15.5 and 16.2 kcal/mol, respectively. The differences between the rigid and partially relaxed models are attributed to the effect of correlated motions of the lattice and the rotating group, which are evident from the structural analysis of the atomic position data as a function of the dihedral angle of the rotator. The displacements of neighboring molecules near the rotary transition states for 1A and 1B can be as large as 2.7 and 1.1 A, respectively. The displacement and oscillation (C2) of interpenetrating phenyl rings from neighboring rotors proximal to the event are significant for both 1A and 1B. In addition, 6-fold (C6) benzene rotations in clathrate 1A were found to be directly correlated to the rotation of the phenylene rotator.     

Ultra-fast rotors for molecular machines and functional materials via halogen bonding: crystals of 1,4-bis(iodoethynyl)bicyclo[2.2.2]octane with distinct gigahertz rotation at two sites

As a point of entry to investigate the potential of halogen-bonding interactions in the construction of functional materials and crystalline molecular machines, samples of 1,4-bis(iodoethynyl)bicyclo[2.2.2]octane (BIBCO) were synthesized and crystallized. Knowing that halogen-bonding interactions are common between electron-rich acetylenic carbons and electron-deficient iodines, it was expected that the BIBCO rotors would be an ideal platform to investigate the formation of a crystalline array of molecular rotors. Variable temperature single crystal X-ray crystallography established the presence of a halogen-bonded network, characterized by lamellarly ordered layers of crystallographically unique BIBCO rotors, which undergo a reversible monoclinic-to-triclinic phase transition at 110 K. In order to elucidate the rotational frequencies and the activation parameters of the BIBCO molecular rotors, variable-temperature (1)H wide-line and (13)C cross-polarization/magic-angle spinning solid-state NMR experiments were performed at temperatures between 27 and 290 K. Analysis of the (1)H spin-lattice relaxation and second moment as a function of temperature revealed two dynamic processes simultaneously present over the entire temperature range studied, with temperature-dependent rotational rates of k(rot) = 5.21 × 10(10) s(-1)·exp(-1.48 kcal·mol(-1)/RT) and k(rot) = 8.00 × 10(10) s(-1)·exp(-2.75 kcal·mol(-1)/RT). Impressively, these correspond to room temperature rotational rates of 4.3 and 0.8 GHz, respectively. Notably, the high-temperature plastic crystalline phase I of bicyclo[2.2.2]octane has a reported activation energy of 1.84 kcal·mol(-1) for rotation about the 1,4 axis, which is 24% larger than E(a) = 1.48 kcal·mol(-1) for the same rotational motion of the fastest BIBCO rotor; yet, the BIBCO rotor has three fewer degrees of translational freedom and two fewer degrees of rotational freedom! Even more so, these rates represent some of the fastest engineered molecular machines, to date. The results of this study highlight the potential of halogen bonding as a valuable construction tool for the design and the synthesis of amphidynamic artificial molecular machines and suggest the potential of modulating properties that depend on the dielectric behavior of crystalline media.     

The role of long-range forces in the determination of translational kinetic energy release. Loss of C4H4+ from the benzene and pyridine cations

Kinetic energy release distributions (KERDs) for the benzene ion fragmenting into C 4H 4 (+) and C 2H 2 have been recorded by double-focusing mass spectrometry in the metastable energy window and by a retarding field experiment up to an energy of 5 eV above the fragmentation threshold. They are compared with those resulting from the HCN loss reaction from the pyridine ion. Both reactions display a similar variation of the kinetic energy release as a function of the internal energy: the average release is smaller than statistically expected, with a further restriction of the phase-space sampling for the C 5H 5N (+) dissociation. Ab initio calculations of the potential-energy profile have been carried out. They reveal a complicated reaction mechanism, the last step of which consists in the dissociation of a weakly bound ion-quadrupole or ion-dipole complex. The KERDs have been analyzed by the maximum entropy method. The fraction of phase space effectively sampled by the pair of fragments has been determined and is similar for both dissociations. Both reactions are constrained by the square root of the released translational energy, epsilon (1/2). This indicates that in the latter stage of the dissociation process, the reaction coordinate is adiabatically decoupled from the bath of the bound degrees of freedom. For the C 6H 6 (+) fragmentation, the analysis of the experimental results strongly suggests that, just as for a spherically symmetric interaction potential, the translational motion is confined to a plane. For the dissociation of the pyridine ion, the main dynamical constraint is also a restriction to a two-dimensional subspace. This dimensionality reduction of the translational phase space is due to the fact that the Hamiltonian of both weakly bound complexes contains a cyclic coordinate.