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.

Missing levels in acoustic resonators

It is shown that the deviations of the experimental statistics of six chaotic acoustic resonators from Wigner-Dyson random matrix theory predictions are explained by a recent model of random missing levels. In these resonatorsa made of aluminum plates a the larger deviations occur in the spectral rigidity (SRs) while the nearest-neighbor distributions (NNDs) are still close to the Wigner surmise. Good fits to the experimental NNDs and SRs are obtained by adjusting only one parameter, which is the fraction of remaining levels of the complete spectra. For two Sinai stadiums, one Sinai stadium without planar symmetry, two triangles, and a sixth of the three-leaf clover shapes, was found that 7%, 4%, 7%, and 2%, respectively, of eigenfrequencies were not detected.

Onset of pi(0) Suppression Studied in Cu plus Cu Collisions at root s(NN)=22.4, 62.4, and 200 GeV

Neutral pion transverse momentum (p(T)) spectra at midrapidity (|y| less than or similar to 0.35) were measured in Cu + Cu collisions at root s(NN) = 22.4, 62.4, and 200 GeV. Relative to pi(0) yields in p + p collisions scaled by the number of inelastic nucleon-nucleon collisions (N(coll)) the pi(0) yields for p(T) greater than or similar to 2 GeV/c in central Cu + Cu collisions are suppressed at 62.4 and 200 GeV whereas an enhancement is observed at 22.4 GeV. A comparison with a jet-quenching model suggests that final state parton energy loss dominates in central Cu + Cu collisions at 62.4 and 200 GeV, while the enhancement at 22.4 GeV is consistent with nuclear modifications in the initial state alone.

Centrality dependence of charged hadron production in deuteron plus gold and nucleon plus gold collisions at root S(NN)=200 GeV

We present transverse momentum (p(T)) spectra of charged hadrons measured in deuteron-gold and nucleon-gold collisions at root s(NN)=200 GeV for four centrality classes. Nucleon-gold collisions were selected by tagging events in which a spectator nucleon was observed in one of two forward rapidity detectors. The spectra and yields were investigated as a function of the number of binary nucleon-nucleon collisions, nu, suffered by deuteron nucleons. A comparison of charged particle yields to those in p+p collisions show that yield per nucleon-nucleon collision saturates with nu for high momentum particles. We also present the charged hadron to neutral pion ratios as a function of p(T).

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.

CMB in a box: Causal structure and the Fourier-Bessel expansion

This paper makes two points. First, we show that the line-of-sight solution to cosmic microwave anisotropies in Fourier space, even though formally defined for arbitrarily large wavelengths, leads to position-space solutions which only depend on the sources of anisotropies inside the past light cone of the observer. This foretold manifestation of causality in position (real) space happens order by order in a series expansion in powers of the visibility gamma = e(-mu), where mu is the optical depth to Thomson scattering. We show that the contributions of order gamma(N) to the cosmic microwave background (CMB) anisotropies are regulated by spacetime window functions which have support only inside the past light cone of the point of observation. Second, we show that the Fourier-Bessel expansion of the physical fields (including the temperature and polarization momenta) is an alternative to the usual Fourier basis as a framework to compute the anisotropies. The viability of the Fourier-Bessel series for treating the CMB is a consequence of the fact that the visibility function becomes exponentially small at redshifts z >> 10(3), effectively cutting off the past light cone and introducing a finite radius inside which initial conditions can affect physical observables measured at our position (x) over right arrow = 0 and time t(0). Hence, for each multipole l there is a discrete tower of momenta k(il) (not a continuum) which can affect physical observables, with the smallest momenta being k(1l) similar to l. The Fourier-Bessel modes take into account precisely the information from the sources of anisotropies that propagates from the initial value surface to the point of observation-no more, no less. We also show that the physical observables (the temperature and polarization maps), and hence the angular power spectra, are unaffected by that choice of basis. This implies that the Fourier-Bessel expansion is the optimal scheme with which one can compute CMB anisotropies.

Deformed Gaussian-orthogonal-ensemble description of small-world networks

The study of spectral behavior of networks has gained enthusiasm over the last few years. In particular, random matrix theory (RMT) concepts have proven to be useful. In discussing transition from regular behavior to fully chaotic behavior it has been found that an extrapolation formula of the Brody type can be used. In the present paper we analyze the regular to chaotic behavior of small world (SW) networks using an extension of the Gaussian orthogonal ensemble. This RMT ensemble, coined the deformed Gaussian orthogonal ensemble (DGOE), supplies a natural foundation of the Brody formula. SW networks follow GOE statistics until a certain range of eigenvalue correlations depending upon the strength of random connections. We show that for these regimes of SW networks where spectral correlations do not follow GOE beyond a certain range, DGOE statistics models the correlations very well. The analysis performed in this paper proves the utility of the DGOE in network physics, as much as it has been useful in other physical systems.

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.

Directed flow at midrapidity in event-by-event hydrodynamics

