Molecular dynamics simulation of thermal transport across a solid/liquid interface created by a meniscus.

Understandings heat transfer across a solid/liquid interface is crucial for establishing novel thermal control pathways in a range of energy applications. One of the major problems raised in this context is the impact of the three-phase contact line between solid, liquid, and gas on heat flux perturbations at the nanoscale. The focus of this research is the thermal transport via nanosized meniscus restricted between two solid walls. The molecular dynamics approach was used to consider different wetting states of the meniscus by varying the interaction potential between atoms of the substrate and the liquid. The influence of the meniscus size on the energy exchange between two solid walls was also studied. It was discovered that possessing a three-phase contact line reduces the interfacial boundary resistance between solid and liquid. Furthermore, the finite element method was employed to connect atomistic simulations with continuum mechanics. We show that the wetting angle and interfacial boundary resistance are essential important parameters for multiscale analysis of thermal engineering issues with precise microscale parametrization.

Theory of length-scale dependent relaxation moduli and stress fluctuations in glass-forming and viscoelastic liquids.

The spatiotemporal correlations of the local stress tensor in supercooled liquids are studied both theoretically and by molecular dynamics simulations of a two-dimensional (2D) polydisperse Lennard-Jones system. Asymptotically exact theoretical equations defining the dynamical structure factor and all components of the stress correlation tensor for low wave-vector q are presented in terms of the generalized (q-dependent) shear and longitudinal relaxation moduli, G(q, t) and K(q, t). We developed a rigorous approach (valid for low q) to calculate K(q, t) in terms of certain bulk correlation functions (for q = 0), the static structure factor S(q), and thermal conductivity κ. The proposed approach takes into account both the thermostatting effect and the effect of polydispersity. The theoretical results for the (q, t)-dependent stress correlation functions are compared with our simulation data, and an excellent agreement is found for qb̄≲0.5 (with b̄ being the mean particle diameter) both above and below the glass transition without any fitting parameters. Our data are consistent with recently predicted (both theoretically and by simulations) long-range correlations of the shear stress quenched in heterogeneous glassy structures.

Simple models for strictly non-ergodic stochastic processes of macroscopic systems.

We investigate simple models for strictly non-ergodic stochastic processes [Formula: see text] (t being the discrete time step) focusing on the expectation value v and the standard deviation [Formula: see text] of the empirical variance [Formula: see text] of finite time series [Formula: see text]. [Formula: see text] is averaged over a fluctuating field [Formula: see text] ([Formula: see text] being the microcell position) characterized by a quenched spatially correlated Gaussian field [Formula: see text]. Due to the quenched [Formula: see text]-field [Formula: see text] becomes a finite constant, [Formula: see text], for large sampling times [Formula: see text]. The volume dependence of the non-ergodicity parameter [Formula: see text] is investigated for different spatial correlations. Models with marginally long-ranged [Formula: see text]-correlations are successfully mapped on shear stress data from simulated amorphous glasses of polydisperse beads.

Relaxation moduli of glass-forming systems: temperature effects and fluctuations.

Equilibrium and dynamical properties of a two-dimensional polydisperse colloidal model system are characterized by means of molecular dynamics (MD) and Monte Carlo (MC) simulations. We employed several methods to prepare quasi-equilibrated systems: in particular, by slow cooling and tempering with MD (method SC-MD), and by tempering with MC dynamics involving swaps of particle diameters (methods Sw-MD, Sw-MC). It is revealed that the Sw-methods are much more efficient for equilibration below the glass transition temperature Tg leading to denser and more rigid systems which show much slower self-diffusion and shear-stress relaxation than their counterparts prepared with the SC-MD method. The shear-stress relaxation modulus G(t) is obtained based on the classical stress-fluctuation relation. We demonstrate that the α-relaxation time τα obtained using a time-temperature superposition of G(t) shows a super-Arrhenius behavior with the VFT temperature T0 well below Tg. We also derive novel rigorous fluctuation relations providing isothermic and adiabatic compression relaxation moduli in the whole time range (including the short-time inertial regime) based on correlation data for thermostatted systems. It is also shown that: (i) the assumption of Gaussian statistics for stress fluctuations leads to accurate predictions of the variances of the fluctuation moduli for both shear (µF) and compression (ηF) at T⪆Tg. (ii) The long-time (quasi-static) isothermic and adiabatic moduli increase on cooling faster than the affine compression modulus ηA, and this leads to a monotonic temperature dependence of ηF which is qualitatively different from µF(T) showing a maximum near Tg.

