Evaluation of data for atmospheric models: master equation/RRKM calculations on the combination reaction, BrO + NO2 --> BrONO2, a conundrum

Experimental data for the title reaction were modeled using master equation (ME)/RRKM methods based on the Multiwell suite of programs. The starting point for the exercise was the empirical fitting provided by the NASA (Sander, S. P.; Finlayson-Pitts, B. J.; Friedl, R. R.; Golden, D. M.; Huie, R. E.; Kolb, C. E.; Kurylo, M. J.; Molina, M. J.; Moortgat, G. K.; Orkin, V. L.; Ravishankara, A. R. Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies, Evaluation Number 15; Jet Propulsion Laboratory: Pasadena, California, 2006)1 and IUPAC (Atkinson, R.; Baulch, D. L.; Cox, R. A.; R. F. Hampson, J.; Kerr, J. A.; Rossi, M. J.; Troe, J. J. Phys. Chem. Ref. Data 2000, 29, 167)2 data evaluation panels, which represents the data in the experimental pressure ranges rather well. Despite the availability of quite reliable parameters for these calculations (molecular vibrational frequencies (Parthiban, S.; Lee, T. J. J. Chem. Phys. 2000, 113, 145)3 and a value (Orlando, J. J.; Tyndall, G. S. J. Phys. Chem. 1996, 100, 19398)4 of the bond dissociation energy, D298(BrO-NO2) = 118 kJ mol-1, corresponding to DeltaH0o = 114.3 kJ mol-1 at 0 K) and the use of RRKM/ME methods, fitting calculations to the reported data or the empirical equations was anything but straightforward. Using these molecular parameters resulted in a discrepancy between the calculations and the database of rate constants of a factor of ca. 4 at, or close to, the low-pressure limit. Agreement between calculation and experiment could be achieved in two ways, either by increasing DeltaH0o to an unrealistically high value (149.3 kJ mol-1) or by increasing DeltaEd, the average energy transferred in a downward collision, to an unusually large value (>5000 cm-1). The discrepancy could also be reduced by making all overall rotations fully active. The system was relatively insensitive to changing the moments of inertia in the transition state to increase the centrifugal effect. The possibility of involvement of BrOONO was tested and cannot account for the difficulties of fitting the data.     

Structural events in the photocycle of green fluorescent protein

Picosecond time-resolved mid-infrared absorption changes of the wild type green fluorescent protein from Aequorea victoria are reported on structural events during the photocycle. Concomitant with rapid H/D transfer following excitation of the neutral A state at 400 nm, a transient signal at 1721/1711 cm(-1) (H/D) developed belonging to protonated glutamate 222, which was definitively assigned using the E222D mutant from the altered proton-transfer kinetics to aspartate in addition to the altered band position and intensity in the spectra. A transient at 1697 cm(-1), assigned to a structural perturbation of glutamine 69, had a H/D kinetic isotope effect of >32, showing the conformational dynamics to be sensitive to the active site H/D vibrations. The kinetic data up to 2 ns after excitation in the 1250-1800 cm(-1) region in D2O provided 10 and 75 ps time constants for the excited-state deuteron transfer and the associated A1* - A1 and A2* - A2 difference spectra and showed the radiative intermediate I state vibrations and the transient accumulation of the long-lived ground-state intermediate I2. Assignments of chromophore modes for the A1, A2, and I2 ground states are proposed on the basis of published model compound studies (Esposito, A. P.; Schellenberg, P.; Parson, W. W.; Reid, P. J. J. Mol. Struct. 2001, 569, 25 and He, X.; Bell, A. F.; Tonge, P. J. J. Phys. Chem. B 2002, 106, 6056). Tentative assignments for the singlet-state intermediates A1*, A2*, and I* are discussed. An unexpected and unassigned band that may be a C=C chromophore vibration was observed in the ground state (1665 cm(-1)) as well as in all photocycle intermediates. Optical dumping of the transient I population was achieved using an additional 532 nm pulse and the directly obtained I2 - I* difference spectrum was highly similar to the I2 - I* photocycle spectrum. The pump-dump-probe spectrum differed from the pump-probe photocycle difference spectrum with respect to the intensity of the phenol 1 mode at 1578 cm(-1), suggesting stronger delocalization of the negative charge onto the phenolic ring of the anionic chromophore in the dumped I2 state. Indication for structural heterogeneity of the chromophore, Glu 222, and the chromophore environment was found in the two parallel proton-transfer reactions and their distinct associated ground- and intermediate-state vibrations. Vibrational spectral markers at 1695 cm(-1) assigned to Gln 69, at 1631 cm(-1) belonging to a C=C mode, and at 1354 cm(-1) belonging to a phenolate vibration further indicated the I2 and I* states to be unrelaxed.     

