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.

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.

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.

Use of correlated potential harmonic basis functions for the description of the (4)He trimer and small clusters

A correlated many-body basis function is used to describe the (4)He trimer and small helium clusters ((4)HeN) with N = 4-9. A realistic helium dimer potential is adopted. The ground state results of the (4)He dimer and trimer are in close agreement with earlier findings. But no evidence is found for the existence of Efimov state in the trimer for the actual (4)He-(4)He interaction. However, decreasing the potential strength we calculate several excited states of the trimer which exhibit Efimov character. We also solve for excited state energies of these clusters which are in good agreement with Monte Carlo hyperspherical description. (C) 2011 American Institute of Physics. [doi:10.1063/1.3583365]

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.

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.

Optical second harmonic generation in the centrosymmetric magnetic semiconductors EuTe and EuSe

The magnetic europium chalcogenide semiconductors EuTe and EuSe are investigated by the spectroscopy of second harmonic generation (SHG) in the vicinity of the optical band gap formed by transitions involving the 4f and 5d electronic orbitals of the magnetic Eu(2+) ions. In these materials with centrosymmetric crystal lattice the electric-dipole SHG process is symmetry forbidden so that no signal is observed in zero magnetic field. Signal appears, however, in applied magnetic field with the SHG intensity being proportional to the square of magnetization. The magnetic field and temperature dependencies of the induced SHG allow us to introduce a type of nonlinear optical susceptibility determined by the magnetic-dipole contribution in combination with a spontaneous or induced magnetization. The experimental results can be described qualitatively by a phenomenological model based on a symmetry analysis and are in good quantitative agreement with microscopic model calculations accounting for details of the electronic energy and spin structure.

Atomistic molecular dynamics study of interface formation: Al on poly(p-phenylene vinylene)

We model interface formation by metal deposition on the conjugated polymer poly-para-phenylene vinylene, studying direct aluminum and layered aluminum-calcium structures Al/PPV and Al/Ca/PPV. To do that we use classical molecular dynamics simulations, checked by ab initio density-functional theory calculations, for selected relevant configurations. We find that Al not only migrates easily into the film, with a strong charge transfer to the neighboring chains, but also promotes rearrangement of the polymer in the interfacial region to the hexagonal structure. On the other hand, our results indicate that a thin Ca layer is sufficient to protect the film and maintain a well-defined metal/polymer interface, and that also a thin Al capping layer may protect the whole from environmental degradation.

Optical absorption and electronic band structure first-principles calculations of alpha-glycine crystals

Light absorption of alpha-glycine crystals grown by slow evaporation at room temperature was measured, indicating a 5.11 +/- 0.02 eV energy band gap. Structural, electronic, and optical absorption properties of alpha-glycine crystals were obtained by first-principles quantum mechanical calculations using density functional theory within the generalized gradient approximation in order to understand this result. To take into account the contribution of core electrons, ultrasoft and norm-conserving pseudopotentials, as well as an all electron approach were considered to compute the electronic density of states and band structure of alpha-glycine crystals. They exhibit three indirect energy band gaps and one direct Gamma-Gamma energy gap around 4.95 eV. The optical absorption related to transitions between the top of the valence band and the bottom of the conduction band involves O 2p valence states and C, O 2p conduction states, with the carboxyl group contributing significantly to the origin of the energy band gap. The calculated optical absorption is highly dependent on the polarization of the incident radiation due to the spatial arrangement of the dipolar glycine molecules; in the case of a polycrystalline sample, the first-principles calculated optical absorption is in good agreement with the measurement when a rigid energy shift is applied.

Scaling properties at freeze-out in relativistic heavy-ion collisions

Identified charged pion, kaon, and proton spectra are used to explore the system size dependence of bulk freeze-out properties in Cu + Cu collisions at root s(NN) = 200 and 62.4 GeV. The data are studied with hydrodynamically motivated blast-wave and statistical model frameworks in order to characterize the freeze-out properties of the system. The dependence of freeze-out parameters on beam energy and collision centrality is discussed. Using the existing results from Au + Au and pp collisions, the dependence of freeze-out parameters on the system size is also explored. This multidimensional systematic study furthers our understanding of the QCD phase diagram revealing the importance of the initial geometrical overlap of the colliding ions. The analysis of Cu + Cu collisions expands the system size dependence studies from Au + Au data with detailed measurements in the smaller system. The systematic trends of the bulk freeze-out properties of charged particles is studied with respect to the total charged particle multiplicity at midrapidity, exploring the influence of initial state effects.

