Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

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.

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.

Nuclear incoherent photoproduction of pi(0) and eta from 4 to 12 GeV

The mechanism of incoherent pi(0) and eta photoproduction from complex nuclei is investigated from 4 to 12 GeV with an extended version of the multicollisional Monte Carlo (MCMC) intranuclear cascade model. The calculations take into account the elementary photoproduction amplitudes via a Regge model and the nuclear effects of photon shadowing, Pauli blocking, and meson-nucleus final-state interactions. The results for pi(0) photoproduction reproduced for the first time the magnitude and energy dependence of the measured rations sigma(gamma A)/sigma(gamma N) for several nuclei (Be, C, Al, Cu, Ag, and Pb) from a Cornell experiment. The results for eta photoproduction fitted the inelastic background in Cornell's yields remarkably well, which is clearly not isotropic as previously considered in Cornell's analysis. With this constraint for the background, the eta -> gamma gamma. decay width was extracted using the Primakoff method, combining Be and Cu data [Gamma(eta ->gamma gamma) = 0.476(62) keV] and using Be data only [Gamma(eta ->gamma gamma) = 0.512(90) keV]; where the errors are only statistical. These results are in sharp contrast (similar to 50-60%) with the value reported by the Cornell group [Gamma(eta ->gamma gamma). = 0.324(46) keV] and in line with the Particle Data Group average of 0.510(26) keV.

Detailed study of the K(L)->pi(0)pi(0)pi(0) Dalitz plot

Using a sample of 68.3x10(6) K(L)->pi(0)pi(0)pi(0) decays collected in 1996-1999 by the KTeV (E832) experiment at Fermilab, we present a detailed study of the K(L)->pi(0)pi(0)pi(0) Dalitz plot density. We report the first observation of interference from K(L)->pi(+)pi(-)pi(0) decays in which pi(+)pi(-) rescatters to pi(0)pi(0) in a final-state interaction. This rescattering effect is described by the Cabibbo-Isidori model, and it depends on the difference in pion scattering lengths between the isospin I=0 and I=2 states, a(0)-a(2). Using the Cabibbo-Isidori model, and fixing (a(0)-a(2))m(pi)(+)=0.268 +/- 0.017 as measured by the CERN-NA48 collaboration, we present the first measurement of the K(L)->pi(0)pi(0)pi(0) quadratic slope parameter that accounts for the rescattering effect: h(000)=(+0.59 +/- 0.20(stat)+/- 0.48(syst)+/- 1.06(ext))x10(-3), where the uncertainties are from data statistics, KTeV systematic errors, and external systematic errors. Fitting for both h(000) and a(0)-a(2), we find h(000)=(-2.09 +/- 0.62(stat)+/- 0.72(syst)+/- 0.28(ext))x10(-3), and m(pi)(+)(a(0)-a(2))=0.215 +/- 0.014(stat)+/- 0.025(syst)+/- 0.006(ext); our value for a(0)-a(2) is consistent with that from NA48.

Search for the rare decay K(L)->pi(0)pi(0)gamma

The KTeV E799 experiment has conducted a search for the rare decay K(L)->pi(0)pi(0)gamma via the topology K(L)->pi(0)pi(0)(D)gamma (where pi(0)(D)->gamma e(+)e(-)). Because of Bose statistics of the pi(0) pair and the real nature of the photon, the K(L)->pi(0)pi(0)gamma decay is restricted to proceed at lowest order by the CP conserving direct emission (DE) of an E2 electric quadrupole photon. The rate of this decay is interesting theoretically since chiral perturbation theory predicts that this process vanishes at level O(p(4)). Therefore, this mode probes chiral perturbation theory at O(p(6)). In this paper we report a determination of an upper limit of 2.43x10(-7) (90% CL) for K(L)->pi(0)pi(0)gamma. This is approximately a factor of 20 lower than previous results.

Final results from the KTeV experiment on the decay KL ->pi(0)gamma gamma

We report on a new measurement of the branching ratio B(K(L) -> pi(0)gamma gamma) using the KTeV detector. We reconstruct 1982 events with an estimated background of 608, that results in B(K(L) -> pi(0)gamma gamma)=(1.29 +/- 0.03(stat) +/- 0.05(syst)) x 10(-6). We also measure the parameter, a(V), which characterizes the strength of vector meson exchange terms in this decay. We find a(V) = -0.31 +/- 0.05(stat) +/- 0.07(syst). These results utilize the full KTeV data set collected from 1997 to 2000 and supersede earlier KTeV measurements of the branching ratio and a(V).

Cross section and double helicity asymmetry for eta mesons and their comparison to pi(0) production in p plus p collisions at root s=200 GeV

Measurements of double-helicity asymmetries in inclusive hadron production in polarized p + p collisions are sensitive to helicity-dependent parton distribution functions, in particular, to the gluon helicity distribution, Delta g. This study focuses on the extraction of the double-helicity asymmetry in eta production ((p) over right arrow + (p) over right arrow -> eta + X), the eta cross section, and the eta/pi(0) cross section ratio. The cross section and ratio measurements provide essential input for the extraction of fragmentation functions that are needed to access the helicity-dependent parton distribution functions.

High p(T) direct photon and pi(0) triggered azimuthal jet correlations and measurement of k(T) for isolated direct photons in p plus p collisions at root s=200 GeV

Correlations of charged hadrons of 1< p(T) < 10 Gev/c with high pT direct photons and pi(0) mesons in the range 5< p(T) < 15 Gev/c are used to study jet fragmentation in the gamma + jet and dijet channels, respectively. The magnitude of the partonic transverse momentum, k(T), is obtained by comparing to a model incorporating a Gaussian kT smearing. The sensitivity of the associated charged hadron spectra to the underlying fragmentation function is tested and the data are compared to calculations using recent global fit results. The shape of the direct photon-associated hadron spectrum as well as its charge asymmetry are found to be consistent with a sample dominated by quark-gluon Compton scattering. No significant evidence of fragmentation photon correlated production is observed within experimental uncertainties.

