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.

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.

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.

Accurate prediction of the Si/SiO(2) interface band offset using the self-consistent ab initio DFT/LDA-1/2 method

We use the density functional theory/local-density approximation (DFT/LDA)-1/2 method [L. G. Ferreira , Phys. Rev. B 78, 125116 (2008)], which attempts to fix the electron self-energy deficiency of DFT/LDA by half-ionizing the whole Bloch band of the crystal, to calculate the band offsets of two Si/SiO(2) interface models. Our results are similar to those obtained with a ""state-of-the-art"" GW approach [R. Shaltaf , Phys. Rev. Lett. 100, 186401 (2008)], with the advantage of being as computationally inexpensive as the usual DFT/LDA. Our band gap and band offset predictions are in excellent agreement with experiments.

Atomistic molecular dynamics study of interface formation: Al on poly(p-phenylene vinylene)

We model interface formation by metal deposition on the conjugated polymer poly-para-phenylene vinylene, studying direct aluminum and layered aluminum-calcium structures Al/PPV and Al/Ca/PPV. To do that we use classical molecular dynamics simulations, checked by ab initio density-functional theory calculations, for selected relevant configurations. We find that Al not only migrates easily into the film, with a strong charge transfer to the neighboring chains, but also promotes rearrangement of the polymer in the interfacial region to the hexagonal structure. On the other hand, our results indicate that a thin Ca layer is sufficient to protect the film and maintain a well-defined metal/polymer interface, and that also a thin Al capping layer may protect the whole from environmental degradation.

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.

Soliton dynamics at an interface between a uniform medium and a nonlinear optical lattice

We study trapping and propagation of a matter-wave soliton through the interface between uniform medium and a nonlinear optical lattice. Different regimes for transmission of a broad and a narrow solitons are investigated. Reflections and transmissions of solitons are predicted as a function of the lattice phase. The existence of a threshold in the amplitude of the nonlinear optical lattice, separating the transmission and reflection regimes, is verified. The localized nonlinear surface state, corresponding to the soliton trapped by the interface, is found. Variational approach predictions are confirmed by numerical simulations for the original Gross-Pitaevskii equation with nonlinear periodic potentials.

Substrate/semiconductor interface effects on the emission efficiency of luminescent polymers

The importance of interface effects for organic devices has long been recognized, but getting detailed knowledge of the extent of such effects remains a major challenge because of the difficulty in distinguishing from bulk effects. This paper addresses the interface effects on the emission efficiency of poly(p-phenylene vinylene) (PPV), by producing layer-by-layer (LBL) films of PPV alternated with dodecylbenzenesulfonate. Films with thickness varying from similar to 15 to 225 nm had the structural defects controlled empirically by converting the films at two temperatures, 110 and 230 degrees C, while the optical properties were characterized by using optical absorption, photoluminescence (PL), and photoluminescence excitation spectra. Blueshifts in the absorption and PL spectra for LBL films with less than 25 bilayers (<40-50 nm) pointed to a larger number of PPV segments with low conjugation degree, regardless of the conversion temperature. For these thin films, the mean free-path for diffusion of photoexcited carriers decreased, and energy transfer may have been hampered owing to the low mobility of the excited carriers. The emission efficiency was then found to depend on the concentration of structural defects, i.e., on the conversion temperature. For thick films with more than 25 bilayers, on the other hand, the PL signal did not depend on the PPV conversion temperature. We also checked that the interface effects were not caused by waveguiding properties of the excited light. Overall, the electronic states at the interface were more localized, and this applied to film thickness of up to 40-50 nm. Because this is a typical film thickness in devices, the implication from the findings here is that interface phenomena should be a primary concern for the design of any organic device. (C) 2011 American Institute of Physics. [doi:10.1063/1.3622143]

Interface roughness in short-period InGaAs/InP superlattices

Electron mobility was studied in lattice-matched short-period InGaAs/InP superlattices as a function of the width of the wells. The decreasing mobility with decreasing well width was shown to occur due to the interface roughness. The roughnesses of InGaAs/InP and GaAs/AlGaAs interfaces were compared. Much smoother InGaAs/InP interfaces resulted in higher electron mobility limited by interface roughness.