An ab initio investigation of the O(3P)-H2(1sigma(g)+) van der Waals well

We report an ab initio study of the van der Waals region of the O(3P)-H2 potential energy surface based on RCCSD(T) calculations with an aug-cc-pVQZ basis supplemented by bond functions. In addition, an open-shell implementation of symmetry-adapted perturbation theory (SAPT) is used to corroborate the RCCSD(T) calculations and to investigate the relative magnitudes of the various contributions to the van der Waals interaction. We also investigate the effect of the spin-orbit coupling on the position and depth of the van der Waals well. We predict the van der Waals minimum to occur in perpendicular geometry, and located at a closer distance than a secondary well in colinear geometry. The potentials obtained in the present study confirm the previous calculations of Alexander [M. H. Alexander, J. Chem. Phys., 1998, 108, 4467], but disagree with the earlier work of Harding and co-workers [Z. Li, V. A. Apkarian and L. B. Harding, J. Chem. Phys., 1997, 106, 942] as well as with recently refitted surfaces of Brand?o and coworkers [J. Brand?o, C. Mogo and B. C. Silva, J. Chem. Phys., 2004, 121, 8861]. Inclusion of spin-orbit coupling reduces the depth of the van der Waals minimum without causing a change in its position.     

Rotationally resolved C2 symmetric conformers of bis-(4-hydroxyphenyl)methane: prototypical examples of excitonic coupling in the S1 and S2 electronic states

Rotationally resolved microwave and ultraviolet spectra of jet-cooled bis-(4-hydroxyphenyl)methane (b4HPM) have been obtained using Fourier-transform microwave and UV laser/molecular beam spectrometers. A recent vibronic level study of b4HPM [Rodrigo, C. P.; Mu?ller, C. W.; Pillsbury, N. R.; James, W. H., III; Plusquellic, D. F.; Zwier, T. S. J. Chem. Phys. 2011, 134, 164312] has assigned two conformers distinguished by the orientation of the in-plane OH groups and has identified two excitonic origins in each conformer. In the present study, the rotationally resolved bands of all four states have been well-fit to asymmetric rotor Hamiltonians. For the lower exciton (S(1)) levels, the transition dipole moment (TDM) orientations are perpendicular to the C(2) symmetry axes and consist of 41(2):59(2) and 34(2):66(2)% a:c hybrid-type character. The S(1) levels are therefore delocalized states of B symmetry and represent the antisymmetric combinations of the zero-order locally excited states of the p-cresol-like chromophores. The TDM polarizations of bands located at ≈132 cm(-1) above the S(1) origins are exclusively b-type and identify them as the upper exciton S(2) origin levels of A symmetry. The TDM orientations and the relative band strengths from the vibronic study have been analyzed within a dipole-dipole coupling model in terms of the localized TDM orientations, μ(loc), on the two chromophores. The out-of-the-ring plane angles of μ(loc) are both near 20° and are similar to results for diphenylmethane [Stearns, J. A.; Pillsbury, N. R.; Douglass, K. O.; Mu?ller, C. W.; Zwier, T. S.; Plusquellic, D. F. J. Chem. Phys. 2008, 129, 224305]. The in-plane angles are, however, rotated by 14 and 18° relative to DPM and, in part, explain the smaller than expected exciton splittings of these two conformers.     

Small-angle neutron scattering studies of charged carboxyl-terminated dendrimers in solutions

Small-angle neutron scattering was used to characterize the solution behavior of charged carboxylic acid terminated "cascade" dendrimers (Z-Cascade/methane [4]/3-oxo-6-oxa-2-azaheptylidyne/3-oxo-2-azaheptylidyne/propanoic acids) of third (G3) and fifth (G5) generations as a function of dendrimer concentration, pH, and ionic strength. An increase in dendrimer concentration leads to a single broad peak in the scattering profile arising from interdendrimer interaction. The dissociation of terminal carboxylate groups also gives rise to an interdendrimer interaction peak, which could be suppressed by the addition of excess salt. The results of contrast matching measurements indicate the accumulation of an excess concentration of tetramethylammonium counterions around the surface of these highly charged particles, and the thickness of these counterions lies somewhere between 4 and 6 A. This conclusion is consistent with our previous potentiometric titration (Zhang, H.; et al. J. Phys. Chem. B 1997, 101, 3494) and capillary electrophoresis (Huang, Q. R.; et al. J. Phys. Chem. B 2000, 104, 898) data.     

