Hitoshi Washizu

Vapor-liquid equilibrium (VLE) data of fragrance components are crucial for product development and separation processes. However, experimentally obtaining these data can often be a high-cost and challenging task. In order to address this issue, simulations of VLE data using molecular dynamics (MD) methods have gained attention, though there are still relatively few studies about the vapor-liquid equilibrium calculations of fragrance components using MD. In this study, we focused on a mixture of d-limonene and 1-pentanol as representative components and conducted MD simulations. The VLE data obtained by varying the molar fraction of d-limonene, including x-y phase diagrams and activity coefficients, showed a high degree of agreement with the experimental data. Additionally, an analysis of the density profiles on a molecular level revealed a slight increase in the concentration of 1-pentanol at the vapor-liquid interface.

Molecular dynamics simulations were conducted using the generalized replica exchange method (gREM) on the 4-cyano-4′-n-alkyl biphenyl (nCB) system with n = 5, 6, 7, and 8, which exhibits a nematic–isotropic (NI) phase transition. Sampling near the phase transition temperature in systems undergoing first-order phase transitions, such as the NI phase transition, is demanding due to the substantial energy gap between the two phases. To address this, gREM, specifically designed for first-order phase transitions, was utilized to enhance sampling near the NI phase transition temperature. Free-energy calculations based on the energy representation (ER) theory were employed to characterize the NI phase transition. ER evaluates the insertion free energy of the nCB molecule for both nematic and isotropic phases, revealing a change in the temperature dependence across the NI phase transition. Further decomposition into energetic and entropic terms quantitatively shows the balance between these contributions at the NI phase transition temperature.

A reactive force field (ReaxFF) molecular dynamics simulation is performed for the sliding of diamond-like carbon (DLC) and yttria-stabilized zirconia (YSZ) under an ethanol gas environment, motivated by the previous experiment of ultralow friction phenomenon (friction fade-out). We observe (i) dissociation of ethanol into ethoxy and hydrogen, both of which simultaneously adsorb on the YSZ surface, and (ii) dissociation of ethanol into ethyl and hydroxy, the former of which forms a bond with another ethanol molecule and the latter of which adsorbs on the DLC surface. Reaction (i) is enhanced by the sliding motion, but occurs even without it, while reaction (ii) only occurs during sliding with a sufficiently high load pressure. The potentials of mean force for the two reactions are also calculated combining the steered MD and Jarzynski equality. It is shown that the activation energies of reactions (i) and (ii) are significantly lowered by the YSZ and DLC surfaces, respectively, as compared to those in a vacuum. The resultant activation energy is higher for reaction (ii) than for reaction (i).

In many physical or biological systems, diffusion can be described by Brownian motions with stochastic diffusion coefficients (DCs). In the present study, we investigate properties of the diffusion with a broad class of stochastic DCs with a novel approach. We show that for a finite time, the propagator is non-Gaussian and heavy-tailed. This means that when the mean square displacements are the same, for a finite time, some of the diffusing particles with stochastic DCs diffuse farther than the particles with deterministic DCs or exhibiting a fractional Brownian motion. We also show that when a stochastic DC is ergodic, the propagator converges to a Gaussian distribution in the long time limit. The speed of convergence is determined by the autocovariance function of the DC.

Organophosphates are well-known as the canonical additives for lubricants. Thus, understanding of the additive behaviour is a key aspect in the design of films on metal surfaces. Different types of phosphates are added to improve their antiwear properties, but the contributions of individual esters to these properties has not been studied using a combination of practical and theoretical approaches. In this study, organophosphates were isolated with high purity and their tribological characteristics were evaluated by using a Bowden-type reciprocating friction tester and a four-ball wear tester. Mono-oleylphosphate had a lower friction than di-oleylphosphate and exhibited excellent antiwear characteristics. Analysis of the sliding surfaces using desorption electrospray ionization-mass spectrometry (DESI-MS) and X-ray photoelectron spectroscopy (XPS) indicated that the film structure could predict the occurrence factor of the tribological characteristics of the oleylphosphates. Then the adsorption energies of the monoester on iron and iron oxide surfaces were higher than those of the diester, as assessed using density functional theory (DFT) calculations, owing to the difference in their chemisorption processes, as confirmed by further DFT analysis. Studies on the reactivity of additives and their interactions with surfaces are important for understanding the tribochemistry of additives.