Fluctuations in the initial geometry of a nucleus-nucleus collision have been recently shown to result in a new type of directed flow (v(1)) that, unlike the usual directed flow, is also present at midrapidity. We compute this new v(1) versus transverse momentum and centrality for Au-Au collisions at RHIC using the hydrodynamic code NeXSPheRIO. We find that the event plane of v(1) is correlated with the angle of the initial dipole of the distribution, as predicted, though with a large dispersion. It is uncorrelated with the reaction plane. Our results are in excellent agreement with results inferred from STAR correlation data.

Hydro-kinetic approach to relativistic heavy ion collisions

We develop a combined hydro-kinetic approach which incorporates a hydrodynamical expansion of the systems formed in A + A collisions and their dynamical decoupling described by escape probabilities. The method corresponds to a generalized relaxation time (tau(rel)) approximation for the Boltzmann equation applied to inhomogeneous expanding systems; at small tau(rel) it also allows one to catch the viscous effects in hadronic component-hadron-resonance gas. We demonstrate how the approximation of sudden freeze-out can be obtained within this dynamical picture of continuous emission and find that hypersurfaces, corresponding to a sharp freeze-out limit, are momentum dependent. The pion m(T) spectra are computed in the developed hydro-kinetic model, and compared with those obtained from ideal hydrodynamics with the Cooper-Frye isothermal prescription. Our results indicate that there does not exist a universal freeze-out temperature for pions with different momenta, and support an earlier decoupling of higher p(T) particles. By performing numerical simulations for various initial conditions and equations of state we identify several characteristic features of the bulk QCD matter evolution preferred in view of the current analysis of heavy ion collisions at RHIC energies.

Quantum fields with noncommutative target spaces

Quantum field theories (QFT's) on noncommutative spacetimes are currently under intensive study. Usually such theories have world sheet noncommutativity. In the present work, instead, we study QFT's with commutative world sheet and noncommutative target space. Such noncommutativity can be interpreted in terms of twisted statistics and is related to earlier work of Oeckl [R. Oeckl, Commun. Math. Phys. 217, 451 (2001)], and others [A. P. Balachandran, G. Mangano, A. Pinzul, and S. Vaidya, Int. J. Mod. Phys. A 21, 3111 (2006); A. P. Balachandran, A. Pinzul, and B. A. Qureshi, Phys. Lett. B 634,434 (2006); A.P. Balachandran, A. Pinzul, B.A. Qureshi, and S. Vaidya, arXiv:hep-th/0608138; A.P. Balachandran, T. R. Govindarajan, G. Mangano, A. Pinzul, B.A. Qureshi, and S. Vaidya, Phys. Rev. D 75, 045009 (2007); A. Pinzul, Int. J. Mod. Phys. A 20, 6268 (2005); G. Fiore and J. Wess, Phys. Rev. D 75, 105022 (2007); Y. Sasai and N. Sasakura, Prog. Theor. Phys. 118, 785 (2007)]. The twisted spectra of their free Hamiltonians has been found earlier by Carmona et al. [J. M. Carmona, J. L. Cortes, J. Gamboa, and F. Mendez, Phys. Lett. B 565, 222 (2003); J. M. Carmona, J. L. Cortes, J. Gamboa, and F. Mendez, J. High Energy Phys. 03 (2003) 058]. We review their derivation and then compute the partition function of one such typical theory. It leads to a deformed blackbody spectrum, which is analyzed in detail. The difference between the usual and the deformed blackbody spectrum appears in the region of high frequencies. Therefore we expect that the deformed blackbody radiation may potentially be used to compute a Greisen-Zatsepin-Kuzmin cutoff which will depend on the noncommutative parameter theta.

Disordered ensembles of random matrices

It is shown that the families of generalized matrix ensembles recently considered which give rise to an orthogonal invariant stable Levy ensemble can be generated by the simple procedure of dividing Gaussian matrices by a random variable. The nonergodicity of this kind of disordered ensembles is investigated. It is shown that the same procedure applied to random graphs gives rise to a family that interpolates between the Erdos-Renyi and the scale free models.

Optical third-harmonic spectroscopy of the magnetic semiconductor EuTe

EuTe possesses the centrosymmetric crystal structure m3m of rocksalt type in which the second-harmonic generation is forbidden in electric dipole approximation but the third-harmonic generation (THG) is allowed. We studied the THG spectra of this material and observed several resonances in the vicinity of the band gap at 2.2-2.5 eV and at higher energies up to 4 eV, which are related to four-photon THG processes. The observed resonances are assigned to specific combinations of electronic transitions between the ground 4f(7) state at the top of the valence band and excited 4f(6)5d(1) states of Eu(2+) ions, which form the lowest energy conduction band. Temperature, magnetic field, and rotational anisotropy studies allowed us to distinguish crystallographic and magnetic-field-induced contributions to the THG. A strong modification of THG intensity for the 2.4 eV band and suppression of the THG for the 3.15 eV band was observed in applied magnetic field. Two main features of the THG spectra were assigned to 5d(t(2g)) and 5d(e(g)) subbands at 2.4 eV and 3.15 eV, respectively. A microscopic quantum-mechanical model of the THG response was developed and its conclusions are in qualitative agreement with the experimental results.