Source author record

J. S. Wettlaufer

J. S. Wettlaufer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.soft physics.ao-ph physics.flu-dyn cond-mat.mtrl-sci physics.geo-ph cond-mat.stat-mech math-ph math.MP astro-ph.EP Biological Physics nlin.CD cond-mat.other math.DS physics.chem-ph physics.data-an

Catalog footprint

What is connected

25works

15topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

Analytical solution of stochastic resonance in the nonadiabatic regime

We generalize stochastic resonance to the nonadiabatic limit by treating the double-well potential using two quadratic potentials. We use a singular perturbation method to determine an approximate analytical solution for the probability density function that asymptotically connects local solutions in boundary layers near the two minima with those in the region of the maximum that separates them. The validity of the analytical solution is confirmed numerically. Free from the constraints of the adiabatic limit, the approach allows us to predict the escape rate from one stable basin to another for systems experiencing a more complex periodic forcing.

preprint2022arXiv

Asymptotic interpretation of the Miles mechanism of wind-wave instability

When wind blows over water, ripples are generated on the water surface. These ripples can be regarded as perturbations of the wind field, which is modelled as a parallel inviscid flow. For a given wavenumber $k$, the perturbed streamfunction of the wind field and the complex phase speed are the eigenfunction and the eigenvalue of the so-called Rayleigh equation in a semi-infinite domain. Because of the small air-water density ratio, $ρ_{\rm{a}}/ρ_{\rm{w}}\equivε<<1$, the wind and the ripples are weakly coupled, and the eigenvalue problem can be solved perturbatively. At the leading order, the eigenvalue is equal to the phase speed $c_0$ of surface waves. At order $ε$, the eigenvalue has a finite imaginary part, which implies growth. Miles (1957) showed that the growth rate is proportional to the square modulus of the leading-order eigenfunction evaluated at the so-called critical level $z=z_c$, where the wind speed is equal to $c_0$ and the waves extract energy from the wind. Here, we construct uniform asymptotic approximations of the leading-order eigenfunction for long waves, which we use to calculate the growth rate as a function of $k$. In the strong wind limit, we find that the fastest growing wave is such that the aerodynamic pressure is in phase with the wave slope. The results are confirmed numerically.

preprint2022arXiv

Biolocomotion and premelting in ice

Biota are found in glaciers, ice sheets and permafrost. Ice bound micro-organisms evolve in a complex mobile environment facilitated or hindered by a range of bulk and surface interactions. When a particle is embedded in a host solid near its bulk melting temperature, a melted film forms at the surface of the particle in a process known as interfacial premelting. Under a temperature gradient, the particle is driven by a thermomolecular pressure gradient toward regions of higher temperatures in a process called thermal regelation. When the host solid is ice and the particles are biota, thriving in their environment requires the development of strategies, such as producing exopolymeric substances (EPS) and antifreeze glycoproteins (AFP) that enhance the interfacial water. Therefore, thermal regelation is enhanced and modified by a process we term {\em bio-enhanced premelting}. Additionally, the motion of bioparticles is influenced by chemical gradients influenced by nutrients within the icy host body. We show how the overall trajectory of bioparticles is controlled by a competition between thermal regelation and directed biolocomotion. By re-casting this class of regelation phenomena in the stochastic framework of active Ornstein-Uhlenbeck dynamics, and using multiple scales analysis, we find that for an attractive (repulsive) nutrient source, that thermal regelation is enhanced (suppressed) by biolocomotion. This phenomena is important in astrobiology, the biosignatures of extremophiles and in terrestrial paleoclimatology.

preprint2022arXiv

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

We use well-resolved numerical simulations to study the combined effects of buoyancy, pressure-driven shear and rotation on the melt rate and morphology of a layer of pure solid overlying its liquid phase in three dimensions at a Rayleigh number $Ra = 1.25 \times 10^5$. During thermal convection we find that the rate of melting of the solid phase varies non-monotonically with the strength of the imposed shear flow. In the absence of rotation, depending on whether buoyancy or shear dominates the flow, we observe either domes or ridges aligned in the direction of the shear flow respectively. Furthermore, we show that the geometry of the phase boundary has important effects on the magnitude and evolution of the heat flux in the liquid layer. In the presence of rotation, the strength of which is characterized by the Rossby number, Ro, we observe that for Ro = O(1) the mean flow in the interior is perpendicular to the direction of the constant horizontal applied pressure gradient. As the magnitude of this pressure gradient increases, the geometry of solid-liquid interface evolves from the voids characteristic of melting by rotating convection, to grooves oriented perpendicular or obliquely to the direction of the pressure gradient.

