Molecular dynamics simulation of the equilibrium liquid–vapor interphase with solidification

The equilibrium structure of the finite, interphase interfacial region that exists between a liquid film and a bulk vapor is resolved by molecular dynamics simulation. Argon systems are considered for a temperature range that extends below the melting point. Physically consistent procedures are developed to define the boundaries between the interphase and the liquid and vapor phases. The procedures involve counting of neighboring molecules and comparing the results with boundary criteria that permit the boundaries to be precisely established. Two-dimensional radial distribution functions at the liquid and vapor boundaries and within the interphase region demonstrate the physical consistency of the boundary criteria and the state of transition within the region. The method developed for interphase boundary definitions can be extended to nonequilibrium systems. Spatial profiles of macroscopic properties across the interphase region are presented. A number of interfacial thermodynamic properties and profile curve-fit parameters are tabulated, including evaporation/condensation coefficients determined from molecular flux statistics. The evaporation/condensation coefficients away from the melting point compare more favorably with transition state theory than those of previous simulations. Near the melting point, transition theory approximations are less valid and the present results differ from the theory. The effects of film substrate wetting on evaporation/condensation coefficients are also presented.     

Molecular dynamics simulation of surfactant effects on ion transport through a liquid–liquid interface between partially miscible liquids

Molecular dynamics (MD) simulations of water + 1-hexanol + NaCl mixtures with and without a surfactant (methanol) were performed to analyze the surfactant's effect on the transport of a sodium ion through the liquid–liquid interface. Without surfactant, the 1-hexanol forms a bilayer at the interface with OH groups directed outward toward the aqueous phase. Addition of the surfactant produces higher concentrations of the surfactant on the aqueous side of the interface without altering the organic bilayer structure. An electrical double-layer is created in both cases as chloride ion concentration is enhanced near the interface and sodium ion concentration is enhanced toward the center of the water phase. A potential of mean force (PMF) was calculated for the transfer of a sodium ion through the interface. Although the surfactant reduced the interfacial tension, the total work required for the ion transfer increased with the addition of the surfactant.     

Molecular simulation of vapor–liquid equilibria of toxic gases

In this paper, we derived the potential parameters for three toxic gases, hydrogen sulfide, phosgene and nitrous oxide, modeled by the effective Stockmayer potential model proposed by Gao et al. [Fluid Phase Equilib. 137 (1997) 87]. The vapor–liquid equilibria (VLE) of these substances have been extensively investigated over a wide range of temperatures by the Gibbs ensemble Monte Carlo (GEMC) technique. The simulated saturated densities and pressures are in good agreement with experimental data. The critical properties obtained by regression of the simulated data also agree well with the experimental values. The present work demonstrates that the effective Stockmayer potential can describe well the toxic gases concerned.     

Why does the intermolecular dynamics of liquid biphenyl so closely resemble that of liquid benzene? Molecular dynamics simulation of the optical-Kerr-effect spectra

The combination of optical-Kerr-effect (OKE) spectroscopy and molecular dynamics simulations has provided us with a newfound ability to delve into the librational dynamics of liquids, revealing, in the process, some surprising commonalities among aromatic liquids. Benzene and biphenyl, for example, have remarkably similar OKE spectra despite marked differences in their shapes, sizes, and moments of inertia--and even more chemically distinct aromatics tend to have noticeable similarities in their spectra. We explore this universality by using a molecular dynamics simulation to investigate the librational dynamics of molten biphenyl and to predict its OKE spectrum, comparing the results with our previous calculations for liquid benzene. We suggest that the impressive level of quantitative agreement between these two liquids is largely a reflection of the fact that librations in these and other aromatic liquids act as torsional oscillations with oscillator frequencies selected from the liquid's librational bands. Since these bands are centered about the librational Einstein frequencies, the quantitative similarities between the liquids are essentially reflections of the near identities of their Einstein frequencies. Why then are the Einstein frequencies themselves so insensitive to molecular details? We show that, for nearly planar molecules, mean-square torques and moments of inertia tend to scale with molecular dimensions in much the same way. We demonstrate that this near cancellation provides both a quantitative explanation of the close relationship between benzene and biphenyl and a more general perspective on the similarities seen in the ultrafast dynamics of aromatic liquids.     

