NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Solvent Effects in Chemical Processes. Water-Assisted Proton Transfer Reaction of Pterin in Aqueous Environment

Pterins are members of a family of heterocyclic compounds present in a wide variety of biological systems and may exist in two forms, corresponding to an acid and a basic tautomer. In this work, the proton transfer reaction between these tautomeric forms was investigated in the gas phase and in aqueous solution. In gas phase, the intramolecular mechanism was carried out for die isolated pterin by quantum mechanical second-order Moller-Plesset Perturbation theory (MP2/aug-cc-pVDZ) calculations and it indicates that the acid form is more stable than the basic form by -1.4 kcal/mol with a barrier of 34.2 kcal/mol with respect to the basic form. In aqueous solution, the role of the water molecules in the proton transfer reaction was analyzed in two separated parts, the direct participation of one water molecule in the reaction path, called water-assisted mechanism, and the complementary participation of the aqueous solvation. The water-assisted mechanism was carried out for one pterin-water cluster by quantum mechanical calculations and it indicates that the acid form is still more stable by -3.3 kcal/mol with a drastic reduction of 70% of the barrier, The bulk solution effect on the intramolecular and water-assisted mechanisms was included by free energy perturbation implemented on Monte Carlo simulations. The bulk water effect is found to be substantial and decisive when the reaction path involves the water-assisted mechanism. In this case, the free energy barrier is only 6.7 kcal/mol and the calculated relative Gibbs free energy for the two tautomers is -11.2 kcal/mol. This value is used to calculate the pK(a) value of 8.2 +/- 0.6 that is in excellent agreement with the experimental result of 7.9.

Dipole polarizability and Rayleigh light scattering by the hydrated electron

Assuming the existence of a confined state of the electron in bulk water the polarizability of the hydrated electron is analyzed. Statistically uncorrelated supermolecular structures composed of seven water molecules (first solvation shell) with an extra electron were extracted from classical Monte Carlo simulation and used in quantum mechanical second-order Moller-Plesset calculations. It is found that the bound excess electron contributes with 274 a.u. to the total dipole polarizability of 345 a.u. for (H(2)O)(7)(-). From the calculated polarizabilities the Rayleigh elastic light scattering properties are inferred and found to considerably enhance activity and light depolarization. (C) 2009 Elsevier B.V. All rights reserved.

The degree of polydispersity in micellar dispersions: A sum rule approach

Statistical properties of a two-dimensional ideal dispersion of polydisperse micelles are derived by analyzing the convergence properties of a sum rule set by mass conservation. Internal micellar degrees of freedom are accounted for by a microscopic model describing small displacements of the constituting amphiphiles with respect to their equilibrium positions. The transfer matrix (TM) method is employed to compute internal micelle partition function. We show that the conditions under which the sum rule is saturated by the largest eigenvalue of the TM determine the value of amphiphile concentration above which the dispersion becomes highly polydisperse and micelle sizes approach a Schultz distribution. (C) 2011 Elsevier B.V. All rights reserved.

On the catalytic hydrogenation of polycyclic aromatic hydrocarbons into less toxic compounds by a facile recoverable catalyst

Here we present the catalytic hydrogenation of polycyclic aromatic hydrocarbons (PAHs) to less toxic mixtures of saturated and partial unsaturated polycyclic hydrocarbons under mild reaction conditions using a magnetically recoverable rhodium catalyst and molecular hydrogen as the exclusive H source. The catalyst is easily recovered after each reaction by placing a permanent magnet on the reactor wall and it can be reused in successive runs without any significant loss of catalytic activity. As an example, anthracene was totally converted into the saturated polycyclic hydrocarbon form (ca. 60%) and the partially hydrogenated form, 1,2,3,4,5,6,7,8-octahydroanthracene (ca. 40%). The catalyst operates in a broad range of temperature and H(2) pressure in both organic and aqueous/organic solutions of anthracene and it also exhibits significant activity at low substrate concentrations (20 ppm). This can be an efficient recycling process for hydrogenation of PAHs present in contaminated fluid waste streams. (C) 2009 Elsevier B.V. All rights reserved.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Polarization and solvatochromic shift of ortho-betaine in water

