Source author record

Luis G. MacDowell

Luis G. MacDowell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Analytical theory for the crossover from retarded to non-retarded interactions between metal plates

The van der Waals force established between two surfaces plays a central role in many phenomena, such as adhesion or friction. However, the dependence of this forces on the distance of separation between plates is very complex. Two widely different non-retarded and retarded regimes are well known, but these have been traditionally studied separately. Much less is known about the important experimentally accesible cross-over regime. In this study, we provide analytical approximations for the van der Waals forces between two plates that interpolates exactly between the short distance and long distance behavior, and provides new insight into the crossover from London to Casimir forces at finite temperature. At short distance, where the behavior is dominated by non-retarded interactions, we work out a very accurate simplified approximation for the Hamaker constant which adopts analytical form for both the Drude and Lorentz models of dielectric response. We apply our analytical expressions for the study of forces between metallic plates, and observe very good agreement with exact results from numerical calculations. Our results show that contributions of interband transitions remain important in the experimentally accessible regime of decades nm for several metals, including gold.

preprint2021arXiv

How ice grows from premelting films and water droplets

Close to the triple point, the surface of ice is covered by a thin liquid layer (so-called quasi-liquid layer) which crucially impacts growth and melting rates. Experimental probes cannot observe the growth processes below this layer, and classical models of growth by vapor deposition do not account for the formation of premelting films. Here, we develop a mesoscopic model of liquid-film mediated ice growth, and identify the various resulting growth regimes. At low saturation, freezing proceeds by terrace spreading, but the motion of the buried solid is conveyed through the liquid to the outer liquid-vapor interface. At higher saturations water droplets condense, a large crater forms below, and freezing proceeds undetectably beneath the droplet. Our approach is a general framework that naturally models freezing close to three phase coexistence and provides a first principle theory of ice growth and melting which may prove useful in the geosciences.

preprint2020arXiv

Rounded Layering Transitions on the Surface of Ice

Understanding the wetting properties of premelting films requires knowledge of the film's equation of state, which is not usually available. Here we calculate the disjoining pressure curve of premelting films, and perform a detailed thermodynamic characterization of premelting behavior on ice. Analysis of the density profiles reveals the signature of weak layering phenomena, from one to two and from two to three water molecular layers. However, disjoining pressure curves, which closely follow expectations from a renormalized mean field liquid state theory, show that there are no layering phase transitions in the thermodynamic sense along the sublimation line. Instead, we find that transitions at mean field level are rounded due to capillary wave fluctuations. We see signatures that true first order layering transitions could arise at low temperatures, for pressures between the metastable line of water/vapor coexistence and the sublimation line. The extrapolation of the disjoining pressure curve above water vapor saturation displays a true first order phase transition from a thin to a thick film consistent with experimental observations.

preprint2020arXiv

Structure and fluctuations of the premelted liquid film of ice at the triple point

In this paper we study the structure of the ice/vapor interface in the neighborhood of the triple point for the TIP4P/2005 model. We probe the fluctuations of the ice/film and film/vapor surfaces that separate the liquid film from the coexisting bulk phases at basal, primary prismatic and secondary prismatic planes. The results are interpreted using a coupled sine Gordon plus Interface Hamiltonian model. At large length-scales, the two bounding surfaces are correlated and behave as a single complex ice/vapor interface. For small length, on the contrary, the ice/film and film/vapor surfaces behave very much like independent ice/water and water/vapor interfaces. The study suggests that the basal facet of the TIP4P/2005 model is smooth, the prismatic facet is close to a roughening transition, and the secondary prismatic facet is rough. For the faceted basal face, our fluctuation analysis allows us to estimate the step free energy in good agreement with experiment. Our results allow for a quantitative characterization of the extent to which the adsorbed quasi-liquid layer behaves as water, and explains experimental observations which reveal similar activation energies for crystals grown in bulk vapor or bulk water.

preprint2020arXiv

Structure and water attachment rates of ice in the atmosphere: role of nitrogen

In this work we perform computer simulations of the ice surface in order to elucidate the role of nitrogen in the crystal growth rates and crystal habits of snow in the atmosphere. In pure water vapor at temperatures typical of ice crystal formation in cirrus clouds, we find that basal and primary prismatic facets exhibit a layer of premelted ice, with thickness in the subnanometer range. For partial pressures of 1 bar, well above the expected values in the troposphere, we find that only small amounts of nitrogen are adsorbed. The adsorption takes place onto the premelted surface, and hardly any nitrogen dissolves within the premelting film. The premelting film thickness does not change either. We quantify the resulting change of the ice/vapor surface tension to be in the hundredth of mN/m and find that the structure of the pristine ice surface is not changed in a significant manner. We perform a trajectory analysis of colliding water molecules, and find that the attachment rates from direct ballistic collision are very close to unity irrespective of the nitrogen pressure. Nitrogen is however at sufficient density to deflect a fraction of trajectories with smaller distance than the mean free path. Our results show explicitly that the reported differences in growth rates measured in pure water vapor and a controlled nitrogen atmosphere are not related to a significant disruption of the ice surface due to nitrogen adsorption. On the contrary, we show clearly from our trajectory analysis that nitrogen slows down the crystal growth rates due to collisions between water molecules with bulk nitrogen gas. This clarifies the long standing controversy of the role of inert gases on crystal growth rates and demonstrates their influence is solely related to the diffusion limited flow of water vapor across the gas phase.

preprint2020arXiv

Surface phase transitions and crystal growth rates of ice in the atmosphere

