eta ->gamma gamma decay width via the Primakoff cross section

Incoherent eta photoproduction in nuclei is evaluated at forward angles within 4 to 9 GeV using a multiple scattering Monte Carlo cascade calculation with full eta-nucleus final-state interactions. The Primakoff, nuclear coherent and nuclear incoherent components of the cross sections fit remarkably well previous measurements for Be and Cu from Cornell, suggesting a destructive interference between the Coulomb and nuclear coherent amplitudes for Cu. The inelastic background of the data is consistently attributed to the nuclear incoherent part, which is clearly not isotropic as previously considered in Cornell's analysis. The respective Primakoff cross sections from Be and Cu give Gamma(eta ->gamma gamma)=0.476(62) keV, where the quoted error is only statistical. This result is consistent with the Particle Data Group average of 0.510(26) keV and in sharp contrast (similar to 50%) with the value of 0.324(46) keV obtained at Cornell.

Diffusion anomaly and dynamic transitions in the Bell-Lavis water model

In this paper we investigate the dynamic properties of the minimal Bell-Lavis (BL) water model and their relation to the thermodynamic anomalies. The BL model is defined on a triangular lattice in which water molecules are represented by particles with three symmetric bonding arms interacting through van der Waals and hydrogen bonds. We have studied the model diffusivity in different regions of the phase diagram through Monte Carlo simulations. Our results show that the model displays a region of anomalous diffusion which lies inside the region of anomalous density, englobed by the line of temperatures of maximum density. Further, we have found that the diffusivity undergoes a dynamic transition which may be classified as fragile-to-strong transition at the critical line only at low pressures. At higher densities, no dynamic transition is seen on crossing the critical line. Thus evidence from this study is that relation of dynamic transitions to criticality may be discarded. (C) 2010 American Institute of Physics. [doi:10.1063/1.3479001]

Delocalized water and fluoride contributions to Dyson orbitals for electron detachment from the hydrated fluoride anion

The experimental vertical electron detachment energy (VEDE) of aqueous fluoride, [F(-)(H(2)O)], is approximately 9.8 eV, but spectral assignment is complicated by interference between F(-) 2p and H(2)O 1b(1) orbitals. The electronic structure of [F(-)(H(2)O)] is analyzed with Monte Carlo and ab initio quantum-mechanical calculations. Electron-propagator calculations in the partial third-order approximation yield a VEDE of 9.4 eV. None of the Dyson orbitals corresponding to valence VEDEs consists primarily of F 2p functions. These results and ground-state atomic charges indicate that the final, neutral state is more appropriately described as [F(-)(H(2)O)(+)] than as [F(H(2)O)]. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3431081]

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.

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

Electron collisions with the CH(2)O-H(2)O complex

We report cross sections for elastic collisions of low-energy electrons with the CH(2)O-H(2)O complex. We employed the Schwinger multichannel method with pseudopotentials in the static-exchange and in the static-exchange-polarization approximations for energies from 0.1 to 20 eV. We considered four different hydrogen-bonded structures for the complex that were generated by classical Monte Carlo simulations. Our aim is to investigate the effect of the water molecule on the pi* shape resonance of formaldehyde. Previous studies reported a pi* shape resonance for CH(2)O at around 1 eV. The resonance positions of the complexes appear at lower energies in all cases due to the mutual polarization between the two molecules. This indicates that the presence of water may favor dissociation by electron impact and may lead to an important effect on strand breaking in wet DNA by electron impact.

Dynamic transitions in a two dimensional associating lattice gas model

Using Monte Carlo simulations we investigate some new aspects of the phase diagram and the behavior of the diffusion coefficient in an associating lattice gas (ALG) model on different regions of the phase diagram. The ALG model combines a two dimensional lattice gas where particles interact through a soft core potential and orientational degrees of freedom. The competition between soft core potential and directional attractive forces results in a high density liquid phase, a low density liquid phase, and a gas phase. Besides anomalies in the behavior of the density with the temperature at constant pressure and of the diffusion coefficient with density at constant temperature are also found. The two liquid phases are separated by a coexistence line that ends in a bicritical point. The low density liquid phase is separated from the gas phase by a coexistence line that ends in tricritical point. The bicritical and tricritical points are linked by a critical lambda-line. The high density liquid phase and the fluid phases are separated by a second critical tau-line. We then investigate how the diffusion coefficient behaves on different regions of the chemical potential-temperature phase diagram. We find that diffusivity undergoes two types of dynamic transitions: a fragile-to-strong transition when the critical lambda-line is crossed by decreasing the temperature at a constant chemical potential; and a strong-to-strong transition when the critical tau-line is crossed by decreasing the temperature at a constant chemical potential.

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]

