Erich A. Müller

Optimizing Water Transport through Graphene-Based Membranes: Insights from Nonequilibrium Molecular Dynamics

Recent experimental results suggest that stacked layers of graphene oxide exhibit strong selective permeability to water. To construe this observation, the transport mechanism of water permeating through a membrane consisting of layered graphene sheets is investigated via nonequilibrium and equilibrium molecular dynamics simulations. The effect of sheet geometry is studied by changing the offset between the entrance and exit slits of the membrane. The simulation results reveal that the permeability is not solely dominated by entrance effects; the path traversed by water molecules has a considerable impact on the permeability. We show that contrary to speculation in the literature, water molecules do not pass through the membrane as a hydrogen-bonded chain; instead, they form well-mixed fluid regions confined between the graphene sheets. The results of the present work are used to provide guidelines for the development of graphene and graphene oxide membranes for desalination and solvent separation.

Surface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water

The test-area (TA) perturbation approach has been gaining popularity as a methodology for the direct computation of the interfacial tension in molecular simulation. Though originally implemented for planar interfaces, the TA approach has also been used to analyze the interfacial properties of curved liquid interfaces. Here, we provide an interpretation of the TA method taking the view that it corresponds to the change in free energy under a transformation of the spatial metric for an affine distortion. By expressing the change in configurational energy of a molecular configuration as a Taylor expansion in the distortion parameter, compact relations are derived for the interfacial tension and its energetic and entropic components for three different geometries: planar, cylindrical, and spherical fluid interfaces. While the tensions of the planar and cylindrical geometries are characterized by first-order changes in the energy, that of the spherical interface depends on second-order contributions. We show that a greater statistical uncertainty is to be expected when calculating the thermodynamic properties of a spherical interface than for the planar and cylindrical cases, and the evaluation of the separate entropic and energetic contributions poses a greater computational challenge than the tension itself. The methodology is employed to determine the vapour-liquid interfacial tension of TIP4P/2005 water at 293 K by molecular dynamics simulation for planar, cylindrical, and spherical geometries. A weak peak in the curvature dependence of the tension is observed in the case of cylindrical threads of condensed liquid at a radius of about 8 Å, below which the tension is found to decrease again. In the case of spherical drops, a marked decrease in the tension from the planar limit is found for radii below ∼ 15 Å; there is no indication of a maximum in the tension with increasing curvature. The vapour-liquid interfacial tension tends towards the planar limit for large system sizes for both the cylindrical and spherical cases. Estimates of the entropic and energetic contributions are also evaluated for the planar and cylindrical geometries and their magnitudes are in line with the expectations of our simple analysis.

Modelling the interfacial behaviour of dilute light-switching surfactant solutions

The direct molecular modelling of an aqueous surfactant system at concentrations below the critical micelle concentration (pre-cmc) conditions is unviable in terms of the presently available computational power. Here, we present an alternative that combines experimental information with tractable simulations to interrogate the surface tension changes with composition and the structural behaviour of surfactants at the water-air interface. The methodology is based on the expression of the surface tension as a function of the surfactant surface excess, both in the experiments and in the simulations, allowing direct comparisons to be made. As a proof-of-concept a coarse-grained model of a light switching non-ionic surfactant bearing a photosensitive azobenzene group is considered at the air-water interface at 298 K. Coarse-grained molecular dynamic simulations are detailed based on the use of the SAFT force field with parameters tuned specifically for this purpose. An excellent agreement is obtained between the simulation predictions and experimental observations; furthermore, the molecular model allows the rationalization of the macroscopic behaviour in terms of the different conformations of the cis and trans surfactants at the surface.

Electron microscopic evidence for a tribologically induced phase transformation as the origin of wear in diamond