Numerical simulation of combustion stability of liquid rocket engine based on chemistry dynamics

Combustion instability of O2/kerosene, O2/kerosene/hydrogen, and O2/kerosene/hydro-gen spray flame is numerically studied. The numerical results of combustion self-oscillation are consistent with the historical experiments. Hydrogen is helpful to stabilizing oxygen/hydrocarbon combustion. High gas injecting velocity of the coaxial injector would increase the combustion stability. Contrary to the former expectation, the most sensitive region for combustion instability is not where the heat releases most intensely but is the low-temperature premixed region near the injectors. According to the simulation, the technology steps, such as adding catalyzer to decrease the reaction activity energy, or improving the injector design to reduce the premixed low temperature region, would improve the combustion stability.     

Molecular simulation of protein adsorption and conformation at gas-liquid,liquid–liquid and solid–liquid interfaces

The adsorption of proteins at surfaces and interfaces is important in a wide range of industries. Understanding and controlling the conformation of adsorbed proteins at surfaces is critical to stability and function in many technological applications including foods and biomedical testing kits or sensors. Studying adsorbed protein conformation is difficult experimentally and so over the past few decades researchers have turned to computer simulation methods to give information at the atomic level on this important area. In this review we summarize some of the significant simulation work over the past four years at both fluid (liquid–liquid and gas–liquid interfaces) and solid–liquid interfaces. Of particular significance is the work on surfactant proteins such as fungal hydrophobins, ranspumin-2 from the túngara frog and the bacteria protein BslA. These have evolved unique structures impart very high surface-active properties to the molecules. A highlight is the elucidation of the clam-shell unhinging mechanism of ranspumin-2 adsorption to the gas–liquid interface that is responsible for its adsorption to and stabilization of the air bubbles in túngara frog foam nests.     

Molecular simulation of the phase behaviour of ternary fluid mixtures: the effect of a third component on vapour–liquid and liquid–liquid coexistence

The Gibbs ensemble algorithm is implemented to determine the vapour–liquid and liquid–liquid phase coexistence of dilute ternary fluid mixtures interacting via a Lennard–Jones potential. Calculations are reported for mixtures with a third component characterised by different intermolecular potential energy parameters. Comparison with binary mixture data indicates that the choice of energy parameter for the third component affects the composition range of vapour–liquid substantially. The addition of a third component lowers the energy of liquid phase while slightly increasing the energy of the vapour phase.     

The influence of annealing on electrochromic properties of Al–B–NiO thin films prepared by sol–gel

Al–B–NiO thin films were prepared using the sol–gel process and deposited on Indium tin oxide (ITO)-coated glass substrates via the dip-coating technique for the purpose of developing high performance electrochromic materials. The influence of the anneal on the structural and electrochromic properties of Al–B–NiO films is reported. Thermogravimetry (TG) and differential thermal analysis (DTA), cyclic voltammetry measurements (CV), UV spectrophotometer, atomic force microscopy (AFM) and X-ray diffraction (XRD) have been used to investigate the structural and electrochromic properties. The thickness of the films was determined by spectrophotometric analysis in 350–1,000 nm wavelength. Results showed that the Al–B–NiO thin films treated at high temperature have both the excellent electrochromic properties and good reversibility. The transmittance change (ΔT) of the film treated at 500 °C reaches still ~50% at the wavelength of 550 nm. The microstructure and the surface morphology were considered to play an important role in the electrochromic properties with different temperatures.     

Molecular dynamics investigations on Lennard–Jones systems near the gas–liquid critical point

