Source author record

Jadran Vrabec

Jadran Vrabec appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2015arXiv

Long range correction for multi-site Lennard-Jones models and planar interfaces

A slab based long range correction approach for multi-site Lennard-Jones models is presented for systems with a planar film geometry that is based on the work by Janecek, J. Phys. Chem. B 110: 6264 (2006). It is efficient because it relies on a center-of-mass cutoff scheme and scales in terms of numerics almost perfectly with the molecule number. For validation, a series of simulations with the two-center Lennard-Jones model fluid, carbon dioxide and cyclohexane is carried out. The results of the present approach, a site-based long range correction and simulations without any long range correction are compared with respect to the saturated liquid density and the surface tension. The present simulation results exhibit only a weak dependence on the cutoff radius, indicating a high accuracy of the implemented long range correction.

preprint2015arXiv

ms2: A molecular simulation tool for thermodynamic properties, new version release

A new version release (2.0) of the molecular simulation tool ms2 [S. Deublein et al., Comput. Phys. Commun. 182 (2011) 2350] is presented. Version 2.0 of ms2 features a hybrid parallelization based on MPI and OpenMP for molecular dynamics simulation to achieve higher scalability. Furthermore, the formalism by Lustig [R. Lustig, Mol. Phys. 110 (2012) 3041] is implemented, allowing for a systematic sampling of Massieu potential derivatives in a single simulation run. Moreover, the Green-Kubo formalism is extended for the sampling of the electric conductivity and the residence time. To remove the restriction of the preceding version to electro-neutral molecules, Ewald summation is implemented to consider ionic long range interactions. Finally, the sampling of the radial distribution function is added.

preprint2014arXiv

ls1 mardyn: The massively parallel molecular dynamics code for large systems

The molecular dynamics simulation code ls1 mardyn is presented. It is a highly scalable code, optimized for massively parallel execution on supercomputing architectures, and currently holds the world record for the largest molecular simulation with over four trillion particles. It enables the application of pair potentials to length and time scales which were previously out of scope for molecular dynamics simulation. With an efficient dynamic load balancing scheme, it delivers high scalability even for challenging heterogeneous configurations. Presently, multi-center rigid potential models based on Lennard-Jones sites, point charges and higher-order polarities are supported. Due to its modular design, ls1 mardyn can be extended to new physical models, methods, and algorithms, allowing future users to tailor it to suit their respective needs. Possible applications include scenarios with complex geometries, e.g. for fluids at interfaces, as well as non-equilibrium molecular dynamics simulation of heat and mass transfer.

preprint2013arXiv

Computational molecular engineering as an emerging technology in process engineering

The present level of development of molecular force field methods is assessed from the point of view of simulation-based engineering, outlining the immediate perspective for further development and highlighting the newly emerging discipline of Computational Molecular Engineering (CME) which makes basic research in soft matter physics fruitful for industrial applications. Within the coming decade, major breakthroughs can be reached if a research focus is placed on processes at interfaces, combining aspects where an increase in the accessible length and time scales due to massively parallel high-performance computing will lead to particularly significant improvements.

preprint2013arXiv

Molecular modelling and simulation of electrolyte solutions, biomolecules, and wetting of component surfaces

Massively-parallel molecular dynamics simulation is applied to systems containing electrolytes, vapour-liquid interfaces, and biomolecules in contact with water-oil interfaces. Novel molecular models of alkali halide salts are presented and employed for the simulation of electrolytes in aqueous solution. The enzymatically catalysed hydroxylation of oleic acid is investigated by molecular dynamics simulation taking the internal degrees of freedom of the macromolecules into account. Thereby, Ewald summation methods are used to compute the long range electrostatic interactions. In systems with a phase boundary, the dispersive interaction, which is modelled by the Lennard-Jones potential here, has a more significant long range contribution than in homogeneous systems. This effect is accounted for by implementing the Janecek cutoff correction scheme. On this basis, the HPC infrastructure at the Steinbuch Centre for Computing was accessed and efficiently used, yielding new insights on the molecular systems under consideration.