preprint2020arXiv

Precursors to Rare Events in Stochastic Resonance

In stochastic resonance, a periodically forced Brownian particle in a double-well potential jumps between minima at rare increments, the prediction of which poses a major theoretical challenge. Here, we use a path-integral method to find a precursor to these transitions by determining the most probable (or "{optimal}") space-time path of a particle. We characterize the optimal path using a direct comparison principle between the Langevin and Hamiltonian dynamical descriptions, allowing us to express the jump condition in terms of the accumulation of noise around the stable periodic path. In consequence, as a system approaches a rare event these fluctuations approach one of the deterministic minimizers, thereby providing a precursor for predicting a stochastic transition. We demonstrate the method numerically, which allows us to determine whether a state is following a stable periodic path or will experience an incipient jump with a high probability. The vast range of systems that exhibit stochastic resonance behavior insures broad relevance of our framework, which allows one to extract precursor fluctuations from data.

preprint2020arXiv

Transient Convective Spin-up Dynamics

We study the formation, longevity and breakdown of convective rings during impulsive spin-up in square and cylindrical containers using direct numerical simulations. The rings, which are axisymmetric alternating regions of up- and down-welling flow that can last for O (100) rotation times, were first demonstrated experimentally and arise due to a balance between Coriolis and viscous effects. We study the formation of these rings in the context of the Greenspan-Howard spin-up process, the disruption of which modifies ring formation and evolution. We show that, unless imprinted by boundary geometry, convective rings can only form when the surface providing buoyancy forcing is a free-slip surface, thereby explaining an apparent disagreement between experimental results in the literature. For Prandtl numbers from 1--5 we find that the longest-lived rings occur for intermediate Prandtl numbers, with a Rossby number dependence. Finally, we find that the constant evaporative heat-flux conditions imposed in the experiments are essential in sustaining the rings and in maintaining the vortices that form in consequence of the ring breakdown.

preprint2019arXiv

Nonadiabatic escape and stochastic resonance

We analyze the fluctuation-driven escape of particles from a metastable state under the influence of a weak periodic force. We develop an asymptotic method to solve the appropriate Fokker-Planck equation with mixed natural and absorbing boundary conditions. The approach uses two boundary layers flanking an interior region; most of the probability is concentrated within the boundary layer near the metastable point of the potential and particles transit the interior region before exiting the domain through the other boundary layer, which is near the unstable maximal point of the potential. The dominant processes in each region are given by approximate time-dependent solutions matched to construct the approximate composite solution, which gives the rate of escape with weak periodic forcing. Using reflection we extend the method to a double well potential influenced by white noise and weak periodic forcing, and thereby derive a two-state stochastic model--the simplest treatment of stochastic resonance theory--in the nonadiabatic limit

preprint2015arXiv

Maximal Stochastic Transport in the Lorenz Equations

We calculate the stochastic upper bounds for the Lorenz equations using an extension of the background method. In analogy with Rayleigh-Bénard convection the upper bounds are for heat transport versus Rayleigh number. As might be expected, the stochastic upper bounds are larger than the deterministic counterpart of \citet{Doering15}, but their variation with noise amplitude exhibits interesting behavior. Below the transition to chaotic dynamics the upper bounds increase monotonically with noise amplitude. However, in the chaotic regime this monotonicity depends on the number of realizations in the ensemble; at a particular Rayleigh number the bound may increase or decrease with noise amplitude. The origin of this behavior is the coupling between the noise and unstable periodic orbits, the degree of which depends on the degree to which the ensemble represents the ergodic set. This is confirmed by examining the close returns plots of the full solutions to the stochastic equations and the numerical convergence of the noise correlations. The numerical convergence of both the ensemble and time averages of the noise correlations is sufficiently slow that it is the limiting aspect of the realization of these bounds. Finally, we note that the full solutions of the stochastic equations demonstrate that the effect of noise is equivalent to the effect of chaos.

preprint2015arXiv

Theory of the sea ice thickness distribution