NVT simulations on Lennard–Jones (L–J) systems near the gas–liquid critical point were performed by a direct approach. As a result, the two necessary conditions for simulating the systems in accordance with the thermodynamic limit were proposed: (i) L/ξ≳20 (L: the box-length, ξ: the correlation length), (ii) the total time of evolution, t_E>500 L–J units, for ξ≈3.5. The proposed conditions are probably very close to the sufficient ones. The influence of finite-size effects on pressure and density of small systems was qualitatively predicted. The prediction was confirmed by the simulations but only for L markedly lower than the length of typical critical wave, 2πξ. For L markedly higher, the evolutions were dominated by an effect called here the instability effect. The effect became negligible just when the condition for L/ξ was fulfilled. The ξ₀′ constant for L–J fluid was estimated from direct measurements of ξ to be 0.27±0.02 (L–J units). The thermodynamic parameters of the critical point, obtained from extrapolation, were in agreement with the results of other authors. The β_C exponent was estimated from minimization for a high range of temperatures to be 0.346. A comparison of the efficiency of NVT and NpT methods was also performed and no distinct differences were noted.     

Molecular dynamics simulation study on zwitterionic structure to maintain the normal conformations of Glutathione

Molecular dynamics simulations were applied to normal conformational Glutathione (GSH) and GSH over zwitterionic and hydrophobic surfaces respectively. Conformational analysis of GSH during the simulation time on RMSD, conformational flexibility and dihedral distribution were performed. The re- sults showed that zwitterionic structure maintains the normal conformations of GSH to a better extent, which should be a first good proof of the hypothesis of "maintain of normal structure".     

Formation of Zn–Al layered double hydroxide thin films intercalated with sulfonated spiropyran

We prepared Zn–Al layered double hydroxide (LDH) thin films intercalated with sulfonated 1,3′,3′-trimethyl-6-nitrospiro[2H-chromene-2,2′-indoline] anions (SP-SO₃ ^?) by immersion of sol–gel derived amorphous Al₂O₃–ZnO thin films in hot water containing SP-SO₃H. Extended interlayer spacing, in comparison to the Zn–Al LDH with carbonate anions, was observed after immersion in distilled water containing SP-SO₃H at 60 °C for 30 min, indicating that we formed Zn–Al LDH films with SP-SO₃ ^? directly on glass substrates. The merocyanine form of SP-SO₃ ^? was shown by UV spectra to have stabilized in the hydroxide layers of LDH.     

Molecular dynamics simulation study of friction and diffusion of a tracer in a Lennard–Jones solvent

The friction and diffusion coefficients of a tracer in a Lennard–Jones (LJ) solvent are evaluated by equilibrium molecular dynamics simulations in a microcanonical ensemble. The solvent molecules interact through a repulsive LJ force each other and the tracer of diameter σ₂ interacts with the solvent molecules through the same repulsive LJ force with a different LJ parameter σ. Positive deviation of the diffusion coefficient D of the tracer from a Stokes–Einstein behavior is observed and the plot of 1/D versus σ₂ shows a linear behavior. It is also observed that the friction coefficient ζ of the tracer varies linearly with σ₂ in accord with the prediction of the Stokes formula but shows a smaller slope than the Stokes prediction. When the values of ratios of sizes between the tracer and solvent molecules are higher than 5 approximately, the behavior of the friction and diffusion coefficients is well described by the Einstein relation D = k _B T/ζ, from which the tracer is considered as a Brownian particle.     

Molecular dynamics study on homonuclear and heteronuclear chains of Lennard–Jones segments

We report molecular dynamics (MD) simulation data for three simulated fluids: a homopolymer with 16 tangent Lennard–Jones (LJ) segments at the reduced temperature of 1.25, an equimolar binary homopolymer fluid with eight tangent LJ segments at 15 state points, and three corresponding copolymers with equimolar segment fraction and varying segment distribution at 15 state points. We find that the compressibility factors and energies do not change as the segment distribution varies in the copolymer example. The simulation data are compared with thermodynamic perturbation theory (TPT1) calculations. The TPT1 compressibility factors compare favorably with the MD data at high reduced temperatures but differ significantly at lower temperatures.     

Electron–ion recombination in radiation tracks in liquid argon: a computer simulation study

