Electron and nuclear dynamics of molecular clusters in ultraintense laser fields. III. Coulomb explosion of deuterium clusters

In this paper we present a theoretical and computational study of the energetics and temporal dynamics of Coulomb explosion of molecular clusters of deuterium (D2)n/2 (n = 480 - 7.6 x 10(4), cluster radius R0 = 13.1 - 70 A) in ultraintense laser fields (laser peak intensity I = 10(15) - 10(20)W cm(-2)). The energetics of Coulomb explosion was inferred from the dependence of the maximal energy EM and the average energy Eav of the product D+ ions on the laser intensity, the laser pulse shape, the cluster radius, and the laser frequency. Electron dynamics of outer cluster ionization and nuclear dynamics of Coulomb explosion were investigated by molecular dynamics simulations. Several distinct laser pulse shape envelopes, involving a rectangular field, a Gaussian field, and a truncated Gaussian field, were employed to determine the validity range of the cluster vertical ionization (CVI) approximation. The CVI predicts that Eav, EM proportional to R0(2) and that the energy distribution is P(E) proportional to E1/2. For a rectangular laser pulse the CVI conditions are satisfied when complete outer ionization is obtained, with the outer ionization time toi being shorter than both the pulse width and the cluster radius doubling time tau2. By increasing toi, due to the increase of R0 or the decrease of I, we have shown that the deviation of Eav from the corresponding CVI value (Eav(CVI)) is (Eav(CVI) - Eav)/Eav(CVI) approximately (toi/2.91tau2)2. The Gaussian pulses trigger outer ionization induced by adiabatic following of the laser field and of the cluster size, providing a pseudo-CVI behavior at sufficiently large laser fields. The energetics manifest the existence of a finite range of CVI size dependence, with the validity range for the applicability of the CVI being R0 < or = (R0)I, with (R0)I representing an intensity dependent boundary radius. Relating electron dynamics of outer ionization to nuclear dynamics for Coulomb explosion induced by a Gaussian pulse, the boundary radius (R0)I and the corresponding ion average energy (Eav)I were inferred from simulations and described in terms of an electrostatic model. Two independent estimates of (R0)I, which involve the cluster size where the CVI relation breaks down and the cluster size for the attainment of complete outer ionization, are in good agreement with each other, as well as with the electrostatic model for cluster barrier suppression. The relation (Eav)I proportional to (R0)I(2) provides the validity range of the pseudo-CVI domain for the cluster sizes and laser intensities, where the energetics of D+ ions produced by Coulomb explosion of (D)n clusters is optimized. The currently available experimental data [Madison et al., Phys. Plasmas 11, 1 (2004)] for the energetics of Coulomb explosion of (D)n clusters (Eav = 5 - 7 keV at I = 2 x 10(18) W cm(-2)), together with our simulation data, lead to the estimates of R0 = 51 - 60 A, which exceed the experimental estimate of R0 = 45 A. The predicted anisotropy of the D+ ion energies in the Coulomb explosion at I = 10(18) W cm(-2) is in accord with experiment. We also explored the laser frequency dependence of the energetics of Coulomb explosion in the range nu = 0.1 - 2.1 fs(-1) (lambda = 3000 - 140 nm), which can be rationalized in terms of the electrostatic model.     

Fragmentation channels of large multicharged clusters

We address unifying features of fragmentation channels driven by long-range Coulomb or pseudo-Coulomb forces in clusters, nuclei, droplets, and optical molasses. We studied the energetics, fragmentation patterns, and dynamics of multicharged (A+)n (n=55, 135, 321) clusters. In Morse clusters the variation of the range of the pair-potential induced changes in the cluster surface energy and in the fissibility parameter X=E(Coulomb)2E(surface). X was varied in the range of X=1-8 for short-range interactions and of X=0.1-1.0 for long-range interactions. Metastable cluster configurations were prepared by vertical ionization of the neutral clusters and by subsequent structural equilibration. The energetics of these metastable ionic clusters was described in terms of the liquid drop model, with the coefficients of the volume and surface energies depending linearly on the Morse band dissociation energy. Molecular-dynamics simulations established two distinct fragmentation patterns of multicharged clusters that involve cluster fission into a small number of large, multicharged clusters for X<1 and Coulomb explosion into a large number of individual ions and small ionic fragments for X>1. The Rayleigh instability limit X=1 separates between spatially anisotropic fission and spatially isotropic Coulomb explosion. Distinct features of the fragmentation energetics and dynamics were unveiled. For fission of n=55 clusters, large kinetic and internal energies of the large fragments are exhibited and the characteristic fragmentation time is approximately 700 fs, while for Coulomb explosion the major energy content of the small fragments involves kinetic energy and the characteristic fragmentation time of approximately 300 fs is shorter. The Rayleigh (X=1) limit, leading to isotropic Coulomb explosion, is transcended by a marked enhancement of the Coulomb energy, which is realized for extremely ionized clusters in ultraintense laser fields, or by a dramatic reduction of the surface energy as is the case for the expansion of optical molasses.     