We use concepts from statistical physics to transform the original evolution equation for the sea ice thickness distribution $g(h)$ due to Thorndike et al., (1975) into a Fokker-Planck like conservation law. The steady solution is $g(h) = {\cal N}(q) h^q \mathrm{e}^{-~ h/H}$, where $q$ and $H$ are expressible in terms of moments over the transition probabilities between thickness categories. The solution exhibits the functional form used in observational fits and shows that for $h \ll 1$, $g(h)$ is controlled by both thermodynamics and mechanics, whereas for $h \gg 1$ only mechanics controls $g(h)$. Finally, we derive the underlying Langevin equation governing the dynamics of the ice thickness $h$, from which we predict the observed $g(h)$. The genericity of our approach provides a framework for studying the geophysical scale structure of the ice pack using methods of broad relevance in statistical mechanics.

preprint2014arXiv

On the Interpretation of Stratonovich Calculus

The Ito-Stratonovich dilemma is revisited from the perspective of the interpretation of Stratonovich calculus using shot noise. Over the long time scales of the displacement of an observable, the principal issue is how to deal with finite/zero autocorrelation of the stochastic noise. The former (non-zero) noise autocorrelation structure preserves the normal chain rule using a mid-point selection scheme, which is the basis Stratonovich calculus, whereas the instantaneous autocorrelation structure of Ito's approach does not. By considering the finite decay of the noise correlations on time scales very short relative to the overall displacement times of the observable, we suggest a generalization of the integral Taylor expansion criterion of Wong and Zakai (1965) for the validity of the Stratonovich approach.

preprint2013arXiv

Can planetesimals form by collisional fusion?

As a test bed for the growth of protoplanetary bodies in a turbulent circumstellar disk we examine the fate of a boulder using direct numerical simulations of particle seeded gas flowing around it. We provide an accurate description of the flow by imposing no-slip and non-penetrating boundary conditions on the boulder surface using the immersed boundary method pioneered by Peskin (2002). Advected by the turbulent disk flow, the dust grains collide with the boulder and we compute the probability density function (PDF) of the normal component of the collisional velocity. Through this examination of the statistics of collisional velocities we test the recently developed concept of collisional fusion which provides a physical basis for a range of collisional velocities exhibiting perfect sticking. A boulder can then grow sufficiently rapidly to settle into a Keplerian orbit on disk evolution time scales.

preprint2013arXiv

Grain boundary melting in ice

We describe an optical scattering study of grain boundary premelting in water ice. Ubiquitous long ranged attractive polarization forces act to suppress grain boundary melting whereas repulsive forces originating in screened Coulomb interactions and classical colligative effects enhance it. The liquid enhancing effects can be manipulated by adding dopant ions to the system. For all measured grain boundaries this leads to increasing premelted film thickness with increasing electrolyte concentration. Although we understand that the interfacial surface charge densities $q_s$ and solute concentrations can potentially dominate the film thickness, we can not directly measure them within a given grain boundary. Therefore, as a framework for interpreting the data we consider two appropriate $q_s$ dependent limits; one is dominated by the colligative effect and one is dominated by electrostatic interactions.

preprint2012arXiv

Modeling the Influence of Antifreeze Proteins on Three-Dimensional Ice Crystal Melt Shapes using a Geometric Approach

The melting of pure axisymmetric ice crystals has been described previously by us within the framework of so-called geometric crystal growth. Nonequilibrium ice crystal shapes evolving in the presence of hyperactive antifreeze proteins (hypAFPs) are experimentally observed to assume ellipsoidal geometries ("lemon" or "rice" shapes). To analyze such shapes we harness the underlying symmetry of hexagonal ice Ih and extend two-dimensional geometric models to three-dimensions to reproduce the experimental dissolution process. The geometrical model developed will be useful as a quantitative test of the mechanisms of interaction between hypAFPs and ice.

preprint2012arXiv

On the existence of stable seasonally varying Arctic sea ice in simple models

Within the framework of lower order thermodynamic theories for the climatic evolution of Arctic sea ice we isolate the conditions required for the existence of stable seasonally-varying solutions, in which ice forms each winter and melts away each summer. This is done by constructing a two-season model from the continuously evolving theory of Eisenman and Wettlaufer (2009) and showing that seasonally-varying states are unstable under constant annual average short-wave radiative forcing. However, dividing the summer season into two intervals (ice covered and ice free) provides sufficient freedom to stabilize seasonal ice. Simple perturbation theory shows that the condition for stability is determined by when the ice vanishes in summer and hence the relative magnitudes of the summer heat flux over the ocean versus over the ice. This scenario is examined within the context of greenhouse gas warming, as a function of which stability conditions are discerned.

