Liquid polymorphism, order-disorder transitions and anomalous behavior: A Monte Carlo study of the Bell-Lavis model for water

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution. (C) 2009 American Institute of Physics. [doi:10.1063/1.3253297]

Photon-hadron jet correlations in p plus p and Au plus Au collisions at s(NN)=200 GeV

We report the observation at the Relativistic Heavy Ion Collider of suppression of back-to-back correlations in the direct photon+jet channel in Au+Au relative to p+p collisions. Two-particle correlations of direct photon triggers with associated hadrons are obtained by statistical subtraction of the decay photon-hadron (gamma-h) background. The initial momentum of the away-side parton is tightly constrained, because the parton-photon pair exactly balance in momentum at leading order in perturbative quantum chromodynamics, making such correlations a powerful probe of the in-medium parton energy loss. The away-side nuclear suppression factor, I(AA), in central Au+Au collisions, is 0.32 +/- 0.12(stat)+/- 0.09(syst) for hadrons of 3 < p(T)(h)< 5 in coincidence with photons of 5 < p(T)(gamma)< 15 GeV/c. The suppression is comparable to that observed for high-p(T) single hadrons and dihadrons. The direct photon associated yields in p+p collisions scale approximately with the momentum balance, z(T)equivalent to p(T)(h)/p(T)(gamma), as expected for a measurement of the away-side parton fragmentation function. We compare to Au+Au collisions for which the momentum balance dependence of the nuclear modification should be sensitive to the path-length dependence of parton energy loss.

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.

Event structure and double helicity asymmetry in jet production from polarized p plus p collisions at root s=200 GeV

We report on the event structure and double helicity asymmetry (A(LL)) of jet production in longitudinally polarized p + p collisions at root s = 200 GeV. Photons and charged particles were measured by the PHENIX experiment at midrapidity vertical bar eta vertical bar < 0.35 with the requirement of a high-momentum (> 2 GeV/c) photon in the event. Event structure, such as multiplicity, p(T) density and thrust in the PHENIX acceptance, were measured and compared with the results from the PYTHIA event generator and the GEANT detector simulation. The shape of jets and the underlying event were well reproduced at this collision energy. For the measurement of jet A(LL), photons and charged particles were clustered with a seed-cone algorithm to obtain the cluster pT sum (p(T)(reco)). The effect of detector response and the underlying events on p(T)(reco) was evaluated with the simulation. The production rate of reconstructed jets is satisfactorily reproduced with the next-to-leading-order and perturbative quantum chromodynamics jet production cross section. For 4< p(T)(reco) < 12 GeV/c with an average beam polarization of < P > = 49% we measured Lambda(LL) = -0.0014 +/- 0.0037(stat) at the lowest p(T)(reco) bin (4-5 GeV= c) and -0.0181 +/- 0.0282(stat) at the highest p(T)(reco) bin (10-12 GeV= c) with a beam polarization scale error of 9.4% and a pT scale error of 10%. Jets in the measured p(T)(reco) range arise primarily from hard-scattered gluons with momentum fraction 0: 02 < x < 0: 3 according to PYTHIA. The measured A(LL) is compared with predictions that assume various Delta G(x) distributions based on the Gluck-Reya-Stratmann-Vogelsang parameterization. The present result imposes the limit -a.1 < integral(0.3)(0.02) dx Delta G(x, mu(2) = GeV2) < 0.4 at 95% confidence level or integral(0.3)(0.002) dx Delta G(x, mu(2) = 1 GeV2) < 0.5 at 99% confidence level.

Particle-species dependent modification of jet-induced correlations in Au plus Au collisions at root s(NN)=200 GeV

Measurements in Au + Au collisions at root s(NN) = 200 GeV of jet correlations for a trigger hadron at intermediate transverse momentum (p(T,trig)) with associated mesons or baryons at lower p(T,assoc) indicate strong modification of the away-side jet. The ratio of jet-associated baryons to mesons increases with centrality and p(T,assoc). For the most central collisions, the ratio is similar to that for inclusive measurements. This trend is incompatible with in-vacuum fragmentation but could be due to jetlike contributions from correlated soft partons, which recombine upon hadronization.

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.

