Analyzing the n→π* Electronic Transition of Formaldehyde in Water. A Sequential Monte Carlo/Time-Dependent Density Functional Theory

The n→π* absorption transition of formaldehyde in water is analyzed using combined and sequential classical Monte Carlo (MC) simulations and quantum mechanics (QM) calculations. MC simulations generate the liquid solute-solvent structures for subsequent QM calculations. Using time-dependent density functional theory in a localized set of gaussian basis functions (TD-DFT/6-311++G(d,p)) calculations are made on statistically relevant configurations to obtain the average solvatochromic shift. All results presented here use the electrostatic embedding of the solvent. The statistically converged average result obtained of 2300 cm-1 is compared to previous theoretical results available. Analysis is made of the effective dipole moment of the hydrogen-bonded shell and how it could be held responsible for the polarization of the solvent molecules in the outer solvation shells.

System-Size Independence of Directed Flow Measured at the BNL Relativistic Heavy-Ion Collider

We measure directed flow (v(1)) for charged particles in Au + Au and Cu + Cu collisions at root s(NN) = 200 and 62.4 GeV, as a function of pseudorapidity (eta), transverse momentum (p(t)), and collision centrality, based on data from the STAR experiment. We find that the directed flow depends on the incident energy but, contrary to all available model implementations, not on the size of the colliding system at a given centrality. We extend the validity of the limiting fragmentation concept to v(1) in different collision systems, and investigate possible explanations for the observed sign change in v(1)(p(t)).

Spin alignment measurements of the K*(0)(892) and phi(1020) vector mesons in heavy ion collisions at root S(NN)=200 GeV

We present the first spin alignment measurements for the K*(0)(892) and phi(1020) vector mesons produced at midrapidity with transverse momenta up to 5 GeV/c at root s(NN) = 200 GeV at RHIC. The diagonal spin-density matrix elements with respect to the reaction plane in Au+Au collisions are rho(00) = 0.32 +/- 0.04 (stat) +/- 0.09 (syst) for the K*(0) (0.8 < p(T) < 5.0 GeV/c) and rho(00) = 0.34 +/- 0.02 (stat) +/- 0.03 (syst) for the phi (0.4 < p(T) < 5.0 GeV/c) and are constant with transverse momentum and collision centrality. The data are consistent with the unpolarized expectation of 1/3 and thus no evidence is found for the transfer of the orbital angular momentum of the colliding system to the vector-meson spins. Spin alignments for K(*0) and phi in Au+Au collisions were also measured with respect to the particle's production plane. The phi result, rho(00) = 0.41 +/- 0.02 (stat) +/- 0.04 (syst), is consistent with that in p+p collisions, rho(00) = 0.39 +/- 0.03 (stat) +/- 0.06 (syst), also measured in this work. The measurements thus constrain the possible size of polarization phenomena in the production dynamics of vector mesons.

Korteveg-de Vries solitons in a cold quark-gluon plasma

The relativistic heavy ion program developed at RHIC and now at LHC motivated a deeper study of the properties of the quark-gluon plasma (QGP) and, in particular, the study of perturbations in this kind of plasma. We are interested on the time evolution of perturbations in the baryon and energy densities. If a localized pulse in baryon density could propagate throughout the QGP for long distances preserving its shape and without loosing localization, this could have interesting consequences for relativistic heavy ion physics and for astrophysics. A mathematical way to prove that this can happen is to derive (under certain conditions) from the hydrodynamical equations of the QGP a Korteveg-de Vries (KdV) equation. The solution of this equation describes the propagation of a KdV soliton. The derivation of the KdV equation depends crucially on the equation of state (EOS) of the QGP. The use of the simple MIT bag model EOS does not lead to KdV solitons. Recently we have developed an EOS for the QGP which includes both perturbative and nonperturbative corrections to the MIT one and is still simple enough to allow for analytical manipulations. With this EOS we were able to derive a KdV equation for the cold QGP.

Identifying Multiquark Hadrons from Heavy Ion Collisions

Identifying hadronic molecular states and/or hadrons with multiquark components either with or without exotic quantum numbers is a long-standing challenge in hadronic physics. We suggest that studying the production of these hadrons in relativistic heavy ion collisions offers a promising resolution to this problem as yields of exotic hadrons are expected to be strongly affected by their structures. Using the coalescence model for hadron production, we find that, compared to the case of a nonexotic hadron with normal quark numbers, the yield of an exotic hadron is typically an order of magnitude smaller when it is a compact multiquark state and a factor of 2 or more larger when it is a loosely bound hadronic molecule. We further find that some of the newly proposed heavy exotic states could be produced and realistically measured in these experiments.

Higher-derivative supersymmetric gauge theory

We study the one-loop low-energy effective action for the higher-derivative superfield gauge theory coupled to chiral matter.

Effective actions at finite temperature

This is a more detailed version of our recent paper where we proposed, from first principles, a direct method for evaluating the exact fermion propagator in the presence of a general background field at finite temperature. This can, in turn, be used to determine the finite temperature effective action for the system. As applications, we discuss the complete one loop finite temperature effective actions for 0+1 dimensional QED as well as for the Schwinger model in detail. These effective actions, which are derived in the real time (closed time path) formalism, generate systematically all the Feynman amplitudes calculated in thermal perturbation theory and also show that the retarded (advanced) amplitudes vanish in these theories. Various other aspects of the problem are also discussed in detail.

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.

The gross theory model for neutrino-nucleus cross-section

