Multiparticle collision dynamics for ferrofluids.

Detailed studies of the intriguing field-dependent dynamics and transport properties of confined flowing ferrofluids require efficient mesoscopic simulation methods that account for fluctuating ferrohydrodynamics. Here, we propose such a new mesoscopic model for the dynamics and flow of ferrofluids, where we couple the multi-particle collision dynamics method as a solver for the fluctuating hydrodynamics equations to the stochastic magnetization dynamics of suspended magnetic nanoparticles. This hybrid model is validated by reproducing the magnetoviscous effect in Poiseuille flow, obtaining the rotational viscosity in quantitative agreement with theoretical predictions. We also illustrate the new method for the benchmark problem of flow around a square cylinder. Interestingly, we observe that the length of the recirculation region is increased, whereas the drag coefficient is decreased in ferrofluids when an external magnetic field is applied compared with the field-free case at the same effective Reynolds number. The presence of thermal fluctuations and the flexibility of this particle-based mesoscopic method provide a promising tool to investigate a broad range of flow phenomena of magnetic fluids, and the method could also serve as an efficient way to simulate solvent effects when colloidal particles are immersed in ferrofluids.

Dynamics of interacting magnetic nanoparticles: effective behavior from competition between Brownian and Néel relaxation.

The intriguing properties of magnetic nanoparticles have sparked a growing number of theoretical studies as well as practical applications. Here, we provide the first comprehensive study of the influence of interactions on the two main relaxation mechanisms: internal (Néel) and Brownian relaxation. While non-interacting magnetic nanoparticles show Debye behavior with an effective relaxation time, many authors use this model also for the interacting case. Since Néel relaxation is typically a thermally activated process on times scales that are many orders of magnitude larger than the underlying micromagnetic times, we use extensive computer simulations employing a Brownian dynamics/Monte-Carlo algorithm to show that dipolar interactions lead to significant deviations from the Debye behavior. We find that Néel and Brownian relaxation can be considered as independent processes for short enough times until dipolar interactions lead to a coupling of these mechanisms, making the interpretation more difficult. We provide mean-field arguments that describe these short and long-time, effective relaxation times well for weak up to moderate interaction strengths. Our findings about the coupling of Brownian and Néel process and the effective relaxation time provide an important theoretical insight that will have also important consequences for the interpretation of magnetic susceptibility measurements and magnetorelaxometry analysis.

Gas-liquid phase equilibrium of a model Langmuir monolayer captured by a multiscale approach.

The gas-liquid expanded phase transition of a Langmuir monolayer happens at very low surface concentrations which makes this phenomenon extremely expensive to explore in finite three-dimensional (3D) atomistic simulations. Starting with a 3D model reference system of amphiphilic surfactants at a 2D vapor-liquid interface, we apply our recently developed approach (Phys. Chem. Chem. Phys., 2018, 20, 16238) and map the entire system to an effective 2D system of surfactant center-of-masses projected onto the interface plane. The coarse-grained interaction potential obtained via a force-matching scheme from the 3D simulations is then used to predict the 2D gas-liquid phase equilibrium of the corresponding Langmuir monolayer. Monte Carlo simulations in the Gibbs ensemble are performed to calculate areal densities, chemical potentials and surface pressures of the gaseous and liquid coexisting phases within the monolayer. We compare these simulations to the results of a 2D density functional approach based on Weeks-Chandler-Anderson perturbation theory. We furthermore use this approach to determine the density profiles across the equilibrium gas-liquid dividing line and the corresponding line tensions.

Entanglement dynamics at flat surfaces: Investigations using multi-chain molecular dynamics and a single-chain slip-spring model.

