BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Marcello Sega

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

ResearcherUnclaimed source record
cond-mat.softBiological Physicsphysics.flu-dyncond-mat.mes-hallphysics.chem-phphysics.comp-phcond-mat.othercond-mat.stat-mechDistributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

15works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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)

preprint2022arXiv

On the accuracy and performance of the lattice Boltzmann method with 64-bit, 32-bit and novel 16-bit number formats

Fluid dynamics simulations with the lattice Boltzmann method (LBM) are very memory-intensive. Alongside reduction in memory footprint, significant performance benefits can be achieved by using FP32 (single) precision compared to FP64 (double) precision, especially on GPUs. Here, we evaluate the possibility to use even FP16 and Posit16 (half) precision for storing fluid populations, while still carrying arithmetic operations in FP32. For this, we first show that the commonly occurring number range in the LBM is a lot smaller than the FP16 number range. Based on this observation, we develop novel 16-bit formats - based on a modified IEEE-754 and on a modified Posit standard - that are specifically tailored to the needs of the LBM. We then carry out an in-depth characterization of LBM accuracy for six different test systems with increasing complexity: Poiseuille flow, Taylor-Green vortices, Karman vortex streets, lid-driven cavity, a microcapsule in shear flow (utilizing the immersed-boundary method) and finally the impact of a raindrop (based on a Volume-of-Fluid approach). We find that the difference in accuracy between FP64 and FP32 is negligible in almost all cases, and that for a large number of cases even 16-bit is sufficient. Finally, we provide a detailed performance analysis of all precision levels on a large number of hardware microarchitectures and show that significant speedup is achieved with mixed FP32/16-bit.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Structure and scaling laws of liquid/vapor interfaces close to the critical point

To reach a deeper understanding of fluid interfaces it is necessary to identify a meaningful coarse-graining length that separates intrinsic fluctuations from capillary ones, given the lack of a proper statistical mechanical definition of the latter. Here, with the help of unsupervised learning techniques, we introduce a new length scale based on the local density of the fluid. This length scale follows a scaling law that diverges more mildly than the bulk correlation length upon approaching the critical point. This allows to distinguish regimes of correlated and uncorrelated capillary waves from that of intrinsic fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capillary interactions between soft capsules protruding through thin fluid films

When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well studied, much is still to be understood about the behaviour of soft particles and the role of their softness during the final stages of film drying. Here, we use our recently-introduced numerical method that couples a fluid described using the lattice Boltzmann approach to a finite element description of deformable objects to investigate the drying process of a film with suspended soft particles. Our measured menisci deformations and lateral capillary forces, which agree well with previous theoretical and experimental works in case of rigid particles, show the deformations become smaller with increasing particles softness, resulting in weaker lateral interaction forces. At large interparticle distances, the force approaches that of rigid particles. Finally, we investigate the time dependent formation of particle clusters at the late stages of the film drying.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Long-range dispersion effects on the water/vapor interface simulated using the most common models

The long-range contribution to dispersion forces is known to have a major impact on the properties of inhomogeneous fluids, and its correct treatment is increasingly recognized as being a necessary requirement to avoid cutoff-related artefacts. Although analytical corrections for quantities like the surface tension are known, these can not take into account the structural changes induced by the long-range contributions. Here, we analyze the interfacial properties of seven popular water models, comparing the results with the cut-off version of the dispersion potential. The differences in surface tension estimates are in all cases found to be less than 2 mN/m.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The importance of chemical potential in the determination of water slip in nanochannels

We investigate the slip properties of water confined in graphite-like nano-channels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle -- but a different microscopic surface potential -- do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kinetic and ion pairing contributions in the dielectric spectra of electrolyte aqueous solutions

Understanding dielectric spectra can reveal important information about the dynamics of solvents and solutes from the dipolar relaxation times down to electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding salt ions to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applied field. So far, this effect has escaped an experimental verification, mainly because of the concomitant appearance of dielectric saturation from which the Hubbard-Onsager decrement cannot be easily separated. Here we develop a novel non-equilibrium molecular dynamics simulation approach to determine this decrement accurately for the first time. Using a thermodynamic consistent all-atom force field we show that for an aqueous solution containing sodium chloride around 4.8 Mol/l, this effect accounts for 12\% of the total dielectric permittivity. The dielectric decrement can be strikingly different if a less accurate force field for the ions is used. Using the widespread GROMOS parameters, we observe in fact an {\it increment} of the dielectric permittivity rather than a decrement. We can show that this increment is caused by ion pairing, introduced by a too low dispersion force, and clarify the microscopic connection between long-living ion pairs and the appearance of specific features in the dielectric spectrum of the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Kinetic dielectric decrement revisited: phenomenology of finite ion concentrations

With the help of a recently developed non-equilibrium approach, we investigate the ionic strength dependence of the Hubbard--Onsager dielectric decrement. We compute the depolarization of water molecules caused by the motion of ions in sodium chloride solutions from the dilute regime (0.035 M) up close to the saturation concentration (4.24 M), and find that the kinetic decrement displays a strong nonmonotonic behavior, in contrast to the prediction of available models. We introduce a phenomenological modification of the Hubbard--Onsager continuum theory, that takes into account the screening due to the ionic cloud at mean field level, and, is able to describe the kinetic decrement at high concentrations including the presence of a pronounced minimum.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mesoscale Structures at Complex Fluid-Fluid Interfaces: a Novel Lattice Boltzmann / Molecular Dynamics Coupling