Frustration and anomalous behavior in the Bell-Lavis model of liquid water

We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively). Both models present liquid-liquid coexistence and several thermodynamic anomalies. This result suggests that anisotropy introduced through orientational variables play no specific role in producing the density anomaly, in agreement with a similar conclusion discussed previously following results for continuous soft core,models. We propose that the presence of liquid anomalies may be related to energetic frustration, a feature common to both models.

Transverse momentum and centrality dependence of dihadron correlations in Au plus Au collisions at root s(NN)=200 GeV: Jet quenching and the response of partonic matter

Azimuthal angle (Delta phi) correlations are presented for charged hadrons from dijets for 0.4 < p(T)< 10 GeV/c in Au+Au collisions at root s(NN)=200 GeV. With increasing p(T), the away-side distribution evolves from a broad and relatively flat shape to a concave shape, then to a convex shape. Comparisons to p+p data suggest that the away-side can be divided into a partially suppressed ""head"" region centered at Delta phi similar to pi and an enhanced ""shoulder"" region centered at Delta phi similar to pi +/- 1.1. The p(T) spectrum for the head region softens toward central collisions, consistent with the onset of jet quenching. The spectral slope for the shoulder region is independent of centrality and trigger p(T), which offers constraints on energy transport mechanisms and suggests that it contains the medium response to energetic jets.

Adiabatic intramolecular movements for water systems

An effective treatment of the intramolecular degrees of freedom is presented for water, where these modes are decoupled from the intermolecular ones, ""adiabatically"" allowing these coordinates to be positioned at their local minimum of the potential energy surface. We perform ab initio Monte Carlo simulations with the configurational energies obtained via density functional theory. We study a water dimer as a prototype system, and even in this simple case the intramolecular relaxations are very important to properly describe properties such as the dipole moment. We show that rigid simulations do not correctly sample the phase space, resulting in an average dipole moment smaller than the one obtained with the adiabatic model, which is closer to the experimental result. (c) 2008 American Institute of Physics.

Experimental studies of di-jet survival and surface emission bias in Au plus Au collisions via angular correlations with respect to back-to-back leading hadrons

We report first results from an analysis based on a new multi-hadron correlation technique, exploring jet-medium interactions and di-jet surface emission bias at the BNL Relativistic Heavy Ion Collider (RHIC). Pairs of back-to-back high-transverse-momentum hadrons are used for triggers to study associated hadron distributions. In contrast with two-and three-particle correlations with a single trigger with similar kinematic selections, the associated hadron distribution of both trigger sides reveals no modification in either relative pseudorapidity Delta eta or relative azimuthal angle Delta phi from d + Au to central Au + Au collisions. We determine associated hadron yields and spectra as well as production rates for such correlated back-to-back triggers to gain additional insights on medium properties.

Long range rapidity correlations and jet production in high energy nuclear collisions

The STAR Collaboration at the Relativistic Heavy Ion Collider presents a systematic study of high-transverse-momentum charged-di-hadron correlations at small azimuthal pair separation Delta phi in d+Au and central Au+Au collisions at s(NN)=200 GeV. Significant correlated yield for pairs with large longitudinal separation Delta eta is observed in central Au+Au collisions, in contrast to d+Au collisions. The associated yield distribution in Delta eta x Delta phi can be decomposed into a narrow jet-like peak at small angular separation which has a similar shape to that found in d+Au collisions, and a component that is narrow in Delta phi and depends only weakly on Delta eta, the ""ridge."" Using two systematically independent determinations of the background normalization and shape, finite ridge yield is found to persist for trigger p(t)>6 GeV/c, indicating that it is correlated with jet production. The transverse-momentum spectrum of hadrons comprising the ridge is found to be similar to that of bulk particle production in the measured range (2 < p(t)< 4 GeV/c).

Ultra-high-energy cosmic rays from centaurus A: jet interaction with gaseous shells

