CASPT2 Study of the Potential Energy Surface of the HSO(2) System

The importance of the HSO(2) system in atmospheric and combustion chemistry has motivated several works dedicated to the study of associated structures and chemical reactions. Nevertheless controversy still exists in connection with the reaction SH + O(2) -> H + SO(2) and also related to the role of the HSOO isomers in the potential energy surface (PES). Here we report high-level ab initio calculation for the electronic ground state of the HSO(2) system. Energetic, geometric, and frequency properties for the major stationary states of the PES are reported at the same level of calculations:,CASPT2/aug-cc-pV(T+d)Z. This study introduces three new stationary points (two saddle points and one minimum). These structures allow the connection of the skewed HSOOs and the HSO(2) minima defining new reaction paths for SH + O(2) -> H + SO(2) and SH + O(2) -> OH + SO. In addition, the location of the HSOO isomers in the reaction pathways have been clarified.

Gamma Knife(A (R)) 3-D Dose Distribution Mapping by Magnetic Resonance Imaging

In medical processes where ionizing radiation is used, dose planning and dose delivery are the key elements to patient safety and treatment success, particularly, when the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Therefore, the radiation dose should be well defined and precisely delivered to the target while minimizing radiation exposure to surrounding normal tissues [1]. Several methods have been proposed to obtain three-dimensional (3-D) dose distribution [2, 3]. In this paper, we propose an alternative method, which can be easily implemented in any stereotactic radiosurgery center with a magnetic resonance imaging (MRI) facility. A phantom with or without scattering centers filled with Fricke gel solution is irradiated with Gamma Knife(A (R)) system at a chosen spot. The phantom can be a replica of a human organ such as head, breast or any other organ. It can even be constructed from a real 3-D MR image of an organ of a patient using a computer-aided construction and irradiated at a specific region corresponding to the tumor position determined by MRI. The spin-lattice relaxation time T (1) of different parts of the irradiated phantom is determined by localized spectroscopy. The T (1)-weighted phantom images are used to correlate the image pixels intensity to the absorbed dose and consequently a 3-D dose distribution with a high resolution is obtained.

Continuum, Discrete, and Explicit Solvation Models for Describing the Low-Lying Absorption Spectrum of the Pterin Acid in Aqueous Environment

The absorption spectrum of the acid form of pterin in water was investigated theoretically. Different procedures using continuum, discrete, and explicit models were used to include the solvation effect on the absorption spectrum, characterized by two bands. The discrete and explicit models used Monte Carlo simulation to generate the liquid structure and time-dependent density functional theory (B3LYP/6-31G+(d)) to obtain the excitation energies. The discrete model failed to give the correct qualitative effect on the second absorption band. The continuum model, in turn, has given a correct qualitative picture and a semiquantitative description. The explicit use of 29 solvent molecules, forming a hydration shell of 6 angstrom, embedded in the electrostatic field of the remaining solvent molecules, gives absorption transitions at 3.67 and 4.59 eV in excellent agreement with the S(0)-S(1) and S(0)-S(2) absorption bands at of 3.66 and 4.59 eV, respectively, that characterize the experimental spectrum of pterin in water environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2371-2377, 2010

Excited state electronic polarization and reappraisal of the n <- pi* emission of acetone in water

Electronic polarization of the acetone molecule in the excited n -> pi* state is considered and its influence on the solvent shift in the emission spectrum is analyzed. Using an iterative procedure the electronic polarizations of both the ground and the excited states are included and compared with previous results obtained with Car-Parrinello dynamics. Analysis of the emission transition obtained using CIS(D)/aug-cc-pVDZ on statistically uncorrelated solute-solvent structures, composed of acetone and twelve explicit water molecules embedded in the electrostatic field of remaining 263 water molecules, corroborates that the solvent effect is mild, calculated here between 80 and 380 cm (1). (c) 2010 Published by Elsevier B.V.

An experimental and theoretical investigation into positron and electron scattering from formaldehyde

We report on measurements of total cross sections (TCSs) for positron scattering from the fundamental organic molecule formaldehyde (CH(2)O). The energy range of these measurements was 0.26-50.3 eV, whereas the energy resolution was similar to 260 meV. To assist us in interpreting these data, Schwinger multichannel level calculations for positron elastic scattering from CH(2)O were also undertaken (0.5-50 eV). These calculations, incorporating an accurate model for the target polarization, are found to be in good qualitative agreement with our measured data. In addition, in order to compare the behaviour of positron and electron scattering from this species, independent atom model-screened additivity rule theoretical electron TCSs, now for energies in the range 1-10 000 eV, are also reported.