High p(T) nonphotonic electron production in p plus p collisions at root s=200 GeV

We present the measurement of nonphotonic electron production at high transverse momentum (p(T) > 2.5 GeV/c) in p + p collisions at root s = 200 GeV using data recorded during 2005 and 2008 by the STAR experiment at the Relativistic Heavy Ion Collider (RHIC). The measured cross sections from the two runs are consistent with each other despite a large difference in photonic background levels due to different detector configurations. We compare the measured nonphotonic electron cross sections with previously published RHIC data and perturbative quantum chromodynamics calculations. Using the relative contributions of B and D mesons to nonphotonic electrons, we determine the integrated cross sections of electrons (e++e-2/2) at 3 GeV/c < p(T) < 10 GeV/c from bottom and charm meson decays to be [(d sigma((B -> e)+(B -> D -> e))/(dy(e))](ye=0) 4.0 +/- 0.5(stat) +/- 1.1(syst) nb and [(d sigma(D -> e))/(dy(e))](ye=0) = 6.2 +/- 0.7(stat) +/- 1.5(syst) nb, respectively.

K*(0) production in Cu plus Cu and Au plus Au collisions at root s(NN)=62.4 GeV and 200 GeV

We report on K*(0) production at midrapidity in Au + Au and Cu + Cu collisions at root s(NN) = 62.4 and 200 GeV collected by the Solenoid Tracker at the Relativistic Heavy Ion Collider detector. The K*(0) is reconstructed via the hadronic decays K*(0) -> K(+)pi(-) and (K*(0)) over bar -> K(+)pi(-). Transverse momentum, p(T), spectra are measured over a range of p(T) extending from 0.2 GeV/c up to 5 GeV/c. The center-of-mass energy and system size dependence of the rapidity density, dN/dy, and the average transverse momentum, < p(T)>, are presented. The measured N(K*(0))/N(K) and N(phi)/N(K*(0)) ratios favor the dominance of rescattering of decay daughters of K*(0) over the hadronic regeneration for the K*(0) production. In the intermediate p(T) region (2.0 < p(T) < 4.0 GeV/c), the elliptic flow parameter, v(2), and the nuclear modification factor, R(CP), agree with the expectations from the quark coalescence model of particle production.

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).

Possible hyperdeformed band in (36)Ar observed in (12)C+(24)Mg elastic scattering

Fifteen strongly oscillating angular distributions of the elastic scattering of (12)C + (24)Mg at energies around the Coulomb barrier (E(c.m). = 10.67-16.00 MeV) are reproduced by adding five Breit-Wigner resonance terms to the l = 2, 4, 6, 7, and 8 elastic S matrix. The nonresonant, background elastic scattering S matrix S(l)(0) is calculated using the Sao Paulo potential. The J = 2, 4, 6, 7, and 8 (h) over bar molecular resonances fit well into a rotational molecular band, together with other higher lying resonances observed in the (16)O + (20)Ne elastic scattering. We propose that the presently observed, largely deformed molecular band corresponds to the hyperdeformed band, which has been found previously in alpha-cluster calculations, as well as in a new Nilsson model calculation. Systematic study of its possible clusterizations predicts the preference of the (12)C + (24)Mg and (16)O + (20)Ne molecular structure, in accordance with our present results.

K/pi Fluctuations at Relativistic Energies

We report K/pi fluctuations from Au+Au collisions at s(NN)=19.6, 62.4, 130, and 200 GeV using the STAR detector at the Relativistic Heavy Ion Collider. K/pi fluctuations in central collisions show little dependence on incident energy and are on the same order as those from NA49 at the Super Proton Synchrotron in central Pb+Pb collisions at s(NN)=12.3 and 17.3 GeV. We report results for the collision centrality dependence of K/pi fluctuations and results for charge-separated fluctuations. We observe that the K/pi fluctuations scale with the charged particle multiplicity density.