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.

Disordered ensembles of random matrices

It is shown that the families of generalized matrix ensembles recently considered which give rise to an orthogonal invariant stable Levy ensemble can be generated by the simple procedure of dividing Gaussian matrices by a random variable. The nonergodicity of this kind of disordered ensembles is investigated. It is shown that the same procedure applied to random graphs gives rise to a family that interpolates between the Erdos-Renyi and the scale free models.

Ionic liquids as recycling solvents for the synthesis of magnetic nanoparticles

This work describes an easy synthesis (one pot) of MFe(2)O(4) (M = Co, Fe, Mn, and Ni) magnetic nanoparticles MNPs by the thermal decomposition of Fe(Acac)(3)/M(Acac)(2) by using BMI center dot NTf(2) (1-n-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) or BMI center dot PF(6) (1-n-butyl-3-methylimidazolium hexafluorophosphate) ionic liquids (ILs) as recycling solvents and oleylamine as the reducing and surface modifier agent. The effects of reaction temperature and reaction time on the features of the magnetic nanomaterials (size and magnetic properties) were investigated. The growth of the MNPs is easily controlled in the IL by adjusting the reaction temperature and time, as inferred from Fe(3)O(4) MNPs obtained at 150 degrees C, 200 degrees C and 250 degrees C with mean diameters of 8, 10 and 15 nm, respectively. However, the thermal decomposition of Fe(Acac)(3) performed in a conventional high boiling point solvent (diphenyl ether, bp 259 degrees C), under a similar Fe to oleylamine molar ratio used in the IL synthesis, does not follow the same growth mechanism and rendered only smaller NPs of 5 nm mean diameter. All MNPs are covered by at least one monolayer of oleylamine making them readily dispersible in non-polar solvents. Besides the influence on the nanoparticles growth, which is important for the preparation of highly crystalline MNPs, the IL was easily recycled and has been used in at least 20 successive syntheses.

Quantum and pseudoclassical descriptions of nonrelativistic spinning particles in noncommutative space

We construct a nonrelativistic wave equation for spinning particles in the noncommutative space (in a sense, a theta modification of the Pauli equation). To this end, we consider the nonrelativistic limit of the theta-modified Dirac equation. To complete the consideration, we present a pseudoclassical model (a la Berezin-Marinov) for the corresponding nonrelativistic particle in the noncommutative space. To justify the latter model, we demonstrate that its quantization leads to the theta-modified Pauli equation. We extract theta-modified interaction between a nonrelativistic spin and a magnetic field from such a Pauli equation and construct a theta modification of the Heisenberg model for two coupled spins placed in an external magnetic field. In the framework of such a model, we calculate the probability transition between two orthogonal Einstein-Podolsky-Rosen states for a pair of spins in an oscillatory magnetic field and show that some of such transitions, which are forbidden in the commutative space, are possible due to the space noncommutativity. This allows us to estimate an upper bound on the noncommutativity parameter.

Theoretical study of one- and two-photon absorption spectra of azoaromatic compounds

In this study, the one- and two-photon absorption spectra of seven azoaromatic compounds (five pseudostilbenes-type and two aminoazobenzenes) were theoretically investigated using the density functional theory combined with the response functions formalism. The equilibrium molecular structure of each compound was obtained at three different levels of theory: Hartree-Fock, density functional theory (DFT), and Moller-Plesset 2. The effect of solvent on the equilibrium structure and the electronic transitions of the compounds were investigated using the polarizable continuum model. For the one-photon absorption, the allowed pi ->pi(*) transition energy showed to be dependent on the molecular structures and the effect of solvent, while the n ->pi(*) and pi ->pi(*)(n) transition energies exhibited only a slight dependence. An inversion between the bands corresponding to the pi ->pi(*) and n ->pi(*) states due to the effect of solvent was observed for the pseudostilbene-type compounds. To characterize the allowed two-photon absorption transitions for azoaromatic compounds, the response functions formalism combined with DFT using the hybrid B3LYP and PBE0 functionals and the long-range corrected CAM-B3LYP functional was employed. The theoretical results support the previous findings based on the three-state model. The model takes into account the ground and two electronic excited states and has already been used to describe and interpret the two-photon absorption spectrum of azoaromatic compounds. The highest energy two-photon allowed transition for the pseudostilbene-type compounds shows to be more effectively affected (similar to 20%) by the torsion of the molecular structure than the lowest allowed transition (similar to 10%). In order to elucidate the effect of the solvent on the two-photon absorption spectra, the lowest allowed two-photon transition (dipolar transition) for each compound was analyzed using a two-state approximation and the polarizable continuum model. The results obtained reveal that the effect of solvent increases drastically the two-photon cross-section of the dipolar transition of the pseudostilbene-type compounds. In general, the features of both one- and two-photon absorption spectra of the azoaromatic compounds are well reproduced by the theoretical calculations.