Critical Properties at the Field-Induced Bose-Einstein Condensation in NiCl(2)-4SC(NH(2))(2)

We report new magnetization measurements on the spin-gap compound NiCl(2)-4SC(NH(2))(2) at the low-field boundary of the magnetic field-induced ordering. The critical density of the magnetization is analyzed in terms of a Bose-Einstein condensation of bosonic quasiparticles. The analysis of the magnetization at the transition leads to the conclusion for the preservation of the U(1) symmetry, as required for Bose-Einstein condensation. The experimental data are well described by quantum Monte Carlo simulations.

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.

Growth of Long Range Forward-Backward Multiplicity Correlations with Centrality in Au plus Au Collisions at root s(NN)=200 GeV

Forward-backward multiplicity correlation strengths have been measured with the STAR detector for Au + Au and p + p collisions at root s(NN) = 200 GeV. Strong short- and long-range correlations (LRC) are seen in central Au + Au collisions. The magnitude of these correlations decrease with decreasing centrality until only short-range correlations are observed in peripheral Au + Au collisions. Both the dual parton model (DPM) and the color glass condensate (CGC) predict the existence of the long-range correlations. In the DPM, the fluctuation in the number of elementary (parton) inelastic collisions produces the LRC. In the CGC, longitudinal color flux tubes generate the LRC. The data are in qualitative agreement with the predictions of the DPM and indicate the presence of multiple parton interactions.

Numerical evidence against a conjecture on the cover time of planar graphs

We investigate a conjecture on the cover times of planar graphs by means of large Monte Carlo simulations. The conjecture states that the cover time tau (G(N)) of a planar graph G(N) of N vertices and maximal degree d is lower bounded by tau (G(N)) >= C(d)N(lnN)(2) with C(d) = (d/4 pi) tan(pi/d), with equality holding for some geometries. We tested this conjecture on the regular honeycomb (d = 3), regular square (d = 4), regular elongated triangular (d = 5), and regular triangular (d = 6) lattices, as well as on the nonregular Union Jack lattice (d(min) = 4, d(max) = 8). Indeed, the Monte Carlo data suggest that the rigorous lower bound may hold as an equality for most of these lattices, with an interesting issue in the case of the Union Jack lattice. The data for the honeycomb lattice, however, violate the bound with the conjectured constant. The empirical probability distribution function of the cover time for the square lattice is also briefly presented, since very little is known about cover time probability distribution functions in general.

Effect of long cyclic exchanges on the magnetic properties of bcc (3)He

Using path-integral Monte Carlo calculations, we have calculated ring exchange frequencies in the bcc phase of solid (3)He for densities from melting to the highest stable density. We evaluate 42 different exchange frequencies from two atoms up to eight atoms and find their Gruneisen exponents. Using a fit to these frequencies, we calculate the contribution to the Curie-Weiss temperature, Theta(CW), and upper critical magnetic field, B(c2), for even longer exchanges using a lattice Monte Carlo procedure. We find that contributions from seven-and eight-particle exchanges make a significant contribution to Theta(CW) and B(c2) at melting density. Comparison with experimental data is given.

Directed Abelian algebras and their application to stochastic models

With each directed acyclic graph (this includes some D-dimensional lattices) one can associate some Abelian algebras that we call directed Abelian algebras (DAAs). On each site of the graph one attaches a generator of the algebra. These algebras depend on several parameters and are semisimple. Using any DAA, one can define a family of Hamiltonians which give the continuous time evolution of a stochastic process. The calculation of the spectra and ground-state wave functions (stationary state probability distributions) is an easy algebraic exercise. If one considers D-dimensional lattices and chooses Hamiltonians linear in the generators, in finite-size scaling the Hamiltonian spectrum is gapless with a critical dynamic exponent z=D. One possible application of the DAA is to sandpile models. In the paper we present this application, considering one- and two-dimensional lattices. In the one-dimensional case, when the DAA conserves the number of particles, the avalanches belong to the random walker universality class (critical exponent sigma(tau)=3/2). We study the local density of particles inside large avalanches, showing a depletion of particles at the source of the avalanche and an enrichment at its end. In two dimensions we did extensive Monte-Carlo simulations and found sigma(tau)=1.780 +/- 0.005.