Interprotein electron transfer in a confined space: uncoupling protein dynamics from electron transfer by sol-gel encapsulation

In this paper, we describe the first observations of photoinitiated interprotein electron transfer (ET) within sol-gels. We have encapsulated three protein-protein complexes, specifically selected because they represent a full range of affinities, are sensitive to different types of dynamic processes, and thus are expected to respond differently to sol-gel encapsulation. The three systems are (i) the [Zn, Fe(3+)L] mixed-metal hemoglobin hybrids, where the alpha(1)-Zn and beta(2)-Fe subunits correspond to a "predocked" protein-protein complex with a crystallographically defined interface (Natan, M. J.; Baxter, W. W.; Kuila, D.; Gingrich, D. J.; Martin, G. S.; Hoffman, B. M. Adv. Chem. Ser. 1991, 228 (Electron-Transfer Inorg., Org., Biol. Syst.), 201-213), (ii) the Zn-cytochrome c peroxidase complex with cytochrome c, [ZnCcP, Fe(3+)Cc], having an intermediate affinity between its partners (Nocek, J. M.; Zhou, J. S.; De Forest, S.; Priyadarshy, S.; Beratan, D. N.; Onuchic, J. N.; Hoffman, B. M. Chem. Rev. 1996, 96, 2459-2489), and (iii) the [Zn-deuteromyoglobin, ferricytochrome b(5)] complex, [ZnDMb, Fe(3+)b(5)], which is loosely bound and highly dynamic (Liang, Z.-X.; Nocek, J.; Huang, K.; Hayes, R. T.; Kurnikov, I. V.; Beratan, D. N.; Hoffman, B. M. J. Am. Chem. Soc. 2002, 124, 6849-6859. Intersubunit ET within the hybrid does not involve second-order processes or subunit rearrangements, and thus is influenced only by perturbations of high-frequency motions coupled to ET. For the latter two complexes, sol-gel encapsulation eliminates second-order processes: protein partners encapsulated as a complex must stay together throughout a photoinitiated ET cycle, while proteins encapsulated alone cannot acquire a partner. It further modulates intracomplex motions of the two partners.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Effect of vertical temperature variation on the oscillatory wetting instability in a fluid Ga-Pb alloy

Oscillatory wetting instabilities driven by capillary-gravitation forces have been explored very recently in the binary fluid Ga-Pb alloy [A. Turchanin, R. Tsekov and W. Freyland, J. Chem. Phys., 2004, 120, 11 171]. This system is characterized by a complete wetting transition at liquid-liquid coexistence. Due to its metallic nature the bulk and interfacial instabilities are strongly coupled via variation of the respective emissivities. In our previous work we have investigated these phenomena at different cooling cycles and at constant temperature inside the miscibility gap. In this study we present for the first time the observations of the oscillatory wetting instabilities also in heating cycles. The interfacial properties of a Ga0.95Pb0.05 alloy at conditions inside the miscibility gap have been investigated by following the second harmonic generation (SHG) intensity changes. Corresponding model calculations of the Pb-rich wetting film instabilities have been performed taking into account the effect of a temperature variation vertical to the bulk sample. The influence of this temperature variation on the occurrence of the oscillations is discussed.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].     

Reaction of OH + NO2: high pressure experiments and falloff analysis

High pressure experiments on the OH + NO2 reaction are presented for 3 different temperatures. At 300 K, experiments in He (p = 2-500 bar) as well as in Ar (p = 2-4 bar) were performed. The rate constants obtained in Ar agree well with values which have been reported earlier by our group (Forster, R.; Frost, M.; Fulle, D.; Hamann, H. F.; Hippler, H.; Schlepegrell, A.; Troe, J. J. Chem. Phys. 1995, 103, 2949. Fulle, D.; Hamann, H. F.; Hippler, H.; Troe, J. J. Chem. Phys. 1998, 108, 5391). In contrast, the rate coefficients determined in He were found to be 15-25% lower than the values given in our earlier publications. Additionally, results for He as bath gas at elevated temperatures (T = 400 K, p = 3-150 bar; T = 600 K, p = 3-150 bar) are reported. The results obtained at elevated pressures are found to be in good agreement with existing literature data. The observed falloff behavior is analyzed in terms of the Troe formalism taking into account two reaction channels: one yielding HNO3 and one yielding HOONO. It is found that the extracted parameters are in agreement with rate constants for vibrational relaxation and isotopic scrambling as well as with experimentally determined branching ratios. Based on our analysis we determine falloff parameters to calculate the rate constant for atmospheric conditions.     