The dynamics of an entangled polymer melt confined in a channel by parallel plates is investigated by Molecular Dynamics (MD) simulations of a detailed, multi-chain model. A primitive path analysis predicts that the density of entanglements remains approximately constant throughout the gap and drops to lower values only in the immediate vicinity of the surface. Based on these observations, we propose a coarse-grained, single-chain slip-spring model with a uniform density of slip-spring anchors and slip-links. The slip-spring model is compared to the Kremer-Grest MD bead-spring model via equilibrium correlation functions of chain orientations. Reasonably good agreement between the single-chain model and the detailed multi-chain model is obtained for chain relaxation dynamics, both away from the surface and for chains whose center of mass positions are at a distance from the surface that is less than the bulk chain radius of gyration, without introducing any additional model parameters. Our results suggest that there is no considerable drop in topological interactions for chains in the vicinity of a single flat surface. We infer from the slip-spring model that the experimental plateau modulus of a confined polymer melt may be different to a corresponding unconfined system even if there is no drop in topological interactions for the confined case.

Effective interaction potentials for model amphiphilic surfactants adsorbed at fluid-fluid interfaces.

Computer simulations are a useful tool to explore the effects of interactions and structure of surfactants on interfacial microstructure and properties. Starting with "molecular-level", three-dimensional reference systems of short amphiphilic surfactants at fluid-fluid interfaces, we here derive effective interaction potentials for the corresponding two-dimensional systems of structureless particles confined to the interface plane. These reference systems are comprised of two immiscible mono atomic fluids (water- and oil-like particles) and nonionic linear amphiphilic surfactants. Our results show that coarse grained interaction potentials are only weakly dependent on surface concentration but their behavior is strongly dependent on surfactant interactions. The coarse grained system preserves the in-plane surfactant center-of-mass pair correlation function at the interface and the results of surface pressure-area isotherms are in a good agreement. This approach can be extended straightforwardly to other types of surfactants at both fluid-fluid and fluid-gas interfaces providing us with an effective pairwise interaction potential for the surfactant monolayer. These effective interactions can be used to explore large-scale self-assembly within the monolayer especially at low surface concentrations where reference simulations are extremely time-consuming.

Transient inhomogeneous flow patterns in supercooled liquids under shear.

Supercooled liquids and other soft glassy systems show characteristic spatial inhomogeneities in their local dynamical properties. Using detailed molecular dynamics simulations, we find that for sufficiently low temperatures and sufficiently high shear rates supercooled liquids also show transient inhomogeneous flow patterns (shear banding) in the start-up of steady shear flow, similar to what has already been observed for many other soft glassy systems. We verify that the onset of transient shear banding coincides quite well with the appearance of a stress overshoot for temperatures in the supercooled regime. We find that the slower bands adapt less well to the imposed deformation and therefore accumulate higher shear stresses compared to the fast bands at comparable local shear rates. Our results also indicate that the shear rates of the fast and slow bands are adjusted such that the local dissipation rate is approximately the same in both bands.

Nuclear Waste Management under Approaching Disaster: A Comparison of Decommissioning Strategies for the German Repository Asse II.

This article compares different strategies for handling low- and medium-level nuclear waste buried in a retired potassium mine in Germany (Asse II) that faces significant risk of uncontrollable brine intrusion and, hence, long-term groundwater contamination. We survey the policy process that has resulted in the identification of three possible so-called decommissioning options: complete backfilling, relocation of the waste to deeper levels in the mine, and retrieval. The selection of a decommissioning strategy must compare expected investment costs with expected social damage costs (economic, environmental, and health damage costs) caused by flooding and subsequent groundwater contamination. We apply a cost minimization approach that accounts for the uncertainty regarding the stability of the rock formation and the risk of an uncontrollable brine intrusion. Since economic and health impacts stretch out into the far future, we examine the impact of different discounting methods and rates. Due to parameter uncertainty, we conduct a sensitivity analysis concerning key assumptions. We find that retrieval, the currently preferred option by policymakers, has the lowest expected social damage costs for low discount rates. However, this advantage is overcompensated by higher expected investment costs. Considering all costs, backfilling is the best option for all discounting scenarios considered.

Polymer Brushes under Shear: Molecular Dynamics Simulations Compared to Experiments.

