Monte Carlo investigations of intermolecular interactions in water-amide mixtures

Monte Carlo simulations of water-amides (amide=fonnamide-FOR, methylfonnamide-NMF and dimethylformamide-DMF) solutions have been carried out in the NpT ensemble at 308 K and 1 atm. The structure and excess enthalpy of the mixtures as a function of the composition have been investigated. The TIP4P model was used for simulating water and six-site models previously optimized in this laboratory were used for simulating the liquid amides. The intermolecular interaction energy was calculated using the classical 6-12 Lennard-Jones potential plus a Coulomb term. The interaction energy between solute and solvent has been partitioned what leads to a better understanding of the behavior of the enthalpy of mixture obtained for the three solutions experimentally. Radial distribution functions for the water-amides correlations permit to explore the intermolecular interactions between the molecules. The results show that three, two and one hydrogen bonds between the water and the amide molecules are formed in the FOR, NMF and DMF-water solutions, respectively. These H-bonds are, respectively, stronger for DMF-water, NMF-water and FOR-water. In the NMF-water solution, the interaction between the methyl group of the NMF and the oxygen of the water plays a role in the stabilization of the aqueous solution quite similar to that of an H-bond in the FOR-water solution. (c) 2005 Elsevier B.V. All rights reserved.

Monte Carlo thermodynamic and structural properties of the TIP4P water model: dependence on the computational conditions

Classical Monte Carlo simulations were carried out on the NPT ensemble at 25°C and 1 atm, aiming to investigate the ability of the TIP4P water model [Jorgensen, Chandrasekhar, Madura, Impey and Klein; J. Chem. Phys., 79 (1983) 926] to reproduce the newest structural picture of liquid water. The results were compared with recent neutron diffraction data [Soper; Bruni and Ricci; J. Chem. Phys., 106 (1997) 247]. The influence of the computational conditions on the thermodynamic and structural results obtained with this model was also analyzed. The findings were compared with the original ones from Jorgensen et al [above-cited reference plus Mol. Phys., 56 (1985) 1381]. It is notice that the thermodynamic results are dependent on the boundary conditions used, whereas the usual radial distribution functions g(O/O(r)) and g(O/H(r)) do not depend on them.

C-H center dot center dot center dot O and N-H center dot center dot center dot O hydrogen bonds in liquid amides investigated by Monte Carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N,N-dimethylformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ... O and N-H ... O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ... O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This Interaction is particularly important in the structure of MF. The intensity of the N-H ... O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. (C) 1997 John Wiley & Sons, Inc.

Theoretical studies of liquids by computer simulations: the water-acetone mixture.

Monte Carlo simulation results for pure liquid acetone and water-acetone mixtures calculated in the isothermal and isobaric (NPT) ensemble at T=298K and p=1.0atm are presented. The TIP4P model was used for water and optimized potential for liquid simulation (OPLS) force field parameters used for acetone. The results obtained for the average configurational energy as a function of the mole fraction are in good accord with experimental data. Energy partitioning and co-ordination numbers results calculated for equimolar water-acetone solution are compared to similar data obtained for other water-organic liquid mixtures. These results show an increase in water-water interaction energy and co-ordination numbers when the interaction between water and organic liquid molecules decrease. Distribution functions for pure liquid acetone and water-acetone mixtures are presented. Dipole-dipole angular correlation functions obtained for pure liquid acetone show a predominance of dimers with parallel alignment of dipole moments. Radial distribution functions from water-acetone interaction show characteristic features of hydrogen bonded liquids. Radial and angular distribution functions for water-water correlation calculated in pure water and in equimolar water-acetone mixture are compared, showing very similar features in both systems. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Application of a new reverse Monte Carlo algorithm to polyatomic molecular systems. I. Liquid water