Tribological testing of a coarse-grained diamond layer, deposited by plasma-enhanced chemical vapor deposition, was performed on a ring-on-ring tribometer with a diamond counterpart. The origin of the wear of diamond and of the low friction coefficient of 0.15 was studied by analyzing the microstructure of worn and unworn regions by transmission and scanning electron microscopy. In the worn regions, the formation of an amorphous carbon layer with a thickness below 100u2009nm is observed. Electron energy loss spectroscopy of the C-K ionization edge reveals the transition from sp3-hybridized C-atoms in crystalline diamond to a high fraction of sp2-hybridized C-atoms in the tribo-induced amorphous C-layer within a transition region of less than 5u2009nm thickness. The mechanically induced phase transformation from diamond to the amorphous phase is found to be highly anisotropic which is clearly seen at a grain boundary, where the thickness of the amorphous layer above the two differently oriented grains abruptly cha...

Quantification of sample thickness and in-concentration of InGaAs quantum wells by transmission measurements in a scanning electron microscope.

High-angle annular dark-field (HAADF) scanning transmission electron microscopy (STEM) images of electron-transparent samples show dominant atomic number (Z-) contrast with a high lateral resolution. HAADF STEM at low electron energies <30 keV is applied in this work for quantitative composition analyses of InGaAs quantum wells. To determine the local composition, normalized experimental image intensities are compared with results of Monte Carlo simulations. For verification of the technique, InGaAs/GaAs quantum-well structures with known In concentration are used. Transmission electron microscopy samples with known thickness are prepared by the focused-ion-beam technique. The method can be extended to other material systems and is particularly promising for the analysis of materials that are sensitive toward knock-on damage.

Bulk Heterojunction Nanomorphology of Fluorenyl Hexa-peri-hexabenzocoronene–Fullerene Blend Films

In this study, the nanomorphology of fluorenyl hexa-peri-hexabenzocoronene:[6,6]-phenyl C61-butyric acid methyl ester (FHBC:PC61BM) absorber layers of organic solar cells was investigated. Different electron microscopical techniques, atomic force microscopy, and grazing incidence wide-angle X-ray scattering were applied for a comprehensive nanomorphology analysis. The development of the nanomorphology upon sample annealing and the associated change of the device performance were investigated. It was shown that the annealing process enhances the phase separation and therefore the bulk heterojunction structure. Due to π-π stacking, the FHBC molecules assemble into columnar stacks, which are already present before annealing. While the nonannealed sample consists of a mixture of homogeneously distributed PC61BM molecules and FHBC stacks with a preferential in-plane stack orientation, crystalline FHBC precipitates occur in the annealed samples. These crystals, which consist of hexagonal arranged FHBC stacks, grow with increased annealing time. They are distributed homogeneously over the whole volume of the absorber layer as revealed by electron tomography. The FHBC stacks, whether in the two phase mixture or in the pure crystalline precipitates, exhibit an edge-on orientation, according to results from grazing incidence wide-angle X-ray scattering (GIWAXS), dark-field transmission electron microscopy (DF TEM) imaging and selective area electron diffraction (SAED). The best solar cell efficiencies were obtained after 20 or 40 s sample annealing. These annealing times induce an optimized degree of phase separation between donor and acceptor material.

Physical insights into the blood–brain barrier translocation mechanisms

The number of individuals suffering from diseases of the central nervous system (CNS) is growing with an aging population. While candidate drugs for many of these diseases are available, most of these pharmaceutical agents cannot reach the brain rendering most of the drug therapies that target the CNS inefficient. The reason is the blood-brain barrier (BBB), a complex and dynamic interface that controls the influx and efflux of substances through a number of different translocation mechanisms. Here, we present these mechanisms providing, also, the necessary background related to the morphology and various characteristics of the BBB. Moreover, we discuss various numerical and simulation approaches used to study the BBB, and possible future directions based on multi-scale methods. We anticipate that this review will motivate multi-disciplinary research on the BBB aiming at the design of effective drug therapies.

Assembly of porous smectic structures formed from interlocking high-symmetry planar nanorings.