preprint2011arXiv

Atomistic simulations of electrolyte solutions and hydrogels with explicit solvent models

Two of the most challenging tasks in molecular simulation consist in capturing the properties of systems with long-range interactions (e.g. electrolyte solutions) as well as systems containing large molecules such as hydrogels. For the development and optimization of molecular force fields and models, a large number of simulation runs have to be evaluated to obtain the sensitivity of the target properties with respect to the model parameters. The present work discusses force field development for electrolytes regarding thermodynamic properties of their aqueous solutions. Furthermore, simulations are conducted for the volume transition of hydrogels in the presence of electrolytes. It is shown that the properties of these complex systems can be captured by molecular simulation.

preprint2011arXiv

Static and dynamic properties of curved vapour-liquid interfaces by massively parallel molecular dynamics simulation

Curved fluid interfaces are investigated on the nanometre length scale by molecular dynamics simulation. Thereby, droplets surrounded by a metastable vapour phase are stabilized in the canonical ensemble. Analogous simulations are conducted for cylindrical menisci separating vapour and liquid phases under confinement in planar nanopores. Regarding the emergence of nanodroplets during nucleation, a non-equilibrium phenomenon, both the non-steady dynamics of condensation processes and stationary quantities related to supersaturated vapours are considered. Results for the truncated and shifted Lennard-Jones fluid and for mixtures of quadrupolar fluids confirm the applicability of the capillarity approximation and the classical nucleation theory.

preprint2011arXiv

The excess equimolar radius of liquid drops

The curvature dependence of the surface tension is related to the excess equimolar radius of liquid drops, i.e., the deviation of the equimolar radius from that defined with the macroscopic capillarity approximation. Based on the Tolman [J. Chem. Phys. 17, 333 (1949)] approach and its interpretation by Nijmeijer et al. [J. Chem. Phys. 96, 565 (1991)], the surface tension of spherical interfaces is analysed in terms of the pressure difference due to curvature. In the present study, the excess equimolar radius, which can be obtained directly from the density profile, is used instead of the Tolman length. Liquid drops of the truncated-shifted Lennard-Jones fluid are investigated by molecular dynamics simulation in the canonical ensemble, with equimolar radii ranging from 4 to 33 times the Lennard-Jones size parameter sigma. In these simulations, the magnitudes of the excess equimolar radius and the Tolman length are shown to be smaller than sigma/2. Other methodical approaches, from which mutually contradicting findings have been reported, are critically discussed, outlining possible sources of inaccuracy.

preprint2010arXiv

Contact angle dependence on the fluid-wall dispersive energy

Vapor-liquid menisci of the truncated and shifted Lennard-Jones fluid between parallel planar walls are investigated by molecular dynamics simulation. Thereby, the characteristic energy of the unlike dispersive interaction between fluid molecules and wall atoms is systematically varied to determine its influence on the contact angle. The temperature is varied as well, covering most of the range between the triple point temperature and the critical temperature of the bulk fluid. The transition between obtuse and acute angles is found to occur at a temperature-independent magnitude of the fluid-wall dispersive interaction energy. On the basis of the present simulation results, fluid-wall interaction potentials can be adjusted to contact angle measurements.

preprint2010arXiv

Development of models and methods for the molecular simulation of large systems and molecules

The most important factor for quantitative results in molecular dynamics simulation are well developed force fields and models. In the present work, the development of new models and the usage of force fields from the literature in large systems are presented. Both tasks lead to time consuming simulations that require massively parallel high performance computing. In the present work, new models for carbon dioxide and cyclohexanolare discussed and a new method for the model development is introduced. Force fields and models for the simulation of PNIPAAm hydrogel in pure water and sodium chloride solution are tested and verified and applied to the simulation of nucleation processes.

preprint2010arXiv

Molecular modeling of hydrogen bonding fluids: Vapor-liquid coexistence and interfacial properties