Reentrant phenomena in soft matter and biosystems have attracted considerable attention because their properties are closely related to high functionality. Here, we report a combined experimental and computational study on the self-assembly and reentrant behavior of a single-component thermotropic smectic liquid crystal toward the realization of dynamically functional materials. We have designed and synthesized a mesogenic molecule consisting of an alicyclic trans,trans -bicyclohexyl mesogen and a polar cyclic carbonate group connected by a flexible tetra(oxyethylene) spacer. The molecule exhibits an unprecedented sequence of layered smectic phases, in the order: smectic A-smectic B-reentrant smectic A. Electron density profiles and large-scale molecular dynamics simulations indicate that competition between the stacking of bicyclohexyl mesogens and the conformational flexibility of tetra(oxyethylene) chains induces this unusual reentrant behavior. Ion-conductive reentrant liquid-crystalline materials have been developed, which undergo the multistep conductivity changes in response to temperature. The reentrant liquid crystals have potential as new mesogenic materials exhibiting switching functions.

In this paper, we classified the types of water in the vicinity of the chitosan polymer and gold plate by applying an electric field of magnitude 1 V Å-1 in various directions at varying temperatures by using molecular dynamics simulation. The three types of water were categorized by analyzing the data through the tetrahedral order method with four water regions separated in the distance from 1 to 6 Å around polymers. The interaction between water molecules and functional groups, such as hydroxyl, ether, and ester, leads to the formation of intermediate and nonfreezing water. Under an electric field, this formation appeared more clearly due to the transformation of liquid water to crystal cubic ice with two structural formations depending on gold plates at a temperature of 300 K. The enhancement of the tetrahedral order of water in cubic ice is related to the existence of a four-fold H-bonded structure and lower ones in the XES experiment.

A key to achieve the accuracy of molecular dynamics (MD) simulation is the set of force fields used to express the atomistic interactions. In particular, the electrostatic interaction remains the main issue for the precise simulation of various ionic soft materials from ionic liquids to their supramolecular compounds. In this study, we test the nonpolarizable force fields of ionic liquids (ILs) and self-assembled ionic liquid crystals (ILCs) for which the intermolecular charge transfer and intramolecular polarization are significant. The self-consistent modeling scheme is adopted to refine the atomic charges of ionic species in a condensed state through the use of density functional theory (DFT) under the periodic boundary condition. The atomic charges of the generalized amber force field (GAFF) are effectively updated to express the electrostatic properties of ionic molecules obtained by the DFT calculation in condensed phase, which improves the prediction accuracy of ionic conductivity with the obtained force field (GAFF-DFT). The derived DFT charges then suggest that the substitution of a hydrophobic liquid-crystalline moiety into IL-based cations enhances the charge localization of ionic groups in the amphiphilic molecules, leading to the amplification of the electrostatic interactions among the hydrophilic/ionic groups in the presence of hydrophobic moieties. In addition, we focus on an ion-conductive pathway hidden in the self-assembled nanostructure. The MD results indicate that the ionic groups of cation and anion interact strongly for keeping the bicontinuous nanosegregation of ionic nanochannel. The partial fractions of hydrophilic/ionic and hydrophobic nanodomains are then quantified with the volume difference from referenced IL systems, while the calculated ionic conductivity decreases in the self-assembled ILCs more than the occupied volume of ionic nanodomains. These analyses suggest that the mobility of ions in the self-assembled ILCs remains quite restricted even with small tetrafluoroborate anions because of strong attractive interaction among ionic moieties.