Ultrasensitive thermal lens spectroscopy of water

A recently developed dual-beam configuration that optimizes the thermal lens technique has been used to obtain the absorption spectrum of pure water from 350 to 528 nm. Our results indicate the minimum linear absorption coefficient smaller than 2 X 10(-5) cm(-1) between 360 and 400 nm. This value is lower than previous literature data, and it is blueshifted. Absorption coefficients as small as 2 X 10(-7) cm(-1) can be measured for water using 1 W of excitation power. A detection limit of similar to 6 X 10(-9) cm(-1) (P=1 W) for CCl(4) was estimated, which represents, to the best of our knowledge, the highest sensitivity obtained in small absorption measurements in liquids. (C) 2009 Optical Society of America

A monohydrate pseudopolymorph of 3,4-dihydroxybenzophenone and the role of water in the crystal assembly of benzophenones

During a polymorphism screening of hydroxybenzophenone derivatives, a monohydrate pseudopolymorph of (3,4-dihydroxyphenyl)(phenyl)methanone, C(13)H(10)O(3)center dot H(2)O, (I), was obtained. Structural relationships and the role of water in crystal assembly were established on the basis of the known anhydrous form [Cox, Kechagias & Kelly (2008). Acta Cryst. B64, 206-216]. The crystal packing of (I) is stabilized by classical intermolecular O-H...O hydrogen bonds, generating a three-dimensional network.

Thermo-solvatochromism in binary mixtures of water and ionic liquids: on the relative importance of solvophobic interactions

The thermo-solvatochromism of 2,6-dibromo-4-[(E)-2-(1-methylpyridinium-4-yl)ethenyl] phenolate, MePMBr(2), has been studied in mixtures of water, W, with ionic liquids, ILs, in the temperature range of 10 to 60 degrees C, where feasible. The objectives of the study were to test the applicability of a recently introduced solvation model, and to assess the relative importance of solute-solvent solvophobic interactions. The ILs were 1-allyl-3-alkylimidazolium chlorides, where the alkyl groups are methyl, 1-butyl, and 1-hexyl, respectively. The equilibrium constants for the interaction of W and the ILs were calculated from density data; they were found to be linearly dependent on N(C), the number of carbon atoms of the alkyl group; van't Hoff equation (log K versus 1/T) applied satisfactorily. Plots of the empirical solvent polarities, E(T) (MePMBr(2)) in kcal mol(-1), versus the mole fraction of water in the binary mixture, chi(w), showed non-linear, i.e., non-ideal behavior. The dependence of E(T) (MePMBr(2)) on chi(w), has been conveniently quantified in terms of solvation by W, IL, and the ""complex"" solvent IL-W. The non-ideal behavior is due to preferential solvation by the IL and, more efficiently, by IL-W. The deviation from linearity increases as a function of increasing N(C) of the IL, and is stronger than that observed for solvation of MePMBr(2) by aqueous 1-propanol, a solvent whose lipophilicity is 12.8 to 52.1 times larger than those of the ILs investigated. The dependence on N(C) is attributed to solute-solvent solvophobic interactions, whose relative contribution to solvation are presumably greater than that in mixtures of water and 1-propanol.

Determination of parathion and carbaryl pesticides in water and food samples using a self assembled monolayer/acetylcholinesterase electrochemical biosensor

An acetylcholinesterase (AchE) based amperometric biosensor was developed by immobilisation of the enzyme onto a self assembled modified gold electrode. Cyclic voltammetric experiments performed with the SAM-AchE biosensor in phosphate buffer solutions ( pH = 7.2) containing acetylthiocholine confirmed the formation of thiocholine and its electrochemical oxidation at E-p = 0.28 V vs Ag/AgCl. An indirect methodology involving the inhibition effect of parathion and carbaryl on the enzymatic reaction was developed and employed to measure both pesticides in spiked natural water and food samples without pre-treatment or pre-concentration steps. Values higher than 91-98.0% in recovery experiments indicated the feasibility of the proposed electroanalytical methodology to quantify both pesticides in water or food samples. HPLC measurements were also performed for comparison and confirmed the values measured amperometrically.