We study the excess free energy due to phase coexistence of fluids by Monte Carlo simulations using successive umbrella sampling in finite L×L×L boxes with periodic boundary conditions. Both the vapor-liquid phase coexistence of a simple Lennard-Jones fluid and the coexistence between A-rich and B-rich phases of a symmetric binary (AB) Lennard-Jones mixture are studied, varying the density ρ in the simple fluid or the relative concentration x(A) of A in the binary mixture, respectively. The character of phase coexistence changes from a spherical droplet (or bubble) of the minority phase (near the coexistence curve) to a cylindrical droplet (or bubble) and finally (in the center of the miscibility gap) to a slablike configuration of two parallel flat interfaces. Extending the analysis of Schrader et al., [Phys. Rev. E 79, 061104 (2009)], we extract the surface free energy γ(R) of both spherical and cylindrical droplets and bubbles in the vapor-liquid case and present evidence that for R→∞ the leading order (Tolman) correction for droplets has sign opposite to the case of bubbles, consistent with the Tolman length being independent on the sign of curvature. For the symmetric binary mixture, the expected nonexistence of the Tolman length is confirmed. In all cases and for a range of radii R relevant for nucleation theory, γ(R) deviates strongly from γ(∞) which can be accounted for by a term of order γ(∞)/γ(R)-1∝R(-2). Our results for the simple Lennard-Jones fluid are also compared to results from density functional theory, and we find qualitative agreement in the behavior of γ(R) as well as in the sign and magnitude of the Tolman length.

We report comprehensive simulations of the critical dynamics of a symmetric binary Lennard-Jones mixture near its consolute point. The self-diffusion coefficient exhibits no detectable anomaly. The data for the shear viscosity and the mutual-diffusion coefficient are fully consistent with the asymptotic power laws and amplitudes predicted by renormalization-group and mode-coupling theories provided finite-size effects and the background contribution to the relevant Onsager coefficient are suitably accounted for. This resolves a controversy raised by recent molecular simulations.

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grand-canonical Monte Carlo (SGMC) and molecular dynamics (MD) methods near a liquid-liquid critical temperature T(c). Choosing equal chemical potentials for the two species, the SGMC switches identities (A-->B-->A) to generate well-equilibrated configurations of the system on the coexistence curve for TT(c). A finite-size scaling analysis of the concentration susceptibility above T(c) and of the order parameter below T(c) is performed, varying the number of particles from N=400 to 12 800. The data are fully compatible with the expected critical exponents of the three-dimensional Ising universality class. The equilibrium configurations from the SGMC runs are used as initial states for microcanonical MD runs, from which transport coefficients are extracted. Self-diffusion coefficients are obtained from the Einstein relation, while the interdiffusion coefficient and the shear viscosity are estimated from Green-Kubo expressions. As expected, the self-diffusion constant does not display a detectable critical anomaly. With appropriate finite-size scaling analysis, we show that the simulation data for the shear viscosity and the mutual diffusion constant are quite consistent both with the theoretically predicted behavior, including the critical exponents and amplitudes, and with the most accurate experimental evidence.

We use inelastic neutron scattering and molecular dynamics (MD) simulation to investigate the chemical short range order (CSRO), visible through prepeaks in the structure factors, and its relation to self diffusion in Al-Ni melts. As a function of composition at 1795 K Ni self diffusion coefficients from experiment and simulation exhibit a non-linear dependence with a pronounced increase on the Al-rich side. This comes along with a change in CSRO with increasing Al content that is related to a more dense packing of the atoms in Ni-rich Al-Ni systems.

Nonequilibrium dynamics in Ising and Ginzburg-Landau models were studied for a nonconserved order parameter that mimics ordering in ferromagnets. The focus was on the understanding of the decay of the two time (t, t(w); t > tw) order-parameter correlation function. For this quantity, a full form has been obtained empirically which, for t ≫ t(w), provides a power-law ∼ (ℓ/ℓ(w))(-λ), ℓ and ℓ(w) being the characteristic lengths at t and tw, respectively. This empirical form was used for a finite-size scaling analysis to obtain the exponent λ in space dimensions d = 2 and 3. Our estimates of λ and understanding of the finite-size effects, for the models considered, provide useful information on the relevance of thermal noise. The values of λ obtained are in good agreement with the predictions of a theory based on Gaussian auxiliary field ansatz.

A combination of experimental techniques and molecular dynamics (MD) computer simulation is used to investigate the diffusion dynamics in Al80Ni20 melts. Experimentally, the self-diffusion coefficient of Ni is measured by the long-capillary (LC) method and by quasielastic neutron scattering. The LC method yields also the interdiffusion coefficient. Whereas the experiments were done in the normal liquid state, the simulations provided the determination of both self-diffusion and interdiffusion constants in the undercooled regime as well. The simulation results show good agreement with the experimental data. In the temperature range 3000 K≥ T ≥ 715 K, the interdiffusion coefficient is larger than the self-diffusion constants. Furthermore the simulation shows that this difference becomes larger in the undercooled regime. This result can be refered to a relatively strong temperature dependence of the thermodynamic factor Φ, which describes the thermodynamic driving force for interdiffusion. The simulations also indicate that the Darken equation is a good approximation, even in the undercooled regime. This implies that dynamic cross correlations play a minor role for the temperature range under consideration.

We present results from a study of finite-size effect in the kinetics of domain growth with conserved order parameter for a critical quench. Our observation of a weak size effect is a significant and surprising result. For diffusive dynamics, appropriate scaling analysis of Monte Carlo results obtained for small systems using a two-dimensional Ising model also shows that the correction to the expected Lifshitz-Slyozov law for the domain growth is very small. The methods used in this work to understand the growth dynamics should find application in other nonequilibrium systems with increasing length scales.

In Ce-doped RMnO 3 the rare earth (R) is partially replaced with cerium ions instead of divalent alkalineearth metals. These compounds are very sensitive to annealing. The resistivity and thermoelectric power ͑TEP, S) of R 0.7 Ce 0.3 MnO 3 ͑RϭLa, Pr, and Nd͒ have been studied as a function of annealing conditions. They show a ferromagnetic to paramagnetic transition, and a resistive peak near T c , marking a metal-insulator ͑MI͒ transition. A pronounced giant magnetoresistance effect is also observed in these compounds, and T c can be changed by annealing in various atmospheres. The (T) curve of oxygen overdoped samples shows double peaks along with a hysteresis loop ͑in the heating and cooling cycles͒, which bears the signature of the first-order nature of the transition. The thermopower of as-prepared samples show a peak near the MI transition, and changes sign from positive to negative at temperatures below the ferromagnetic transition. The magnitude of S at low temperatures is enhanced to an abnormally large value. With increasing of oxygen content, several characteristic features develop in the temperature dependence of the thermopower which is consistent with the resistive behavior. At high temperatures, the resistivity and thermopower follow the predictions of the Emin-Holstein theory of adiabatic polaron hopping.