Density functional study of the pressure tensor for inhomogeneous Lennard—Jones fluids

Based on classical density functional theory,an expression of the pressure tensor for inhomogeneous fluids is presented.This takes into account greater correlation between particles,especially for systems that are geometrically confined or involve an interface.The density and pressure components of Lennard-Jones fluids confined in hard and softened nano-cavities are calculated.A comparison between the results of this work and IK expression suggests that the agreement depends on temperature.The interfacial tension for hard sphere fluids agrees well with the Monte Carlo result when the bulk density is not too large.The results of the solid-fluid interfacial tension for Lennard-Jones fluids demonstrate that different types of external potentials modulate the interfacial tension in different manners.     

Influence of precursor feeding rate on vapor–liquid–solid nanowire growth

The yield of Ge nanowires (NWs) synthesized using the vapor–liquid–solid (VLS) method was discovered to be highly sensitive to the rate of precursor feeding. When other parameters were fixed, fast filling of precursors yielded nearly 100% Ge NWs with regard to the growth seeds. By contrast, slow feeding produced nearly no or very low yield of Ge NWs. The dramatic difference was attributed to a layer of Ge coating on the surface of growth seeds. The coating formed at relatively low precursor pressures as a result of the imbalance in the VLS process. The results shed new light on the VLS mechanism in general. Electronic Supplementary Material  The online version of this article () contains supplementary material, which is available to authorized users.     

Loci of extrema of thermodynamic response functions for the Lennard–Jones fluid

Lines of extrema along isotherms and isobars for the residual isochoric heat capacity, the residual isobaric heat capacity, the isobaric thermal expansivity, and the isothermal compressibility of the Lennard–Jones fluid have been studied from popular equations of state due to Johnson et al. [Mol. Phys. 78, 591 (1993)], Kolafa and Nezbeda [Fluid Phase Equilib. 100, 1 (1994)], and Mecke et al. [Int. J. Thermophys. 17, 391 (1996)]. On depicting such loci in the pressure–temperature plane, the characteristic behaviour of thermodynamic response functions in the ideal-gas limit (at high enough temperatures or low enough pressures), the close-packed-fluid limit (at low enough temperatures or high enough pressures) as well as in the liquid, critical, and supercritical regions is identified. The present analysis is informative itself, but it also stimulates further work in order to tackle more complicated cases of study including associated fluids.     

Vapour–liquid phase equilibrium and surface tension of fully flexible Lennard–Jones chains

We report measurements of the vapour–liquid coexistence densities and surface tension of fully flexible Lennard–Jones chain molecules ranging in length from 4 to 60 beads. We demonstrate that the surface tension for all chain lengths collapses to a single master curve when plotted according to the universal parachor correlation. We find a universal parachor exponent 3.79 ± 0.05 for conditions close to the critical point, with a deviation observed for the longest chains far below the critical point.     

Comprehensive study of the vapour–liquid coexistence of the truncated and shifted Lennard–Jones fluid including planar and spherical interface properties

Vapour–liquid equilibria of the Lennard–Jones potential, truncated and shifted at 2.5σ, are studied using molecular dynamics simulations, an attractive option for studying inhomogeneous systems. Comprehensive simulation data are reported for three cases: no interface, a planar interface, and a spherical interface between the coexisting phases, covering a wide range of temperatures. Spherical droplets are also studied for a range of radii between 5 and 16σ. The size dependence of the surface tension, based on the Irving–Kirkwood pressure tensor, and other properties is quantified for spherical interfaces. All simulation results are correlated with a consistent set of empirical equations. A comparison with the results of other authors as well as with experimental data for noble gases and methane is also presented.     

Ultrafast solid–liquid–vapor phase change of a gold film induced by pico- to femtosecond lasers