A simulation method is proposed to model electron–ion recombination in radiation tracks in liquid argon at 87 K. The method is applied to calculate the electron escape probability in clusters of up to 20 pairs of electrons and cations that represent a fragment of the track. The results reproduce the basic features of the track recombination in liquid argon observed in experiment.     

Molecular dynamics simulation of liquid CH2Cl2 with 3×3 and 5×5 site-site interactions

A molecular dynamics simulation of liquid CH₂Cl₂ is compared with the far infrared spectrum at the same state point (293K, 1 bar). Two representations of the force field are used, a 3×3 and 5×5 site-site interaction consisting of Lennard-Jones and charge terms. The far infra-red spectrum shows unambiguously that the 5×5 representation is more realistic in the sense that it reproduces the observed spectrum more closely.     

Study on the microstructure and liquid–solid correlation of Al–Mg alloys

Based on the available structural models and theories of electrical resistivity (ER) of liquid alloys, the structure and the liquid–solid correlation of Al _(100-x) Mg_x (x = 0, 10, 20, 30, 40, 50) alloys have been qualitatively studied by measuring the ER during the heating/cooling process using the direct-current (DC) four-probe method, as well as by characterizing the solidification morphology and testing the hardness. The result shows that the ER of Al–Mg alloys increases with the increasing temperature and the Mg content; thermal state and history have an effect on the solidification structure and properties: the ER of Al–Mg alloys exhibits a lag phenomenon of structure change during the heating/cooling process. A higher heating/cooling rate contributes to the more obvious relaxation effect of ER and the more uniform structure. Furthermore, higher pouring temperature (PT) leads the melts and solidification structure to be more homogeneous, which increases the hardness.     

Investigation of immiscible alloy system of Al–Sn thin films as anodes for lithium ion batteries

The Al–Sn, which is immiscible alloy, film was prepared by e-beam deposition to explore the possibility as anode material for lithium ion batteries for the first time. The film has a complex structure with tiny Sn particles dispersed homogeneously in the Al active matrix. The diffusion coefficients of Li⁺ in these Al–Sn alloy films were determined to be 2.1–3.2 × 10⁻⁸ cm²/s by linear sweep voltammetry. The film electrode with high Al content (Al–33wt%Sn) delivered a high initial discharge capacity of 972.8 mA h g⁻¹, while the film electrode with high Sn content (Al–64wt%Sn) with an initial discharge capacity of 552 mA h g⁻¹ showed good cycle performance indicated by retaining a capacity of about 381 mA h g⁻¹ after 60 cycles. Our preliminary results demonstrate that Al–Sn immiscible alloy is a potential candidate for anodic material of lithium ion batteries.     

Influence of the nickel content on the electrocatalytic activity of thin nanostructured Co–Te–Ni–O films

The influence of nickel addition in Co–Te–O catalytic films, obtained by vacuum co-evaporation of Co, Ni, and TeO₂ on electrocatalytic activity toward oxygen reactions in alkaline media has been investigated. Bifunctional gas-diffusion oxygen electrodes were prepared by direct deposition of catalyst films on gas-diffusion membranes consisting of hydrophobized carbon blacks. The method used allows the deposition of nanostructured films consisting of intertwined nanowires with high surface area. Thus, obtained electrodes with different atomic ratio R _(Co+Ni)/Te of the catalyst, fresh and thermally treated at 100 °C temperatures were electrochemically tested by means of cyclic voltammetry and steady-state voltammetry. It has been shown that the partial replacement of Co with about 30 at.% Ni leads to the increase in the film catalytic activity toward oxygen evolution reaction.     

Lyotropic liquid crystal guest–host material and anisotropic thin films for optical applications

We propose a new approach for the production of thin film optical functional materials. The method is based on molecular design whereby two different types of lyotropic liquid crystals (LC), lyotropic LC based on columnar supramolecules and water-soluble rod-like polymer molecules are mixed. The resulting lyotropic guest–host system allows production of optical retardation films with tunable optical anisotropy controlled by composition of the guest–host system. Coatable retarders can be used in modern liquid crystal displays and TVs for optical compensation and enhancement of the LCD's performance.