Hyper-Rayleigh scattering with picosecond pulse trains

We propose a method for measuring hyper-Rayleigh scattering employing pulse trains produced by a Q-switched and mode-locked Nd:YAG laser. The use of the entire pulse train under the Q-switch envelope avoids the need of any device to scan the irradiance, as is usually done with nanosecond and femtosecond single-pulse lasers. To verify the feasibility of the technique, we performed measurements in different solutions of para-nitroaniline and compared the results with those obtained with nanosecond pulses. In both cases, the agreement with the hyperpolarizability values reported in the literature is about the same, but the measurements carried out with pulse trains are at least 20 times faster. Besides the advantage of acquisition speed, the use of pulse trains also allows the instantaneous inspection of slow luminescence contributions arising from multiphoton absorption. (C) 2008 Optical Society of America.

INFERENCE FOR EIGENVALUES AND EIGENVECTORS OF GAUSSIAN SYMMETRIC MATRICES

This article presents maximum likelihood estimators (MLEs) and log-likelihood ratio (LLR) tests for the eigenvalues and eigenvectors of Gaussian random symmetric matrices of arbitrary dimension, where the observations are independent repeated samples from one or two populations. These inference problems are relevant in the analysis of diffusion tensor imaging data and polarized cosmic background radiation data, where the observations are, respectively, 3 x 3 and 2 x 2 symmetric positive definite matrices. The parameter sets involved in the inference problems for eigenvalues and eigenvectors are subsets of Euclidean space that are either affine subspaces, embedded submanifolds that are invariant under orthogonal transformations or polyhedral convex cones. We show that for a class of sets that includes the ones considered in this paper, the MLEs of the mean parameter do not depend on the covariance parameters if and only if the covariance structure is orthogonally invariant. Closed-form expressions for the MLEs and the associated LLRs are derived for this covariance structure.

A note on the likelihood and moments of the skew-normal distribution

In this paper an alternative approach to the one in Henze (1986) is proposed for deriving the odd moments of the skew-normal distribution considered in Azzalini (1985). The approach is based on a Pascal type triangle, which seems to greatly simplify moments computation. Moreover, it is shown that the likelihood equation for estimating the asymmetry parameter in such model is generated as orthogonal functions to the sample vector. As a consequence, conditions for a unique solution of the likelihood equation are established, which seem to hold in more general setting.

The impact of the subtropical South Atlantic SST on South American precipitation

The Community Climate Model (CCM3) from the National Center for Atmospheric Research (NCAR) is used to investigate the effect of the South Atlantic sea surface temperature (SST) anomalies on interannual to decadal variability of South American precipitation. Two ensembles composed of multidecadal simulations forced with monthly SST data from the Hadley Centre for the period 1949 to 2001 are analysed. A statistical treatment based on signal-to-noise ratio and Empirical Orthogonal Functions (EOF) is applied to the ensembles in order to reduce the internal variability among the integrations. The ensemble treatment shows a spatial and temporal dependence of reproducibility. High degree of reproducibility is found in the tropics while the extratropics is apparently less reproducible. Austral autumn (MAM) and spring (SON) precipitation appears to be more reproducible over the South America-South Atlantic region than the summer (DJF) and winter (JJA) rainfall. While the Inter-tropical Convergence Zone (ITCZ) region is dominated by external variance, the South Atlantic Convergence Zone (SACZ) over South America is predominantly determined by internal variance, which makes it a difficult phenomenon to predict. Alternatively, the SACZ over western South Atlantic appears to be more sensitive to the subtropical SST anomalies than over the continent. An attempt is made to separate the atmospheric response forced by the South Atlantic SST anomalies from that associated with the El Nino - Southern Oscillation (ENSO). Results show that both the South Atlantic and Pacific SSTs modulate the intensity and position of the SACZ during DJF. Particularly, the subtropical South Atlantic SSTs are more important than ENSO in determining the position of the SACZ over the southeast Brazilian coast during DJF. On the other hand, the ENSO signal seems to influence the intensity of the SACZ not only in DJF but especially its oceanic branch during MAM. Both local and remote influences, however, are confounded by the large internal variance in the region. During MAM and JJA, the South Atlantic SST anomalies affect the magnitude and the meridional displacement of the ITCZ. In JJA, the ENSO has relatively little influence on the interannual variability of the simulated rainfall. During SON, however, the ENSO seems to counteract the effect of the subtropical South Atlantic SST variations on convection over South America.