Intracluster ion-molecule reactions of Ti+ with C2H5OH and CF3CH2OH clusters: influence of fluorine substituents on chemical reactivity

A laser ablation-molecular beam/reflectron time-of-flight mass spectrometric technique was used to investigate the ion-molecule reactions that proceed within Ti+(ROH)n (R = C2H5, CF3CH2) heterocluster ions. The mass spectra exhibit a major sequence of cluster ions with the formula Ti+(OR)m(ROH)n (m = 1, 2), which is attributed to sequential insertions of Ti+ into the O-H bond of C2H5OH or CF3CH2OH molecules within the heteroclusters, followed by H eliminations. The TiO+ and TiOH+ ions produced from the reactions of Ti+ with C2H5OH are interpreted as arising from insertion of Ti+ into the C-O bond, followed by C2H5 and C2H6 eliminations, respectively. When Ti+ reacted with CF3CH2OH, by contrast, considerable contributions from TiFOH+, TiF2+, and TiF2OH+ ions were observed in the mass spectrum of the reaction products, indicating that F and OH abstractions are the dominant product channels. Ab initio calculations of the complex of Ti+ with 2,2,2-trifluoroethanol show that the minimum energy structure is that in which Ti+ is attached to the O atom and one of the three F atoms of 2,2,2-trifluoroethanol, forming a five-membered ring. Isotope-labeling experiments additionally show that the chemical reactivity of heterocluster ions is greatly influenced by the presence of fluorine substituents and cluster size. The reaction energetics and formation mechanisms of the observed heterocluster ions are discussed.     

Wavelength-dependent Coulomb explosion in carbon disulphide (CS2) clusters: generation of energetic multiply charged carbon and sulphur ions

A gigawatt laser-induced Coulomb explosion has been observed in carbon disulphide (CS(2)) clusters generating energetic, multiply charged [C](m+) (m = 1-4) and [S](n+) (n = 1-6) atomic ions of carbon and sulphur. The Coulomb explosion shows wavelength dependence. Comparison of these results with our earlier work shows that the polarizability and dipole moment might help in energy absorption from the laser field but they are not mandatory conditions for this low-intensity Coulomb explosion. The results show that in a field of 10(9) W/cm(2), absorption of 266 and 355 nm laser radiation by CS(2) clusters leads to multiphoton dissociation/ionization whereas at 532 nm the whole cluster explodes generating multiply charged atomic ions.     

Ionization and dissociation of CH3I in intense laser field

The ionization-dissociation of methyl iodide in intense laser field has been studied using a reflection time-of-flight mass spectrometry (RTOF-MS), at a laser intensity of < or =6.6x10(14) W/cm(2), lambda=798 nm, and a pulse width of 180 fs. With the high resolution of RTOF-MS, the fragment ions with the same M/z but from different dissociation channels are resolved in the mass spectra, and the kinetic energy releases (KERs) of the fragment ions such as I(q+) (q=1-6), CH(m) (+) (m=0-3), C(2+), and C(3+) are measured. It is found that the KERs of the fragment ions are independent of the laser intensity. The fragments CH(3) (+) and I(+) with very low KERs (<1 eV for CH(3) (+) and <0.07 eV for I(+)) are assigned to be produced by the multiphoton dissociation of CH(3)I(+). For the fragments CH(3) (+) and I(+) from CH(3)I(2+), they are produced by the Coulomb explosion of CH(3)I(2+) with the interaction from the covalent force of the remaining valence electrons. The split of the KER of the fragments produced from CH(3)I(2+) dissociation is observed experimentally and explained with the energy split of I(+)((3)P(2)) and I(+)((3)P(0,1)). The dissociation CH(3)I(3+)-->CH(3) (+)+I(2+) is caused by Coulomb explosion. The valid charge distance R(c) between I(2+) and CH(3) (+), at which enhanced ionization of methyl iodide occurs, is obtained to be 3.7 A by the measurements of the KERs of the fragments CH(3) (+) and I(2+). For the CH(3)I(n+) (n> or =3), the KERs of the fragment ions CH(3) (p+) and I(q+) are attributed to the Coulomb repulsion between CH(3) (p+) and I(q+) from R(c) approximately 3.7 A. The dissociation of the fragment CH(3) (+) is also discussed. By the enhanced ionization mechanism and using the measured KER of I(q+), all the possible Coulomb explosion channels are identified. By comparing the abundance of fragment ions in mass spectrum, it is found that the asymmetric dissociation channels with more charges on iodine, q>p, are the dominant channels.     