Betaine dyes are known to show very large transition energy shifts in different solvents. The ortho-betaine molecule - a simple two-ring prototype of the E-T(30) Reichardt dye - has been investigated theoretically from a combined statistical and quantum mechanics approach. Using sequential Monte Carlo (MC) simulations and MP2/cc-pVDZ calculations the in-water dipole moment of ortho-betaine is obtained as 12.30 +/- 0.05 D. This result shows a considerable increase of 75% compared to the in-vacuum dipole moment. For comparison, the use of a polarizable continuum model using the same MP2/cc-pVDZ leads to an in-water dipole moment of 11.6 D, in good agreement. This large polarization is incorporated in the classical potential for another MC simulation to generate solute-solvent configurations and to obtain the contribution of the polarization effect in the solvatochromic shift. Using statistically uncorrelated configurations and supermolecular INDO/CIS calculations, including the solute and, explicitly, 230 solvent water molecules, the statistically converged calculated shift is obtained here as 6360 cm(-1), in good agreement with the experimental result of 7550 cm(-1). (c) 2007 Elsevier B.V. All rights reserved.

On the interaction of the anthraquinone barbaloin with negatively charged DMPG bilayers

Barbaloin is a bioactive glycosilated 1,8-dihydroxyanthraquinone present in several exudates from plants, Such as Aloe vera, which are used for cosmetic or food purposes. It has been shown that barbaloin interacts with DMPG (dimyristoylphosphatidylglycerol) model membranes, altering the bilayer structure (Alves, D. S.; Perez-Fons, L.; Estepa, A.; Micol, V. Biochem. Pharm. 2004, 68, 549). Considering that ESR (electron spin resonance) of spin labels is one of the best techniques to monitor structural properties at the molecular level, the alterations caused by the anthraquinone barbaloin on phospholipid bilayers will be discussed here via the ESR signal of phospholipid spin probes intercalated into the membranes. In DMPG at high ionic strength (10 mM Hepes pH 7.4 + 100 mM NaCl), a system that presents a gel-fluid transition around 23 degrees C, 20 mol % barbaloin turns the gel phase more rigid, does not alter much the fluid phase packing, but makes the lipid thermal transition less sharp. However, in a low-salt DMPG dispersion (10 mM Hepes pH 7.4 + 2 mM NaCl), which presents a rather complex gel-fluid thermal transition (Lamy-Freund, M. T.; Riske, K. A. Chem. Phys. Lipids 2003, 122, 19), barbaloin strongly affects bilayer structural properties, both in the gel and fluid phases, extending the transition region to much higher temperature values. The position of barbaloin in DMPG bilayers will be discussed on the basis of ESR results, in parallel with data from sample viscosity, DSC (differential scanning calorimetry), and SAXS (small-angle X-ray scattering).

Self-assembling of phenothiazine compounds investigated by small-angle X-ray scattering and electron paramagnetic resonance spectroscopy

Small-angle X-ray scattering (SAXS) and electron paramagnetic resonance (EPR) have been carried out to investigate the structure of the self-aggregates of two phenothiazine drugs, chlorpromazine (CPZ) and trifluoperazine (TFP), in aqueous solution. In the SAXS studies, drug solutions of 20 and 60 mM, at pH 4.0 and 7.0, were investigated and the best data fittings were achieved assuming several different particle form factors with a homogeneous electron density distribution in respect to the water environment. Because of the limitation of scattering intensity in the q range above 0.15 angstrom(-1), precise determination of the aggregate shape was not possible and all of the tested models for ellipsoids, cylinders, or parallelepipeds fitted the experimental data equally well. The SAXS data allows inferring, however, that CPZ molecules might self-assemble in a basis set of an orthorhombic cell, remaining as nanocrystallites in solution. Such nanocrystals are composed of a small number of unit cells (up to 10, in c-direction), with CPZ aggregation numbers of 60-80. EPR spectra of 5- and 16-doxyl stearic acids bound to the aggregates were analyzed through simulation, and the dynamic and magnetic parameters were obtained. The phenothiazine concentration in EPR experiments was in the range of 5-60 mM. Critical aggregation concentration of TFP is lower than that for CPZ, consistent with a higher hydrophobicity of TFP. At acidic pH 4.0 a significant residual motion of the nitroxide relative to the aggregate is observed, and the EPR spectra and corresponding parameters are similar to those reported for aqueous surfactant micelles. However, at pH 6.5 a significant motional restriction is observed, and the nitroxide rotational correlation times correlate very well with those estimated for the whole aggregated particle from SAXS data. This implies that the aggregate is densely packed at this pH and that the nitroxide is tightly bound to it producing a strongly immobilized EPR spectrum. Besides that, at pH 6.5 the differences in motional restriction observed between 5- and 16-DSA are small, which is different from that observed for aqueous surfactant micelles.