Ultra-high-energy cosmic rays (UHECRs), with energies above similar to 6 x 10(19) eV, seem to show a weak correlation with the distribution of matter relatively near to us in the universe. It has earlier been proposed that UHECRs could be accelerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen A. We show that UHECR production at a spatially intermediate location about 15 kpc northeast from the nucleus, where the jet emerging from the nucleus is observed to strike a large star-forming shell of gas, is a plausible alternative. A relativistic jet is capable of accelerating lower energy heavy seed cosmic rays (CRs) to UHECRs on timescales comparable to the time it takes the jet to pierce the large gaseous cloud. In this model, many CRs arising from a starburst, with a composition enhanced in heavy elements near the knee region around PeV, are boosted to ultra-high energies by the relativistic shock of a newly oriented jet. This model matches the overall spectrum shown by the Auger data and also makes a prediction for the chemical composition as a function of particle energy. We thus predict an observable anisotropy in the composition at high energy in the sense that lighter nuclei should preferentially be seen toward the general direction of Cen A. Taking into consideration the magnetic field models for the Galactic disk and a Galactic magnetic wind, this scenario may resolve the discrepancy between HiRes and Auger results concerning the chemical composition of UHECRs.

Ultrasensitive thermal lens spectroscopy of water

A recently developed dual-beam configuration that optimizes the thermal lens technique has been used to obtain the absorption spectrum of pure water from 350 to 528 nm. Our results indicate the minimum linear absorption coefficient smaller than 2 X 10(-5) cm(-1) between 360 and 400 nm. This value is lower than previous literature data, and it is blueshifted. Absorption coefficients as small as 2 X 10(-7) cm(-1) can be measured for water using 1 W of excitation power. A detection limit of similar to 6 X 10(-9) cm(-1) (P=1 W) for CCl(4) was estimated, which represents, to the best of our knowledge, the highest sensitivity obtained in small absorption measurements in liquids. (C) 2009 Optical Society of America

A monohydrate pseudopolymorph of 3,4-dihydroxybenzophenone and the role of water in the crystal assembly of benzophenones

During a polymorphism screening of hydroxybenzophenone derivatives, a monohydrate pseudopolymorph of (3,4-dihydroxyphenyl)(phenyl)methanone, C(13)H(10)O(3)center dot H(2)O, (I), was obtained. Structural relationships and the role of water in crystal assembly were established on the basis of the known anhydrous form [Cox, Kechagias & Kelly (2008). Acta Cryst. B64, 206-216]. The crystal packing of (I) is stabilized by classical intermolecular O-H...O hydrogen bonds, generating a three-dimensional network.

Thermo-solvatochromism in binary mixtures of water and ionic liquids: on the relative importance of solvophobic interactions

The thermo-solvatochromism of 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), has been studied in mixtures of water, W, with ionic liquids, ILs, in the temperature range of 10 to 60 degrees C, where feasible. The objectives of the study were to test the applicability of a recently introduced solvation model, and to assess the relative importance of solute-solvent solvophobic interactions. The ILs were 1-allyl-3-alkylimidazolium chlorides, where the alkyl groups are methyl, 1-butyl, and 1-hexyl, respectively. The equilibrium constants for the interaction of W and the ILs were calculated from density data; they were found to be linearly dependent on N(C), the number of carbon atoms of the alkyl group; van't Hoff equation (log K versus 1/T) applied satisfactorily. Plots of the empirical solvent polarities, E(T) (MePMBr(2)) in kcal mol(-1), versus the mole fraction of water in the binary mixture, chi(w), showed non-linear, i.e., non-ideal behavior. The dependence of E(T) (MePMBr(2)) on chi(w), has been conveniently quantified in terms of solvation by W, IL, and the ""complex"" solvent IL-W. The non-ideal behavior is due to preferential solvation by the IL and, more efficiently, by IL-W. The deviation from linearity increases as a function of increasing N(C) of the IL, and is stronger than that observed for solvation of MePMBr(2) by aqueous 1-propanol, a solvent whose lipophilicity is 12.8 to 52.1 times larger than those of the ILs investigated. The dependence on N(C) is attributed to solute-solvent solvophobic interactions, whose relative contribution to solvation are presumably greater than that in mixtures of water and 1-propanol.