Structure and electronic properties of hydrated mesityl oxide: a sequential quantum mechanics/molecular mechanics approach

The hydration of mesityl oxide (MOx) was investigated through a sequential quantum mechanics/molecular mechanics approach. Emphasis was placed on the analysis of the role played by water in the MOx syn-anti equilibrium and the electronic absorption spectrum. Results for the structure of the MOx-water solution, free energy of solvation and polarization effects are also reported. Our main conclusion was that in gas-phase and in low-polarity solvents, the MOx exists dominantly in syn-form and in aqueous solution in anti-form. This conclusion was supported by Gibbs free energy calculations in gas phase and in-water by quantum mechanical calculations with polarizable continuum model and thermodynamic perturbation theory in Monte Carlo simulations using a polarized MOx model. The consideration of the in-water polarization of the MOx is very important to correctly describe the solute-solvent electrostatic interaction. Our best estimate for the shift of the pi-pi* transition energy of MOx, when it changes from gas-phase to water solvent, shows a red-shift of -2,520 +/- 90 cm(-1), which is only 110 cm(-1) (0.014 eV) below the experimental extrapolation of -2,410 +/- 90 cm(-1). This red-shift of around -2,500 cm(-1) can be divided in two distinct and opposite contributions. One contribution is related to the syn -> anti conformational change leading to a blue-shift of similar to 1,700 cm(-1). Other contribution is the solvent effect on the electronic structure of the MOx leading to a red-shift of around -4,200 cm(-1). Additionally, this red-shift caused by the solvent effect on the electronic structure can by composed by approximately 60 % due to the electrostatic bulk effect, 10 % due to the explicit inclusion of the hydrogen-bonded water molecules and 30 % due to the explicit inclusion of the nearest water molecules.

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

Tailoring Molecular Architectures with Cobalt Tetrasulfonated Phthalocyanine: Immobilization in Layer-by-Layer Films and Sensing Applications

In the field of organic thin films, manipulation at the nanoscale can be obtained by immobilization of different materials on platforms designed to enhance a specific property via the layer-by-layer technique. In this paper we describe the fabrication of nanostructured films containing cobalt tetrasulfonated phthalocyanine (CoTsPc) obtained through the layer-by-layer architecture and assembled with linear poly(allylamine hydrochloride) (PAH) and poly(amidoamine) dendrimer (PAMAM) polyelectrolytes. Film growth was monitored by UV-vis spectroscopy following the Q band of CoTsPc and revealed a linear growth for both systems. Fourier transform infrared (FTIR) spectroscopy showed that the driving force keeping the structure of the films was achieved upon interactions of CoTsPc sulfonic groups with protonated amine groups present in the positive polyelectrolyte. A comprehensive SPR investigation on film growth reproduced the deposition process dynamically and provided an estimation of the thicknesses of the layers. Both FTIR and SPR techniques suggested a preferential orientation of the Pc ring parallel to the substrate. The electrical conductivity of the PAH films deposited on interdigitated electrodes was found to be very sensitive to water vapor. These results point to the development of a phthalocyanine-based humidity sensor obtained from a simple thin film deposition technique, whose ability to tailor molecular organization was crucial to achieve high sensitivity.

Molecular organization and doping in poly(2-methoxyaniline)/Ni(dmit)(2) films obtained with the Langmuir-Blodgett technique

The control of the properties of materials at the molecular level is pursued for many applications, especially those associated with nanostructures. In this paper, we show that the coordination compound [Ni(dmit)(2)], where (dmit) is the 1,3-dithiole-2-thione-4,5-dithiolate ligand, can induce doping of poly(2-methoxyaniline) (POMA) in molecularly ordered Langmuir and Langmuir-Blodgett (LB) films. Doping was associated with interactions between the components and the compression of the Langmuir film at the air-water interface, according to polarization-modulated infrared reflection-absorption spectroscopy (PM-IRRAS) data. Taking these results together with in situ UV-Vis absorption measurements, we could identify the molecular groups involved in the interaction, including the way they were reoriented upon film compression. The Langmuir films were sufficiently stable to be transferred as Y-type LB films, while the hybrid POMA/[Ni(dmit)(2)] films remain doped in the solid state. As expected, the molecular charges affected the film morphology, as observed from combined atomic and electric force microscopy measurements. In summary, with adequate spectroscopy and microscopy tools we characterized molecular-level interactions, which may allow one to design molecular electronic devices with controlled electrical properties.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Effect of alkyl chain length of alcohols on nematic uniaxial-to-biaxial phase transitions in a potassium laurate/alcohol/K2SO4/water lyotropic mixture