With climate modeling predicting a raise of at least 2 C by year 2100, the fate of ice has become a serious concern, but we still do not understand how ice grows (or melts). In the atmosphere, crystal growth rates of basal and prismatic facets exhibit an enigmatic temperature dependence, and crossover up to three times in a range between 0 and -40 C. Here we use large scale computer simulations to characterize the ice surface and identify a sequence of novel phase transitions on the main facets of ice crystallites. Unexpectedly, we find that as temperature is increased, the crystal surface transforms from a disordered phase with proliferation of steps, to a smooth phase with small step density. This causes the anomalous increase of step free energies and provides the long sought explanation for the enigmatic crossover of snow crystal growth rates found in the atmosphere.

preprint2016arXiv

Premelting-Induced Smoothening of the Ice-Vapor Interface

We perform computer simulations of the quasiliquid layer of ice formed at the ice-vapor interface close to the ice Ih-liquid-vapor triple point of water. Our study shows that the two distinct surfaces bounding the film behave at small wavelengths as atomically rough and independent ice-water and water-vapor interfaces. For long wavelengths, however, the two surfaces couple, large scale parallel fluctuations are inhibited, and the ice-vapor interface becomes smooth. Our results could help explain the complex morphology of ice crystallites.

preprint2014arXiv

A study of the ice-water interface using the TIP4P/2005 water model

In this work we study the ice-water interface under coexistence conditions by means of molecular simulations using the TIP4P/2005 water model. Following the methodology proposed by Hoyt and co-workers [J. J. Hoyt, M. Asta and A. Karma, Phys. Rev. Lett., 86, 5530, (2001)] we measure the interfacial free energy of ice with liquid water by analysing the spectrum of capillary fluctuations of the interface. We get an orientationally averaged interfacial free energy of 27(2) mN/m, in good agreement with a recent estimate obtained from simulation data of the size of critical clusters [E. Sanz, C. Vega, J. R. Espinosa, R. Caballero-Bernal, J. L. F. Abascal and C. Valeriani, JACS, 135, 15008, (2013)]. We also estimate the interfacial free energy of different planes and obtain 27(2), 28(2)and 28(2) mN/m for the basal, the primary prismatic and the secondary prismatic planes respectively. Finally, we inspect the structure of the interface and find that its thickness is of approximately 4-5 molecular diameters. Moreover, we find that when the basal plane is exposed to the fluid the interface alternates regions of cubic ice with regions of hexagonal ice.

preprint2014arXiv

Computer simulation study of surface wave dynamics at the crystal--melt interface

We study, by means of computer simulations, the crystal-melt interface of three different systems: hard-spheres, Lennard Jones and the TIP4P/2005 water model. In particular, we focus on the dynamics of surface waves. We observe that the processes involved in the relaxation of surface waves are characterized by distinct time scales: a slow one related to the continuous recrystallization and melting, that is governed by capillary forces; and a fast one which we suggest to be due to a combination of processes that quickly cause small perturbations to the shape of the interface (like e. g. Rayleigh waves, subdiffusion, or attachment/detachment of particles to/from the crystal). The relaxation of surface waves becomes dominated by the slow process as the wavelength increases. Moreover, we see that the slow relaxation is not influenced by the details of the microscopic dynamics. In a time scale characteristic for the diffusion of the liquid phase, the relaxation dynamics of the crystal-melt interface of water is around one order of magnitude slower than that of Lennard Jones or hard spheres, which we ascribe to the presence of orientational degrees of freedom in the water molecule. Finally, we estimate the rate of crystal growth from our analysis of the capillary wave dynamics and compare it with previous simulation studies and with experiments for the case of water.

preprint2014arXiv

Damped reaction field method and the accelerated convergence of the real space Ewald summation

In this paper we study a general theoretical framework which allows to approximate the real space Ewald sum by means of effective force shifted screened potentials, together with a self term. Using this strategy it is possible to generalize the reaction field method, as a means to approximate the real space Ewald sum. We show that this method exhibits faster convergence of the Coulomb energy than several schemes proposed recently in the literature while enjoying a much more sound and clear electrostatic significance. In terms of the damping parameter of the screened potential, we are able to identify two clearly distinct regimes of convergence. Firstly, a reaction field regime corresponding to the limit of small screening, where effective pair potentials converge faster than the Ewald sum. Secondly, an Ewald regime, where the plain real space Ewald sum converges faster. Tuning the screening parameter for optimal convergence occurs essentially at the crossover. The implication is that effective pair potentials are an alternative to the Ewald sum only in those cases where optimization of the convergence error is not possible.

preprint2014arXiv

Disjoining Pressure and the Film-Height-Dependent Surface Tension of Thin Liquid Films: New Insight from Capillary Wave Fluctuations

In this paper we review simulation and experimental studies of thermal capillary wave fluctuations as an ideal means for probing the underlying disjoining pressure and surface tensions, and more generally, fine details of the Interfacial Hamiltonian Model. We discuss recent simulation results that reveal a film-height-dependent surface tension not accounted for in the classical Interfacial Hamiltonian Model. We show how this observation may be explained bottom-up from sound principles of statistical thermodynamics and discuss some of its implications.

preprint2013arXiv

Capillary Fluctuations and Film-Height-Dependent Surface Tension of an Adsorbed Liquid Film

Our understanding of both structure and dynamics of adsorbed liquids heavily relies on the capillary wave Hamiltonian, but a thorough test of this model is still lacking. Here we study the capillary wave fluctuations of a liquid film with short range forces adsorbed on a solid exhibiting van der Waals interactions. We show for the first time that the measured capillary wave spectrum right above the first order wetting transition provides an interface potential consistent with independent calculations from thermodynamic integration. However, the surface tension exhibits an oscillatory film thick dependence which reveals a hitherto unnoticed capillary wave broadening mechanism beyond mere interfacial displacements.