Using a new reverse Monte Carlo algorithm, we present simulations that reproduce very well several structural and thermodynamic properties of liquid water. Both Monte Carlo, molecular dynamics simulations and experimental radial distribution functions used as input are accurately reproduced using a small number of molecules and no external constraints. Ad hoc energy and hydrogen bond analysis show the physical consistency and limitations of the generated RMC configurations. (C) 2001 American Institute of Physics.

X-ray diffraction investigation of Co(II) ions in borosilicate glasses

The radial distribution functions of cobalt glasses with 6%, 8%, and 14% CoO are compared with those of suitable cobalt-free borosilicate matrices leading to difference distribution curves representative of the cobalt structural arrangement. Analysis of the curves indicates that cobalt ions are surrounded by approximately four oxygen neighbors at the distance expected for fourfold coordination. © 1986 American Institute of Physics.

C - H ⋯ 0 and N - H ⋯ 0 hydrogen bonds in liquid amides investigated by monte carlo simulation

Monte Carlo simulations of liquid formamide, N-methylformamide (MF), and N, N-dimethytformamide (DMF) have been performed in the isothermal and isobaric ensemble at 298 K and 1 atm, aiming to investigate the C-H ⋯ O and N-H ⋯ O hydrogen bonds. The interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term on a rigid six-site molecular model with the potential parameters being optimized in this work. Theoretical values obtained for heat of vaporization and liquid densities are in good agreement with the experimental data. The radial distribution function [RDF, g(r)] obtained compare well with R-X diffraction data available. The RDF and molecular mechanics (MM2) minimization show that the C-H ⋯ O interaction has a significant role in the structure of the three liquids. These results are supported by ab initio calculations. This interaction is particularly important in the structure of MF. The intensity of the N - H ⋯ O hydrogen bond is greater in the MF than formamide. This could explain some anomalous properties verified in MF. © 1997 John Wiley & Sons, Inc.

Study of water and dimethylformamide interaction by computer simulation

Monte Carlo simulations of water-dimethylformamide (DMF) mixtures were performed in the isothermal and isobaric ensemble at 298.15 K and 1 atm. The intermolecular interaction energy was calculated using the classical 6-12 Lennard-Jones pairwise potential plus a Coulomb term. The TIP4P model was used for simulating water molecules, and a six-site model previously optimised by us was used to represent DMF. The potential energy for the water-DMF interaction was obtained via standard geometric combining rules using the original potential parameters for the pure liquids. The radial distribution functions calculated for water-DMF mixtures show well characterised hydrogen bonds between the oxygen site of DMF and hydrogen of water. A structureless correlation curve was observed for the interaction between the hydrogen site of the carbonyl group and the oxygen site of water. Hydration effects on the stabilisation of the DMF molecule in aqueous solution have been investigated using statistical perturbation theory. The results show that energetic changes involved in the hydration process are not strong enough to stabilise another configuration of DMF than the planar one.

Effective pair potentials and structure of water

Classical Monte Carlo calculations have been performed in order to investigate the ability of the TIP4P, SPC, and SPCE water models to reproduce the structural features of liquid water. The simulations were carried out in the NPT ensemble at 4 thermodynamic conditions. The results are compared with recent neutron diffraction data. Essentially, the three models capture equally well the thermodynamic and structural features of water. Although they were parametrized to reproduce the water properties at ambient conditions, the agreement with the experimental pair correlation functions was even better at supercritical conditions. This is because the effective pair potentials have some difficulty to reproduce cooperative interactions, like hydrogen bonds. These interactions are less effective at supercritical conditions, where the liquid behaves roughly like a gas.

Hydrogen bonding donation of N-methylformamide with dimethylsulfoxide and water

20% N-methylformamide (NMF) mixtures with water and with dimethylsulfoxide (DMSO) have been studied. A comparison between the hydrogen bonding (H-bond) donation of N-methylformamide with both solvents in the mixtures is presented. Results of radial distribution functions, pair distribution energies, molecular dipole moment correlation, and geometry of the H-bonded species in each case are shown. The results indicate that the NMF-solvent H-bond is significantly stronger with DMSO than with water. The solvation shell is best organized in the DMSO mixture than in the aqueous one. © 2013 Elsevier B.V. All rights reserved.

