Universität Stuttgart
Permanent URI for this communityhttps://elib.uni-stuttgart.de/handle/11682/1
Browse
40 results
Search Results
Item Open Access Efficient algorithms for electrostatic interactions including dielectric contrasts(2013) Arnold, Axel; Breitsprecher, Konrad; Fahrenberger, Florian; Kesselheim, Stefan; Lenz, Olaf; Holm, ChristianCoarse grained models of soft matter are usually combined with implicit solvent models that take the electrostatic polarizability into account via a dielectric background. In biophysical or nanoscale simulations that include water, this constant can vary greatly within the system. Performing molecular dynamics or other simulations that need compute exact electrostatic interactions between charges in those systems is computationally demanding. We review here several algorithms developped by us that perform exactly this task. For planar dielectric surfaces in partial periodic boundary conditions, the arising image charges can be either treated with the MMM2D algorithm in a very efficient and accurate way, or with the ELC term that enables the user to use his favorite 3D periodic Coulomb solver . Arbitrarily shaped interfaces can be dealt with using induced surface charges with the ICC algorithm. Finally, the local electrostatics algorithm MEMD (Maxwell Equations Molecular Dynamics) allows even to employ a smoothly varying dielectric constant in the systems. We introduce the concepts of these three algorithms, and an extension for the inclusion of boundaries that are to be held fixed at constant potential (metal conditions). For each method, we present a showcase application to highlight the importance of dielectric interfaces.Item Open Access A brief review of capillary number and its use in capillary desaturation curves(2022) Guo, Hu; Song, Kaoping; Hilfer, R.Capillary number, understood as the ratio of viscous force to capillary force, is one of the most important parameters in enhanced oil recovery (EOR). It continues to attract the interest of scientists and engineers, because the nature and quantification of macroscopic capillary forces remain controversial. At least 41 different capillary numbers have been collected here from the literature. The ratio of viscous and capillary force enters crucially into capillary desaturation experiments. Although the ratio is length scale dependent, not all definitions of capillary number depend on length scale, indicating potential inconsistencies between various applications and publications. Recently, new numbers have appeared and the subject continues to be actively discussed. Therefore, a short review seems appropriate and pertinent.Item Open Access Machine learning-driven investigation of the structure and dynamics of the BMIM-BF4 room temperature ionic liquid(2024) Zills, Fabian; Schäfer, Moritz René; Tovey, Samuel; Kästner, Johannes; Holm, ChristianRoom-temperature ionic liquids are an exciting group of materials with the potential to revolutionize energy storage. Due to their chemical structure and means of interaction, they are challenging to study computationally. Classical descriptions of their inter- and intra-molecular interactions require time intensive parametrization of force-fields which is prone to assumptions. While ab initio molecular dynamics approaches can capture all necessary interactions, they are too slow to achieve the time and length scales required. In this work, we take a step towards addressing these challenges by applying state-of-the-art machine-learned potentials to the simulation of 1-butyl-3-methylimidazolium tetrafluoroborate. We demonstrate a learning-on-the-fly procedure to train machine-learned potentials from single-point density functional theory calculations before performing production molecular dynamics simulations. Obtained structural and dynamical properties are in good agreement with computational and experimental references. Furthermore, our results show that hybrid machine-learned potentials can contribute to an improved prediction accuracy by mitigating the inherent shortsightedness of the models. Given that room-temperature ionic liquids necessitate long simulations to address their slow dynamics, achieving an optimal balance between accuracy and computational cost becomes imperative. To facilitate further investigation of these materials, we have made our IPSuite-based training and simulation workflow publicly accessible, enabling easy replication or adaptation to similar systems.Item Open Access Percolativity of porous media(2022) Hilfer, Rudolf; Hauskrecht, J.Connectivity and connectedness are nonadditive geometric functionals on the set of pore scale structures. They determine transport of mass, volume or momentum in porous media, because without connectivity there cannot be transport. Percolativity of porous media is introduced here as a geometric descriptor of connectivity, that can be computed from the pore scale and persists to the macroscale through a suitable upscaling limit. It is a measure that combines local percolation probabilities with a probability density of ratios of eigenvalues of the tensor of local percolating directions. Percolativity enters directly into generalized effective medium approximations. Predictions from these generalized effective medium approximations are found to be compatible with apparently anisotropic Archie correlations observed in experiment.Item Open Access On extremal domains and codomains for convolution of distributions and fractional calculus(2022) Kleiner, T.; Hilfer, R.It is proved that the class of c-closed distribution spaces contains extremal domains and codomains to make convolution of distributions a well-defined bilinear mapping. The distribution spaces are systematically endowed with topologies and bornologies that make convolution hypocontinuous whenever defined. Largest modules and smallest algebras for convolution semigroups are constructed along the same lines. The fact that extremal domains and codomains for convolution exist within this class of spaces is fundamentally related to quantale theory. The quantale theoretic residual formed from two c-closed spaces is characterized as the largest c-closed subspace of the corresponding space of convolutors. The theory is applied to obtain maximal distributional domains for fractional integrals and derivatives, for fractional Laplacians, Riesz potentials and for the Hilbert transform. Further, maximal joint domains for families of these operators are obtained such that their composition laws are preserved.Item Open Access Water structuring induces nonuniversal hydration repulsion between polar surfaces : quantitative comparison between molecular simulations, theory, and experiments(2024) Schlaich, Alexander; Daldrop, Jan O.; Kowalik, Bartosz; Kanduč, Matej; Schneck, Emanuel; Netz, Roland R.Polar surfaces in water typically repel each other at close separations, even if they are charge-neutral. This so-called hydration repulsion balances the van der Waals attraction and gives rise to a stable nanometric water layer between the polar surfaces. The resulting hydration water layer is crucial for the properties of concentrated suspensions of lipid membranes and hydrophilic particles in biology and technology, but its origin is unclear. It has been suggested that surface-induced molecular water structuring is responsible for the hydration repulsion, but a quantitative proof of this water-structuring hypothesis is missing. To gain an understanding of the mechanism causing hydration repulsion, we perform molecular simulations of different planar polar surfaces in water. Our simulated hydration forces between phospholipid bilayers agree perfectly with experiments, validating the simulation model and methods. For the comparison with theory, it is important to split the simulated total surface interaction force into a direct contribution from surface-surface molecular interactions and an indirect water-mediated contribution. We find the indirect hydration force and the structural water-ordering profiles from the simulations to be in perfect agreement with the predictions from theoretical models that account for the surface-induced water ordering, which strongly supports the water-structuring hypothesis for the hydration force. However, the comparison between the simulations for polar surfaces with different headgroup architectures reveals significantly different decay lengths of the indirect water-mediated hydration-force, which for laterally homogeneous water structuring would imply different bulk-water properties. We conclude that laterally inhomogeneous water ordering, induced by laterally inhomogeneous surface structures, shapes the hydration repulsion between polar surfaces in a decisive manner. Thus, the indirect water-mediated part of the hydration repulsion is caused by surface-induced water structuring but is surface-specific and thus nonuniversal.Item Open Access PDADMAC/PSS oligoelectrolyte multilayers : internal structure and hydration properties at early growth stages from atomistic simulations(2020) Sánchez, Pedro A.; Vögele, Martin; Smiatek, Jens; Qiao, Baofu; Sega, Marcello; Holm, ChristianWe analyze the internal structure and hydration properties of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) oligoelectrolyte multilayers at early stages of their layer-by-layer growth process. Our study is based on large-scale molecular dynamics simulations with atomistic resolution that we presented recently [Sánchez et al., Soft Matter 2019, 15, 9437], in which we produced the first four deposition cycles of a multilayer obtained by alternate exposure of a flat silica substrate to aqueous electrolyte solutions of such polymers at 0.1M of NaCl. In contrast to any previous work, here we perform a local structural analysis that allows us to determine the dependence of the multilayer properties on the distance to the substrate. We prove that the large accumulation of water and ions next to the substrate observed in previous overall measurements actually decreases the degree of intrinsic charge compensation, but this remains as the main mechanism within the interface region. We show that the range of influence of the substrate reaches approximately 3 nm, whereas the structure of the outer region is rather independent from the position. This detailed characterization is essential for the development of accurate mesoscale models able to reach length and time scales of technological interest.Item Open Access Permeability estimation of regular porous structures : a benchmark for comparison of methods(2021) Wagner, Arndt; Eggenweiler, Elissa; Weinhardt, Felix; Trivedi, Zubin; Krach, David; Lohrmann, Christoph; Jain, Kartik; Karadimitriou, Nikolaos; Bringedal, Carina; Voland, Paul; Holm, Christian; Class, Holger; Steeb, Holger; Rybak, IrynaThe intrinsic permeability is a crucial parameter to characterise and quantify fluid flow through porous media. However, this parameter is typically uncertain, even if the geometry of the pore structure is available. In this paper, we perform a comparative study of experimental, semi-analytical and numerical methods to calculate the permeability of a regular porous structure. In particular, we use the Kozeny-Carman relation, different homogenisation approaches (3D, 2D, very thin porous media and pseudo 2D/3D), pore-scale simulations (lattice Boltzmann method, Smoothed Particle Hydrodynamics and finite-element method) and pore-scale experiments (microfluidics). A conceptual design of a periodic porous structure with regularly positioned solid cylinders is set up as a benchmark problem and treated with all considered methods. The results are discussed with regard to the individual strengths and limitations of the used methods. The applicable homogenisation approaches as well as all considered pore-scale models prove their ability to predict the permeability of the benchmark problem. The underestimation obtained by the microfluidic experiments is analysed in detail using the lattice Boltzmann method, which makes it possible to quantify the influence of experimental setup restrictions.Item Open Access Insights into Hildebrand solubility parameters : contributions from cohesive energies or electrophilicity densities?(2023) Miranda‐Quintana, Ramón Alain; Chen, Lexin; Smiatek, JensWe introduce certain concepts and expressions from conceptual density functional theory (DFT) to study the properties of the Hildebrand solubility parameter. The original form of the Hildebrand solubility parameter is used to qualitatively estimate solubilities for various apolar and aprotic substances and solvents and is based on the square root of the cohesive energy density. Our results show that a revised expression allows the replacement of cohesive energy densities by electrophilicity densities, which are numerically accessible by simple DFT calculations. As an extension, the reformulated expression provides a deeper interpretation of the main contributions and, in particular, emphasizes the importance of charge transfer mechanisms. All calculated values of the Hildebrand parameters for a large number of common solvents are compared with experimental values and show good agreement for non‐ or moderately polar aprotic solvents in agreement with the original formulation of the Hildebrand solubility parameters. The observed deviations for more polar and protic solvents define robust limits from the original formulation which remain valid. Likewise, we show that the use of machine learning methods leads to only slightly better predictability.Item Open Access Theoretical insights into specific ion effects and strong‐weak acid‐base rules for ions in solution : deriving the law of matching solvent affinities from first principles(2020) Miranda‐Quintana, Ramón Alain; Smiatek, JensWe present a detailed study of specific ion effects, volcano plots and the law of matching solvent affinities by means of a conceptual density functional theory (DFT) approach. Our results highlight that specific ion effects and the corresponding implications on the solvation energy are mainly due to differences in the electric chemical potentials and chemical hardnesses of the ions and the solvent. Our approach can be further used to identify reliable criteria for the validity of the law of matching solvent affinities. Basic expressions are derived, which allow us to study the limiting conditions for this empirical observation with regard to matching chemical reactivity indices. Moreover, we show that chaotropic and kosmotropic concepts and their implications for the stability of ion pairs are directly related to a generalized strong and weak acids and bases (SWAB) principle for ions in solution, which is also applicable to rationalize the shape of volcano plots for different solvents. In contrast to previous assumptions, all empirical findings can be explained by the properties of local solvent‐ion complexes which dominate the specific global behavior of ion pairs in solution.