Melting, vaporization and resolidification processes of thin gold film irradiated by a femtosecond pulse laser are studied numerically. The nonequilibrium heat transfer in electrons and lattice is described using a two-temperature model. The solid–liquid interfacial velocity, as well as elevated melting temperature and depressed solidification temperature, is obtained by considering the interfacial energy balance and nucleation dynamics. An iterative procedure based on energy balance and gas kinetics law to track the location of liquid–vapor interface is utilized to obtain the material removal by vaporization. The effect of surface heat loss by thermal radiation was discussed. The influences of laser fluence and duration on the evaporation process are studied. Results show that higher laser fluence and shorter laser pulse width lead to higher interfacial temperature, deeper melting and ablation depths.     

Spontaneous translation of nanodroplet over a heterogeneous surface due to thermal cycles: role of solid–liquid interfacial fluctuations

ABSTRACT

We study the molecular-scale features of the solid surface that result in the spontaneous motion of a nanodroplet due to the periodic variation of temperature. We first employ a thermodynamic model to predict the variation of solid–fluid interfacial properties that can result in the above motion. The model identifies a composite (surface couple) made of two surfaces that are characterised by a large difference between the entropic parts of the solid–liquid interfacial free energies. In order to understand the molecular-scale features of the two surfaces that may form a surface couple, we performed grand canonical Monte Carlo simulations of Lennard Jones fluid and crystalline surfaces made of Lennard Jones-like atoms. We then used the cumulant expansions of the perturbation formulas to divide the interfacial entropy into two parts: The one that is directly affected by the solid–fluid attraction (direct part), and the other (indirect part) that is indirectly affected by the solid–fluid attraction via the alteration of interfacial fluctuations. Our results indicate that two surfaces form a surface couple if the differences between their chemical natures lead to large differences in the indirect part of the interfacial entropy, while the direct part remains relatively unaffected.     

A critical behavior of the Lennard–Jones dimeric fluid in two-dimensions.: A Monte Carlo study

Monte Carlo simulation with hyper-parallel tempering, the histogram reweighting and finite-size scaling techniques are used to explore the critical behavior of Lennard-Jones fluids of dimers in the two-dimensional system. The dimers are built of two distinct segments. The critical parameters for selected types of dimers are estimated. We analyze the influence of the bond length and differences in the chemical nature of segments on the critical parameters. The results can be used to predict the critical parameters of dimers whose molecular parameters are known.     

Li adsorption on a Fullerene–Single wall carbon nanotube hybrid system: Density functional theory approach

In this study, we investigate Li adsorption mechanisms on the C₆₀-SWCNT hybrid system using density functional theory. It is found that the Li adsorption energy of the C₆₀-SWCNT hybrid system is increased in comparison to that of the pure SWCNT. The Li adsorption energy ranges from −1.917 eV to −2.642 eV for the single-Li adsorbed system and from −2.351 eV to −2.636 eV for the double-Li adsorbed system. It is also found that the adsorption energy becomes similar at most positions throughout the structure. In addition, the Li adsorption energy of 31-Li system is calculated to be −1.863 eV, which is significantly lower than the Li–Li binding energy (−1.030 eV). These results infer that Li atoms will be adsorbed on the space 1) between C₆₀ and C₆₀; 2) between SWCNT and C₆₀; 3) the rest of the space (e.g. between SWCNTs), rather than form Li clusters. As more Li atoms are adsorbed onto the C₆₀-SWCNT hybrid system due to such improved Li adsorption capability, the metallic character of the system is enhanced, which is confirmed via the band structure and electronic density of states.     

Transition region width of nanowire hetero- and pn-junctions grown using vapor–liquid–solid processes

The transition region width of nanowire heterojunctions and pn-junctions grown using vapor–liquid–solid (VLS) processes has been modeled. With two constituents or dopants I and II, the achievable width or abruptness of the junctions is attributed to the residual I atom/molecule stored in the liquid droplet at the onset of introducing II to grow the junction, and the stored I atom/molecule consumption into the subsequently grown crystal layers. The model yields satisfactory quantitative fits to a set of available Si-Ge junction data. Moreover, the model provides a satisfactory explanation to the relative junction width or abruptness differences between elemental and compound semiconductor junction cases, as well as a guideline for achieving the most desirable pn-junction widths. PACS 81.07.-b; 64.75.Jk; 61.46.Km     

Density functional theory study on the structure and properties of Si3N4 clusters