Surfaces coated with polymer brushes in a good solvent are known to exhibit excellent tribological properties. We have performed coarse-grained equilibrium and nonequilibrium molecular dynamics (MD) simulations to investigate dextran polymer brushes in an aqueous environment in molecular detail. In a first step, we determined simulation parameters and units by matching experimental results for a single dextran chain. Analyzing this model when applied to a multichain system, density profiles of end-tethered polymer brushes obtained from equilibrium MD simulations compare very well with expectations based on self-consistent field theory. Simulation results were further validated against and correlated with available experimental results. The simulated compression curves (normal force as a function of surface separation) compare successfully with results obtained with a surface forces apparatus. Shear stress (friction) obtained via nonequilibrium MD is contrasted with nanoscale friction studies employing colloidal-probe lateral force microscopy. We find good agreement in the hydrodynamic regime and explain the observed leveling-off of the friction forces in the boundary regime by means of an effective polymer-wall attraction.

Modelling the rheology of anisotropic particles adsorbed on a two-dimensional fluid interface.

We present a general approach based on nonequilibrium thermodynamics for bridging the gap between a well-defined microscopic model and the macroscopic rheology of particle-stabilised interfaces. Our approach is illustrated by starting with a microscopic model of hard ellipsoids confined to a planar surface, which is intended to simply represent a particle-stabilised fluid-fluid interface. More complex microscopic models can be readily handled using the methods outlined in this paper. From the aforementioned microscopic starting point, we obtain the macroscopic, constitutive equations using a combination of systematic coarse-graining, computer experiments and Hamiltonian dynamics. Exemplary numerical solutions of the constitutive equations are given for a variety of experimentally relevant flow situations to explore the rheological behaviour of our model. In particular, we calculate the shear and dilatational moduli of the interface over a wide range of surface coverages, ranging from the dilute isotropic regime, to the concentrated nematic regime.

Influence of inherent structure shear stress of supercooled liquids on their shear moduli.

Configurations of supercooled liquids residing in their local potential minimum (i.e., in their inherent structure, IS) were found to support a non-zero shear stress. This IS stress was attributed to the constraint to the energy minimization imposed by boundary conditions, which keep size and shape of the simulation cell fixed. In this paper, we further investigate the influence of these boundary conditions on the IS stress. We investigate its importance for the computation of the low frequency shear modulus of a glass obtaining a consistent picture for the low- and high frequency shear moduli over the full temperature range. Hence, we find that the IS stress corresponds to a non-thermal contribution to the fluctuation term in the Born-Green expression. This leads to an unphysical divergence of the moduli in the low temperature limit if no proper correction for this term is applied. Furthermore, we clarify the IS stress dependence on the system size and put its origin on a more formal basis.

Soft modes and nonaffine rearrangements in the inherent structures of supercooled liquids.

We find that the hierarchical organization of the potential energy landscape in a model supercooled liquid can be related to a change in the spatial distribution of soft normal modes. For groups of nearby minima, between which fast relaxation processes typically occur, the localization of the soft modes is very similar. The spatial distribution of soft regions changes, instead, for minima between which transitions relevant to structural relaxation occur. This may be the reason why the soft modes are able to predict spatial heterogeneities in the dynamics. Nevertheless, the very softest modes are only weakly correlated with dynamical heterogeneities and instead show higher statistical overlap with regions in the local minima that would undergo nonaffine rearrangements if subjected to a shear deformation. This feature of the supercooled liquid is reminiscent of the behavior of nonaffine deformations in amorphous solids, where the very softest modes identify the loci of plastic instabilities.

Adsorption energies of poly(ethylene oxide)-based surfactants and nanoparticles on an air-water surface.

The self-assembly of polymer-based surfactants and nanoparticles on fluid-fluid interfaces is central to many applications, including dispersion stabilization, creation of novel 2D materials, and surface patterning. Very often these processes involve compressing interfacial monolayers of particles or polymers to obtain a desired material microstructure. At high surface pressures, however, even highly interfacially active objects can desorb from the interface. Methods of directly measuring the energy which keeps the polymer or particles bound to the interface (adsorption/desorption energies) are therefore of high interest for these processes. Moreover, though a geometric description linking adsorption energy and wetting properties through the definition of a contact angle can be established for rigid nano- or microparticles, such a description breaks down for deformable or aggregating objects. Here, we demonstrate a technique to quantify desorption energies directly, by comparing surface pressure-density compression measurements using a Wilhelmy plate and a custom-microfabricated deflection tensiometer. We focus on poly(ethylene oxide)-based polymers and nanoparticles. For PEO-based homo- and copolymers, the adsorption energy of PEO chains scales linearly with molecular weight and can be tuned by changing the subphase composition. Moreover, the desorption surface pressure of PEO-stabilized nanoparticles corresponds to the saturation surface pressure for spontaneously adsorbed monolayers, yielding trapping energies of ∼10(3) k(B)T.