Lyotropic liquid crystalline quaternary mixtures of potassium laurate (KL), potassium sulphate (K2SO4)/alcohol (n-OH)/water, with the alcohols having different numbers of carbon atoms in the alkyl chain (n), from 1-octanol to 1-hexadecanol, were investigated by optical techniques (optical microscopy and laser conoscopy). The biaxial nematic phase domain is present in a window of values of n = n(KL) +/- 2, where n(KL) = 11 is the number of carbon atoms in the alkyl chain of KL. The biaxial phase domain became smaller and the uniaxial-to-biaxial phase transition temperatures shifted to relatively higher temperatures upon going from 1-nonanol to 1-tridecanol. Moreover, compared with other lyotropic mixtures these new mixtures present high birefringence values, which we expect to be related to the micellar shape anisotropy. Our results are interpreted assuming that alcohol molecules tend to segregate in the micelles in a way that depends on the relative value of n with respect to nKL. The larger the value of n, the more alcohol molecules tend to be located in the curved parts of the micelle, favoring the uniaxial nematic calamitic phase with respect to the biaxial and uniaxial discotic nematic phases.

An environmentally friendly analytical procedure for nickel determination by atomic and molecular spectrometry after cloud point extraction in different samples

Cloud point extraction (CPE) was employed for separation and preconcentration prior to the determination of nickel by graphite furnace atomic absorption spectrometry (GFAAS), flame atomic absorption spectrometry (FAAS) or UV-Vis spectrophotometry. Di-2-pyridyl ketone salicyloylhydrazone (DPKSH) was used for the first time as a complexing agent in CPE. The nickel complex was extracted from the aqueous phase using the Triton X-114 surfactant. Under optimized conditions, limits of detection obtained with GFAAS, FAAS and UV-Vis spectrophotometry were 0.14, 0.76 and 1.5 mu g L-1, respectively. The extraction was quantitative and the enrichment factor was estimated to be 27. The method was applied to natural waters, hemodialysis concentrates, urine and honey samples. Accuracy was evaluated by analysis of the NIST 1643e Water standard reference material.

Epidemiological, parasitological and molecular aspects of Giardia duodenalis infection in children attending public daycare centers in southeastern Brazil

The purpose of this study was to determine the prevalence, associated risk factors and genotype of Giardia duodenalis infection in children attending public daycare centers in the city of Araguari, state of Minas Gerais, Brazil. Fecal samples were collected from 245 children aged 0-5 years, and questionnaires were asked about sociodemographic and hygiene-related characteristics. At the daycare centers where children tested positive, fecal samples were collected from the staff handling food, and from family members and domestic animals. Positive samples were analyzed at the dehydrogenase glutamate (gdh) locus to determine the genotype. The prevalence of G. duodenalis was 51.8%, and drinking unfiltered and unboiled water (OR 2.12, CI 1.26-3.69, p<0.001) and washing hands only with water (OR 2.14, Cl 1.19-4.04, p<0.001) were related risk factors. No association was found between test-positive children anti their family members, domestic animals and food handlers. An analysis of the sequences of 30 samples revealed that they all belonged to genotype B. (C) 2012 Royal Society of Tropical Medicine and Hygiene. Published by Elsevier Ltd. All rights reserved.

Ontogenetic and evolutionary change of external morphology of the neotropical shrimp potimirim (Holthuis, 1954) explained by a molecular phylogeny of the genus