Azimuthal Anisotropy of pi(0) Production in Au plus Au Collisions at root s(NN)=200 GeV: Path-Length Dependence of Jet Quenching and the Role of Initial Geometry

We have measured the azimuthal anisotropy of pi(0) production for 1 < p(T) < 18 GeV/c for Au + Au collisions at root s(NN) = 200 GeV. The observed anisotropy shows a gradual decrease for 3 less than or similar to p(T) less than or similar to 7-10 GeV/c, but remains positive beyond 10 GeV/c. The magnitude of this anisotropy is underpredicted, up to at least similar to 10 GeV/c, by current perturbative QCD (PQCD) energy-loss model calculations. An estimate of the increase in anisotropy expected from initial-geometry modification due to gluon saturation effects and fluctuations is insufficient to account for this discrepancy. Calculations that implement a path-length dependence steeper than what is implied by current PQCD energy-loss models show reasonable agreement with the data.

High-p(T) pi(0) production with respect to the reaction plane in Au plus Au collisions at s(NN)=200 GeV

Measurements of the azimuthal anisotropy of high-p(T) neutral pion (pi(0)) production in Au+Au collisions at s(NN)=200 GeV by the PHENIX experiment are presented. The data included in this article were collected during the 2004 Relativistic Heavy Ion Collider running period and represent approximately an order of magnitude increase in the number of analyzed events relative to previously published results. Azimuthal angle distributions of pi(0) mesons detected in the PHENIX electromagnetic calorimeters are measured relative to the reaction plane determined event-by-event using the forward and backward beam-beam counters. Amplitudes of the second Fourier component (v(2)) of the angular distributions are presented as a function of pi(0) transverse momentum (p(T)) for different bins in collision centrality. Measured reaction plane dependent pi(0) yields are used to determine the azimuthal dependence of the pi(0) suppression as a function of p(T), R(AA)(Delta phi,p(T)). A jet-quenching motivated geometric analysis is presented that attempts to simultaneously describe the centrality dependence and reaction plane angle dependence of the pi(0) suppression in terms of the path lengths of hypothetical parent partons in the medium. This set of results allows for a detailed examination of the influence of geometry in the collision region and of the interplay between collective flow and jet-quenching effects along the azimuthal axis.

Gluon-Spin Contribution to the Proton Spin from the Double-Helicity Asymmetry in Inclusive pi(0) Production in Polarized p plus p Collisions at root s=200 GeV

The double helicity asymmetry in neutral pion production for p(T) = 1 to 12 GeV/c was measured with the PHENIX experiment to access the gluon-spin contribution, Delta G, to the proton spin. Measured asymmetries are consistent with zero, and at a theory scale of mu 2 = 4 GeV(2) a next to leading order QCD analysis gives Delta G([0.02,0.3]) = 0.2, with a constraint of -0.7 < Delta G([0.02,0.3]) < 0.5 at Delta chi(2) = 9 (similar to 3 sigma) for the sampled gluon momentum fraction (x) range, 0.02 to 0.3. The results are obtained using predictions for the measured asymmetries generated from four representative fits to polarized deep inelastic scattering data. We also consider the dependence of the Delta G constraint on the choice of the theoretical scale, a dominant uncertainty in these predictions.

Inclusive cross section and double helicity asymmetry for pi(0) production in p plus p collisions at root s=62.4 GeV

The PHENIX experiment presents results from the RHIC 2006 run with polarized p + p collisions at root s = 62.4 GeV, for inclusive pi(0) production at midrapidity. Unpolarized cross section results are measured for transverse momenta p(T) = 0.5 to 7 GeV/c. Next-to-leading order perturbative quantum chromodynamics calculations are compared with the data, and while the calculations are consistent with the measurements, next-to-leading logarithmic corrections improve the agreement. Double helicity asymmetries A(LL) are presented for p(T) = 1 to 4 GeV/c and probe the higher range of Bjorken x of the gluon (x(g)) with better statistical precision than our previous measurements at root s = 200 GeV. These measurements are sensitive to the gluon polarization in the proton for 0.06 < x(g) < 0.4.

Combined Monte Carlo and quantum mechanics study of the solvatochromism of phenol in water. The origin of the blue shift of the lowest pi-pi* transition

A combined and sequential use of Monte Carlo simulations and quantum mechanical calculations is made to analyze the spectral shift of the lowest pi-pi* transition of phenol in water. The solute polarization is included using electrostatic embedded calculations at the MP2/aug-cc-pVDZ level giving a dipole moment of 2.25 D, corresponding to an increase of 76% compared to the calculated gas-phase value. Using statistically uncorrelated configurations sampled from the MC simulation,first-principle size-extensive calculations are performed to obtain the solvatochromic shift. Analysis is then made of the origin of the blue shift. Results both at the optimized geometry and in room-temperature liquid water show that hydrogen bonds of water with phenol promote a red shift when phenol is the proton-donor and a blue shift when phenol is the proton-acceptor. In the case of the optimized clusters the calculated shifts are in very good agreement with results obtained from mass-selected free jet expansion experiments. In the liquid case the contribution of the solute-solvent hydrogen bonds partially cancels and the total shift obtained is dominated by the contribution of the outer solvent water molecules. Our best result, including both inner and outer water molecules, is 570 +/- 35 cm(-1), in very good agreement with the small experimental shift of 460 cm(-1) for the absorption maximum.