A microscopic view of substitution reactions solvated by ionic liquids

The solvation effect of the ionic liquid 1-N-butyl-3-methylimidazolium hexafluorophosphate on nucleophilic substitution reactions of halides toward the aliphatic carbon of methyl p-nitrobenzenesulfonate (pNBS) was investigated by computer simulations. The calculations were performed by using a hybrid quantum-mechanical/molecular-mechanical (QM/MM) methodology. A semiempirical Hamiltonian was first parametrized on the basis of comparison with ab initio calculations for Cl(-) and Br(-) reaction with pNBS at gas phase. In condensed phase, free energy profiles were obtained for both reactions. The calculated reaction barriers are in agreement with experiment. The structure of species solvated by the ionic liquid was followed along the reaction progress from the reagents, through the transition state, to the final products. The simulations indicate that this substitution reaction in the ionic liquid is slower than in nonpolar molecular solvents proper to significant stabilization of the halide anion by the ionic liquid in comparison with the transition state with delocalized charge. Solute-solvent interactions in the first solvation shell contain several hydrogen bonds that are formed or broken in response to charge density variation along the reaction coordinate. The detailed structural analysis can be used to rationalize the design of new ionic liquids with tailored solvation properties. (c) 2008 American Institute of Physics.

Ionic mobility of the solvated proton and acid-base titration in a four-compartment capillary electrophoresis system

Although H(+) and OH(-) are the most common ions in aqueous media, they are not usually observable in capillary electrophoresis (CE) experiments, because of the extensive use of buffer solutions as the background electrolyte. In the present work, we introduce CE equipment designed to allow the determination of such ions in a similar fashion as any other ion. Basically, it consists of a four-compartment piece of equipment for electrolysis-separated experiments (D. P. de Jesus et at, Anal. Chem., 2005, 77, 607). In such a system, the ends of the capillary are placed in two reservoirs, which are connected to two other reservoirs through electrolyte-filled tubes. The electrodes of the high-voltage power source are positioned in these reservoirs. Thus, the electrolysis products are kept away from the inputs of the capillary. The detection was provided by two capacitively coupled contactless conductivity detectors (CD), each one positioned about 11 cm from the end of the capillary. Two applications were demonstrated: titration-like procedures for nanolitre samples and mobility measurements. Strong and weak acids (pK(a) < 5), pure or mixtures, could be titrated. The analytical curve is linear from 50 mu M up to 10 mM of total dissociable hydrogen (r = 0.99899 for n =10) in 10-nL samples. By including D(2)O in the running electrolyte, we could demonstrate how to measure the mixed proton/deuteron mobility. When H(2)O/D(2)O (9 : 1 v/v) was used as the solvent, the mobility was 289.6 +/- 0.5 x 10(-5) cm(2) V(-1) s(-1). Due to the fast conversion of the species, this value is related to the overall behaviour of all isotopologues and isotopomers of the Zundel and Eigen structures, as well as the Stokesian mobility of proton and deuteron. The effect of neutral (o-phenanthroline) and negatively charged (chloroacetate) bases and aprotic solvent (DMSO) over the H(+) mobility was also demonstrated.

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.

Understanding solvation

The effects of solvents on different chemical phenomena, including reactivity, spectroscopic data, and swelling of biopolymers can be rationalized by use of solvatochromic probes, substances whose UV-vis spectra, absorption, or emission are sensitive to the properties of the medium. Thermo-solvatochromism refers to the effect of temperature on solvatochromism. The study of both phenomena sheds light on the relative importance of the factors that contribute to solvation, namely, properties of the probe, those of the solvent (acidity, basicity, dipolarity/polarizability, and lipophilicity), and the temperature. Solvation in binary solvent mixtures is complex because of ""preferential solvation"" of the probe by some component of the mixture. A recently introduced solvent exchange model is based on the presence in the binary solvent mixture of the organic component (molecular solvent or ionic liquid), S, water, W, and a 1:1 hydrogen-bonded species (S-W). Solvation by the latter is more efficient than by its precursor solvents, due to probe-solvent hydrogen-bonding and hydrophobic interactions; dimethyl sulfoxide (DMSO)-W is an exception. Solvatochromic data are employed in order to explain apparently disconnected phenomena, namely, medium effect on the pH-independent hydrolysis of esters, (1)H NMR data of water-ionic liquid (IL) mixtures, and the swelling of cellulose.