The Landau free energy of hard ellipses obtained from microscopic simulations.

Systems of two-dimensional hard ellipses of varying aspect ratios and packing fractions are studied by Monte Carlo simulations in the generalised canonical ensemble. From this microscopic model, we extract a coarse-grained macroscopic Landau-de Gennes free energy as a function of packing fraction and orientational order parameter. We separate the free energy into the ideal orientational entropy of non-interacting two-dimensional spins and an excess free energy associated with excluded volume interactions. We further explore the isotropic-nematic phase transition using our empirical expression for the free energy and find that the nature of the phase transition is continuous for the aspect ratios we studied.

Nonequilibrium response of magnetic nanoparticles to time-varying magnetic fields: Contributions from Brownian and Néel processes.

Many technical and biomedical applications of magnetic nanoparticles rely on their response to time-varying magnetic fields. While well-established models exist for either immobile or thermally blocked nanoparticles, the intermediate regime where Brownian as well as Néel relaxation occur at the same time is less well explored. Here, we use an efficient model that allows us to study the nonlinear dynamics of individual magnetic nanoparticles in response to different time-varying magnetic fields over a broad range of field parameters, taking into account both relaxation mechanisms. We provide quasiexact solutions for the longitudinal dynamics as well as approximate formulas from dynamic mean-field theory. Our results are relevant, e.g., for magnetorelaxometry, magnetic fluid hyperthermia, and magnetic particle imaging. For these example applications, we show that the ratio of characteristic Brownian to Néel relaxation time can have a profound impact on characteristic response quantities, especially at large field strengths.

Chondrosarcoma evaluation using hematein-based x-ray staining and high-resolution 3D micro-CT: a feasibility study.

BACKGROUND: Chondrosarcomas are rare malignant bone tumors diagnosed by analyzing radiological images and histology of tissue biopsies and evaluating features such as matrix calcification, cortical destruction, trabecular penetration, and tumor cell entrapment. METHODS: We retrospectively analyzed 16 cartilaginous tumor tissue samples from three patients (51-, 54-, and 70-year-old) diagnosed with a dedifferentiated chondrosarcoma at the femur, a moderately differentiated chondrosarcoma in the pelvis, and a predominantly moderately differentiated chondrosarcoma at the scapula, respectively. We combined a hematein-based x-ray staining with high-resolution three-dimensional (3D) microscopic x-ray computed tomography (micro-CT) for nondestructive 3D tumor assessment and tumor margin evaluation. RESULTS: We detected trabecular entrapment on 3D micro-CT images and followed bone destruction throughout the volume. In addition to staining cell nuclei, hematein-based staining also improved the visualization of the tumor matrix, allowing for the distinction between the tumor and the bone marrow cavity. The hematein-based staining did not interfere with further conventional histology. There was a 5.97 ± 7.17% difference between the relative tumor area measured using micro-CT and histopathology (p = 0.806) (Pearson correlation coefficient r = 0.92, p = 0.009). Signal intensity in the tumor matrix (4.85 ± 2.94) was significantly higher in the stained samples compared to the unstained counterparts (1.92 ± 0.11, p = 0.002). CONCLUSIONS: Using nondestructive 3D micro-CT, the simultaneous visualization of radiological and histopathological features is feasible. RELEVANCE STATEMENT: 3D micro-CT data supports modern radiological and histopathological investigations of human bone tumor specimens. It has the potential for being an integrative part of clinical preoperative diagnostics. KEY POINTS: • Matrix calcifications are a relevant diagnostic feature of bone tumors. • Micro-CT detects all clinically diagnostic relevant features of x-ray-stained chondrosarcoma. • Micro-CT has the potential to be an integrative part of clinical diagnostics.