A phylogenetic analysis of a fragment of the mitochondrial gene 16S was used to test the monophyletic status of Potimirim. Existing doubts on the taxonomic status of brasiliana (once P glabra) and P potimirim (once P mexicana) were clarified. Potimirim mexicana and P potimirim are distinct species according to molecular data and appendix masculina morphology. A new species (Potimirim sp. 1) from Puerto Rico was revealed with molecular data, and it is evolutionarily related to P potimirim and P mexicana according to our analysis. We found out three distinct species under the name P glabra. Then, we recommend the application of the name P glabra for the populations of the Pacific slope of Central America and revalidation of P brasiliana for the Brazilian ones. The need for a new name to those "P glabra" of the Caribbean is highlighted, and it was provisionally referred as Potimirim sp. 2. The ontogenetic (juveniles to adults) development of the appendix masculina of P brasiliana was observed and compared to the other species of Potimirim (adults). In the light of our phylogenetic hypothesis, we postulate a pattern of character addition for the evolution of the appendix masculina of Potimirim. This hypothesis is plausible for two key reasons. First. Potimirim is a monophyletic group according to our hypothesis. Second, the shape of appendix masculina found in adults of P. americana is similar and comparable to those found in the earliest juvenile stages of P brasiliana, a derived species according to our phylogeny (P americana, ((P mexicana, Potimirim sp. 1. P potimirim), (P glabra, (brasiliana, Potimirim sp. 2)))). As so, the basal P americana retain the ancestral morphological state of the appendix masculina when compared to the other species of Potimirim. In our interpretation the ontogeny of the appendix masculina recapitulated the proposed phylogeny, giving further support to it.

Molecular Dynamics Investigations of PRODAN in a DLPC Bilayer

Molecular dynamics computer simulations have been performed to identify preferred positions of the fluorescent probe PRODAN in a fully hydrated DLPC bilayer in the fluid phase. In addition to the intramolecular charge-transfer first vertical excited state, we considered different charge distributions for the electronic ground state of the PRODAN molecule by distinct atomic charge models corresponding to the probe molecule in vacuum as well as polarized in a weak and a strong dielectric solvent (cyclohexane and water). Independent on the charge distribution model of PRODAN, we observed a preferential orientation of this molecule in the bilayer with the dimethylamino group pointing toward the membrane's center and the carbonyl oxygen toward the membrane's interface. However, changing the charge distribution model of PRODAN, independent of its initial position in the equilibrated DLPC membrane, we observed different preferential positions. For the ground state representation without polarization and the in-cyclohexane polarization, the probe maintains its position close to the membrane's center. Considering the in-water polarization model, the probe approaches more of the polar headgroup region of the bilayer, with a strong structural correlation with the choline group, exposing its oxygen atom to water molecules. PRODAN's representation of the first vertical excited state with the in-water polarization also approaches the polar region of the membrane with the oxygen atom exposed to the bilayer's hydration shell. However, this model presents a stronger structural correlation with the phosphate groups than the ground state. Therefore, we conclude that the orientation of the PRODAN molecule inside the DLPC membrane is well-defined, but its position is very sensitive to the effect of the medium polarization included here by different models for the atomic charge distribution of the probe.

Self-Assembly of Poly(4-vynil-N-alkylpyridinium) Bromide onto Silica: Effect of Side-Chain Length on Structural Aspects at a Molecular Level