Elucidating the state of interfacial water, especially the hydrogen-bond configurations, is considered to be key for a better understanding of the functions of polymers that are exhibited in the presence of water. Here, an analysis in this direction is conducted for two water-insoluble biocompatible polymers, poly(2-methoxyethyl acrylate) and cyclic(poly(2-methoxyethyl acrylate)), and a non-biocompatible polymer, poly(n-butyl acrylate), by measuring their IR spectra under humidified conditions and by carrying out theoretical calculations on model complex systems. It is found that the OH stretching bands of water are decomposed into four components, and while the higher-frequency components (with peaks at ∼3610 and ∼3540 cm-1) behave in parallel with the C═O and C-O-C stretching and CH deformation bands of the polymers, the lower-frequency components (with peaks at ∼3430 and ∼3260 cm-1) become pronounced to a greater extent with increasing humidity. From the theoretical calculations, it is shown that the OH stretching frequency that is distributed from ∼3650 to ∼3200 cm-1 is correlated to the hydrogen-bond configurations and is mainly controlled by the electric field that is sensed by the vibrating H atom. By combining these observed and calculated results, the configurations of water at the interface of the polymers are discussed.

This study focuses on designing solid lubricant particles by combining graphene and iron nanoparticles (namely, graphene-iron (GI) particles) and carrying out studies for behaviors of their lubrication for the iron contact by molecular dynamics simulations. By the annealing process of melting and cooling iron, we can create the lubricant particle, where the iron nanoparticle tightly holds the graphene sheet. In the sliding friction investigations, it is found that the influences of orientation of the graphene sheets inside the contact, size and configuration of the GI particles, and lubrication with the bare iron nanoparticles on friction are strong at low pressure and very slight at high pressure. The GI particles provide stability of the friction coefficient over a wide range of pressure; however, it strongly increases with pressure in the lubrication behaviors by the bare iron particles due to the deformation of the particles. The iron contact in the presence of the GI particles can achieve the ultralow values of the friction coefficient from 0.009 to 0.042. The contact surfaces are not nearly damaged (slightly elastic deformation) with the pressure up to 2.0 GPa. From the comparisons between the results in this study and previous reports, the GI particles have better lubrication than graphene coated on a surface and well stabilize under pressure compared to the different lubricant nanoparticles. The main reason for this is due to the contributions of graphene, besides reduction of the contact area resulted from the configuration of the nanoparticle, which promotes sliding and sharing of the pressure, preventing collision between the lubricant particles.

In materials science, water plays an important part, especially at the molecular level. It shows various properties when sorbed onto surfaces of polymers. The structure of the molecular water ensemble in the vicinity of the polymers is under discussion. In this study, we used molecular dynamics methods to analyze the structure of water in the vicinity of the polymer polyrotaxane (PR), composed of α-cyclodextrins (α-CDs), a poly(ethylene glycol) (PEG) axial chain, and α-lipoic acid linkers, at various temperatures. The distribution of water around the functional groups, hydrogen bond network, and tetrahedral order were analyzed to classify the various types of water around the polymer. We found that the tetrahedral order of water had a strained relationship from the XES experiment. Four water regions were separated from each other in the vicinity of 1 to 5 Å around PR. The intermediate and non-freezing water were formed due to the interaction between water molecules and the functional groups, such as hydroxyl, ether, and ester.

The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in the gREM, allowing for the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyzed with different system sizes and compared with that of the temperature replica-exchange method (tREM). It was observed that gREM is capable of sampling configurations at sufficient replica-exchange acceptance ratios even around the NI transition temperature. A bimodal distribution of the order parameter at the transition region was found, which is in agreement with the mean-field theory. In contrast, tREM is ineffective around the transition temperature owing to the potential energy gap between the nematic and isotropic phases.

Self-assembled ionic liquid crystals can transport water and ions via the periodic nanochannels, and these materials are promising candidates as water treatment membranes. Molecular insights on the water transport process are, however, less investigated because of computational difficulties of ionic soft matters and the self-assembly. Here we report specific behavior of water molecules in the nanochannels by using the self-consistent modeling combining density functional theory and molecular dynamics and the large-scale molecular dynamics calculation. The simulations clearly provide the one-dimensional (1D) and 3D-interconnected nanochannels of self-assembled columnar and bicontinuous structures, respectively, with the precise mesoscale order observed by x-ray diffraction measurement. Water molecules are then confined inside the nanochannels with the formation of hydrogen bonding network. The quantitative analyses of free energetics and anisotropic diffusivity reveal that, the mesoscale geometry of 1D nanodomain profits the nature of water transport via advantages of dissolution and diffusion mechanisms inside the ionic nanochannels.