Electron magnetic resonance and density functional theory study of room temperature X-irradiated β-D-fructose single crystals

Stable free radical formation in fructose single crystals X-irradiated at room temperature was investigated using Q-band electron paramagnetic resonance (EPR), electron nuclear double resonance (ENDOR), and ENDOR induced EPR (EIE) techniques. ENDOR angular variations in the three main crystallographic planes allowed an unambiguous determination of 12 proton HFC tensors. From the EIE studies, these hyperfine interactions were assigned to six different radical species, labeled F1-F6. Two of the radicals (F1 and F2) were studied previously by Vanhaelewyn et al. [Vanhaelewyn, G. C. A. M.; Pauwels, E.; Callens, F. J.; Waroquier, M.; Sagstuen, E.; Matthys, P. J. Phys. Chem. A 2006, 110, 2147.] and Tarpan et al. [Tarpan, M. A.; Vrielinck, H.; De Cooman, H.; Callens, F. J. J. Phys. Chem. A 2009, 113, 7994.]. The other four radicals are reported here for the first time and periodic density functional theory (DFT) calculations were used to aid their structural identification. For the radical F3 a C3 carbon centered radical with a carbonyl group at the C4 position is proposed. The close similarity in HFC tensors suggests that F4 and F5 originate from the same type of radical stabilized in two slightly different conformations. For these radicals a C2 carbon centered radical model with a carbonyl group situated at the C3 position is proposed. A rather exotic C2 centered radical model is proposed for F6.     

Excluded volume effect for large and small solutes in water

The cavitation effect, i.e., the process of the creation of a void of excluded volume in bulk solvent (a cavity), is considered. The cavitation free energy is treated in terms of the information theory (IT) approach [Hummer, G.; Garde, S.; Garcia, A. E.; Paulaitis, M. E.; Pratt, L. R. J. Phys. Chem. B 1998, 102, 10469]. The binomial cell model suggested earlier is applied as the IT default distribution p(m) for the number m of solute (water) particles occupying a cavity of given size and shape. In the present work, this model is extended to cover the entire range of cavity size between small ordinary molecular solutes and bulky biomolecular structures. The resulting distribution consists of two binomial peaks responsible for producing the free energy contributions, which are proportional respectively to the volume and to the surface area of a cavity. The surface peak dominates in the large cavity limit, when the two peaks are well separated. The volume effects become decisive in the opposite limit of small cavities, when the two peaks reduce to a single-peak distribution as considered in our earlier work. With a proper interpolation procedure connecting these two regimes, the MC simulation results for model spherical solutes with radii increasing up to R = 10 A [Huang, D. H.; Geissler, P. L.; Chandler, D. J. Phys. Chem. B 2001, 105, 6704] are well reproduced. The large cavity limit conforms to macroscopic properties of bulk water solvent, such as surface tension, isothermal compressibility and Tolman length. The computations are extended to include nonspherical solutes (hydrocarbons C1-C6).     

Breakdown of the pseudopotential approximation for magnetizabilities and electric multipole moments. II. The importance of gauge invariance for large-core semi-local pseudopotentials

In a previous paper we pointed out that core contributions to the static magnetizability are non-negligible, and can therefore lead to erroneous results within the pseudopotential approximation [P. Schwerdtfeger, B. Assadollahzadeh, U. Rohrmann, R. Scha?fer, and J. R. Cheeseman, J. Chem. Phys. 134, 204102 (2011)]. In a recent paper van Wu?llen showed that additional terms arising from the gauge-invariant condition for the semi-local part of a pseudopotential operator can lead to non-negligible contributions to the magnetizability tensor, which are sensitive to the gauge origin and basis set chosen [C. van Wu?llen, J. Chem. Phys. 136, 114110 (2012)]. These terms were neglected in previous calculations as they were assumed to be small. In this paper we analyze the importance of the gauge-dependent semi-local pseudopotential correction term in detail for AuF and clusters of Sn showing that it leads indeed to very large corrections to the paramagnetic term for large-core pseudopotentials. Without this correction the results become very sensitive to the basis set applied. This now resolves some of the unusual large paramagnetic contributions reported before for both AuF and Sn(2).     