Velocity map imaging and theoretical study of the Coulomb explosion of CH3I under intense femtosecond IR pulses

The Coulomb explosion of CH(3)I in an intense (10-100 TW cm(-2)), ultrashort (50 fs) and nonresonant (804 nm) laser field has been studied experimentally and justified theoretically. Ion images have been recorded using the velocity map imaging (VMI) technique for different singly and multiply charged ion fragments, CH(3)(p+) (p = 1) and I(q+) (q ≤ 3), arising from different Coulomb explosion channels. The fragment kinetic energy distributions obtained from the measured images for these ion fragments show significantly lower energies than those expected considering only Coulomb repulsion forces. The experimental results have been rationalized in terms of one-dimensional wave packet calculations on ab initio potential energy curves of the different multiply charged species. The calculations reveal the existence of a potential energy barrier due to a bound minimum in the potential energy curve of the CH(3)I(2+) species and a strong stabilization with respect to the pure Coulombic repulsion for the higher charged CH(3)I(n+) (n = 3, 4) species.     

Tuneable intramolecular intermetallic interactions as a new tool for programming linear heterometallic 4f-4f complexes

Statistical mechanics predicts that the design of pure organized heteropolymetallic chains of metal ions bound to linear receptors depends on controlled deviations from the mixing rule DeltaE(MiMj) = 1/2 (DeltaE(MiMi) + DeltaE(MjMj)), whereby DeltaE(MiMj) is the intramolecular intermetallic interaction between neighboring metal i and metal j along the receptor. A thorough investigation of linear polymetallic trivalent lanthanide triple-stranded helicates shows that such deviations are amplified by an increase in the nuclearity of the final complexes and are thus easily evidenced in the tetranuclear heterobimetallic helicates [La(4-y)Lu(y)(L6)3](12+) (y = 0-4). The chemical and physical origins of this unprecedented behavior are discussed together with its practical consequences for programming pure heteropolymetallic 4f-4f complexes.     

The anomalous shape of the cross section for the formation of SF3+ fragment ions produced by electron impact on SF6 revisited

The partial ionization cross section for the formation of SF(3) (+) fragment ions following electron impact on SF(6) is known to have a pronounced structure in the cross section curve slightly above 40 eV. We used the mass-analyzed ion kinetic energy (MIKE) scan technique to demonstrate the presence of a channel contributing to the SF(3) (+) partial ionization cross section that we attribute to the Coulomb explosion of doubly charged metastable SF(4) (2+) ions into two singly charged ions SF(3) (+) and F(+), with a threshold energy of about 45.5 eV. Thus the observed unusual shape of the SF(3) (+) partial ionization cross section is the result of two contributions, (i) the direct formation of SF(3) (+) fragment ions via dissociative ionization of SF(6) with a threshold energy of 22 eV and (ii) the Coulomb explosion of metastable SF(4) (2+) ions with a threshold energy of about 45.5 eV. A detailed analysis of the MIKE spectrum reveals an average kinetic energy release of about 5 eV in the Coulomb explosion of the SF(4) (2+) ions with evidence of a second channel corresponding to an average kinetic energy release of about 1.1 eV.     

Density functional theory study of 11-atom germanium clusters: effect of electron count on cluster geometry

Density functional theory (DFT) at the hybrid B3LYP level has been applied to the germanium clusters Ge(11)(z) (z = -6, -4, -2, 0, +2, +4, +6) starting from eight different initial configurations. The global minimum within the Ge(11)(2-) set is an elongated pentacapped trigonal prism distorted from D(3)(h) to C(2v) symmetry. However, the much more spherical edge-coalesced icosahedron, also of C(2v) symmetry, expected by the Wade-Mingos rules for a 2n + 2 skeletal electron system and found experimentally in B(11)H(11)(2-) and isoelectronic carboranes, is of only slightly higher energy (+5.2 kcal/mol). Even more elongated D(3)(h) pentacapped trigonal prisms are the global minima for the electron-rich structures Ge(11)(4-) and Ge(11)(6-). For Ge(11)(4-) the C(5v) 5-capped pentagonal antiprism analogous to the dicarbollide ligand C(2)B(9)H(11)(2-) is of significantly higher energy (approximately 28 kcal/mol) than the D(3h) global minimum. The C(2v) edge-coalesced icosahedron is also the global minimum for the electron-poor Ge(11) similar to its occurrence in experimentally known 11-vertex "isocloso" metallaboranes of the type (eta(6)-arene)RuB(10)H(10). The lowest energy polyhedral structures computed for the more hypoelectronic Ge(11)(4+) and Ge(11)(6+) clusters are very similar to those found experimentally for the isoelectronic ions E(11)(7-) (E = Ga, In, Tl) and Tl(9)Au(2)(9-) in intermetallics in the case of Ge(11)(4+) and Ge(11)(6+), respectively. These DFT studies predict an interesting D(5h) centered pentagonal prismatic structure for Ge(11)(2+) and isoelectronic metal clusters.     