The nuclear gross theory, originally formulated by Takahashi and Yamada (1969 Prog. Theor. Phys. 41 1470) for the beta-decay, is applied to the electronic-neutrino nucleus reactions, employing a more realistic description of the energetics of the Gamow-Teller resonances. The model parameters are gauged from the most recent experimental data, both for beta(-)-decay and electron capture, separately for even-even, even-odd, odd-odd and odd-even nuclei. The numerical estimates for neutrino-nucleus cross-sections agree fairly well with previous evaluations done within the framework of microscopic models. The formalism presented here can be extended to the heavy nuclei mass region, where weak processes are quite relevant, which is of astrophysical interest because of its applications in supernova explosive nucleosynthesis.

Approximation to density functional theory for the calculation of band gaps of semiconductors

The local-density approximation (LDA) together with the half occupation (transitionstate) is notoriously successful in the calculation of atomic ionization potentials. When it comes to extended systems, such as a semiconductor infinite system, it has been very difficult to find a way to half ionize because the hole tends to be infinitely extended (a Bloch wave). The answer to this problem lies in the LDA formalism itself. One proves that the half occupation is equivalent to introducing the hole self-energy (electrostatic and exchange correlation) into the Schrodinger equation. The argument then becomes simple: The eigenvalue minus the self-energy has to be minimized because the atom has a minimal energy. Then one simply proves that the hole is localized, not infinitely extended, because it must have maximal self-energy. Then one also arrives at an equation similar to the self- interaction correction equation, but corrected for the removal of just 1/2 electron. Applied to the calculation of band gaps and effective masses, we use the self- energy calculated in atoms and attain a precision similar to that of GW, but with the great advantage that it requires no more computational effort than standard LDA.

Measurement of the Bottom Quark Contribution to Nonphotonic Electron Production in p plus p Collisions at root s=200 GeV

The contribution of B meson decays to nonphotonic electrons, which are mainly produced by the semileptonic decays of heavy-flavor mesons, in p + p collisions at root s = 200 GeV has been measured using azimuthal correlations between nonphotonic electrons and hadrons. The extracted B decay contribution is approximately 50% at a transverse momentum of p(T) >= 5 GeV/c. These measurements constrain the nuclear modification factor for electrons from B and D meson decays. The result indicates that B meson production in heavy ion collisions is also suppressed at high p(T).

Parton energy loss in heavy-ion collisions via direct-photon and charged-particle azimuthal correlations

Charged-particle spectra associated with direct photon (gamma(dir)) and pi(0) are measured in p + p and Au + Au collisions at center-of-mass energy root(S)(NN) = 200 GeV with the STAR detector at the Relativistic Heavy Ion Collider. A shower-shape analysis is used to partially discriminate between gamma(dir) and pi(0). Assuming no associated charged particles in the gamma(dir) direction ( near side) and small contribution from fragmentation photons (gamma(frag)), the associated charged-particle yields opposite to gamma(dir) (away side) are extracted. In central Au + Au collisions, the charged-particle yields at midrapidity (vertical bar eta vertical bar < 1) and high transverse momentum (3 < (assoc)(PT) < 16 GeV/c) associated with gamma(dir) and pi(0) (vertical bar eta vertical bar < 0.9, 8 < (trig)(PT) < 16 GeV/c) are suppressed by a factor of 3-5 compared with p + p collisions. The observed suppression of the associated charged particles is similar for gamma(dir) and pi(0) and independent of the gamma(dir) energy within uncertainties. These measurements indicate that, in the kinematic range covered and within our current experimental uncertainties, the parton energy loss shows no sensitivity to the parton initial energy, path length, or color charge.

Observation of charge-dependent azimuthal correlations and possible local strong parity violation in heavy-ion collisions

Parity (P)-odd domains, corresponding to nontrivial topological solutions of the QCD vacuum, might be created during relativistic heavy-ion collisions. These domains are predicted to lead to charge separation of quarks along the orbital momentum of the system created in noncentral collisions. To study this effect, we investigate a three-particle mixed-harmonics azimuthal correlator which is a P-even observable, but directly sensitive to the charge-separation effect. We report measurements of this observable using the STAR detector in Au + Au and Cu + Cu collisions at root s(NN) = 200 and 62 GeV. The results are presented as a function of collision centrality, particle separation in rapidity, and particle transverse momentum. A signal consistent with several of the theoretical expectations is detected in all four data sets. We compare our results to the predictions of existing event generators and discuss in detail possible contributions from other effects that are not related to P violation.

(In)stability of D-dimensional black holes in Gauss-Bonnet theory

We make an extensive study of evolution of gravitational perturbations of D-dimensional black holes in Gauss-Bonnet theory. There is an instability at higher multipoles l and large Gauss-Bonnet coupling alpha for D = 5, 6, which is stabilized at higher D. Although a small negative gap of the effective potential for the scalar type of gravitational perturbations exists for higher D and whatever alpha, it does not lead to any instability.

Nonlinear waves in a quark gluon plasma

We study the propagation of perturbations in the quark gluon plasma. This subject has been addressed in other works and in most of the theoretical descriptions of this phenomenon the hydrodynamic equations have been linearized for simplicity. We propose an alternative approach, also based on hydrodynamics but taking into account the nonlinear terms of the equations. We show that these terms may lead to localized waves or even solitons. We use a simple equation of state for the QGP and expand the hydrodynamic equations around equilibrium configurations. The resulting differential equations describe the propagation of perturbations in the energy density. We solve them numerically and find that localized perturbations can propagate for long distances in the plasma. Under certain conditions our solutions mimic the propagation of Korteweg-de Vries solitons.