Significance The formation of low-density porous structures is currently a topic of significant interest due to the advantageous electrical, optical, and chemical properties of the materials. We have observed the formation of unique highly open liquid-crystalline smectic phases formed by the interlocking of highly symmetrical planar nanorings. In particular, we demonstrate the relationship between the size of the internal cavity of the rings and their symmetry on the formation of stable liquid-crystalline phases with free volumes of up to 95%. Materials comprising porous structures, often in the form of interconnected concave cavities, are typically assembled from convex molecular building blocks. The use of nanoparticles with a characteristic nonconvex shape provides a promising strategy to create new porous materials, an approach that has been recently used with cagelike molecules to form remarkable liquids with “scrabbled” porous cavities. Nonconvex mesogenic building blocks can be engineered to form unique self-assembled open structures with tunable porosity and long-range order that is intermediate between that of isotropic liquids and of crystalline solids. Here we propose the design of highly open liquid-crystalline structures from rigid nanorings with ellipsoidal and polygonal geometry. By exploiting the entropic ordering characteristics of athermal colloidal particles, we demonstrate that high-symmetry nonconvex rings with large internal cavities interlock within a 2D layered structure leading to the formation of distinctive liquid-crystalline smectic phases. We show that these smectic phases possess uniquely high free volumes of up to ∼95%, a value significantly larger than the 50% that is typically achievable with smectic phases formed by more conventional convex rod- or disklike mesogenic particles.

In Silico Determination of Gas Permeabilities by Non-Equilibrium Molecular Dynamics: CO2 and He through PIM-1

We study the permeation dynamics of helium and carbon dioxide through an atomistically detailed model of a polymer of intrinsic microporosity, PIM-1, via non-equilibrium molecular dynamics (NEMD) simulations. This work presents the first explicit molecular modeling of gas permeation through a high free-volume polymer sample, and it demonstrates how permeability and solubility can be obtained coherently from a single simulation. Solubilities in particular can be obtained to a very high degree of confidence and within experimental inaccuracies. Furthermore, the simulations make it possible to obtain very specific information on the diffusion dynamics of penetrant molecules and yield detailed maps of gas occupancy, which are akin to a digital tomographic scan of the polymer network. In addition to determining permeability and solubility directly from NEMD simulations, the results shed light on the permeation mechanism of the penetrant gases, suggesting that the relative openness of the microporous topology promotes the anomalous diffusion of penetrant gases, which entails a deviation from the pore hopping mechanism usually observed in gas diffusion in polymers.

Water droplet excess free energy determined by cluster mitosis using guided molecular dynamics

Atmospheric aerosols play a vital role in affecting climate by influencing the properties and lifetimes of clouds and precipitation. Understanding the underlying microscopic mechanisms involved in the nucleation of aerosol droplets from the vapour phase is therefore of great interest. One key thermodynamic quantity in nucleation is the excess free energy of cluster formation relative to that of the saturated vapour. In our current study, the excess free energy is extracted for clusters of pure water modelled with the TIP4P/2005 intermolecular potential using a method based on nonequilibrium molecular dynamics and the Jarzynski relation. The change in free energy associated with the mitosis or division of a cluster of N water molecules into two N/2 sub-clusters is evaluated. This methodology is an extension of the disassembly procedure used recently to calculate the excess free energy of argon clusters [H. Y. Tang and I. J. Ford, Phys. Rev. E 91, 023308 (2015)]. Our findings are compared to the corresponding excess free energies obtained from classical nucleation theory (CNT) as well as internally consistent classical theory (ICCT). The values of the excess free energy that we obtain with the mitosis method are consistent with CNT for large cluster sizes but for the smallest clusters, the results tend towards ICCT; for intermediate sized clusters, we obtain values between the ICCT and CNT predictions. Furthermore, the curvature-dependent surface tension which can be obtained by regarding the clusters as spherical droplets of bulk density is found to be a monotonically increasing function of cluster size for the studied range. The data are compared to other values reported in the literature, agreeing qualitatively with some but disagreeing with the values determined by Joswiak et al. [J. Phys. Chem. Lett. 4, 4267 (2013)] using a biased mitosis approach; an assessment of the differences is the main motivation for our current study.