Kinetic energy distribution of multiply charged ions in Coulomb explosion of Xe clusters

We report on the calculations of kinetic energy distribution (KED) functions of multiply charged, high-energy ions in Coulomb explosion (CE) of an assembly of elemental Xe(n) clusters (average size (n) = 200-2171) driven by ultra-intense, near-infrared, Gaussian laser fields (peak intensities 10(15) - 4?×?10(16) W cm(-2), pulse lengths 65-230 fs). In this cluster size and pulse parameter domain, outer ionization is incomplete∕vertical, incomplete∕nonvertical, or complete∕nonvertical, with CE occurring in the presence of nanoplasma electrons. The KEDs were obtained from double averaging of single-trajectory molecular dynamics simulation ion kinetic energies. The KEDs were doubly averaged over a log-normal cluster size distribution and over the laser intensity distribution of a spatial Gaussian beam, which constitutes either a two-dimensional (2D) or a three-dimensional (3D) profile, with the 3D profile (when the cluster beam radius is larger than the Rayleigh length) usually being experimentally realized. The general features of the doubly averaged KEDs manifest the smearing out of the structure corresponding to the distribution of ion charges, a marked increase of the KEDs at very low energies due to the contribution from the persistent nanoplasma, a distortion of the KEDs and of the average energies toward lower energy values, and the appearance of long low-intensity high-energy tails caused by the admixture of contributions from large clusters by size averaging. The doubly averaged simulation results account reasonably well (within 30%) for the experimental data for the cluster-size dependence of the CE energetics and for its dependence on the laser pulse parameters, as well as for the anisotropy in the angular distribution of the energies of the Xe(q+) ions. Possible applications of this computational study include a control of the ion kinetic energies by the choice of the laser intensity profile (2D∕3D) in the laser-cluster interaction volume.     

On the viability of small endohedral hydrocarbon cage complexes: X@C4H4, X@C8H8, X@C8H14, X@C10H16, X@C12H12, AND X@C16H16

Small hydrocarbon complexes (X@cage) incorporating cage-centered endohedral atoms and ions (X = H(+), H, He, Ne, Ar, Li(0,+), Be(0,+,2+), Na(0,+), Mg(0,+,2+)) have been studied at the B3LYP/6-31G(d) hybrid HF/DFT level of theory. No tetrahedrane (C(4)H(4), T(d)()) endohedral complexes are minima, not even with the very small hydrogen atom or beryllium dication. Cubane (C(8)H(8), O(h)()) and bicyclo[2.2.2]octane (C(8)H(14), D(3)(h)()) minima are limited to encapsulating species smaller than Ne and Na(+). Despite its intermediate size, adamantane (C(10)H(16), T(d)()) can enclose a wide variety of endohedral atoms and ions including H, He, Ne, Li(0,+), Be(0,+,2+), Na(0,+), and Mg(2+). In contrast, the truncated tetrahedrane (C(12)H(12), T(d)()) encapsulates fewer species, while the D(4)(d)() symmetric C(16)H(16) hydrocarbon cage (see Table of Contents graphic) encapsulates all but the larger Be, Mg, and Mg(+) species. The host cages have more compact geometries when metal atoms, rather than cations, are inside. This is due to electron donation from the endohedral metals into C-C bonding and C-H antibonding cage molecular orbitals. The relative stabilities of endohedral minima are evaluated by comparing their energies (E(endo)) to the sum of their isolated components (E(inc) = E(endo) - E(cage) - E(x)) and to their exohedral isomer energies (E(isom) = E(endo) - E(exo)). Although exohedral binding is preferred to endohedral encapsulation without exception (i.e., E(isom) is always exothermic), Be(2+)@C(10)H(16) (T(d)(); -235.5 kcal/mol), Li(+)@C(12)H(12) (T(d)(); 50.2 kcal/mol), Be(2+)@C(12)H(12) (T(d)(); -181.2 kcal/mol), Mg(2+)@C(12)H(12) (T(d)(); -45.0 kcal/mol), Li(+)@C(16)H(16) (D(4)(d)(); 13.3 kcal/mol), Be(+)@C(16)H(16) (C(4)(v)(); 31.8 kcal/mol), Be(2+)@C(16)H(16) (D(4)(d)(); -239.2 kcal/mol), and Mg(2+)@C(16)H(16) (D(4)(d)(); -37.7 kcal/mol) are relatively stable as compared to experimentally known He@C(20)H(20) (I(h)()), which has an E(inc) = 37.9 kcal/mol and E(isom) = -35.4 kcal/mol. Overall, endohedral cage complexes with low parent cage strain energies, large cage internal cavity volumes, and a small, highly charged guest species are the most viable synthetic targets.     

