Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.

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.

Transverse momentum dependence of eta meson suppression in Au plus Au collisions at root s(NN)=200 GeV

New measurements by the PHENIX experiment at the Relativistic Heavy Ion Collider for. production at midrapidity as a function of transverse momentum ((PT)) and collision centrality in root s(NN) = 200 GeV Au + Au and p + p collisions are presented. They indicate nuclear modification factors (R(AA)) which are similar in both magnitude and trend to those found in earlier pi(0) measurements. Linear fits to R(AA) as a function of (PT) in 5-20 GeV/c show that the slope is consistent with zero within two standard deviations at all centralities, although a slow rise cannot be excluded. Having different statistical and systematic uncertainties, the pi(0) and eta measurements are complementary at high (PT); thus, along with the extended (PT) range of these data they can provide additional constraints for theoretical modeling and the extraction of transport properties.

Transverse momentum dependence of J/psi polarization at midrapidity in p plus p collisions at root s=200 GeV

We report the measurement of the transverse momentum dependence of inclusive J/psi polarization in p + p collisions at root s = 200 GeV performed by the PHENIX Experiment at the Relativistic Heavy Ion Collider. The J/psi polarization is studied in the helicity, Gottfried-Jackson, and Collins-Soper frames for p(T) < 5 GeV/c and vertical bar y vertical bar < 0.35. The polarization in the helicity and Gottfried-Jackson frames is consistent with zero for all transverse momenta, with a slight (1.8 sigma) trend towards longitudinal polarization for transverse momenta above 2 GeV/c. No conclusion is allowed due to the limited acceptance in the Collins-Soper frame and the uncertainties of the current data. The results are compared to observations for other collision systems and center of mass energies and to different quarkonia production models.

Double-helicity dependence of jet properties from dihadrons in longitudinally polarized p plus p collisions at root s=200 GeV

It has been postulated that partonic orbital angular momentum can lead to a significant double-helicity dependence in the net transverse momentum of Drell-Yan dileptons produced in longitudinally polarized p + p collisions. Analogous effects are also expected for dijet production. If confirmed by experiment, this hypothesis, which is based on semiclassical arguments, could lead to a new approach for studying the contributions of orbital angular momentum to the proton spin. We report the first measurement of the double-helicity dependence of the dijet transverse momentum in longitudinally polarized p + p collisions at root s = 200 GeV from data taken by the PHENIX experiment in 2005 and 2006. The analysis deduces the transverse momentum of the dijet from the widths of the near-and far-side peaks in the azimuthal correlation of the dihadrons. When averaged over the transverse momentum of the triggered particle, the difference of the root mean square of the dijet transverse momentum between like-and unlike-helicity collisions is found to be -37 +/- 88(stat) +/- 14(sys)t MeV/c.

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.

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.

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

Mn-doped cubic BN as an atomiclike memory device: A density functional study

We investigate the electronic properties of Mn(B) substitutional doping in cubic boron nitride (BN), for different charge states, using density functional theory (DFT) calculations. We show that the neutral Mn has a nonmagnetic ground state (S=0). Upon charge injection, it is unambiguously shown that the Mn(B)(-) has a high-spin configuration with a strong, localized magnetic moment of 5 mu(Bohr). We developed a simple model, parameterized by the DFT results, that allows us to interpret the rules played by the crystal-field and exchange-correlation splitting in the magnetization process.

High density flux of Co nanoparticles produced by a simple gas aggregation apparatus

Gas aggregation is a well known method used to produce clusters of different materials with good size control, reduced dispersion, and precise stoichiometry. The cost of these systems is relatively high and they are generally dedicated apparatuses. Furthermore, the usual sample production speed of these systems is not as fast as physical vapor deposition devices posing a problem when thick samples are needed. In this paper we describe the development of a multipurpose gas aggregation system constructed as an adaptation to a magnetron sputtering system. The cost of this adaptation is negligible and its installation and operation are both remarkably simple. The gas flow for flux in the range of 60-130 SCCM (SCCM denotes cubic centimeter per minute at STP) is able to completely collimate all the sputtered material, producing spherical nanoparticles. Co nanoparticles were produced and characterized using electron microscopy techniques and Rutherford back-scattering analysis. The size of the particles is around 10 nm with around 75 nm/min of deposition rate at the center of a Gaussian profile nanoparticle beam.

Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids

One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory [Perdew-Burke-Ernzerhof (PBE) GGA] and a recently proposed modification designed specifically for solids (PBEsol) are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules, and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints stemming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.

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.

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.

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.

Charged and strange hadron elliptic flow in Cu plus Cu collisions at root s(NN)=62.4 and 200 GeV

We present the results of an elliptic flow, v(2), analysis of Cu + Cu collisions recorded with the solenoidal tracker detector (STAR) at the BNL Relativistic Heavy Ion Collider at root s(NN) = 62.4 and 200 GeV. Elliptic flow as a function of transverse momentum, v(2)(p(T)), is reported for different collision centralities for charged hadrons h(+/-) and strangeness-ontaining hadrons K(S)(0), Lambda, Xi, and phi in the midrapidity region vertical bar eta vertical bar < 1.0. Significant reduction in systematic uncertainty of the measurement due to nonflow effects has been achieved by correlating particles at midrapidity, vertical bar eta vertical bar < 1.0, with those at forward rapidity, 2.5 < vertical bar eta vertical bar < 4.0. We also present azimuthal correlations in p + p collisions at root s = 200 GeV to help in estimating nonflow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au + Au collisions at root s(NN) = 200 GeV. We observe that v(2)(p(T)) of strange hadrons has similar scaling properties as were first observed in Au + Au collisions, that is, (i) at low transverse momenta, p(T) < 2 GeV/c, v(2) scales with transverse kinetic energy, m(T) - m, and (ii) at intermediate p(T), 2 < p(T) < 4 GeV/c, it scales with the number of constituent quarks, n(q.) We have found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of v(2)(p(T)) for K(S)(0) and Lambda. Eccentricity scaled v(2) values, v(2)/epsilon, are larger in more central collisions, suggesting stronger collective flow develops in more central collisions. The comparison with Au + Au collisions, which go further in density, shows that v(2)/epsilon depends on the system size, that is, the number of participants N(part). This indicates that the ideal hydrodynamic limit is not reached in Cu + Cu collisions, presumably because the assumption of thermalization is not attained.