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.

Delocalized water and fluoride contributions to Dyson orbitals for electron detachment from the hydrated fluoride anion

The experimental vertical electron detachment energy (VEDE) of aqueous fluoride, [F(-)(H(2)O)], is approximately 9.8 eV, but spectral assignment is complicated by interference between F(-) 2p and H(2)O 1b(1) orbitals. The electronic structure of [F(-)(H(2)O)] is analyzed with Monte Carlo and ab initio quantum-mechanical calculations. Electron-propagator calculations in the partial third-order approximation yield a VEDE of 9.4 eV. None of the Dyson orbitals corresponding to valence VEDEs consists primarily of F 2p functions. These results and ground-state atomic charges indicate that the final, neutral state is more appropriately described as [F(-)(H(2)O)(+)] than as [F(H(2)O)]. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3431081]

Hyperpolarizabilities of the methanol molecule: A CCSD calculation including vibrational corrections

In this work we present the results for hyperpolarizabilities of the methanol molecule including vibrational corrections and electron correlation effects at the CCSD level. Comparisons to random phase approximation results previously reported show that the electron correlation is in general important for both electronic contribution and vibrational corrections. The role played by the anharmonicities on the calculations of the vibrational corrections has also been analyzed and the obtained results indicate that the anharmonic terms are important for the dc-Pockels and dc-Kerr effects. For the other nonlinear optical properties studied the double-harmonic approximation is found to be suitable. Comparison to available experimental result in gas phase for the dc-second harmonic generation second hyperpolarizability shows a very good agreement with the electronic contribution calculated here while our total value is 14% larger than the experimental value.

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.

Dynamic polarizability, Cauchy moments, and the optical absorption spectrum of liquid water: A sequential molecular dynamics/quantum mechanical approach

The dynamic polarizability and optical absorption spectrum of liquid water in the 6-15 eV energy range are investigated by a sequential molecular dynamics (MD)/quantum mechanical approach. The MD simulations are based on a polarizable model for liquid water. Calculation of electronic properties relies on time-dependent density functional and equation-of-motion coupled-cluster theories. Results for the dynamic polarizability, Cauchy moments, S(-2), S(-4), S(-6), and dielectric properties of liquid water are reported. The theoretical predictions for the optical absorption spectrum of liquid water are in good agreement with experimental information.

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.

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: A sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP/6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP/6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7 +/- 0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm. (C) 2008 American Institute of Physics.

Cold nuclear matter effects on J/psi production as constrained by deuteron-gold measurements at root S(NN)=200 GeV

We present a new analysis of J/psi production yields in deuteron-gold collisions at root s(NN) =200 GeV using data taken from the PHENIX experiment in 2003 and previously published in S. S. Adler [Phys. Rev. Lett 96, 012304 (2006)]. The high statistics proton-proton J/psi data taken in 2005 are used to improve the baseline measurement and thus construct updated cold nuclear matter modification factors (R(dAu)). A suppression of J/psi in cold nuclear matter is observed as one goes forward in rapidity (in the deuteron-going direction), corresponding to a region more sensitive to initial-state low-x gluons in the gold nucleus. The measured nuclear modification factors are compared to theoretical calculations of nuclear shadowing to which a J/psi (or precursor) breakup cross section is added. Breakup cross sections of sigma(breakup)=2.8(-1.4)(+1.7) (2.2(-1.5)(+1.6)) mb are obtained by fitting these calculations to the data using two different models of nuclear shadowing. These breakup cross-section values are consistent within large uncertainties with the 4.2 +/- 0.5 mb determined at lower collision energies. Projecting this range of cold nuclear matter effects to copper-copper and gold-gold collisions reveals that the current constraints are not sufficient to firmly quantify the additional hot nuclear matter effect.

Electronic properties of liquid ammonia: A sequential molecular dynamics/quantum mechanics approach

The electronic properties of liquid ammonia are investigated by a sequential molecular dynamics/quantum mechanics approach. Quantum mechanics calculations for the liquid phase are based on a reparametrized hybrid exchange-correlation functional that reproduces the electronic properties of ammonia clusters [(NH(3))(n); n=1-5]. For these small clusters, electron binding energies based on Green's function or electron propagator theory, coupled cluster with single, double, and perturbative triple excitations, and density functional theory (DFT) are compared. Reparametrized DFT results for the dipole moment, electron binding energies, and electronic density of states of liquid ammonia are reported. The calculated average dipole moment of liquid ammonia (2.05 +/- 0.09 D) corresponds to an increase of 27% compared to the gas phase value and it is 0.23 D above a prediction based on a polarizable model of liquid ammonia [Deng , J. Chem. Phys. 100, 7590 (1994)]. Our estimate for the ionization potential of liquid ammonia is 9.74 +/- 0.73 eV, which is approximately 1.0 eV below the gas phase value for the isolated molecule. The theoretical vertical electron affinity of liquid ammonia is predicted as 0.16 +/- 0.22 eV, in good agreement with the experimental result for the location of the bottom of the conduction band (-V(0)=0.2 eV). Vertical ionization potentials and electron affinities correlate with the total dipole moment of ammonia aggregates. (c) 2008 American Institute of Physics.

