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.

Liquid polymorphism, order-disorder transitions and anomalous behavior: A Monte Carlo study of the Bell-Lavis model for water

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution. (C) 2009 American Institute of Physics. [doi:10.1063/1.3253297]

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.

Suppression Pattern of Neutral Pions at High Transverse Momentum in Au plus Au Collisions at root S(NN)=200 GeV and Constraints on Medium Transport Coefficients

For Au + Au collisions at 200 GeV, we measure neutral pion production with good statistics for transverse momentum, p(T), up to 20 GeV/c. A fivefold suppression is found, which is essentially constant for 5 < p(T) < 20 GeV/c. Experimental uncertainties are small enough to constrain any model-dependent parametrization for the transport coefficient of the medium, e. g., <(q) over cap > in the parton quenching model. The spectral shape is similar for all collision classes, and the suppression does not saturate in Au + Au collisions.

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.

Quantitative constraints on the transport properties of hot partonic matter from semi-inclusive single high transverse momentum pion suppression in Au plus Au collisions at root S(NN)=200 GeV

The PHENIX experiment has measured the suppression of semi-inclusive single high-transverse-momentum pi(0)'s in Au+Au collisions at root s(NN) = 200 GeV. The present understanding of this suppression is in terms of energy loss of the parent (fragmenting) parton in a dense color-charge medium. We have performed a quantitative comparison between various parton energy-loss models and our experimental data. The statistical point-to-point uncorrelated as well as correlated systematic uncertainties are taken into account in the comparison. We detail this methodology and the resulting constraint on the model parameters, such as the initial color-charge density dN(g)/dy, the medium transport coefficient <(q) over cap >, or the initial energy-loss parameter epsilon(0). We find that high-transverse-momentum pi(0) suppression in Au+Au collisions has sufficient precision to constrain these model-dependent parameters at the +/- 20-25% (one standard deviation) level. These constraints include only the experimental uncertainties, and further studies are needed to compute the corresponding theoretical uncertainties.

Frustration and anomalous behavior in the Bell-Lavis model of liquid water

We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively). Both models present liquid-liquid coexistence and several thermodynamic anomalies. This result suggests that anisotropy introduced through orientational variables play no specific role in producing the density anomaly, in agreement with a similar conclusion discussed previously following results for continuous soft core,models. We propose that the presence of liquid anomalies may be related to energetic frustration, a feature common to both models.

Thermal Transport and Strong Mass Renormalization in NiCl(2)-4SC(NH(2))(2)

Several quantum paramagnets exhibit magnetic-field-induced quantum phase transitions to an anti-ferromagnetic state that exists for H(c1) <= H <= H(c2). For some of these compounds, there is a significant asymmetry between the low-and high-field transitions. We present specific heat and thermal conductivity measurements in NiCl(2)-4SC(NH(2))(2), together with calculations which show that the asymmetry is caused by a strong mass renormalization due to quantum fluctuations for H <= H(c1) that are absent for H >= H(c2). We argue that the enigmatic lack of asymmetry in thermal conductivity is due to a concomitant renormalization of the impurity scattering.

Transport in disordered two-dimensional topological insulators

The transport properties of the ""inverted"" semiconductor HgTe-based quantum well, recently shown to be a two-dimensional topological insulator, are studied experimentally in the diffusive regime. Nonlocal transport measurements are performed in the absence of magnetic field, and a large signal due to the edge states is observed. This shows that the edge states can propagate over a long distance, similar to 1 mm, and therefore, there is no difference between local and nonlocal electrical measurements in a topological insulator. In the presence of an in-plane magnetic field a strong decrease of the local resistance and complete suppression of the nonlocal resistance is observed. We attribute this behavior to an in-plane magnetic-field-induced transition from the topological insulator state to a conventional bulk metal state.

Nonlinear transport and oscillating magnetoresistance in double quantum wells

The nonlinear regime of low-temperature magnetoresistance of double quantum wells in the region of magnetic fields below 1 T is studied both experimentally and theoretically. The observed inversion of the magnetointersubband oscillation peaks with increasing electric current and splitting of these peaks are described by the theory based on the kinetic equation for the isotropic nonequilibrium part of electron distribution function. The inelastic-scattering time of electrons is determined from the current dependence of the inversion field.

Electronic, structural, and transport properties of Ni-doped graphene nanoribbons

We have investigated the electronic and transport properties of zigzag Ni-adsorbed graphene nanoribbons (Ni/GNRs) using ab initio calculations. We find that the Ni adatoms lying along the edge of zigzag GNRs represent the energetically most stable configuration, with an energy difference of approximately 0.3 eV when compared to the adsorption in the middle of the ribbon. The carbon atoms at the ribbon edges still present nonzero magnetic moments as in the pristine GNR even though there is a quenching by a factor of almost five in the value of the local magnetic moments at the C atoms bonded to the Ni. This quenching decays relatively fast and at approximately 9 A from the Ni adsorption site the magnetic moments have already values close to the pristine ribbon. At the opposite edge and at the central carbon atoms the changes in the magnetic moments are negligible. The energetic preference for the antiparallel alignment between the magnetization at the opposite edges of the ribbon is still maintained upon Ni adsorption. We find many Ni d-related states within an energy window of 1 eV above and below the Fermi energy, which gives rise to a spin-dependent charge transport. These results suggest the possibility of manufacturing spin devices based on GNRs doped with Ni atoms.

Linear and nonlinear transport in a small charge-tunable open quantum ring

We experimentally study the Aharonov-Bohm-conductance oscillations under external gate voltage in a semiconductor quantum ring with a radius of 80 nm. We find that, in the linear regime, the resistance-oscillation plot in the voltage-magnetic-field plane corresponds to the quantum ring energy spectra. The chessboard pattern assembled by resistance diamonds, while loading the ring, is attributed to a short electron lifetime in the open configuration, which agrees with calculations within the single-particle model. Remarkably, the application of a small dc current allows observing strong deviations in the oscillation plot from this pattern accompanied by a magnetic-field symmetry break. We relate such behavior to the higher-order-conductance coefficients determined by electron-electron interactions in the nonlinear regime.

Transport properties of single vacancies in nanotubes

We present density of states and electronic transport calculations of single vacancies in carbon nanotubes. We confirm that the defect reconstructs into a pentagon and a nonagon, following the removal of a single carbon atom. This leads to the formation of a dangling bond. Finally, we demonstrate that care must be taken when calculating the density of states of impurities in one-dimensional systems in general. Traditional treatments of these systems using periodic boundary conditions leads to the formation of minigaps even in the limit of large unit cells.