Characterization and spectroscopic analysis of phenol-ethanol hydrogen bonded clusters

Different hydrogen bonded clusters involving phenol and ethanol are studied theoretically using MP2/aug-cc-pVDZ. Nine different 1: 1 clusters are obtained and analyzed according to their stability and spectroscopic properties. Different isomeric forms of ethanol are considered. Attention is also devoted to the spectral shift of the characteristic pi -> pi* transition of phenol. Using TDHF, CIS, CIS(D) and TDB3LYP in aug-cc-pVDZ basis set, all results agree that a red shift is obtained when phenol is the hydrogen donor and a blue shift is obtained in the opposite case. These results are used to rationalize the red shift observed for phenol in liquid ethanol. (C) 2010 Elsevier B.V. All rights reserved.

Determination of Very Slow Diffusion of Paramagnetic Ions by NMR Spectroscopy

In this paper, we propose a new method of measuring the very slow paramagnetic ion diffusion coefficient using a commercial high-resolution spectrometer. If there are distinct paramagnetic ions influencing the hydrogen nuclear magnetic relaxation time differently, their diffusion coefficients can be measured separately. A cylindrical phantom filled with Fricke xylenol gel solution and irradiated with gamma rays was used to validate the method. The Fricke xylenol gel solution was prepared with 270 Bloom porcine gelatin, the phantom was irradiated with gamma rays originated from a (60)Co source and a high-resolution 200 MHz nuclear magnetic resonance (NMR) spectrometer was used to obtain the phantom (1)H profile in the presence of a linear magnetic field gradient. By observing the temporal evolution of the phantom NMR profile, an apparent ferric ion diffusion coefficient of 0.50 mu m(2)/ms due to ferric ions diffusion was obtained. In any medical process where the ionizing radiation is used, the dose planning and the dose delivery are the key elements for the patient safety and success of treatment. These points become even more important in modern conformal radio therapy techniques, such as stereotactic radiosurgery, where the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Several methods have been proposed to obtain the three-dimensional (3-D) dose distribution. Recently, we proposed an alternative method for the 3-D radiation dose mapping, where the ionizing radiation modifies the local relative concentration of Fe(2+)/Fe(3+) in a phantom containing Fricke gel and this variation is associated to the MR image intensity. The smearing of the intensity gradient is proportional to the diffusion coefficient of the Fe(3+) and Fe(2+) in the phantom. There are several methods for measurement of the ionic diffusion using NMR, however, they are applicable when the diffusion is not very slow.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

Solvent effects on the electronic absorption spectrum of camphor using continuum, discrete or explicit approaches

We address the effect of solvation on the lowest electronic excitation energy of camphor. The solvents considered represent a large variation in-solvent polarity. We consider three conceptually different ways of accounting for the solvent using either an implicit, a discrete or an explicit solvation model. The solvatochromic shifts in polar solvents are found to be in good agreement with the experimental data for all three solvent models. However, both the implicit and discrete solvation models are less successful in predicting solvatochromic shifts for solvents of low polarity. The results presented suggest the importance of using explicit solvent molecules in the case of nonpolar solvents. (C) 2009 Elsevier B.V. All rights reserved.

Electronic spectroscopy of biomolecules in solution: fluorescein dianion in water

A combined and sequential Monte Carlo-quantum mechanics methodology is used to describe the electronic absorption spectrum of the fluorescein dianion in water. Different sets of 100 statistically relevant configurations composed of the solute and several solvent molecules are sampled from the Monte Carlo simulation for a posteriori quantum mechanical calculations of the spectra. In the largest case the configurations are composed of fluorescein and 90 explicit water molecules embedded in the electrostatic field of all remaining water molecules within a distance of 11.3 angstrom. These configurations include 305 atoms and 842 valence electrons, justifying the use of a semi-empirical approach. The electronic spectrum is then calculated using the INDO/CIS method. The solvatochromic shift of fluorescein in water, compared with in isolation, is calculated using the discrete and explicit solvent models. The use of electrostatically embedded explicit water molecules, in INDO/CIS calculations, gives a good description of the spectral shift of the fluorescein dianion in aqueous environment. The results are verified to converge both statistically and with respect to the number of explicit solvent molecules used.

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.

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.