Marked difference in the electronic structure of cyanide-ligated ferric protoglobins and myoglobin due to heme ruffling

Electron paramagnetic resonance experiments reveal a significant difference between the principal g values (and hence ligand-field parameters) of the ferric cyanide-ligated form of different variants of the protoglobin of Methanosarcina acetivorans (MaPgb) and of horse heart myoglobin (hhMb). The largest principal g value of the ferric cyanide-ligated MaPgb variants is found to be significantly lower than for any of the other globins reported so far. This is at least partially caused by the strong heme distortions as proven by the determination of the hyperfine interaction of the heme nitrogens and mesoprotons. Furthermore, the experiments confirm recent theoretical predictions [Forti, F.; Boechi, L., Bikiel, D., Martí, M.A.; Nardini, M.; Bolognesi, M.; Viappiani, C.; Estrin, D.; Luque, F. J. J. Phys. Chem. B2011, 115, 13771-13780] that Phe(G8)145 plays a crucial role in the ligand modulation in MaPgb. Finally, the influence of the N-terminal 20 amino-acid chain on the heme pocket in these protoglobins is also proven.     

Complex Langmuir-Blodgett films of SiO2 and ZnO nanoparticles with advantageous optical and photocatalytical properties

Multifunctional Langmuir-Blodgett (LB) films were fabricated on the surface of glass substrates using sol-gel derived ZnO and SiO2 particles. ZnO particles of 6 and 110 nm diameter were synthesized according to the methods of Meulenkamp and Seelig et al. (Meulenkamp, E. A. J. Phys. Chem. B 1998, 102, 5566; Seelig, E. W.; Tang, B.; Yamilov, A.; Cao, H.; Chang, R. P. H. Mater. Chem. Phys. 2003, 80, 257). Silica particles of 37 and 96 nm were prepared by the Stober method (Stober, W.; Fink, A.; Bohn, E. J. Colloid Interface Sci. 1968, 26, 62). Alternate deposition of monoparticulate Langmuir films of SiO2 and ZnO nanoparticles provided complex (six- and nine-layered) LB films with both antireflective and photocatalytic properties. The LB films were investigated with scanning electron microscopy (morphology and structure) and UV-vis spectroscopy (optical properties and stability). The photocatalytic activity was measured by immersing the UV-irradiated films into an aqueous solution of Methyl Orange and following the photodegradation of the dye by optical spectroscopy. Adding ZnO particles to the silica films slightly lowered the antireflection property but ensured strong photocatalytic activity. Both the photocatalytic activity and antireflection properties were proved to be sensitive to the sequence of the silica and ZnO layers, with optimum properties in the case of nine-layered films with a repeated (SiO2-ZnO-ZnO) structure.     

Solid-state 23Na NMR study of sodium lariat ether receptors exhibiting cation-pi interactions

Noncovalent cation-pi interactions are important in a variety of supramolecular and biochemical systems. We present a 23Na solid-state nuclear magnetic resonance (SSNMR) study of two sodium lariat ether complexes, 1 and 2, in which a sodium cation interacts with an indolyl group that models the side chain of tryptophan. Sodium-23 SSNMR spectra of magic-angle spinning (MAS) and stationary powdered samples have been acquired at three magnetic field strengths (9.4, 11.75, 21.1 T) and analyzed to provide key information on the sodium electric field gradient and chemical shift (CS) tensors which are representative of the cation-pi binding environment. Triple-quantum MAS NMR spectra acquired at 21.1 T clearly reveal two crystallographically distinct sites in both 1 and 2. The quadrupolar coupling constants, CQ(23Na), range from 2.92 +/- 0.05 MHz for site A of 1 to 3.33 +/- 0.05 MHz for site B of 2; these values are somewhat larger than those reported previously by Wong et al. (Wong, A.; Whitehead, R. D.; Gan, Z.; Wu, G. J. Phys. Chem. A 2004, 108, 10551) for NaBPh4, but very similar to the values obtained for sodium metallocenes by Willans and Schurko (Willans, M. J.; Schurko, R. W. J. Phys. Chem. B 2003, 107, 5144). We conclude from the 21.1 T data that the spans of the sodium CS tensors are less than 20 ppm for 1 and 2 and that the largest components of the EFG and CS tensors are non-coincident. Quantum chemical calculations of the NMR parameters substantiate the experimental findings and provide additional insight into the dependence of CQ(23Na) on the proximity of the indole ring to Na+. Taken together, this work has provided novel information on the NMR interaction tensors characteristic of a sodium cation interacting with a biologically important arene.     