General relations to obtain the time-dependent heat capacity from isothermal simulations.

It is well-known that time-dependent correlation functions related to temperature and energy can crucially depend on the thermostatting mechanism used in computer simulations of molecular systems. We argue, however, that linear response functions must be considered as universal properties of physical systems. This implies that the classical fluctuation equation for the transient heat capacity, cv(t), is not applicable to the thermostatted molecular dynamics (apart from long enough times). To improve on this point, we derive a number of exact general expressions for the frequency-dependent heat capacity in terms of energy correlation functions, valid for the Nosé-Hoover and some other thermostats. We also establish a general relation between auto- and cross correlation functions of energy and temperature. Recommendations on how to use these relations to maximize the numerical precision are provided. It is demonstrated that our approach allows us to obtain cv(t) for a supercooled liquid system with high precision and over many decades in time reflecting all pertinent relaxation processes.

Fluctuations of non-ergodic stochastic processes.

We investigate the standard deviation [Formula: see text] of the variance [Formula: see text] of time series [Formula: see text] measured over a finite sampling time [Formula: see text] focusing on non-ergodic systems where independent "configurations" c get trapped in meta-basins of a generalized phase space. It is thus relevant in which order averages over the configurations c and over time series k of a configuration c are performed. Three variances of [Formula: see text] must be distinguished: the total variance [Formula: see text] and its contributions [Formula: see text], the typical internal variance within the meta-basins, and [Formula: see text], characterizing the dispersion between the different basins. We discuss simplifications for physical systems where the stochastic variable x(t) is due to a density field averaged over a large system volume V. The relations are illustrated for the shear-stress fluctuations in quenched elastic networks and low-temperature glasses formed by polydisperse particles and free-standing polymer films. The different statistics of [Formula: see text] and [Formula: see text] are manifested by their different system-size dependences.

Ensemble fluctuations matter for variances of macroscopic variables.

Extending recent work on stress fluctuations in complex fluids and amorphous solids we describe in general terms the ensemble average [Formula: see text] and the standard deviation [Formula: see text] of the variance [Formula: see text] of time series [Formula: see text] of a stochastic process x(t) measured over a finite sampling time [Formula: see text]. Assuming a stationary, Gaussian and ergodic process, [Formula: see text] is given by a functional [Formula: see text] of the autocorrelation function h(t). [Formula: see text] is shown to become large and similar to [Formula: see text] if [Formula: see text] corresponds to a fast relaxation process. Albeit [Formula: see text] does not hold in general for non-ergodic systems, the deviations for common systems with many microstates are merely finite-size corrections. Various issues are illustrated for shear-stress fluctuations in simple coarse-grained model systems.

Composition fluctuations in polydisperse liquids: Glasslike effects well above the glass transition.

We study a two-dimensional glass-forming system of slightly polydisperse (LJ) particles using molecular dynamics simulations and demonstrate that in the liquid regime (well above the vitrification temperature) this model shows a number of features typical of the glass transition: (i) the relation between compressibility and structure factor S(q) is strongly violated; (ii) the dynamical structure factor S(q,t) at low q shows a two-step relaxation; (iii) the time-dependent heat capacity c_{v}(t) shows a long-time power-law tail. We show that these phenomena can be rationalized with the idea of composition fluctuations and provide a quantitative theory for the effects (i) and (ii). It implies that such effects must be inherent in all polydisperse colloidal models, including binary LJ mixtures.

Long-range stress correlations in viscoelastic and glass-forming fluids.

A simple and rigorous approach to obtain stress correlations in viscoelastic liquids (including supercooled liquid and equilibrium amorphous systems) is proposed. The long-range dynamical correlations of local shear stress are calculated and analyzed in 2-dimensional space. It is established how the long-range character of the stress correlations gradually emerges as the relevant dynamical correlation length l grows in time. The correlation range l is defined by momentum propagation due to acoustic waves and vorticity diffusion which are the basic mechanisms for transmission of shear stress perturbations. We obtain the general expression defining the time- and distance-dependent stress correlation tensor in terms of material functions (generalized relaxation moduli). The effect of liquid compressibility is quantitatively analyzed; it is shown to be important at large distances and/or short times. The revealed long-range stress correlation effect is shown to be dynamical in nature and unconnected with static structural correlations in liquids (correlation length ξs). Our approach is based on the assumption that ξs is small enough as reflected in weak wave-number dependencies of the generalized relaxation moduli. We provide a simple physical picture connecting the elucidated long-range fluctuation effect with anisotropic correlations of the (transient) inherent stress field, and discuss its implications.