preprint2012arXiv

Trends, noise and reentrant long-term persistence in Arctic sea ice

We examine the long-term correlations and multifractal properties of daily satellite retrievals of Arctic sea ice albedo and extent, for periods of $\sim$ 23 years and 32 years respectively. The approach harnesses a recent development called Multifractal Temporally Weighted Detrended Fluctuation Analysis (MF-TWDFA), which exploits the intuition that points closer in time are more likely to be related than distant points. In both data sets we extract multiple crossover times, as characterized by generalized Hurst exponents, ranging from synoptic to decadal. The method goes beyond treatments that assume a single decay scale process, such as a first-order autoregression, which cannot be justifiably fit to these observations. Importantly, the strength of the seasonal cycle "masks" long term correlations on time scales beyond seasonal. When removing the seasonal cycle from the original record, the ice extent data exhibits white noise behavior from seasonal to bi-seasonal time scales, whereas the clear fingerprints of the short (weather) and long ($\sim$ 7 and 9 year) time scales remain, the latter associated with the recent decay in the ice cover. Therefore, long term persistence is reentrant beyond the seasonal scale and it is not possible to distinguish whether a given ice extent minimum/maximum will be followed by a minimum/maximum that is larger or smaller in magnitude.

preprint2011arXiv

A fundamental measure theory for the sticky hard sphere fluid

We construct a density functional theory (DFT) for the sticky hard sphere (SHS) fluid which, like Rosenfeld's fundamental measure theory (FMT) for the hard sphere fluid [Phys. Rev. Lett. {\bf 63}, 980 (1989)], is based on a set of weighted densities and an exact result from scaled particle theory (SPT). It is demonstrated that the excess free energy density of the inhomogeneous SHS fluid $Φ_{\text{SHS}}$ is uniquely defined when (a) it is solely a function of the weighted densities from Kierlik and Rosinberg's version of FMT [Phys. Rev. A {\bf 42}, 3382 (1990)], (b) it satisfies the SPT differential equation, and (c) it yields any given direct correlation function (DCF) from the class of generalized Percus-Yevick closures introduced by Gazzillo and Giacometti [J. Chem. Phys. {\bf 120}, 4742 (2004)]. The resulting DFT is shown to be in very good agreement with simulation data. In particular, this FMT yields the correct contact value of the density profiles with no adjustable parameters. Rather than requiring higher order DCFs, such as perturbative DFTs, our SHS FMT produces them. Interestingly, although equivalent to Kierlik and Rosinberg's FMT in the case of hard spheres, the set of weighted densities used for Rosenfeld's original FMT is insufficient for constructing a DFT which yields the SHS DCF.

preprint2011arXiv

A Low Order Theory of Arctic Sea Ice Stability

We analyze the stability of a low-order coupled sea ice and climate model and extract the essential physics governing the time scales of response as a function of greenhouse gas forcing. Under present climate conditions the stability is controlled by longwave radiation driven heat conduction. However, as greenhouse gas forcing increases and the ice cover decays, the destabilizing influence of ice-albedo feedback acts on equal footing with longwave stabilization. Both are seasonally out of phase and as the system warms towards a seasonal ice state these effects, which underlie the bifurcations between climate states, combine exhibiting a "slowing down" to extend the intrinsic relaxation time scale from ~ 2 yr to 5 yr.

preprint2011arXiv

Decadal to seasonal variability of Arctic sea ice albedo

A controlling factor in the seasonal and climatological evolution of the sea ice cover is its albedo $α$. Here we analyze Arctic data from the Advanced Very High Resolution Radiometer (AVHRR) Polar Pathfinder and assess the seasonality and variability of broadband albedo from a 23 year daily record. We produce a histogram of daily albedo over ice covered regions in which the principal albedo transitions are seen; high albedo in late winter and spring, the onset of snow melt and melt pond formation in the summer, and fall freeze up. The bimodal late summer distribution demonstrates the combination of the poleward progression of the onset of melt with the coexistence of perennial bare ice with melt ponds and open water, which then merge to a broad peak at $α\gtrsim $ 0.5. We find the interannual variability to be dominated by the low end of the $α$ distribution, highlighting the controlling influence of the ice thickness distribution and large-scale ice edge dynamics. The statistics obtained provide a simple framework for model studies of albedo parameterizations and sensitivities.

