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.

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.

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.

Charged hadron multiplicity fluctuations in Au+Au and Cu+Cu collisions from root S(NN) = 22.5 to 200 GeV

A comprehensive survey of event-by-event fluctuations of charged hadron multiplicity in relativistic heavy ions is presented. The survey covers Au+Au collisions at s(NN)=62.4 and 200 GeV, and Cu+Cu collisions at s(NN)=22.5,62.4, and 200 GeV. Fluctuations are measured as a function of collision centrality, transverse momentum range, and charge sign. After correcting for nondynamical fluctuations due to fluctuations in the collision geometry within a centrality bin, the remaining dynamical fluctuations expressed as the variance normalized by the mean tend to decrease with increasing centrality. The dynamical fluctuations are consistent with or below the expectation from a superposition of participant nucleon-nucleon collisions based upon p+p data, indicating that this dataset does not exhibit evidence of critical behavior in terms of the compressibility of the system. A comparison of the data with a model where hadrons are independently emitted from a number of hadron clusters suggests that the mean number of hadrons per cluster is small in heavy ion collisions.

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.

1/N expansion in noncommutative quantum mechanics

We study the 1/N expansion in noncommutative quantum mechanics for the anharmonic and Coulombian potentials. The expansion for the anharmonic oscillator presented good convergence properties, but for the Coulombian potential, we found a divergent large N expansion when using the usual noncommutative generalization of the potential. We proposed a modified version of the noncommutative Coulombian potential which provides a well-behaved 1/N expansion.

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.

On the mixing property for a class of states of relativistic quantum fields

Let omega be a factor state on the quasilocal algebra A of observables generated by a relativistic quantum field, which, in addition, satisfies certain regularity conditions [satisfied by ground states and the recently constructed thermal states of the P(phi)(2) theory]. We prove that there exist space- and time-translation invariant states, some of which are arbitrarily close to omega in the weak * topology, for which the time evolution is weakly asymptotically Abelian. (C) 2010 American Institute of Physics. [doi: 10.1063/1.3372623]

Noncommutativity due to spin

Using the Berezin-Marinov pseudoclassical formulation of the spin particle we propose a classical model of spin noncommutativity. In the nonrelativistic case, the Poisson brackets between the coordinates are proportional to the spin angular momentum. The quantization of the model leads to the noncommutativity with mixed spatial and spin degrees of freedom. A modified Pauli equation, describing a spin half particle in an external electromagnetic field is obtained. We show that nonlocality caused by the spin noncommutativity depends on the spin of the particle; for spin zero, nonlocality does not appear, for spin half, Delta x Delta y >= theta(2)/2, etc. In the relativistic case the noncommutative Dirac equation was derived. For that we introduce a new star product. The advantage of our model is that in spite of the presence of noncommutativity and nonlocality, it is Lorentz invariant. Also, in the quasiclassical approximation it gives noncommutativity with a nilpotent parameter.

Hard thermal loops in static external fields

We examine, in the imaginary-time formalism, the high temperature behavior of n-point thermal loops in static Yang-Mills and gravitational fields. We show that in this regime, any hard thermal loop gives the same leading contribution as the one obtained by evaluating the loop integral at zero external energies and momenta.

Quantum reflection: The invisible quantum barrier

We construct an invisible quantum barrier which represents the phenomenon of quantum reflection using available data on atom-wall and Bose-Einstein-condensate-wall reflection. We use the Abel equation to invert the data. The resulting invisible quantum barrier is double valued in both axes. We study this invisible barrier in the case of atom and Bose-Einstein condensate (BEC) reflection from a solid silicon surface. A time-dependent, one-spatial-dimension Gross-Pitaevskii equation is solved for the BEC case. We found that the BEC behaves very similarly to the single atom except for size effects, which manifest themselves in a maximum in the reflectivity at small distances from the wall. The effect of the atom-atom interaction on the BEC reflection and correspondingly on the invisible barrier is found to be appreciable at low velocities and comparable to the finite-size effect. The trapping of an ultracold atom or BEC between two walls is discussed.

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.

Magnetoelectric effects in an organometallic quantum magnet

Metal-organic materials constitute a new field in which to search for ferroelectricity and coupling between electricity and magnetism. We observe a magnetic field-induced change in the electric polarization, Delta P(H), that reaches 50 mu C/m(2) in single crystals of NiCl(2)-4SC(NH(2))(2) (DTN). DTN forms a tetragonal structure that breaks inversion symmetry with the electrically polar thiourea molecules [SC(NH(2))] all tilted in the same direction along the c axis. The field H induces canted antiferromagnetism of the Ni S = 1 spins between 2 and 12 T and our measurements show that the electric polarization increases monotonically in this range, saturating above 12 T. By modeling the microscopic origin of this magnetoelectric effect, we find that the leading contribution to Delta P comes from the change in the crystal electric field, with a smaller contribution from magnetic exchange striction. The finite value of Delta P induced by magnetostriction results from the polar nature of the thiourea molecules bonded to the Ni atoms, and it is amplified by the softness of these organic molecules.

Hyperfine interactions in silicon quantum dots

A fundamental interaction for electrons is their hyperfine interaction (HFI) with nuclear spins. HFI is well characterized in free atoms and molecules, and is crucial for purposes from chemical identification of atoms to trapped ion quantum computing. However, electron wave functions near atomic sites, therefore HFI, are often not accurately known in solids. Here we perform an all-electron calculation for conduction electrons in silicon and obtain reliable information on HFI. We verify the outstanding quantum spin coherence in Si, which is critical for fault-tolerant solid state quantum computing.

Thermally activated intersubband scattering and oscillating magnetoresistance in quantum wells

Experimental studies of magnetoresistance in high-mobility wide quantum wells reveal oscillations which appear with an increase in temperature to 10 K and whose period is close to that of Shubnikov-de Haas oscillations. The observed phenomenon is identified as magnetointersubband oscillations caused by the scattering of electrons between two occupied subbands and the third subband which becomes occupied as a result of thermal activation. These small-period oscillations are less sensitive to thermal suppression than the large-period magnetointersubband oscillations caused by the scattering between the first and the second subbands. Theoretical study, based on consideration of electron scattering near the edge of the third subband, gives a reasonable explanation of our experimental findings.