<title>Abstract</title>The optimal method of the polymer Materials Informatics (MI) has not been developed because the amorphous nature of the higher-order structure affects these properties. We have now tried to develop the polymer MI’s descriptor of the higher-order structure using persistent homology as the topological method. We have experimentally studied the influence of the MD simulation cell size as the higher-order structure of the polymer on its electrical properties important for a soft material sensor or actuator device. The all-atom MD simulation of the polymer has been calculated and the obtained atomic coordinate has been analyzed by the persistent homology. The change in the higher-order structure by different cell size simulations affects the dielectric constant, although these changes are not described by a radial distribution function (RDF). On the other hand, using the 2nd order persistent diagram (PD), it was found that when the cell size is small, the island-shaped distribution become smoother as the cell size increased. There is the same tendency for the condition of change in the monomer ratio, the polymer chain length or temperature. As a result, the persistent homology may express the higher-order structure generated by the MD simulation as a descriptor of the polymer MI.

<title>Abstract</title> The study focuses on monitoring influence of the alumina coatings on friction and stability of the microscale iron contacts by the smoothed particle hydrodynamics. The obtained results show a better stability and a higher value of the friction coefficient of the coated surface compared to the uncoated one. This study also supposes that in concern of stability of the surface the coating should be done with only the substrate surface accompanied with the roughness of the coating layer. The proportion of the coated particles is found to be strongly resulting in the friction properties and the roughness of the coating layer also slightly results in those. The surface reaches the most stability at the proportion of around 70%.

<title>Abstract</title> This paper investigates the friction and friction heat of the micronscale iron under the influences such as the velocity of the slider and temperature of the substrate by using the smoothed particle hydrodynamics simulations. It is found that in the velocity range of 10–100 m/s, the sliding velocity–friction coefficient relationship well complies with the fitted exponent or hyperbolic tangent function, and the friction coefficient approaches a stable value of 0.3 at around the velocity of 50 m/s after a rapidly increasing situation. The steady friction coefficient maintains over the temperature range of 200–400 K. The friction heat is detailed analyzed versus the sliding time. The sliding time–system temperature relationship is well fitted by the sigmoidal functions, except the interfacial particle layers. The layer causing friction shows the highest steady temperature and largest temperature rise. The increment between the initial temperatures of the slider and the substrate strongly results in the temperature rise while it does not affect the configuration of the sliding time–system temperature curves.

<p>The paper focuses on examining agreement of the adaptive smoothed particle hydrodynamics (ASPH) in the investigation of the sliding friction of silica at micronscale throughout observation of several friction characteristics. It is found that the ASPH approach well presents the friction of micronscale silica due to agreement of the friction coefficient and the applied load-friction coefficient relationship between the present results and the previously experimental reports. The shape of the particle modeled in the ASPH almost does not effect on the detected results for the hard system due to the very slight variation of the particles during the sliding. However, the variation of the particles can explain for the discrepancy between the stick time and the slip time and the unsharp change between the stick state and the slip one. The study is also extended for the contacts of the two sinusoidal rough surfaces and finds that the friction coefficient is almost independent of the wavelength while it linearly increases with the amplitude.</p>

A large-scale parallel computing model based on the SPH (Smoothed Particle Hydrodynamics) method was developed for dry shear friction between elastic-plastic solids. Our main purpose is to elucidate the mechanism of the frictional wear and heat generation between the asperities on the interface in meso-scale. In our model, the elastic-plastic motion is expressed by the general SPH method for the solid. The frictional interaction between the asperities is also added as a two-body interaction between SPH particles. In thus general SPH model, the isotropic weighting function is used and the SPH particle is the isotropic sphere. However, the surface roughness of the real metal interfaces is very small with respect to the size of the frictional surface area. The aspect ratio between the sliding direction and the thickness direction of the asperities is very large. For the modified model, we propose a method using coordinate transformation, in which the unit length between the shearing direction and the direction perpendicular to the shear plane is different. We call it the “Disk-like SPH method”. We show the validity of thus model and think it as the effective coarsening scheme to clear the wear and heat generation at the real surface.