Self-assembly of poly(4-vynil-N-alkyl)pyridinium bromide with alkyl side chains of 2, 5, 7, 10, or 16 carbons from ethanolic solutions onto flat silica surfaces was studied by means of ellipsometry, atomic force microscopy (AFM), contact angle measurements, and sum-frequency generation (SFG) vibrational spectroscopy in the CH3 and CH2 stretch region. Ab initio quantum-chemical calculations on the N-alkylpyridinium side-group with restricted Hartree-Fock (RHF) method and 6-311G (d,p) basis set were C one to estimate the charge distribution along the pyridinium ring and the alkyl side-chain. SFG results showed that longer side chains promote the disorientation of the alkyl groups at the surface, corroborating with the contact angle values. AFM images revealed film homogeneity, regardless the alkyl side group. However, after 24 h contact with water, ringlike structures appeared on the film surfaces, when the polycation alkyl side chain had 7 or less carbons, and as the alkyl chain increased to 10 or 16 carbons, the films dewetted because the hydrophobic interactions prevailed over the electrostatic interactions between the pyridinium charged groups and the negatively charged SiO2 surface. Under acid conditions (HCl 0.1 mol.L-1), the film mean thickness values decreased up to 50% of original values when the alkyl side chains were ethyl or pentyl groups due to ion-pair disruption, but for longer groups they remained unchanged. Quantum-chemical optimization and Mulliken electron population showed that (i) from C2 to C15 the positive charge at the headgroup (HG) decreased 0.025, while the charge at combined HG + alpha-CH2 increased 0.037; and (ii) for C6 or longer, the alkyl side group presents a tilt in the geometry, moving away from the plane. Such effects summed up over the whole polymer chain give support to suggest that when the side chains are longer than 7 carbons, the hydrophobic interaction decreases film stability and increases acid resistance.

Applying many-body expansion aiming to make ab in initio molecular dynamics more efficient

In this present work we present a methodology that aims to apply the many-body expansion to decrease the computational cost of ab initio molecular dynamics, keeping acceptable accuracy on the results. We implemented this methodology in a program which we called ManBo. In the many-body expansion approach, we partitioned the total energy E of the system in contributions of one body, two bodies, three bodies, etc., until the contribution of the Nth body [1-3]: E = E1 + E2 + E3 + …EN. The E1 term is the sum of the internal energy of the molecules; the term E2 is the energy due to interaction between all pairs of molecules; E3 is the energy due to interaction between all trios of molecules; and so on. In Manbo we chose to truncate the expansion in the contribution of two or three bodies, both for the calculation of the energy and for the calculation of the atomic forces. In order to partially include the many-body interactions neglected when we truncate the expansion, we can include an electrostatic embedding in the electronic structure calculations, instead of considering the monomers, pairs and trios as isolated molecules in space. In simulations we made we chose to simulate water molecules, and use the Gaussian 09 as external program to calculate the atomic forces and energy of the system, as well as reference program for analyzing the accuracy of the results obtained with the ManBo. The results show that the use of the many-body expansion seems to be an interesting approach for reducing the still prohibitive computational cost of ab initio molecular dynamics. The errors introduced on atomic forces in applying such methodology are very small. The inclusion of an embedding electrostatic seems to be a good solution for improving the results with only a small increase in simulation time. As we increase the level of calculation, the simulation time of ManBo tends to largely decrease in relation to a conventional BOMD simulation of Gaussian, due to better scalability of the methodology presented. References [1] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 3, 46 (2007). [2] E. E. Dahlke and D. G. Truhlar; J. Chem. Theory Comput., 4, 1 (2008). [3] R. Rivelino, P. Chaudhuri and S. Canuto; J. Chem. Phys., 118, 10593 (2003).

Atomistic Study of Water Confined in Silica.

In this work, we have used a combined of atomistic simulation methods to explore the effects of confinement of water molecules between silica surfaces. Firstly, the mechanical properties of water severe confined (~3A) between two silica alpha-quartz was determined based on first principles calculations within the density functional theory (DFT). Simulated annealing methods were employed due to the complex potential energry surface, and the difficulties to avoid local minima. Our results suggest that much of the stiffness of the material (46%) remains, even after the insertion of a water monolayer in the silica. Secondly, in order to access typical time scales for confined systems, classical molecular dynamics was used to determine the dynamical properties of water confined in silica cylindrical pores, with diameters varying from 10 to 40A. in this case we have varied the passivation of the silica surface, from 13% to 100% of SiOH, and the other terminations being SiOH2 and SiOH3, the distribution of the different terminations was obtained with a Monte Carlo simulation. The simulations indicates a lowering of the diffusion coefficientes as the diameter decreases, due to the structuration of hydrogen bonds of water molecules; we have also obtained the density profiles of the confined water and the interfacial tension.