Argon predissociation spectroscopy of the OH-center dot H2O and Cl-center dot H2O complexes in the 1000-1900 cm(-1) region: Intramolecular bending transitions and the search for the shared-proton fundamental in the hydroxide monohydrate

We present argon predissociation vibrational spectra of the OH-.H2O and Cl-.H2O complexes in the 1000-1900 cm(-1) energy range, far below the OH stretching region reported in previous studies. This extension allows us to explore the fundamental transitions of the intramolecular bending vibrations associated with the water molecule, as well as that of the shared proton inferred from previous assignments of overtones in the higher energy region. Although the water bending fundamental in the Cl-.H2O spectrum is in very good agreement with expectations, the OH-.H2O spectrum is quite different than anticipated, being dominated by a strong feature at 1090 cm(-1). New full-diniensionality calculations of the OH-.H2O vibrational level structure using diffusion Monte Carlo and the VSCF/CI methods indicate this band arises from excitation of the shared proton.

Extracting quantitative structural parameters for disordered polymers from neutron scattering data

The organization of non-crystalline polymeric materials at a local level, namely on a spatial scale between a few and 100 a, is still unclear in many respects. The determination of the local structure in terms of the configuration and conformation of the polymer chain and of the packing characteristics of the chain in the bulk material represents a challenging problem. Data from wide-angle diffraction experiments are very difficult to interpret due to the very large amount of information that they carry, that is the large number of correlations present in the diffraction patterns.We describe new approaches that permit a detailed analysis of the complex neutron diffraction patterns characterizing polymer melts and glasses. The coupling of different computer modelling strategies with neutron scattering data over a wide Q range allows the extraction of detailed quantitative information on the structural arrangements of the materials of interest. Proceeding from modelling routes as diverse as force field calculations, single-chain modelling and reverse Monte Carlo, we show the successes and pitfalls of each approach in describing model systems, which illustrate the need to attack the data analysis problem simultaneously from several fronts.

Numerical analysis of an incremented statistical sampling procedure in MCNP

MCNP has stood so far as one of the main Monte Carlo radiation transport codes. Its use, as any other Monte Carlo based code, has increased as computers perform calculations faster and become more affordable along time. However, the use of Monte Carlo method to tally events in volumes which represent a small fraction of the whole system may turn to be unfeasible, if a straight analogue transport procedure (no use of variance reduction techniques) is employed and precise results are demanded. Calculations of reaction rates in activation foils placed in critical systems turn to be one of the mentioned cases. The present work takes advantage of the fixed source representation from MCNP to perform the above mentioned task in a more effective sampling way (characterizing neutron population in the vicinity of the tallying region and using it in a geometric reduced coupled simulation). An extended analysis of source dependent parameters is studied in order to understand their influence on simulation performance and on validity of results. Although discrepant results have been observed for small enveloping regions, the procedure presents itself as very efficient, giving adequate and precise results in shorter times than the standard analogue procedure. (C) 2007 Elsevier Ltd. All rights reserved.

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.

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.

A dedicated setup for the measurement of the electron transport parameters in gases at large electric fields

Electron transport parameters are important in several areas ranging from particle detectors to plasma-assisted processing reactors. Nevertheless, especially at high fields strengths and for complex gases, relatively few data are published. A dedicated setup has been developed to measure the electron drift velocity and the first Townsend coefficient in parallel plate geometry. An RPC-like cell has been adopted to reach high field strengths without the risk of destructive sparks. The validation data obtained with pure Nitrogen will be presented and compared to a selection of the available literature and to calculations performed with Magboltz 2 version 8.6. The new data collected in pure Isobutane will then be discussed. This is the first time the electron drift velocity in pure Isobutane is measured well into the saturation region. Good agreement is found with expectations from Magboltz. (C) 2009 Elsevier B.V. All rights reserved.

Tightened Lieb-Oxford Bound for Systems of Fixed Particle Number

The Lieb-Oxford bound is a constraint upon approximate exchange-correlation functionals. We explore a nonempirical tightening of that bound in both universal and electron number-dependent form. The test functional is PBE. Regarding both atomization energies (slightly worsened) and bond lengths (slightly improved), we find the PBE functional to be remarkably insensitive to the value of the Lieb-Oxford bound. This both rationalizes the use of the original Lieb-Oxford constant in PBE and suggests that enhancement factors more sensitive to sharpened constraints await discovery.

Empirical analysis of the Lieb-Oxford bound in ions and molecules

Universal properties of the Coulomb interaction energy apply to all many-electron systems. Bounds on the exchange-correlation energy, in particular, are important for the construction of improved density functionals. Here we investigate one such universal property-the Lieb-Oxford lower bound-for ionic and molecular systems. In recent work [J Chem Phys 127, 054106 (2007)], we observed that for atoms and electron liquids this bound may be substantially tightened. Calculations for a few ions and molecules suggested the same tendency, but were not conclusive due to the small number of systems considered. Here we extend that analysis to many different families of ions and molecules, and find that for these, too, the bound can be empirically tightened by a similar margin as for atoms and electron liquids. Tightening the Lieb-Oxford bound will have consequences for the performance of various approximate exchange-correlation functionals. (C) 2008 Wiley Periodicals Inc.

SPECTRUM OF SCREENING MASSES FOR SU(2) GAUGE THEORY: UNIVERSALITY CLASS

In this work we study the spectrum of the lowest screening masses for Yang-Mills theories on the lattice. We used the SU(2) gauge group in (3 + 1) dmensions. We adopted the multiple exponential method and the so-called ""variational"" method, in order to detect possible excited states. The calculations were done near the critical temperature of the confinement-deconfinement phase transition. We obtained values for the ratios of the screening masses consistent with predictions from universality arguments. A Monte Carlo evolution of the screening masses in the gauge theory confirms the validity of the predictions.

Solvatochromic Shift of the pi -> pi* Transition in All-Trans, Cis-13, Cis-11, Cis-9, and Cis-7 Retinal Isomers Induced by Water and Methanol

The solvatochromic shift of the lowest singlet it pi -> pi* electronic transition in the all-trans, cis-13, cis-11, cis-9, and cis-7 retinal isomers were computed under the influence of water, methanol, and benzene solvents. Excitation energies were calculated in gas phase and in solution. The calculations in solution were performed considering the sequential Monte Carlo (MC) /Quantum Mechanical approach. The MC simulations were performed considering the full retinal isomer molecules and 900 water molecules, 900 methanol, or 400 benzene ones. The OPLS/AA parametrization was chosen for retinal, methanol, and benzene molecules and the SPC model was used for water one. From the MC calculations 100 independent configurations were selected, with 100 solvent molecules in thermodynamical equilibrium at T = 298.15 K. Average point-charges were obtained from those independent configurations for water, methanol, and benzene solvent. TDDFT and CASSCF//CASPT2 methodologies were used to compute the vertical excitation energy of the retinal isomers in different environment. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110: 2076-2087, 2010