Stepwise hapticity changes in sequential one-electron redox reactions of indenyl-molybdenum complexes: combined electrochemical, ESR, X-ray, and theoretical studies.

Reduction of the dication [(eta5-Ind)(Cp)Mo[P(OMe)3]2]2+ (1(2+)) and oxidation of the neutral complex (eta3-Ind)(Cp)Mo[P(OMe)3]2 (1) proceed through a one-electron intermediate, 1+. The structures of 1(2+) and 1 have been determined by X-ray diffraction studies, which show the slip-fold distortion angle, Omega, of the indenyl ring increasing from 4.1 degrees in 1(2+) to 21.7 degrees in 1. Cyclic voltammetry and bulk electrolysis were employed to define the thermodynamics and heterogeneous charge-transfer kinetics of reactions 1(2+) + e(-) <==> 1+ and 1+ + e(-) <==> 1: DeltaE1/2 = 113 mV in CH3CN and 219 mV in CH2Cl2/0.1 M [NBu4][PF6]; k(s) = 0.4 cm x s(-1) for 1(2+)/1+ couple, 1.0 cm x s(-1) for 1+/1 couple in CH3CN. ESR spectra of 1+ displayed a surprisingly large hyperfine splitting (7.4 x 10(-4) x cm(-1)) from a single 1H nucleus, and spectra of the partially deuterated indenyl analogue confirmed assignment of a(H) to the H2 proton of the indenyl ring. The related eta5 18-electron complexes [(eta5-Ind)(Cp)Mo(dppe)]2+ (2(2+)) (dppe = diphenylphosphinoethane) and (eta5-Ind)(Cp)Mo(CN)2 (3) may also be reduced in two successive one-electron steps; ESR spectra of the radicals 2+ and 3- showed a similarly large a(H2) (8.7 x 10(-4) and 6.4 x 10(-4) x cm(-1), respectively). Molecular orbital calculations (density functional theory, DFT, and extended Hückel, EH) predict metal-indenyl bonding in 1+ that is approximately midway between that of the eta5 and eta3 hapticities (e.g., Omega = 11.4 degrees ). DFT results show that the large value of a(H2) arises from polarization of the indenyl-H2 by both inner-sphere orbitals and the singly occupied molecular orbital (SOMO) of 1+. The measured ks values are consistent with only minor inner-sphere reorganizational energies being necessary for the electron-transfer reactions, showing that a full eta5/eta3 hapticity change may require only small inner-sphere reorganization energies when concomitant with a pair of stepwise one-electron-transfer processes. The indenyl ligand in 1+ is best described as donating approximately four pi-electrons to Mo by combining a traditional eta3 linkage with two "half-strength" Mo-C bonds.     

Visualizing hydrogen atoms migrating in acetylene dication by time-resolved three-body and four-body Coulomb explosion imaging

The visualization of ultrafast isomerization of deuterated acetylene dication (C(2)D(2)(2+)) is demonstrated by time-resolved Coulomb explosion imaging with sub-10 fs intense laser pulses (9 fs, 0.13 PW cm(-2), 800 nm). The Coulomb explosion imaging monitoring the three-body explosion process, C(2)D(2)(3+)→ D(+) + C(+) + CD(+), as a function of the delay between the pump and probe pulses revealed that the migration of a deuterium atom proceeds in a recurrent manner; One of the deuterium atoms first shifts from one carbon site to the other in a short timescale (～90 fs), and then migrates back to the original carbon site by 280 fs, in competition with the molecular dissociation. Correlated motion of the two deuterium atoms associated with the hydrogen migration and structural deformation to non-planar geometry are identified by the time-resolved four-body Coulomb explosion imaging, C(2)D(2)(4+)→ D(+) + C(+) + C(+) + D(+).     

DFT calculations of 57Fe Mössbauer isomer shifts and quadrupole splittings for iron complexes in polar dielectric media: applications to methane monooxygenase and ribonucleotide reductase