Complex fluid-fluid interfaces featuring mesoscale structures with adsorbed particles are key components of newly designed materials which are continuously enriching the field of soft matter. Simulation tools which are able to cope with the different scales characterizing these systems are fundamental requirements for efficient theoretical investigations. In this paper we present a novel simulation method, based on the approach of Ahlrichs and Dünweg [Ahlrichs and Dünweg, Int. J. Mod. Phys. C, 1998, 9, 1429], that couples the "Shan-Chen" multicomponent Lattice Boltzmann technique to off-lattice molecular dynamics to simulate efficiently complex fluid-fluid interfaces. We demonstrate how this approach can be used to study a wide class of challenging problems. Several examples are given, with an accent on bicontinuous phases formation in polyelectrolyte solutions and ferrofluid emulsions. We also show that the introduction of solvation free energies in the particle-fluid interaction unveils the hidden, multiscale nature of the particle-fluid coupling, allowing to treat symmetrically (and interchangeably) the on-lattice and off-lattice components of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Regularization of the slip length divergence in water nanoflows by inhomogeneities at the Angstrom scale

We performed non-equilibrium Molecular Dynamics simulations of water flow in nano-channels with the aim of discriminating {\it static} from {\it dynamic} contributions of the solid surface to the slip length of the molecular flow. We show that the regularization of the slip length divergence at high shear rates, formerly attributed to the wall dynamics, is controlled instead by its static properties. Surprisingly, we find that atomic displacements at the Angstrom scale are sufficient to remove the divergence of the slip length and realize the no-slip condition. Since surface thermal fluctuations at room temperature are enough to generate these displacements, we argue that the no-slip condition for water can be achieved also for ideal surfaces, which do not present any surface roughness.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Identifying interfacial molecules in nonplanar interfaces: the generalized ITIM algorithm

We present a generalized version of the ITIM algorithm for the identification of interfacial molecules, which is able to treat arbitrarily shaped interfaces. The algorithm exploits the similarities between the concept of probe sphere used in ITIM and the circumsphere criterion used in the alpha-shapes approach, and can be regarded either as a reference-frame independent version of the former, or as an extended version of the latter that includes the atomic excluded volume. The new algorithm is applied to compute the intrinsic orientational order parameters of water around a DPC and a cholic acid micelle in aqueous environment, and to the identification of solvent-reachable sites in four model structures for soot. The additional algorithm introduced for the calculation of intrinsic density profiles in arbitrary geometries proved to be extremely useful also for planar interfaces, as it allows to solve the paradox of smeared intrinsic profiles far from the interface

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Influence of pore dielectric boundaries on the translocation barrier of DNA

We investigate the impact of dielectric boundary forces on the translocation process of charged rigid DNA segments through solid neutral nanopores. We assess the electrostatic contribution to the translocation free energy barrier of a model DNA segment by evaluating the potential of mean force in absence and presence of polarization effects by means of coarse-grained molecular dynamics simulations. The effect of induced polarization charges has been taken into account by employing ICC*, a recently developed algorithm that can efficiently compute induced polarization charges induced on suitably discretized dielectric boundaries. Since water has a higher dielectric constant than the pore walls, polarization effects repel charged objects in the vicinity of the interface, with the effect of significantly increasing the free energy barrier. Another investigated side effect is the change of the counterion distribution around the charged polymer in presence of the induced pore charges. Furthermore we investigate the influence of adding salt to the solution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Puckering Free Energy of Pyranoses: an NMR and Metadynamics--Umbrella Sampling Investigation

We present the results of a combined metadynamics--umbrella sampling investigation of the puckered conformers of pyranoses described using the gromos 45a4 force field. The free energy landscape of Cremer--Pople puckering coordinates has been calculated for the whole series of alpha and beta aldohexoses, showing that the current force field parameters fail in reproducing proper puckering free energy differences between chair conformers. We suggest a modification to the gromos 45a4 parameter set which improves considerably the agreement of simulation results with theoretical and experimental estimates of puckering free energies. We also report on the experimental measurement of altrose conformers populations by means of NMR spectroscopy, which show good agreement with the predictions of current theoretical models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The ICC* Algorithm: A fast way to include dielectric boundary effects into molecular dynamics simulations

We employ a fast and accurate algorithm to treat dielectric interfaces within molecular dynamics simulations and demonstrate the importance of dielectric boundary forces (DBFs) in two systems of interests in soft-condensed matter science. We investigate a salt solution confined to a slit pore, and a model of a DNA fragment translocating thorugh a narrow pore.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mesoscopic simulations of the counterion-induced electroosmotic flow - a comparative study

We present mesoscopic simulations of the counterion-induced electroosmotic flow in different electrostatic coupling regimes. Two simulation methods are compared, Dissipative Particle Dynamics (DPD) and coupled Lattice-Boltzmann/Molecular Dynamics (LB/MD). A general mapping scheme to match DPD to LB/MD is developed. For the weak-coupling regime, analytic expressions for the flow profiles in the presence of partial-slip as well as no-slip boundary conditions are derived from the Poisson-Boltzmann and Stokes equations, which are in good agreement with the numerical results. The influence of electrofriction and partial slip on the flow profiles is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the Calculation of Puckering Free Energy Surfaces

Cremer-Pople puckering coordinates appear to be the natural candidate variables to explore the conformational space of cyclic compounds, and in literature different parametrizations have been used to this end. However, while every parametrization is equivalent in identifying conformations, it is not obvious that they can also act as proper collective variables for the exploration of the puckered conformations free energy surface. It is shown that only the polar parametrization is fit to produce an unbiased estimate of the free energy landscape. As an example, the case of a six-membered ring, glucuronic acid, is presented, showing the artefacts that are generated when a wrong parametrization is used.

0 citations0 reviewsNo code linkedOpen access