Molecular dynamics (MD) simulations are highly attractive for studying the influence of interfacial effects, such as the enrichment of components, on the mass transfer through the interface. In a recent work, we have presented a steady-state MD simulation method for investigating this phenomenon and tested it using model mixtures with and without interfacial enrichment. The present study extends this work by introducing a nonstationary MD simulation method. A rectangular simulation box that contains a mixture of two components 1 + 2 with a vapor phase in the middle and two liquid phases on both sides is used. Starting from a vapor−liquid equilibrium state, a non-stationary molar flux of component 2 is induced by inserting particles of component 2 into the center of the vapor phase in a pulse-like manner. During the isothermal relaxation process, particles of component 2 pass through the vapor phase, cross the vapor−liquid interface, and enter the liquid phase. The system thereby relaxes into a new vapor−liquid equilibrium state. During the relaxation process, spatially resolved responses for the component densities, fluxes, and pressure are sampled. To reduce the noise and provide measures for the uncertainty of the observables, a set of replicas of simulations is carried out. The new simulation method was applied to study mass transfer in two binary Lennard-Jones mixtures: one that exhibits a strong enrichment of the low-boiling component 2 at the vapor− liquid interface and one that shows no enrichment. Even though both mixtures have similar transport coefficients in the bulk phases, the results for mass transfer differ significantly, indicating that the interfacial enrichment influences the mass transfer.