To predict the isomer shifts of Fe complexes in different oxidation and spin states more accurately, we have performed linear regression between the measured isomer shifts (delta(exp)) and DFT (PW91 potential with all-electron triple-zeta plus polarization basis sets) calculated electron densities at Fe nuclei [rho(0)] for the Fe(2+,2.5+) and Fe(2.5+,3+,3.5+,4+) complexes separately. The geometries and electronic structures of all complexes in the training sets are optimized within the conductor like screening (COSMO) solvation model. Based on the linear correlation equation delta(exp) = alpha[rho(0) - 11884.0] + C, the best fitting for 17 Fe(2+,2.5+) complexes (totally 31 Fe sites) yields alpha = -0.405 +/- 0.042 and C = 0.735 +/- 0.047 mm s(-1). The correlation coefficient is r = -0.876 with a standard deviation of SD = 0.075 mm s(-1). In contrast, the linear fitting for 19 Fe(2.5+,3+,3.5+,4+) complexes (totally 30 Fe sites) yields alpha = -0.393 +/- 0.030 and C = 0.435 +/- 0.014 mm s(-1), with the correlation coefficient r = -0.929 and a standard deviation SD = 0.077 mm s(-1). We provide a physical rationale for separating the Fe(2+,2.5+) fit from the Fe(2.5+,3+,3.5+,4+) fit, which also is clearly justified on a statistical empirical basis. Quadrupole splittings have also been calculated for these systems. The correlation between the calculated (DeltaE(Q(cal))) and experimental (DeltaE(Q(exp))) quadrupole splittings based on |DeltaE(Q(exp))| = A |DeltaE(Q(cal))| + B yields slope A, which is almost the ideal value 1.0 (A = 1.002 +/- 0.030) and intercept B almost zero (B = 0.033 +/- 0.068 mm s(-1)). Further calculations on the reduced diferrous and oxidized diferric active sites of class-I ribonucleotide reductase (RNR) and the hydroxylase component of methane monooxygenase (MMOH), and on a mixed-valent [(tpb)Fe3+(mu-O)(mu-CH3CO2)Fe4+(Me3[9]aneN3)]2+ (S = 3/2) complex and its corresponding diferric state have been performed. Calculated results are in very good agreement with the experimental data.     

Density functional theory-based prediction of some aqueous-phase chemistry of superheavy element 111. Roentgenium(I) is the "softest" metal ion

A previous approach (Hancock, R. D.; Bartolotti, L. J. Inorg. Chem. 2005, 44, 7175) using DFT calculations to predict log K1 (formation constant) values for complexes of NH3 in aqueous solution was used to examine the solution chemistry of Rg(I) (element 111), which is a congener of Cu(I), Ag(I), and Au(I) in Group 1B. Rg(I) has as its most stable presently known isotope a t(1/2) of 3.6 s, so that its solution chemistry is not easily accessible. LFER (Linear free energy relationships) were established between DeltaE(g) calculated by DFT for the formation of monoamine complexes from the aquo ions in the gas phase, and DeltaG(aq) for the formation of the corresponding complexes in aqueous solution. For M2+, M3+, and M4+ ions, the gas-phase reaction was [M(H2O)6]n+(g) + NH3(g) = [M(H2O)5NH3]n+(g) + H2O(g) (1), while for M+ ions, the reaction was [M(H2O)2]+(g) + NH3(g) = [M(H2O)NH3]+(g) + H2O(g) (2). A value for DeltaG(aq) and for DeltaE for the formation of M = Cu2+ in reaction 1, not obtained previously, was calculated by DFT and shown to correlate well with the LFER obtained previously for other M2+ ions, supporting the LFER approach used here. The simpler use of DeltaE values instead of DeltaG(aq) values calculated by DFT for formation of monoamine complexes in the gas phase leads to LFER as good as the DeltaG-based correlations. Values of DeltaE were calculated by DFT to construct LFER with M+ = H+, and the Group 1B metal ions Cu+, Ag+, Au+, and Rg+, and with L = NH3, H2S, and PH3 in reaction 3: [M(H2O)2]+(g) + L(g) = [M(H2O)L]+g) + H2O(g) (3). Correlations involving DeltaE calculated by DMol3 for H+, Cu+, Ag+, and Au+ could reliably be used to construct LFER and estimate unknown log K1 values for Rg(I) complexes of NH3, PH3, and H2S calculated using the ADF (Amsterdam Density Functional) code. Log K1 values for Rg(I) complexes are predicted that suggest the Rg(I) ion to be a very strong Lewis acid that is extremely "soft" in the Pearson hard and soft acids and bases sense.     