Calculations of vibrational frequencies, Raman activities and degrees of depolarization for complexes involving water, methanol and ethanol

Raman activities and degrees of depolarization are reported for 14 complexes involving methanol, ethanol and water using the MP2/aug-cc-pVDZ model. For ethanol both trans and gauche isomers are considered. The red-shifts of the OH stretching and the blue shifts of the bending tau(CO-OH) mode were analyzed for the proton-donor molecules upon hydrogen bond. The shift of the nu(CO) stretching mode of the alcohol molecules are also analyzed and found to be specific giving characterization of the amphoteric relation, being positive for the proton-acceptor and negative for the proton-donor molecule. (c) 2008 Elsevier B.V. All rights reserved.

Conformational behavior of different possible ways of oligoglycine formation in a solvent-free environment

The possible ways for glycine oligopeptide formation in gas phase, both in the extended P-strand like conformation and folded 2(7)-ribbon like conformations are analyzed using quantum chemical calculations. We focus on the sequential formation of peptide bond through upgradation of the immediate lower order molecule and observe the consequences in other related processes like oligoglycine formation through simultaneous peptide linkage of n glycine monomers and interchange of molecular conformation through peptide linkage. A comparison is made between the structures and binding energies obtained for both conformers. All binding energies are increased by the zero-point energy contribution. The role of electron correlation effects is briefly analyzed. The folded 2(7)-ribbon-like conformations in vacuo are found to be more stable in comparison to the extended structure. (c) 2007 Elsevier B.V. All rights reserved.

Electronic and magnetic properties of Ti and Fe on graphene

The structural, electronic and magnetic properties of Fe and Ti atomic wires and the complete covering when adsorbed on graphene are presented through ab initio calculations based on density functional theory. The most stable configurations are investigated for Fe and Ti in different concentrations adsorbed on the graphene surface, and the corresponding binding energies are calculated. The results show a tendency of the Ti atoms to cover uniformly the graphene surface, whereas the Fe atoms form clusters. The adsorption of the transition metal on the graphene surface changes significantly the electronic density of states near the graphene Fermi region. In all arrangements studied, a charge transfer is observed from the adsorbed species to the graphene surface due to the high hybridizations between the systems.

Ab initio study of SO2 molecules interacting with pristine and transition metal covered fullerenes as a possible route for nanofilters

In this work, a systematic study of SO2 molecules interacting with pristine and transition metal (TM) covered C-60 is presented by means of first principles calculations. It is observed that the SO2 molecule interacts weakly with the pristine C-60 fullerene, although the resulting interaction is largely increased when the C-60 structure is covered with Fe, Mn, or Ti atoms and the SO2 Molecules are bounded through the TM atoms. The number of bounded SO2 molecules per TM atoms, in addition to the elevated binding energies per molecules, allows us to conclude that such composites can be used as a template for efficient devices to remove SO2 molecules or, alternatively, as SO2 gas sensor.

Recoverable rhodium nanoparticles: Synthesis, characterization and catalytic performance in hydrogenation reactions

We here report the first magnetically recoverable Rh(0) nanoparticle-supported catalyst with extraordinary recovery and recycling properties. Magnetic separation has been suggested as a very promising technique to improve recovery of metal-based catalysts in liquid-phase batch reactions. The separation method is significantly simple, as it does not require filtration, decantation, centrifugation, or any other separation technique thereby, overcoming traditional time- and solvent-consuming procedures. Our new magnetically separable catalytic system, comprised of Rh nanoparticles immobilized on silica-coated magnetite nanoparticles, is highly active and could be reused for up to 20 times for hydrogenation of cyclohexene (180,000 mol/mol(Rh)) and benzene (11,550 mol/mol(Rh) under mild conditions. (c) 2007 Elsevier B. V. All fights reserved.

Driving forces for the adsorption of cyclopentene on InP(001)

In this work the interaction of cyclopentene with a set of InP(001) surfaces is investigated by means of the density functional theory. We propose a simple approach for evaluating the surface strain and based on it we have found a linear relation between bond and strain energies and the adsorption energy. Our results also indicate that the higher the bond energy, the more disperse the charge distribution is around the adsorption site associated to the high occupied state, a key feature that characterizes the adsorption process. Different adsorption coverages are used to evaluate the proposed equation. Our results suggest that the proposed approach might be extended to other systems where the interaction of the semiconductor surface and the molecule is restricted to first neighbor sites. (C) 2011 Elsevier B.V. All rights reserved.