Field- and concentration-dependent relaxation of magnetic nanoparticles and optimality conditions for magnetic fluid hyperthermia.

The field-dependent relaxation dynamics of suspended magnetic nanoparticles continues to present a fascinating topic of basic science that at the same time is highly relevant for several technological and biomedical applications. Renewed interest in the intriguing behavior of magnetic nanoparticles in response to external fields has at least in parts be driven by rapid advances in magnetic fluid hyperthermia research. Although a wealth of experimental, theoretical, and simulation studies have been performed in this field in recent years, several contradictory findings have so far prevented the emergence of a consistent picture. Here, we present a dynamic mean-field theory together with comprehensive computer simulations of a microscopic model system to systematically discuss the influence of several key parameters on the relaxation dynamics, such as steric and dipolar interactions, the external magnetic field strength and frequency, as well as the ratio of Brownian and Néel relaxation time. We also discuss the specific and intrinsic loss power as measures of the efficiency of magnetic fluid heating and discuss optimality conditions in terms of fluid and field parameters. Our results are helpful to reconcile contradictory findings in the literature and provide an important step towards a more consistent understanding. In addition, our findings also help to select experimental conditions that optimize magnetic fluid heating applications.

Deformation of inherent structures to detect long-range correlations in supercooled liquids.

We propose deformations of inherent structures as a suitable tool for detecting structural changes underlying the onset of cooperativity in supercooled liquids. The non-affine displacement (NAD) field resulting from the applied deformation shows characteristic differences between the high temperature liquid and supercooled state, which are typically observed in dynamic quantities. The average magnitude of the NAD is very sensitive to temperature changes in the supercooled regime and is found to be strongly correlated with the inherent structure energy. In addition, the NAD field is characterized by a correlation length that increases upon lowering the temperature towards the supercooled regime.

Simulating the flow of interacting ferrofluids with multiparticle collision dynamics.

Ferrofluid flow is fascinating since its fluid properties can conveniently be manipulated by external magnetic fields. Novel applications in micro- and nanofluidics as well as in biomedicine have renewed the interest in the flow of colloidal magnetic nanoparticles with a focus on small-scale behavior. Traditional flow simulations of ferrofluids, however, often use simplified constitutive models and do not include fluctuations that are relevant at small scales. Here we address these challenges by proposing a hybrid scheme that combines the multiparticle collision dynamics method for modeling hydrodynamics with Brownian dynamics simulations of a reliable kinetic model describing the microstructure, magnetization dynamics, and resulting stresses. Since both multiparticle collision dynamics and Brownian dynamics are mesoscopic methods that naturally include fluctuations, this hybrid scheme presents a promising alternative to more traditional approaches, also because of the flexibility to model different geometries and modifying the constitutive model. The scheme was tested in several ways. Poiseuille flow was simulated for various model parameters and effective viscosities were determined from the resulting flow profiles. The effective, field-dependent viscosities are found to be in very good agreement with theoretical predictions. We also study Stokes' second flow problem for ferrofluids. For weak amplitudes and low frequencies of the oscillating plate, we find that the velocity profiles are well described by the result for Newtonian fluids at the corresponding, field-dependent viscosity. Furthermore, the time-dependent profiles of the nonequilibrum magnetization component are well approximated by their steady-state values in stationary shear when evaluated with the instantaneous local shear rate. Finally, we also apply our scheme to simulate ferrofluid shear flow over a rough surface. We find characteristic differences in the nonequilibrium magnetization components when the external field is oriented in flow and in a gradient direction.

Probing a critical length scale at the glass transition.

We give evidence of a clear structural signature of the glass transition, in terms of a static correlation length with the same dependence on the system size, which is typical of critical phenomena. Our approach is to introduce an external, static perturbation to extract the structural information from the system's response. In particular, we consider the transformation behavior of the local minima of the underlying potential energy landscape (inherent structures), under a static deformation. The finite-size scaling analysis of our numerical results indicate that the correlation length diverges at a temperature Tc, below the temperatures where the system can be equilibrated. Our numerical results are consistent with random first order theory, which predicts such a divergence with a critical exponent ν=2/3 at the Kauzmann temperature, where the extrapolated configurational entropy vanishes.