preprint2011arXiv

Density functional theory for Baxter's sticky hard spheres in confinement

It has recently been shown that a free energy for Baxter's sticky hard sphere fluid is uniquely defined within the framework of fundamental measure theory (FMT) for the inhomogeneous hard sphere fluid, provided that it obeys scaled-particle theory and the Percus-Yevick (PY) result for the direct correlation function [Hansen-Goos and Wettlaufer, J. Chem. Phys. {\bf 134}, 014506 (2011)]. Here, combining weighted densities from common versions of FMT with a new vectorial weighted density, we derive a regularization of the divergences of the associated strongly confined limit. Moreover, the simple free energy that emerges is exact in the zero-dimensional limit, leaves the underlying equation of state unaffected, and yields a direct correlation function distinct from the PY expression. Comparison with simulation data for both the bulk pair correlation function and the density profiles in confinement shows that the new theory is significantly more accurate than the PY-based results. Finally, the resulting free energy is applicable to a glass of adhesive hard spheres.

preprint2011arXiv

Particle-scale structure in frozen colloidal suspensions from small angle x-ray scattering

During directional solidification of the solvent in a colloidal suspension, the colloidal particles segregate from the growing solid, forming high-particle-density regions with structure on a hierarchy of length scales ranging from that of the particle-scale packing to the large-scale spacing between these regions. Previous work has concentrated mostly on the medium- to large-length scale structure, as it is the most accessible and thought to be more technologically relevant. However, the packing of the colloids at the particle-scale is an important component not only in theoretical descriptions of the segregation process, but also to the utility of freeze-cast materials for new applications. Here we present the results of experiments in which we investigated this structure across a wide range of length scales using a combination of small angle x-ray scattering and direct optical imaging. As expected, during freezing the particles were concentrated into regions between ice dendrites forming a microscopic pattern of high- and low-particle-density regions. X-ray scattering indicates that the particles in the high density regions were so closely packed as to be touching. However, the arrangement of the particles does not conform to that predicted by standard inter-particle pair potentials, suggesting that the particle packing induced by freezing differs from that formed during equilibrium densification processes.

preprint2010arXiv

Accretion in Protoplanetary Disks by Collisional Fusion

The formation of a solar system is believed to have followed a multi-stage process around a protostar. Whipple first noted that planetesimal growth by particle agglomeration is strongly influenced by gas drag; there is a "bottleneck" at the meter scale with such bodies rapidly spiraling into the central star, whereas much smaller or larger particles do not. Thus, successful planetary accretion requires rapid planetesimal growth to km scale. A commonly accepted picture is that for collisional velocities $V_c$ above a certain threshold collisional velocity, ${V_{th}} \sim$ 0.1-10 cm s$^{-1}$, particle agglomeration is not possible; elastic rebound overcomes attractive surface and intermolecular forces. However, if perfect sticking is assumed for all collisions the bottleneck can be overcome by rapid planetesimal growth. While previous work has dealt explicitly with the influences of collisional pressures and the possibility of particle fracture or penetration, the basic role of the phase behavior of matter--phase diagrams, amorphs and polymorphs--has been neglected. Here it is demonstrated that novel aspects of surface phase transitions provide a physical basis for efficient sticking through collisional melting or amphorph-/polymorphization and fusion to extend the collisional velocity range of primary accretion to $ΔV_c \sim$ 1-100 m s$^{-1}$, which bound both turbulent RMS speeds and the velocity differences between boulder sized and small grains $\sim$ 1-50 m s$^{-1}$. Thus, as inspiraling meter sized bodies collide with smaller particles in this high velocity collisional fusion regime they grow rapidly to km scales and hence settle into stable Keplerian orbits in $\sim$ 10$^5$ years before photoevaporative wind clears the disk of source material.

preprint2010arXiv

Dynamics of colloidal particles in ice

We use X-ray Photon Correlation Spectroscopy (XPCS) to probe the dynamics of colloidal particles in polycrystalline ice. During freezing, the dendritic ice morphology and rejection of particles from the ice created regions of high-particle-density, where some of the colloids were forced into contact and formed disordered aggregates. We find that the particles in these high density regions underwent ballistic motion coupled with both stretched and compressed exponential decays of the intensity autocorrelation function, and that the particles' characteristic velocity increased with temperature. We explain this behavior in terms of ice grain boundary migration.

preprint2010arXiv