A post-Hartree-Fock model of intermolecular interactions: inclusion of higher-order corrections

We have previously demonstrated that the dipole moment of the exchange hole can be used to derive intermolecular C(6) dispersion coefficients [J. Chem. Phys. 122, 154104 (2005)]. This was subsequently the basis for a novel post-Hartree-Fock model of intermolecular interactions [J. Chem. Phys. 123, 024101 (2005)]. In the present work, the model is extended to include higher-order dispersion coefficients C(8) and C(10). The extended model performs very well for prediction of intermonomer separations and binding energies of 45 van der Waals complexes. In particular, it performs twice as well as basis-set extrapolated MP2 theory for dispersion-bound complexes, with minimal computational cost.     

Existence of a size-dependent diffusivity maximum for uncharged solutes in water and its implications

Recent studies suggest that there exists a size-dependent diffusivity maximum in binary mixtures interacting via Lennard-Jones potential when the size of one of the two components is varied (Ghorai, P. K.; Yashonath, S. J. Phys. Chem., 2005, 109, 5824). We discuss in the present paper the importance of the existence of a size-dependent maximum for an uncharged solute in liquid or amorphous solid water and its relation to the ionic conductivity maximum in water. We report molecular dynamics investigations into the size dependence of the self-diffusivity, D, of the uncharged solutes in water at low temperatures (30 K) with immobile as well as mobile water. We find that a maximum in self-diffusivity exists as a function of the size of solute diffusing within water at low temperatures but not at high temperatures. This is due to the relatively weak interactions between the solute and the water compared to the kinetic energy at room temperature. Previously, we have shown that a similar maximum exists for guests sorbed in zeolites and is known as the levitation effect (LE). Thus, it appears that the existence of a size-dependent maximum is universal and extends from zeolites to simple liquids to solvents of polyatomic species. We examine the implications of this for the size-dependent maximum in ionic conductivity in polar solvents known for over a hundred years. These results support the view that the size-dependent maximum seen for ions in water has its origin in the LE (see Ghorai, P. Kr.; Yashonath, S.; Lynden-Bell, R. M. J. Phys. Chem. 2005, 109, 8120).     

Grand canonical Monte Carlo simulations of the 129Xe NMR line shapes of xenon adsorbed in (+/-)-[Co(en)3]Cl3

The 129Xe NMR line shapes of xenon adsorbed in the nanochannels of the (+/-)-[Co(en)3]Cl3 ionic crystal have been calculated by grand canonical Monte Carlo (GCMC) simulations. The results of our GCMC simulations illustrate their utility in predicting 129Xe NMR chemical shifts in systems containing a transition metal. In particular, the nanochannels of (+/-)-[Co(en)3]Cl3 provide a simple, yet interesting, model system that serves as a building block toward understanding xenon chemical shifts in more complex porous materials containing transition metals. Using only the Xe-C and Xe-H potentials and shielding response functions derived from the Xe@CH4 van der Waals complex to model the interior of the channel, the GCMC simulations correctly predict the 129Xe NMR line shapes observed experimentally (Ueda, T.; Eguchi, T.; Nakamura, N.; Wasylishen, R. E. J. Phys. Chem. B 2003, 107, 180-185). At low xenon loading, the simulated 129Xe NMR line shape is axially symmetric with chemical-shift tensor components delta(parallel) = 379 ppm and delta(perpendicular) = 274 ppm. Although the simulated isotropic chemical shift, delta(iso) = 309 ppm, is overestimated, the anisotropy of the chemical-shift tensor is correctly predicted. The simulations provide an explanation for the observed trend in the 129Xe NMR line shapes as a function of the overhead xenon pressure: delta(perpendicular) increased from 274 to 292 ppm, while delta(parallel) changed by only 3 ppm over the entire xenon loading range. The overestimation of the isotropic chemical shifts is explained based upon the results of quantum mechanical 129Xe shielding calculations of xenon interacting with an isolated (+/-)-[Co(en)3]Cl3 molecule. The xenon chemical shift is shown to be reduced by about 12% going from the Xe@[Co(en)3]Cl3 van der Waals complex to the Xe@C2H6 fragment.