DSpace Collection:
http://hdl.handle.net/2117/395
Sat, 01 Nov 2014 04:20:28 GMT2014-11-01T04:20:28Zwebmaster.bupc@upc.eduUniversitat Politècnica de Catalunya. Servei de Biblioteques i DocumentaciónoPolyaniline emeraldine salt in the amorphous solid state: polaron versus bipolaron
http://hdl.handle.net/2117/24516
Title: Polyaniline emeraldine salt in the amorphous solid state: polaron versus bipolaron
Authors: Canales Gabriel, Manel; Torras Costa, Juan; Fabregat Jove, Georgina; Meneguzzi, Alvaro; Alemán Llansó, Carlos
Abstract: The polaronic and bipolaronic forms of polyaniline emeraldine salt (PAni-ES) in the amorphous solid state have been simulated using classical molecular dynamics (MD) and hybrid quantum mechanical/molecular mechanical-molecular dynamics (QM/MM-MD) approaches. It should be remarked that the electronic state of PAni-ES has been theoretically investigated in the gas phase, solution phase, and crystalline state, but this is the first study in the amorphous solid state, which is the most typical for this conducting polymer. MD simulations were carried out using force-field parametrizations explicitly developed for polaronic and bipolaronic models. QM/MM-MD calculations were performed using a quantum mechanical zone defined by four repeat units. In addition of the structural and electronic characteristics of the two forms of PAni-ES, MD and QM/MM-MD simulations indicate that the bipolaronic is the most stable state of amorphous PAni-ES. Complementary studies have been carried out using different experimental techniques. Although the morphology and topography of doped and undoped PAni are very similar, comparison of their UV–vis spectra supports the preference toward the bipolaronic form of PAni-ES.Thu, 30 Oct 2014 15:28:04 GMThttp://hdl.handle.net/2117/245162014-10-30T15:28:04ZCanales Gabriel, Manel; Torras Costa, Juan; Fabregat Jove, Georgina; Meneguzzi, Alvaro; Alemán Llansó, CarlosnoThe polaronic and bipolaronic forms of polyaniline emeraldine salt (PAni-ES) in the amorphous solid state have been simulated using classical molecular dynamics (MD) and hybrid quantum mechanical/molecular mechanical-molecular dynamics (QM/MM-MD) approaches. It should be remarked that the electronic state of PAni-ES has been theoretically investigated in the gas phase, solution phase, and crystalline state, but this is the first study in the amorphous solid state, which is the most typical for this conducting polymer. MD simulations were carried out using force-field parametrizations explicitly developed for polaronic and bipolaronic models. QM/MM-MD calculations were performed using a quantum mechanical zone defined by four repeat units. In addition of the structural and electronic characteristics of the two forms of PAni-ES, MD and QM/MM-MD simulations indicate that the bipolaronic is the most stable state of amorphous PAni-ES. Complementary studies have been carried out using different experimental techniques. Although the morphology and topography of doped and undoped PAni are very similar, comparison of their UV–vis spectra supports the preference toward the bipolaronic form of PAni-ES.Mass and size dependence of single ion dynamics in molten monohalides
http://hdl.handle.net/2117/24493
Title: Mass and size dependence of single ion dynamics in molten monohalides
Authors: Alcaraz Sendra, Olga; Trullàs Simó, Joaquim
Abstract: This work is concerned with four molten monohalides with different ionic radii ratios (RbCl, NaI, AgCl, and CuCl) and ideal isotopic systems of these salts with different ionic mass ratios. The velocity autocorrelation functions of the two ionic species in each melt have been studied by both a theoretical approximation and molecular dynamics simulations. It is found that their main features may be qualitatively predicted by considering suitable combinations of the second and fourth frequency moments of their spectra. The analysis of these two parameters allows us to determine how the structure (strongly dependent on the ionic size difference) and the ionic masses contribute to the shape of the velocity autocorrelation functions. The results show that the averaged microscopic motion of the small ions is mainly determined by the first neighboring shell of unlike ions, whereas the nearest shell of like ions also affects the dynamics of the large ions. This effect is more pronounced as the size difference is greater. Furthermore, it is concluded that the size differences encourage the rattling motion of the large ions, whereas the mass difference encourages the backscattering and oscillations of the velocity autocorrelation function of the light ions. A simple rule is derived to determine the interplay between these two effects. Comparison between the mass and nearest distance ratios enables the prediction as to which species will experience a more pronounced backscatteringmotion. The size difference effects prevail in the hydrodynamics regime and the self-diffusion coefficient of the small ions is higher than that of the large ones. The difference between the self-diffusion coefficient increases as the size differences increases.Mon, 27 Oct 2014 18:54:09 GMThttp://hdl.handle.net/2117/244932014-10-27T18:54:09ZAlcaraz Sendra, Olga; Trullàs Simó, JoaquimnoBackscattering
Self diffusion
Hydrodynamics
Kinematics
molecular dynamicsThis work is concerned with four molten monohalides with different ionic radii ratios (RbCl, NaI, AgCl, and CuCl) and ideal isotopic systems of these salts with different ionic mass ratios. The velocity autocorrelation functions of the two ionic species in each melt have been studied by both a theoretical approximation and molecular dynamics simulations. It is found that their main features may be qualitatively predicted by considering suitable combinations of the second and fourth frequency moments of their spectra. The analysis of these two parameters allows us to determine how the structure (strongly dependent on the ionic size difference) and the ionic masses contribute to the shape of the velocity autocorrelation functions. The results show that the averaged microscopic motion of the small ions is mainly determined by the first neighboring shell of unlike ions, whereas the nearest shell of like ions also affects the dynamics of the large ions. This effect is more pronounced as the size difference is greater. Furthermore, it is concluded that the size differences encourage the rattling motion of the large ions, whereas the mass difference encourages the backscattering and oscillations of the velocity autocorrelation function of the light ions. A simple rule is derived to determine the interplay between these two effects. Comparison between the mass and nearest distance ratios enables the prediction as to which species will experience a more pronounced backscatteringmotion. The size difference effects prevail in the hydrodynamics regime and the self-diffusion coefficient of the small ions is higher than that of the large ones. The difference between the self-diffusion coefficient increases as the size differences increases.The bridge functions of molten salts
http://hdl.handle.net/2117/24492
Title: The bridge functions of molten salts
Authors: Tasseven, Çetin; Enrique González, Luis; Silbert, Moises; Alcaraz Sendra, Olga; Trullàs Simó, Joaquim
Abstract: The bridge functions of molten NaCl and AgI near melting were obtained by using the model potentials of Born–Huggins–Mayer for NaCl and Vashishta–Rahman for AgI. The calculations of the bridge functions involved molecular dynamics simulations, the extension of the procedure originally proposed by Poll et al. [Phys. Rev. A 37, 1672 (1988)] and the numerical solution of the Ornstein–Zernike equations. The calculated bridge functions do not conform with the universality ansatz. They also differ from the bridge functions obtained for model electrolytes. Following the results obtained for electrolytes and those of this work for molten salts, it is conjectured that the universality ansatz for the bridge functions does not apply for systems whose attractive interactions play a decisive role in their structural ordering.Mon, 27 Oct 2014 18:30:12 GMThttp://hdl.handle.net/2117/244922014-10-27T18:30:12ZTasseven, Çetin; Enrique González, Luis; Silbert, Moises; Alcaraz Sendra, Olga; Trullàs Simó, JoaquimnoElectrolytes
Molten salts
Molecular dynamics
Numerical solutionsThe bridge functions of molten NaCl and AgI near melting were obtained by using the model potentials of Born–Huggins–Mayer for NaCl and Vashishta–Rahman for AgI. The calculations of the bridge functions involved molecular dynamics simulations, the extension of the procedure originally proposed by Poll et al. [Phys. Rev. A 37, 1672 (1988)] and the numerical solution of the Ornstein–Zernike equations. The calculated bridge functions do not conform with the universality ansatz. They also differ from the bridge functions obtained for model electrolytes. Following the results obtained for electrolytes and those of this work for molten salts, it is conjectured that the universality ansatz for the bridge functions does not apply for systems whose attractive interactions play a decisive role in their structural ordering.Space-dependent self-diffusion processes in molten copper halides: a molecular dynamics study
http://hdl.handle.net/2117/24489
Title: Space-dependent self-diffusion processes in molten copper halides: a molecular dynamics study
Authors: Alcaraz Sendra, Olga; Trullàs Simó, Joaquim
Abstract: This work is concerned with single ion dynamics in molten copper halides (CuI and CuCl) which exhibit fast ionic conduction before melting. The self-dynamic structure factor of the two ionic species in each melt have been calculated by molecular dynamics simulations and the corresponding effective wavelength-dependent self-diffusion coefficients have been studied. The results have been compared with those obtained for molten alkali halides (KCl and RbCl).Mon, 27 Oct 2014 15:29:18 GMThttp://hdl.handle.net/2117/244892014-10-27T15:29:18ZAlcaraz Sendra, Olga; Trullàs Simó, JoaquimnoMolecular dynamics
Copper
Self diffusion
Ionic conductionThis work is concerned with single ion dynamics in molten copper halides (CuI and CuCl) which exhibit fast ionic conduction before melting. The self-dynamic structure factor of the two ionic species in each melt have been calculated by molecular dynamics simulations and the corresponding effective wavelength-dependent self-diffusion coefficients have been studied. The results have been compared with those obtained for molten alkali halides (KCl and RbCl).The force matching approach to multiscale simulations: merits, shortcomings, and future perspectives
http://hdl.handle.net/2117/23978
Title: The force matching approach to multiscale simulations: merits, shortcomings, and future perspectives
Authors: Masia, Marco; Guàrdia Manuel, Elvira; Nicolini, Paolo
Abstract: Among the various approaches to multiscale simulations, in recent years, force matching has been known for a quick growth. The method is based on a least-square fit of reference properties obtained from simulations at a certain scale, to parameterize the force field for coarser-grained scale simulations. Its advantage with respect to conventional schemes used for parameterizing force fields, lies in that only physically accessible configurations are sampled, and that the number of reference data per configuration is large. In this perspective article, we discuss some recent findings on the tailoring of the objective function, on the choice of the empirical potential, and on the way to improve the quality of the reference calculations. We present pros and cons of the algorithm, and we propose a road map to future developments. (C) 2014 Wiley Periodicals, Inc.Wed, 03 Sep 2014 15:14:58 GMThttp://hdl.handle.net/2117/239782014-09-03T15:14:58ZMasia, Marco; Guàrdia Manuel, Elvira; Nicolini, Paolonoforce matching, multiscale simulations, force field parameterization, POTENTIALS, FIELDS, PARAMETRIZATIONAmong the various approaches to multiscale simulations, in recent years, force matching has been known for a quick growth. The method is based on a least-square fit of reference properties obtained from simulations at a certain scale, to parameterize the force field for coarser-grained scale simulations. Its advantage with respect to conventional schemes used for parameterizing force fields, lies in that only physically accessible configurations are sampled, and that the number of reference data per configuration is large. In this perspective article, we discuss some recent findings on the tailoring of the objective function, on the choice of the empirical potential, and on the way to improve the quality of the reference calculations. We present pros and cons of the algorithm, and we propose a road map to future developments. (C) 2014 Wiley Periodicals, Inc.Atomic monolayer deposition on the surface of nanotube mechanical resonators
http://hdl.handle.net/2117/23221
Title: Atomic monolayer deposition on the surface of nanotube mechanical resonators
Authors: Tavernarakis, Alexandros; Chaste, Julien; Eichler, Alexander; Ceballo, Gustavo; Gordillo Bargueño, Maria Carmen; Boronat Medico, Jordi; Bachtold, Adrian
Abstract: We study monolayers of noble gas atoms (Xe, Kr, Ar, and Ne) deposited on individual ultraclean suspended nanotubes. For this, we record the resonance frequency of the mechanical motion of the nanotube, since it provides a direct measure of the coverage. The latter is the number of adsorbed atoms divided by the number of the carbon atoms of the suspended nanotube. Monolayers form when the temperature is lowered in a constant pressure of noble gas atoms. The coverage of Xe monolayers remains constant at 1/6 over a large temperature range. This finding reveals that Xe monolayers are solid phases with a triangular atomic arrangement, and are commensurate with the underlying carbon nanotube. By comparing our measurements to theoretical calculations, we identify the phases of Ar and Ne monolayers as fluids, and we tentatively describe Kr monolayers as solid phases. These results underscore that mechanical resonators made from single nanotubes are excellent probes for surface science.Fri, 13 Jun 2014 17:24:45 GMThttp://hdl.handle.net/2117/232212014-06-13T17:24:45ZTavernarakis, Alexandros; Chaste, Julien; Eichler, Alexander; Ceballo, Gustavo; Gordillo Bargueño, Maria Carmen; Boronat Medico, Jordi; Bachtold, AdriannoWe study monolayers of noble gas atoms (Xe, Kr, Ar, and Ne) deposited on individual ultraclean suspended nanotubes. For this, we record the resonance frequency of the mechanical motion of the nanotube, since it provides a direct measure of the coverage. The latter is the number of adsorbed atoms divided by the number of the carbon atoms of the suspended nanotube. Monolayers form when the temperature is lowered in a constant pressure of noble gas atoms. The coverage of Xe monolayers remains constant at 1/6 over a large temperature range. This finding reveals that Xe monolayers are solid phases with a triangular atomic arrangement, and are commensurate with the underlying carbon nanotube. By comparing our measurements to theoretical calculations, we identify the phases of Ar and Ne monolayers as fluids, and we tentatively describe Kr monolayers as solid phases. These results underscore that mechanical resonators made from single nanotubes are excellent probes for surface science.We study monolayers of noble gas atoms (Xe, Kr, Ar, and Ne) deposited on individual ultraclean suspended nanotubes. For this, we record the resonance frequency of the mechanical motion of the nanotube, since it provides a direct measure of the coverage. The latter is the number of adsorbed atoms divided by the number of the carbon atoms of the suspended nanotube. Monolayers form when the temperature is lowered in a constant pressure of noble gas atoms. The coverage of Xe monolayers remains constant at 1/6 over a large temperature range. This finding reveals that Xe monolayers are solid phases with a triangular atomic arrangement, and are commensurate with the underlying carbon nanotube. By comparing our measurements to theoretical calculations, we identify the phases of Ar and Ne monolayers as fluids, and we tentatively describe Kr monolayers as solid phases. These results underscore that mechanical resonators made from single nanotubes are excellent probes for surface science.Dynamical aspects of intermolecular proton transfer in liquid water and low-density amorphous ices
http://hdl.handle.net/2117/23020
Title: Dynamical aspects of intermolecular proton transfer in liquid water and low-density amorphous ices
Authors: Tahat, Amani; Martí Rabassa, Jordi
Abstract: The microscopic dynamics of an excess proton in water and in low-density amorphous ices has been studied by means of a series of molecular dynamics simulations. Interaction of water with the proton species was modelled using a multistate empirical valence bond Hamiltonian model. The analysis of the effects of low temperatures on proton diffusion and transfer rates has been considered for a temperature range between 100 and 298 K at the constant density of 1 g cm -3 . We observed a marked slowdown of proton transfer rates at low temperatures, but some episodes are still seen at 100 K. In a similar fashion, mobility of the lone proton gets significantly reduced when temperature decreases below 273 K. The proton transfer in low-density amorphous ice is an activated process with energy barriers between 1–10 kJ/mol depending of the temperature range considered and eventually showing Arrhenius-like behavior. Spectroscopic data indicated the survival of both Zundel and Eigen structures along the whole temperature range, revealed by significant spectral frequency shifts.Tue, 20 May 2014 15:34:20 GMThttp://hdl.handle.net/2117/230202014-05-20T15:34:20ZTahat, Amani; Martí Rabassa, Jordinomicroscopic dynamics, excess proton in water, low-density amorphous ices, proton transferThe microscopic dynamics of an excess proton in water and in low-density amorphous ices has been studied by means of a series of molecular dynamics simulations. Interaction of water with the proton species was modelled using a multistate empirical valence bond Hamiltonian model. The analysis of the effects of low temperatures on proton diffusion and transfer rates has been considered for a temperature range between 100 and 298 K at the constant density of 1 g cm -3 . We observed a marked slowdown of proton transfer rates at low temperatures, but some episodes are still seen at 100 K. In a similar fashion, mobility of the lone proton gets significantly reduced when temperature decreases below 273 K. The proton transfer in low-density amorphous ice is an activated process with energy barriers between 1–10 kJ/mol depending of the temperature range considered and eventually showing Arrhenius-like behavior. Spectroscopic data indicated the survival of both Zundel and Eigen structures along the whole temperature range, revealed by significant spectral frequency shifts.Molecular dynamics simulation study of methanesulfonic acid
http://hdl.handle.net/2117/22794
Title: Molecular dynamics simulation study of methanesulfonic acid
Authors: Canales Gabriel, Manel; Alemán Llansó, Carlos
Abstract: A molecular dynamics simulation study of methanesulfonic acid has been carried out using a reliable force field in a large range of temperatures. Several thermodynamic, structural, and dynamical properties have been calculated and compared with the available experimental data. The density, the shear viscosity, the heat of vaporization, and the melting temperature results, calculated from this force field, are in a good agreement with the experimental data. Analysis of the influence of the hydrogen bonds in structural and dynamical properties has also been performed. The continuous and interrupted methodologies to compute hydrogen bonding lifetimes have been applied. The interrupted hydrogen bond lifetimes values are consistent with the diffusion and viscosity coefficients. The activation energies of the self-diffusion, the reorientational motions, and the hydrogen bonding lifetimes are coincident.Wed, 30 Apr 2014 15:55:44 GMThttp://hdl.handle.net/2117/227942014-04-30T15:55:44ZCanales Gabriel, Manel; Alemán Llansó, CarlosnoEngineering controlled terms: Activation energy, Computer simulation, Molecular dynamics, Organic acidsA molecular dynamics simulation study of methanesulfonic acid has been carried out using a reliable force field in a large range of temperatures. Several thermodynamic, structural, and dynamical properties have been calculated and compared with the available experimental data. The density, the shear viscosity, the heat of vaporization, and the melting temperature results, calculated from this force field, are in a good agreement with the experimental data. Analysis of the influence of the hydrogen bonds in structural and dynamical properties has also been performed. The continuous and interrupted methodologies to compute hydrogen bonding lifetimes have been applied. The interrupted hydrogen bond lifetimes values are consistent with the diffusion and viscosity coefficients. The activation energies of the self-diffusion, the reorientational motions, and the hydrogen bonding lifetimes are coincident.Spin-polarized hydrogen adsorbed on the surface of superfluid He-4
http://hdl.handle.net/2117/22078
Title: Spin-polarized hydrogen adsorbed on the surface of superfluid He-4
Authors: Marín, J. M.; Vranješ Markic, Leandra; Boronat Medico, Jordi
Abstract: The experimental realization of a thin layer of spin-polarized hydrogen H double down arrow adsorbed on top of the surface of superfluid He-4 provides one of the best examples of a stable, nearly two-dimensional(2D) quantum Bose gas. We report a theoretical study of this system using quantum Monte Carlo methods in the limit of zero temperature. Using the full Hamiltonian of the system, composed of a superfluid He-4 slab and the adsorbed H double down arrow layer, we calculate the main properties of its ground state using accurate models for the pair interatomic potentials. Comparing the results for the layer with the ones obtained for a strictly 2D setup, we analyze the departure from the 2D character when the density increases. Only when the coverage is rather small the use of a purely 2D model is justified. The condensate fraction of the layer is significantly larger than in 2D at the same surface density, being as large as 60% at the largest coverage studied. (c) 2013 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4843375]Fri, 14 Mar 2014 15:12:24 GMThttp://hdl.handle.net/2117/220782014-03-14T15:12:24ZMarín, J. M.; Vranješ Markic, Leandra; Boronat Medico, JordinoBOSE-EINSTEIN CONDENSATION, ATOMIC-HYDROGEN, LIQUID-HELIUM, ALIGNED HYDROGEN, QUANTUM SYSTEMS, GAS, ADSORPTION, SCATTERING, HEThe experimental realization of a thin layer of spin-polarized hydrogen H double down arrow adsorbed on top of the surface of superfluid He-4 provides one of the best examples of a stable, nearly two-dimensional(2D) quantum Bose gas. We report a theoretical study of this system using quantum Monte Carlo methods in the limit of zero temperature. Using the full Hamiltonian of the system, composed of a superfluid He-4 slab and the adsorbed H double down arrow layer, we calculate the main properties of its ground state using accurate models for the pair interatomic potentials. Comparing the results for the layer with the ones obtained for a strictly 2D setup, we analyze the departure from the 2D character when the density increases. Only when the coverage is rather small the use of a purely 2D model is justified. The condensate fraction of the layer is significantly larger than in 2D at the same surface density, being as large as 60% at the largest coverage studied. (c) 2013 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4843375]Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes
http://hdl.handle.net/2117/21998
Title: Diffusion and spectroscopy of water and lipids in fully hydrated dimyristoylphosphatidylcholine bilayer membranes
Authors: Yang, Jing; Calero Borrallo, Carles; Martí Rabassa, Jordi
Abstract: Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10-5 cm2/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10-8 cm2/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2–8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a “wagging of the tails” frequency around 30 cm-1, which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.Tue, 11 Mar 2014 15:30:00 GMThttp://hdl.handle.net/2117/219982014-03-11T15:30:00ZYang, Jing; Calero Borrallo, Carles; Martí Rabassa, Jordinolipid bilayer membrane, molecular dynamics, dynamics (diffusion, spectroscopy)Microscopic structure and dynamics of water and lipids in a fully hydrated dimyristoylphosphatidylcholine phospholipid lipid bilayer membrane in the liquid-crystalline phase have been analyzed with all-atom molecular dynamics simulations based on the recently parameterized CHARMM36 force field. The diffusive dynamics of the membrane lipids and of its hydration water, their reorientational motions as well as their corresponding spectral densities, related to the absorption of radiation, have been considered for the first time using the present force field. In addition, structural properties such as density and pressure profiles, a deuterium-order parameter, surface tension, and the extent of water penetration in the membrane have been analyzed. Molecular self-diffusion, reorientational motions, and spectral densities of atomic species reveal a variety of time scales playing a role in membrane dynamics. The mechanisms of lipid motion strongly depend on the time scale considered, from fast ballistic translation at the scale of picoseconds (effective diffusion coefficients of the order of 10-5 cm2/s) to diffusive flow of a few lipids forming nanodomains at the scale of hundreds of nanoseconds (diffusion coefficients of the order of 10-8 cm2/s). In the intermediate regime of sub-diffusion, collisions with nearest neighbors prevent the lipids to achieve full diffusion. Lipid reorientations along selected directions agree well with reported nuclear magnetic resonance data and indicate two different time scales, one about 1 ns and a second one in the range of 2–8 ns. We associated the two time scales of reorientational motions with angular distributions of selected vectors. Calculated spectral densities corresponding to lipid and water reveal an overall good qualitative agreement with Fourier transform infrared spectroscopy experiments. Our simulations indicate a blue-shift of the low frequency spectral bands of hydration water as a result of its interaction with lipids. We have thoroughly analyzed the physical meaning of all spectral features from lipid atomic sites and correlated them with experimental data. Our findings include a “wagging of the tails” frequency around 30 cm-1, which essentially corresponds to motions of the tail-group along the instantaneous plane formed by the two lipid tails, i.e., in-plane oscillations are clearly of bigger importance than those along the normal-to-the plane direction.Pattern recognition and data mining software based on artificial neural networks applied to proton transfer in aqueous environments
http://hdl.handle.net/2117/21794
Title: Pattern recognition and data mining software based on artificial neural networks applied to proton transfer in aqueous environments
Authors: Tahat, Amani; Martí Rabassa, Jordi; Khwaldeh, Ali; Tahat, Kaher
Abstract: In computational physics proton transfer phenomena could be viewed as pattern classification problems based on a set of input features allowing to classify the proton motion into two categories: transfer‘occurred’and transfer‘not occurred’. The goal of this paper is to evaluate the use of artificial neural networks in the classification of proton transfer events, based on the feed-forward back propagation neural network, used as a classifier to distinguish between the two transfer cases. In this paper, we use a new developed data mining and pattern recognition tool for automating, controlling, and drawing charts of the output data of an Empirical Valence Bond existing code. The study analyzes the need of pattern recognition in aqueous proton transfer processes and how the learning approach in error back propagation (multilayer perceptron algorithms) could be satisfactorily employed in the present case. We present a tool for pattern recognition and validate the code including a real physical case study. The results of applying the artificial neural networks methodology to crowd patterns based upon selected physical properties (e.g., temperature, density) show the abilities of the network to learn proton transfer patterns corresponding to properties of the aqueous environments, which is in turn proved to be fully compatible with previous proton transfer studies.Thu, 27 Feb 2014 14:47:41 GMThttp://hdl.handle.net/2117/217942014-02-27T14:47:41ZTahat, Amani; Martí Rabassa, Jordi; Khwaldeh, Ali; Tahat, Kahernoproton transfer, pattern recognition, neural networks, water environments, chart pattern, data mining, artificial neural network, empirical valence bondIn computational physics proton transfer phenomena could be viewed as pattern classification problems based on a set of input features allowing to classify the proton motion into two categories: transfer‘occurred’and transfer‘not occurred’. The goal of this paper is to evaluate the use of artificial neural networks in the classification of proton transfer events, based on the feed-forward back propagation neural network, used as a classifier to distinguish between the two transfer cases. In this paper, we use a new developed data mining and pattern recognition tool for automating, controlling, and drawing charts of the output data of an Empirical Valence Bond existing code. The study analyzes the need of pattern recognition in aqueous proton transfer processes and how the learning approach in error back propagation (multilayer perceptron algorithms) could be satisfactorily employed in the present case. We present a tool for pattern recognition and validate the code including a real physical case study. The results of applying the artificial neural networks methodology to crowd patterns based upon selected physical properties (e.g., temperature, density) show the abilities of the network to learn proton transfer patterns corresponding to properties of the aqueous environments, which is in turn proved to be fully compatible with previous proton transfer studies.Possible superfluidity of molecular hydrogen in a two-dimensional crystal phase of sodium
http://hdl.handle.net/2117/21500
Title: Possible superfluidity of molecular hydrogen in a two-dimensional crystal phase of sodium
Authors: Cazorla Silva, Claudio; Boronat Medico, Jordi
Abstract: We theoretically investigate the ground-state properties of a molecular para-hydrogen (p-H 2 ) film in which crystallization is energetically frustrated by embedding sodium (Na) atoms periodically distributed in a triangular lattice. In order to fully deal with the quantum nature of p-H 2 molecules, we employ the diffusion Monte Carlo method and realistic semiempirical pairwise potentials describing the interactions between H 2 -H 2 and Na-H 2 species. In particular, we calculate the energetic, structural, and superfluid properties of two-dimensional Na-H 2 systems within a narrow density interval around equilibrium at zero temperature. In contrast to previous computational studies considering other alkali metal species such as rubidium and potassium, we find that the p-H 2 ground state is a liquid with a significantly large superfluid fraction of ρ s /ρ=0.29(2) . The appearance of p-H 2 superfluid response is due to the fact that the interactions between Na atoms and H 2 molecules are less attractive than between H 2 molecules. This induces a considerable reduction of the hydrogen density which favors the stabilization of the liquid phase.Mon, 10 Feb 2014 16:17:42 GMThttp://hdl.handle.net/2117/215002014-02-10T16:17:42ZCazorla Silva, Claudio; Boronat Medico, JordinoMolecular hydrogen, two-dimensional, crystal phase of sodiumWe theoretically investigate the ground-state properties of a molecular para-hydrogen (p-H 2 ) film in which crystallization is energetically frustrated by embedding sodium (Na) atoms periodically distributed in a triangular lattice. In order to fully deal with the quantum nature of p-H 2 molecules, we employ the diffusion Monte Carlo method and realistic semiempirical pairwise potentials describing the interactions between H 2 -H 2 and Na-H 2 species. In particular, we calculate the energetic, structural, and superfluid properties of two-dimensional Na-H 2 systems within a narrow density interval around equilibrium at zero temperature. In contrast to previous computational studies considering other alkali metal species such as rubidium and potassium, we find that the p-H 2 ground state is a liquid with a significantly large superfluid fraction of ρ s /ρ=0.29(2) . The appearance of p-H 2 superfluid response is due to the fact that the interactions between Na atoms and H 2 molecules are less attractive than between H 2 molecules. This induces a considerable reduction of the hydrogen density which favors the stabilization of the liquid phase.Shortcomings of the standard Lennard-Jones dispersion term in water models, studied with force matching
http://hdl.handle.net/2117/21290
Title: Shortcomings of the standard Lennard-Jones dispersion term in water models, studied with force matching
Authors: Nicolini, Paolo; Guàrdia Manuel, Elvira; Masia, Marco
Abstract: In this work, ab initio parametrization of water force field is used to get insights into the functional form of empirical potentials to properly model the physics underlying dispersion interactions. We exploited the force matching algorithm to fit the interaction forces obtained with dispersion corrected density functional theory based molecular dynamics simulations. We found that the standard Lennard-Jones interaction potentials poorly reproduce the attractive character of dispersion forces. This drawback can be resolved by accounting for the distinctive short range behavior of dispersion interactions, multiplying the r −6 term by a damping function. We propose two novel parametrizations of the force field using different damping functions. Structural and dynamical properties of the new models are computed and compared with the ones obtained from the non-damped force field, showing an improved agreement with reference first principle calculations.Mon, 20 Jan 2014 17:45:49 GMThttp://hdl.handle.net/2117/212902014-01-20T17:45:49ZNicolini, Paolo; Guàrdia Manuel, Elvira; Masia, MarconoDispersion interaction, Dispersion-corrected density functional, Dynamical properties, Empirical potentials, First principle calculations, Force-matching algorithm, Lennard-Jones interaction potential, Molecular dynamics simulationsIn this work, ab initio parametrization of water force field is used to get insights into the functional form of empirical potentials to properly model the physics underlying dispersion interactions. We exploited the force matching algorithm to fit the interaction forces obtained with dispersion corrected density functional theory based molecular dynamics simulations. We found that the standard Lennard-Jones interaction potentials poorly reproduce the attractive character of dispersion forces. This drawback can be resolved by accounting for the distinctive short range behavior of dispersion interactions, multiplying the r −6 term by a damping function. We propose two novel parametrizations of the force field using different damping functions. Structural and dynamical properties of the new models are computed and compared with the ones obtained from the non-damped force field, showing an improved agreement with reference first principle calculations.Zero-temperature phase diagram of D2 physisorbed on graphane
http://hdl.handle.net/2117/20842
Title: Zero-temperature phase diagram of D2 physisorbed on graphane
Authors: Carbonell Coronado, Carmen; De Soto Borrera, Feliciano Carlos; Cazorla Silva, Claudio; Boronat Medico, Jordi; Gordillo Bargueño, Maria Carmen
Abstract: We determined the zero-temperature phase diagram of D2 physisorbed on graphane using the diffusion Monte Carlo method. The substrate used was C-graphane, an allotropic form of the compound that has been experimentally obtained through hydrogenation of graphene. We found that the ground state is the δ phase, a commensurate structure observed experimentally when D2 is adsorbed on graphite, and not the registered structure characteristic of H2 on the same substrate.Thu, 28 Nov 2013 17:56:18 GMThttp://hdl.handle.net/2117/208422013-11-28T17:56:18ZCarbonell Coronado, Carmen; De Soto Borrera, Feliciano Carlos; Cazorla Silva, Claudio; Boronat Medico, Jordi; Gordillo Bargueño, Maria CarmennoAllotropic forms, Commensurate structure, Diffusion Monte Carlo method, Graphane, Physisorbed, Structure characteristic, Zero temperaturesWe determined the zero-temperature phase diagram of D2 physisorbed on graphane using the diffusion Monte Carlo method. The substrate used was C-graphane, an allotropic form of the compound that has been experimentally obtained through hydrogenation of graphene. We found that the ground state is the δ phase, a commensurate structure observed experimentally when D2 is adsorbed on graphite, and not the registered structure characteristic of H2 on the same substrate.Characterization of the methane–graphene hydrophobic interaction in aqueous solution from ab initio simulations
http://hdl.handle.net/2117/20817
Title: Characterization of the methane–graphene hydrophobic interaction in aqueous solution from ab initio simulations
Authors: Calero Borrallo, Carles; Martí Rabassa, Jordi; Guàrdia Manuel, Elvira; Masia, Marco
Abstract: In this article, the interaction between a methane molecule and a graphene plane in liquid water has been characterized employing DFT-based free energy Molecular Dynamics calculations. This system represents a good model to understand the generic interaction between a small hydrophobic solute (methane molecule) and an extense hydrophobic surface (graphene plane). The structural and dynamical properties of graphene and methane hydration water are analyzed and found to be closely related to the main features of the potential of mean force. The results could be used in coarse-grained models to take into account the effect of the hydrophobic interaction in realistic systems relevant to experiment.Wed, 27 Nov 2013 14:39:00 GMThttp://hdl.handle.net/2117/208172013-11-27T14:39:00ZCalero Borrallo, Carles; Martí Rabassa, Jordi; Guàrdia Manuel, Elvira; Masia, Marconomethane-graphene interactions, hydrophobic forces, ab-initio simulationIn this article, the interaction between a methane molecule and a graphene plane in liquid water has been characterized employing DFT-based free energy Molecular Dynamics calculations. This system represents a good model to understand the generic interaction between a small hydrophobic solute (methane molecule) and an extense hydrophobic surface (graphene plane). The structural and dynamical properties of graphene and methane hydration water are analyzed and found to be closely related to the main features of the potential of mean force. The results could be used in coarse-grained models to take into account the effect of the hydrophobic interaction in realistic systems relevant to experiment.