Identified charged hadron production in p plus p collisions at root s=200 and 62.4 GeV

Transverse momentum distributions and yields for pi(+/-), K(+/-), p, and (p) over bar in p + p collisions at root s = 200 and 62.4 GeV at midrapidity are measured by the PHENIX experiment at the Relativistic Heavy Ion Collider (RHIC). These data provide important baseline spectra for comparisons with identified particle spectra in heavy ion collisions at RHIC. We present the inverse slope parameter T(inv), mean transverse momentum < p(T)>, and yield per unit rapidity dN/dy at each energy, and compare them to other measurements at different root s in p + p and p + (p) over bar collisions. We also present the scaling properties such as m(T) scaling and x(T) scaling on the p(T) spectra between different energies. To discuss the mechanism of the particle production in p + p collisions, the measured spectra are compared to next-to-leading-order or next-to-leading-logarithmic perturbative quantum chromodynamics calculations.

Continuum-continuum coupling and polarization potentials for weakly bound systems

We investigate the influence of couplings among continuum states in collisions of weakly bound nuclei. For this purpose, we compare cross sections for complete fusion, breakup, and elastic scattering evaluated by continuum discretized coupled channel (CDCC) calculations, including and not including these couplings. In our study, we discuss this influence in terms of the polarization potentials that reproduces the elastic wave function of the coupled channel method in single channel calculations. We find that the inclusion of couplings among continuum states renders the real part of the polarization potential more repulsive, whereas it leads to weaker absorption to the breakup channel. We show that the noninclusion of continuum-continuum couplings in CDCC calculations may lead to qualitative and quantitative wrong conclusions.

KINETICS VERSUS HYDRODYNAMICS: GENERALIZATION OF LANDAU/COOPER-FRYE PRESCRIPTION FOR FREEZE-OUT

The problem of spectra formation in hydrodynamic approach to A + A collisions is considered within the Boltzmann equations. It is shown analytically and illustrated by numerical calculations that the particle momentum spectra can be presented in the Cooper-R-ye form despite freeze-out is not sharp and has the finite temporal width. The latter is equal to the inverse of the particle collision rate at points (t(sigma) (r, p), r) of the maximal emission at a fixed momentum p. The set of these points forms the hypersurfaces t(sigma)(r,p) which strongly depend on the values of p and typically do not enclose completely the initially dense matter. This is an important difference from the standard Cooper-Frye prescription (CFp), with a common freeze-out hypersurface for all p, that affects significantly the predicted spectra. Also, the well known problem of CFp as for negative contributions to the spectra from non-space-like parts of the freeze-out hypersurface is naturally eliminated in this improved prescription.

Kinetic energy sum spectra in nonmesonic weak decay of hypernuclei

We evaluate the coincidence spectra in the nonmesonic weak decay (NMWD) Lambda N -> nN of Lambda hypernuclei (4)(Lambda)He, (5)(Lambda)He, (12)(Lambda)C, (16)(Lambda)O, and (28)(Lambda)Si, as a function of the sum of kinetic energies E(nN)=E(n)+E(N) for N=n,p. The strangeness-changing transition potential is described by the one-meson-exchange model, with commonly used parametrization. Two versions of the independent-particle shell model (IPSM) are employed to account for the nuclear structure of the final residual nuclei. They are as follows: (a) IPSM-a, where no correlation, except for the Pauli principle, is taken into account and (b) IPSM-b, where the highly excited hole states are considered to be quasistationary and are described by Breit-Wigner distributions, whose widths are estimated from the experimental data. All np and nn spectra exhibit a series of peaks in the energy interval 110 MeV < E(nN)< 170 MeV, one for each occupied shell-model state. Within the IPSM-a, and because of the recoil effect, each peak covers an energy interval proportional to A(-1) , going from congruent to 4 MeV for (28)(Lambda)Si to congruent to 40 MeV for (4)(Lambda)He. Such a description could be pretty fair for the light (4)(Lambda)He and (5)(Lambda)He hypernuclei. For the remaining, heavier, hypernuclei it is very important, however, to consider as well the spreading in strength of the deep-hole states and bring into play the IPSM-b approach. Notwithstanding the nuclear model that is employed the results depend only very weakly on the details of the dynamics involved in the decay process proper. We propose that the IPSM is the appropriate lowest-order approximation for the theoretical calculations of the of kinetic energy sum spectra in the NMWD. It is in comparison to this picture that one should appraise the effects of the final-state interactions and of the two-nucleon-induced decay mode.