Monte Carlo simulation of water-tetrahydrofuran mixtures

Monte Carlo simulations of water-tetrahydrofuran (THF) mixtures were performed in the isothermal-isobaric ensemble (NPT) at T = 298 K and p = 1 atm. The interaction energy was calculated using the TIP4P model for water and a five-site united atom representation for the THF molecule. The potential energy surfaces for water-THF interactions were obtained by using combining rules and the original potential functions used for pure liquids. Theoretical values obtained for the average interaction energy as a function of concentration are in good agreement with available experimental data. Results from the partitioning of the total interaction energy into water-water, water-THF and THF-THF contributions are presented. These results are useful to distinguish between the quantitative contributions of these molecular interactions to the energetic behavior of the water-THF mixing process. The radial distribution functions for HW-OTHF and OW-OTHF site-site interactions show the salient features of hydrogen-bonded liquids. Comparison of the average number of water-water complexes interacting through hydrogen bonding in water-THF and water-methanol mixtures shows an enhancement of the water-water coordination number in a THF rich environment. © 1995.

Network reconfiguration and loss allocation in a deregulated environment of distribution systems

Low flexibility and reliability in the operation of radial distribution networks make those systems be constructed with extra equipment as sectionalising switches in order to reconfigure the network, so the operation quality of the network can be improved. Thus, sectionalising switches are used for fault isolation and for configuration management (reconfiguration). Moreover, distribution systems are being impacted by the increasing insertion of distributed generators. Hence, distributed generation became one of the relevant parameters in the evaluation of systems reconfiguration. Distributed generation may affect distribution networks operation in various ways, causing noticeable impacts depending on its location. Thus, the loss allocation problem becomes more important considering the possibility of open access to the distribution networks. In this work, a graphic simulator for distribution networks with reconfiguration and loss allocation functions, is presented. Reconfiguration problem is solved through a heuristic methodology, using a robust power flow algorithm based on the current summation backward-forward technique, considering distributed generation. Four different loss allocation methods (Zbus, Direct Loss Coefficient, Substitution and Marginal Loss Coefficient) are implemented and compared. Results for a 32-bus medium voltage distribution network, are presented and discussed.

Short-term expansion planning of radial electrical distribution systems using mixed-integer linear programming

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Reconfiguration of radial electric power distribution system via a scatter search algorithm

In this paper, to solve the reconfiguration problem of radial distribution systems a scatter search, which is a metaheuristic-based algorithm, is proposed. In the codification process of this algorithm a structure called node-depth representation is used. It then, via the operators and from the electrical power system point of view, results finding only radial topologies. In order to show the effectiveness, usefulness, and the efficiency of the proposed method, a commonly used test system, 135-bus, and a practical system, a part of Sao Paulo state's distribution network, 7052 bus, are conducted. Results confirm the efficiency of the proposed algorithm that can find high quality solutions satisfying all the physical and operational constraints of the problem.

Comparative study between RBF and radial-PPS neural networks

The study of function approximation is motivated by the human limitation and inability to register and manipulate with exact precision the behavior variations of the physical nature of a phenomenon. These variations are referred to as signals or signal functions. Many real world problem can be formulated as function approximation problems and from the viewpoint of artificial neural networks these can be seen as the problem of searching for a mapping that establishes a relationship from an input space to an output space through a process of network learning. Several paradigms of artificial neural networks (ANN) exist. Here we will be investigated a comparative of the ANN study of RBF with radial Polynomial Power of Sigmoids (PPS) in function approximation problems. Radial PPS are functions generated by linear combination of powers of sigmoids functions. The main objective of this paper is to show the advantages of the use of the radial PPS functions in relationship traditional RBF, through adaptive training and ridge regression techniques.