A major challenge for molecular modeling consists in optimizing the unlike interaction potentials. In many cases, combination rules are generally suboptimal when accurate predictions of properties like the mixture vapor pressure are needed. However, the well known Lorentz-Berthelot rule performs quite well and can be used as a starting point. If more accurate results are required, it is advisable to adjust the dispersive interaction energy parameter. In the present study, mixture properties are explored for binary systems containing hydrogen bonding components. Furthermore, vapor-liquid interface cluster criteria and contact angles are discussed and remarks on computational details are given. Finally, a sterically accurate generic model for benzyl alcohol is introduced and evaluated.

preprint2010arXiv

Molecular simulation of fluid dynamics on the nanoscale

Molecular dynamics simulation is applied to Poiseuille flow of liquid methane in planar graphite channels, covering channel diameters between 3 and 135 nm. On this length scale, a transition is found between the regime where local ordering induced by the wall dominates the entire system and larger channel diameters where the influence of boundary slip is still present, but of a more limited extent. The validity of Darcy's law for pressure-driven flow through porous media is not affected by the transition between these regimes.

preprint2010arXiv

Research on the behavior of liquid fluids atop superhydrophobic gas-bubbled surfaces

Superhydrophobic surfaces play an important role in the development of new product coatings such as cars, but also in mechanical engineering, especially design of turbines and compressors. Thus a vital part of the design of these surfaces is the computational simulation of such with a special interest on variation of shape and size of minor pits grooved into plane surfaces. In the present work, the dependence of the contact angle on the fluid-wall dispersive energy is determined by molecular simulation and static as well as dynamic properties of unpolar fluids in contact with extremely rough surfaces are obtained.

preprint2010arXiv

The air pressure effect on the homogeneous nucleation of carbon dioxide by molecular simulation

Vapour-liquid equilibria (VLE) and the influence of an inert carrier gas on homogeneous vapour to liquid nucleation are investigated by molecular simulation for quaternary mixtures of carbon dioxide, nitrogen, oxygen, and argon. Canonical ensemble molecular dynamics simulation using the Yasuoka-Matsumoto method is applied to nucleation in supersaturated vapours that contain more carbon dioxide than in the saturated state at the dew line. Established molecular models are employed that are known to accurately reproduce the VLE of the pure fluids as well as their binary and ternary mixtures. On the basis of these models, also the quaternary VLE properties of the bulk fluid are determined with the Grand Equilibrium method. Simulation results for the carrier gas influence on the nucleation rate are compared with the classical nucleation theory (CNT) considering the "pressure effect" [Phys. Rev. Lett. 101: 125703 (2008)]. It is found that the presence of air as a carrier gas decreases the nucleation rate only slightly and, in particular, to a significantly lower extent than predicted by CNT. The nucleation rate of carbon dioxide is generally underestimated by CNT, leading to a deviation between one and two orders of magnitude for pure carbon dioxide in the vicinity of the spinodal line and up to three orders of magnitude in presence of air as a carrier gas. Furthermore, CNT predicts a temperature dependence of the nucleation rate in the spinodal limit, which cannot be confirmed by molecular simulation.

preprint2009arXiv

Grand canonical steady-state simulation of nucleation

Grand canonical molecular dynamics (GCMD) is applied to the nucleation process in a metastable phase near the spinodal, where nucleation occurs almost instantaneously and is limited to a very short time interval. With a variant of Maxwell's demon, proposed by McDonald [Am. J. Phys. 31: 31 (1963)], all nuclei exceeding a specified size are removed. In such a steady-state simulation, the nucleation process is sampled over an arbitrary timespan and all properties of the metastable state, including the nucleation rate, can be obtained with an increased precision. As an example, a series of GCMD simulations with McDonald's demon is carried out for homogeneous vapor to liquid nucleation of the truncated-shifted Lennard-Jones (tsLJ) fluid, covering the entire relevant temperature range. The results are in agreement with direct non-equilibrium MD simulation in the canonical ensemble. It is confirmed for supersaturated vapors of the tsLJ fluid that the classical nucleation theory underpredicts the nucleation rate by two orders of magnitude.