Optimal and hysteretic fluxes in alloy solidification: Variational principles and chimney spacing

We take a numerical approach to analyze the mechanisms controlling the spacing of chimneys -- channels devoid of solid -- in two-dimensional mushy layers formed by solidifying a binary alloy. Chimneys are the principal conduits through which buoyancy effects transport material out of the mushy layer and into the liquid from which it formed. Experiments show a coarsening of chimney spacing and we pursue the hypothesis that this observation is a consequence of a variational principle: the chimney spacing adjusts to optimize material transport and hence maximize the rate of removal of potential energy stored in the mushy layer. The optimal solute flux increases approximately linearly with the mushy layer Rayleigh number. However, for spacings below a critical value the chimneys collapse and solute fluxes cease, revealing a hysteresis between chimney convection and no flow.

preprint2010arXiv

Streaks to Rings to Vortex Grids: Generic Patterns in Transient Convective Spin-Up

We observe the transient formation of a ringed pattern state during spin-up of an evaporating fluid on a time scale of order a few Ekman spin-up times. The ringed state is probed using infrared thermometry and particle image velocimetry and it is demonstrated to be a consequence of the transient balance between Coriolis and viscous forces which dominate inertia, each of which are extracted from the measured velocity field. The breakdown of the ringed state is quantified in terms of the antiphasing of these force components which drives a Kelvin-Helmholtz instability and we show that the resulting vortex grid spacing scales with the ring wavelength. This is the fundamental route to quasi-two dimensional turbulent vortex flow and thus may have implications in astrophysics and geophysics wherein rotating convection is ubiquitous. sics.

preprint2010arXiv

Theory of ice premelting in porous media

Premelting describes the confluence of phenomena that are responsible for the stable existence of the liquid phase of matter in the solid region of its bulk phase diagram. Here we develop a theoretical description of the premelting of water ice contained in a porous matrix, made of a material with a melting temperature substantially larger than ice itself, to predict the amount of liquid water in the matrix at temperatures below its bulk freezing point. Our theory combines the interfacial premelting of ice in contact with the matrix, grain boundary melting in the ice, and impurity and curvature induced premelting, the latter occurring in regions which force the ice-liquid interface into a high curvature configuration. These regions are typically found at points where the matrix surface is concave, along contact lines of a grain boundary with the matrix, and in liquid veins. Both interfacial premelting and curvature induced premelting depend on the concentration of impurities in the liquid, which, due to the small segregation coefficient of impurities in ice are treated as homogeneously distributed in the premelted liquid. Our principal result is an equation for the fraction of liquid in the porous medium as a function of the undercooling, which embodies the combined effects of interfacial premelting, curvature induced premelting, and impurities. The result is analyzed in detail and applied to a range of experimentally relevant settings.

J. S. Wettlaufer

What is connected

Connect this record

See the researcher in context

Building this map preview

25 published item(s)

Analytical solution of stochastic resonance in the nonadiabatic regime

Asymptotic interpretation of the Miles mechanism of wind-wave instability

Biolocomotion and premelting in ice

The combined effects of buoyancy, rotation, and shear on phase boundary evolution

Precursors to Rare Events in Stochastic Resonance

Transient Convective Spin-up Dynamics

Nonadiabatic escape and stochastic resonance

Maximal Stochastic Transport in the Lorenz Equations

Theory of the sea ice thickness distribution

On the Interpretation of Stratonovich Calculus

Can planetesimals form by collisional fusion?

Grain boundary melting in ice

Modeling the Influence of Antifreeze Proteins on Three-Dimensional Ice Crystal Melt Shapes using a Geometric Approach

On the existence of stable seasonally varying Arctic sea ice in simple models

Trends, noise and reentrant long-term persistence in Arctic sea ice

A fundamental measure theory for the sticky hard sphere fluid

A Low Order Theory of Arctic Sea Ice Stability

Decadal to seasonal variability of Arctic sea ice albedo

Density functional theory for Baxter's sticky hard spheres in confinement

Particle-scale structure in frozen colloidal suspensions from small angle x-ray scattering

Accretion in Protoplanetary Disks by Collisional Fusion

Dynamics of colloidal particles in ice

Optimal and hysteretic fluxes in alloy solidification: Variational principles and chimney spacing

Streaks to Rings to Vortex Grids: Generic Patterns in Transient Convective Spin-Up

Theory of ice premelting in porous media