The hybrid density functional theory B3LYP with basis sets 6-31G＊ has been used to study on the equilibrium geometries and electronic structures of possible isomers of Si3N4 clusters. 24 possible isomers are obtained. The most stable isomer of Si3N4 is a 3D structure with 7 Si-N bonds and 2 N-N bonds that could beformed by 3 quadrangles. The bond properties of the most stable isomer was analyzed by using natural bond orbital method （NBO）, the results suggest that the charges on Si and N atoms in Si-N bonds are quite large, so theinteraction of N-Si atoms in Si3N4 cluster is of strongly electric interaction. The primary IR and Raman vibrational frequency located at 1033.40 cm^-1, 473.63 cm^-1 respectively. The polarizabilities and hyperpolarizabilities of the most stable isomer are also analyzed.     

Mechanism from particle compaction to fluidization of liquid–solid two-phase flow

A new model of particle yield stress including cohesive strength is proposed, which considers the friction and cohesive strength between particles. A calculation method for the fluidization process of liquid–solid two-phase flow in compact packing state is given, and the simulation and experimental studies of fluidization process are carried out by taking the sand–water two-phase flow in the jet dredging system as an example, and the calculation method is verified.     

Effect of solid-state characteristics on the critical parameters of the vapor–liquid phase transition

Database for the critical point parameters of almost all metals (including transition metals) and semiconductors is used to derive a number of empirical expressions to relate these parameters to the heat of evaporation, the normal density, and the isothermal bulk modulus of these substances in a solid state under normal conditions. The database is obtained using the thermodynamic model proposed earlier.     

Density functional theory study of Mg_2Ni_n(n= 1–8) clusters

The density functional theory B3PW91 with LANL2DZ basis sets has been used to study the possible geometries of Mg2Nin(n = 1–8) clusters. For the lowest energy structures of the clusters, stabilities, electronic properties, and natural bond orbital(NBO) are calculated and discussed. The results show that the doped Mg atoms reduce the stabilities of pure Ni clusters. The Mg2Ni2, Mg2Ni4, and Mg2Ni6clusters are more stable than neighboring clusters. The system appears magic number characteristics. In addition, the hybridization phenomenon occurs, owing to the interaction of Mg and Ni. The result of charge transfer is that Ni atom is negative and the Mg atom is positive. We also conclude that the 3p and 4d orbitals of the Ni atom have an effect on the stabilities of the clusters.     

Flows in Complex Networks: Theory,Algorithms, and Application to Lennard–Jones Cluster Rearrangement

A set of analytical and computational tools based on transition path theory (TPT) is proposed to analyze flows in complex networks. Specifically, TPT is used to study the statistical properties of the reactive trajectories by which transitions occur between specific groups of nodes on the network. Sampling tools are built upon the outputs of TPT that allow to generate these reactive trajectories directly, or even transition paths that travel from one group of nodes to the other without making any detour and carry the same probability current as the reactive trajectories. These objects permit to characterize the mechanism of the transitions, for example by quantifying the width of the tubes by which these transitions occur, the location and distribution of their dynamical bottlenecks, etc. These tools are applied to a network modeling the dynamics of the Lennard–Jones cluster with 38 atoms ( \(\mathrm{LJ}_{38}\) ) and used to understand the mechanism by which this cluster rearranges itself between its two most likely states at various temperatures.     

The Mayer Series of the Lennard–Jones Gas: Improved Bounds for the Convergence Radius

We provide a lower bound for the convergence radius of the Mayer series of the Lennard–Jones gas which strongly improves on the classical bound obtained by Penrose and Ruelle 1963. To obtain this result we use an alternative estimate recently proposed by Morais et al. (J Stat Phys 2014) for a restricted class of stable and tempered pair potentials (namely those which can be written as the sum of a non-negative potential plus an absolutely integrable and stable potential) combined with a method developed by Locatelli and Schoen (J Glob Optim 22, 175–190 2002) for establishing a lower bound for the minimal interatomic distance between particles interacting via a Morse potential in a cluster of minimum-energy configurations.