Luminescent piano-stool complexes incorporating 1-(4-cyanophenyl)imidazole: synthesis, spectral, and structural studies

Three novel luminescent piano-stool arene ruthenium complexes of general formula [(eta(6)-arene)RuCl(2)(CPI)] (eta(6)-arene = benzene, 1, p-cymene, 2, and hexamethylbenzene, 3; CPI=1-(4-cyanophenyl)imidazole were prepared. The molecular structures of 2 and 3 were determined crystallographically. Reaction of 1-3 with EPh(3) (E = P, As, or Sb) and N-N donor bases such as 2,2'-bipyridine and 1,10-phenanthroline afforded cationic mononuclear complexes of general formula [(eta(6)-arene)RuCl(CPI)(EPh(3))](+) (eta(6)-arene = C(6)H(6), E = P (1a), E = As (1b), E = Sb(1c); eta(6)-arene = C(10)H(14), E = P (2a), E = As (2b), E = Sb (2c); eta(6)-arene = C(6)Me(6), E = P (3a), E = As (3b), E = Sb (3c)) and [(eta(6)-arene)Ru(N-N)(CPI)](2+) (eta(6)-arene = C(6)H(6), N-N = bipy (1d), N-N = phen (1e); eta(6)-arene = C(10)H(14), N-N = bipy (2d), N-N = phen (2e); eta(6)-arene = C(6)Me(6), N-N = bipy (3d), N-N = phen (3e)). Molecular structures of 1a and 2a were also confirmed by X-ray crystallography. Structural studies of the complexes 2, 3, 1a, and 2a supported coordination of CPI through the imidazole nitrogen and the presence of a pendant nitrile group. Structural data also revealed stabilization of crystal packing in the complexes 2, 3, and 2a by C-H...X (X = Cl, F) type inter- and intramolecular interactions and in complex 1a by pi-pi stacking. Moreover, neutral homonuclear bimetallic complexes 2f,g were prepared by using complex 2 as a metallo-ligand, where CPI acts as a bridge between two metal centers. Emission spectra of the mononuclear complexes [(eta(6)-arene)RuCl(2)(CPI)] and its derivatives exhibited intense luminescence when excited in the metal to ligand charge-transfer band.     

Palladium-arene interactions in catalytic intermediates: an experimental and theoretical investigation of the soft rearrangement between eta1 and eta2 coordination modes

A series of dichloro-bridged arylbicycloheptylpalladium complexes have been synthesized and characterized by means of NMR spectroscopy. The compound [(C16H19)PdCl]2*CH2Cl2 with ortho and para methyl substituents at the arene has been characterized by means of X-ray diffraction techniques. The C(ipso) atom of the arene lies almost at the fourth planar coordination site of the metal [Pd-C(ipso) = 2.22(1) A (average)], and due to the arene's tilting, the substituted C(ortho) atom is relatively close to the metal atom [2.54(1) A (average)]. The coordinated C(ipso)-C(ortho) linkage, in a seemingly dihapto coordination, is anti with respect to the CH2 bridge of the bicycloheptyl unit. Variable-temperature NMR experiments for the para-substituted dimer 9 reveal restricted rotation of the two aryl groups about the corresponding C-C(ipso) bonds (DeltaE < or =17 kcal x mol(-1)). DFT-B3LYP calculations have been carried out on the known and similar monomer (phenylbicycloheptenyl)Pd(PPh3)I (4) and its related substituted derivatives. The essential results are as follows: (i) The potential energy surface for twisting the phenyl ring away from the symmetric eta(1) coordination in 4 is very flat (DeltaE < or = 1 kcal x mol(-1)) whereas an Atoms in Molecules analysis excludes the existence of an actual Pd-C(ortho) bond in the seemingly eta2-type conformer. (ii) Complete rotation of the unsubstituted phenyl ring is not facile but feasible. A significant strain affects the transition-state structure featuring a Pd-HC(aryl) agostic-type bond. The calculated destabilization of 10.3 kcal x mol(-1), with respect to the ground state, can be compared to the experimental barrier of the dimer 9. (iii) Various methyl-substituted derivatives of 4 have been optimized, and their structural and energetic trends are discussed. An almost ideal eta1 coordination is shown by the anti conformer of the C(ortho)-substituted complex due steric effects. For all of the other cases, a slipped eta2 coordination may be described. As a general conclusion, the unsaturated metal center receives pi electron density of the arene mainly through its C(ipso) atom. The effect may be slightly improved if the C(ortho) atom also gets closer to the metal, but in no case, does the slipped eta2 coordination seem to be crucial for the stability of the system.     

Anisotropic bulletlike emission of terminal ethynyl fragment ions: ionization of ethynylbenzene-d under intense femtosecond laser fields

The authors investigated Coulomb explosions of ethynylbenzenes under intense femtosecond laser fields. Deuteration on the edge of the triple bond gave information about specific fragment emissions and the contribution of hydrogen migration. Some fragments not resulting from migration were emitted in the direction of laser polarization. These were ethynyl fragment ions (D(+), CD(+), C(2)D(+), and C(3)D(+)). Although two bonds have to be cleaved to produce C(3)D(+), the rigid character of the triple bond was maintained in the Coulomb explosion process. In contrast, fragment ions, which are formed after single or double hydrogen migration, showed isotropic emissions with distinct kinetic energies. The character of the substituents has been found to hold even under strong laser light fields where violent fragmentation took place. The ethynyl parts were emitted like bullets from the molecular frame of ethynylbenzene despite the explosion into pieces of the main body of benzene ring.     

On the mechanism of (PCP)Ir-catalyzed acceptorless dehydrogenation of alkanes: a combined computational and experimental study

Pincer complexes of the type ((R)PCP)IrH(2), where ((R)PCP)Ir is [eta(3)-2,6-(R(2)PCH(2))(2)C(6)H(3)]Ir, are the most effective catalysts reported to date for the "acceptorless" dehydrogenation of alkanes to yield alkenes and free H(2). We calculate (DFT/B3LYP) that associative (A) reactions of ((Me)PCP)IrH(2) with model linear (propane, n-PrH) and cyclic (cyclohexane, CyH) alkanes may proceed via classical Ir(V) and nonclassical Ir(III)(eta(2)-H(2)) intermediates. A dissociative (D) pathway proceeds via initial loss of H(2), followed by C-H addition to ((Me)PCP)Ir. Although a slightly higher energy barrier (DeltaE(+ +)) is computed for the D pathway, the calculated free-energy barrier (DeltaG(+ +)) for the D pathway is significantly lower than that of the A pathway. Under standard thermodynamic conditions (STP), C-H addition via the D pathway has DeltaG(o)(+ +) = 36.3 kcal/mol for CyH (35.1 kcal/mol for n-PrH). However, acceptorless dehydrogenation of alkanes is thermodynamically impossible at STP. At conditions under which acceptorless dehydrogenation is thermodynamically possible (for example, T = 150 degrees C and P(H)2 = 1.0 x 10(-7) atm), DeltaG(+ +) for C-H addition to ((Me)PCP)Ir (plus a molecule of free H(2)) is very low (17.5 kcal/mol for CyH, 16.7 kcal/mol for n-PrH). Under these conditions, the rate-determining step for the D pathway is the loss of H(2) from ((Me)PCP)IrH(2) with DeltaG(D)(+ +) approximately DeltaH(D)(+ +) = 27.2 kcal/mol. For CyH, the calculated DeltaG(o)(+ +) for C-H addition to ((Me)PCP)IrH(2) on the A pathway is 35.2 kcal/mol (32.7 kcal/mol for n-PrH). At catalytic conditions, the calculated free energies of C-H addition are 31.3 and 33.7 kcal/mol for CyH and n-PrH addition, respectively. Elimination of H(2) from the resulting "seven-coordinate" Ir-species must proceed with an activation enthalpy at least as large as the enthalpy change of the elimination step itself (DeltaH approximately 11-13 kcal/mol), and with a small entropy of activation. The free energy of activation for H(2) elimination (DeltaG(A)(+ +)) is hence found to be greater than ca. 36 kcal/mol for both CyH and n-PrH under catalytic conditions. The overall free-energy barrier of the A pathway is calculated to be higher than that of the D pathway by ca. 9 kcal/mol. Reversible C-H(D) addition to ((R)PCP)IrH(2) is predicted to lead to H/D exchange, because the barriers for hydride scrambling are extremely low in the "seven-coordinate" polyhydrides. In agreement with calculation, H/D exchange is observed experimentally for several deuteriohydrocarbons with the following order of rates: C(6)D(6) > mesitylene-d(12) > n-decane-d(22) > cyclohexane-d(12). Because H/D exchange in cyclohexane-d(12) solution is not observed even after 1 week at 180 degrees C, we estimate that the experimental barrier to cyclohexane C-D addition is greater than 36.4 kcal/mol. This value is considerably greater than the experimental barrier for the full catalytic dehydrogenation cycle for cycloalkanes (ca. 31 kcal/mol). Thus, the experimental evidence, in agreement with calculation, strongly indicates that the A pathway is not kinetically viable as a segment of the "acceptorless" dehydrogenation cycle.