Bubble raft

Last updated November 09, 2022

A bubble raft is an array of bubbles. It demonstrates materials' microstructural and atomic length-scale behavior by modelling the {111} plane of a close-packed crystal. A material's observable and measurable mechanical properties strongly depend on its atomic and microstructural configuration and characteristics. This fact is intentionally ignored in continuum mechanics, which assumes a material to have no underlying microstructure and be uniform and semi-infinite throughout.

Bubble rafts assemble bubbles on a water surface, often with the help of amphiphilic soaps. These assembled bubbles act like atoms, diffusing, slipping, ripening, straining, and otherwise deforming in a way that models the behavior of the {111} plane of a close-packed crystal. The ideal (lowest energy) state of the assembly would undoubtedly be a perfectly regular single crystal, but just as in metals, the bubbles often form defects, grain boundaries, and multiple crystals.

History of bubble rafts

The concept of bubble raft modelling was first presented in 1947 by Nobel Laureate Sir William Lawrence Bragg and John Nye of Cambridge University's Cavendish Laboratory in Proceedings of the Royal Society A.^[1] Legend claims that Bragg conceived of bubble raft models while pouring oil into his lawn mower. He noticed that bubbles on the surface of the oil assembled into rafts resembling the {111} plane of close-packed crystals.^[2] Nye and Bragg later presented a method of generating and controlling bubbles on the surface of a glycerine-water-oleic acid-triethanolamine solution, in assemblies of 100,000 or more sub-millimeter sized bubbles. In their paper, they go on at length about the microstructural phenomena observed in bubble rafts and hypothesized in metals.^[1]

Dynamics

Avalanches of rupturing bubbles can give rise to self-organized criticality, similar to the Abelian sandpile model.

Bubble rafts exhibit complex dynamics, as illustrated in the video. This is triggered by rupture of a first bubble, driven by thermal fluctuations^[3] and a cascade of subsequent bursting bubbles, which can give rise to self-organized criticality, and a power-law distribution of avalanches.^[4]

Relation to crystal lattices

In deforming a crystal lattice, one changes the energy and the interatomic potential felt by the atoms of the lattice. This interatomic potential is popularly (and mostly qualitatively) modeled using the Lennard-Jones potential, which consists of a balance between attractive and repulsive forces between atoms.

The "atoms" in Bubble Rafts also exhibit such attractive and repulsive forces:^[2]

$U(\rho )=-\pi R^{4}\rho _{solution}g\left({\frac {\mathrm {B} }{\alpha }}\right)^{2}{\mathit {A}}K_{0}(\alpha \rho )+{\begin{cases}0~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~\rho \geq \ 2\\\pi R^{4}\rho _{solution}g\left({\frac {(2-\rho )^{2}}{\alpha ^{2}}}\right)~~~\rho \leq \ 2\end{cases}}$

The portion of the equation to the left of the plus sign is the attractive force, and the portion to the right represents the repulsive force.

$U(\rho )$ is the interbubble potential

$R$ is the average bubble radius

$\rho _{solution}$ is the density of the solution from which the bubbles are formed

$g$ is the gravitational constant

$\rho$ is the ratio of the distance between bubbles to the bubble radius

$\mathrm {B}$ is the radius of ring contact

$\alpha$ is the ratio R/a of the bubble radius to the Laplace constant a, where

$a^{2}={\frac {T}{\rho _{solution}g}}$

$T$ is the surface tension

$A$ is a constant dependent upon the boundary conditions of the calculation

$K_{0}$ is a zeroth-order modified Bessel function of the second kind.^[2]

Bubble rafts can display numerous phenomena seen in the crystal lattice. This includes such things as point defects (vacancies, substitutional impurities, interstitial atoms), edge dislocations and grains. A screw dislocation can't be modeled in a 2D bubble raft because it extends outside the plane. It is even possible to replicate some microstructure treats such as annealing. The annealing process is simulated by stirring the bubble raft. This anneals out the dislocations (recovery) and promotes recrystallization.

Related Research Articles

A centripetal force is a force that makes a body follow a curved path. Its direction is always orthogonal to the motion of the body and towards the fixed point of the instantaneous center of curvature of the path. Isaac Newton described it as "a force by which bodies are drawn or impelled, or in any way tend, towards a point as to a centre". In Newtonian mechanics, gravity provides the centripetal force causing astronomical orbits.

<span class="mw-page-title-main">Crystal structure</span> Ordered arrangement of atoms, ions, or molecules in a crystalline material

In crystallography, crystal structure is a description of the ordered arrangement of atoms, ions or molecules in a crystalline material. Ordered structures occur from the intrinsic nature of the constituent particles to form symmetric patterns that repeat along the principal directions of three-dimensional space in matter.

Electrical resistivity is a fundamental property of a material that measures how strongly it resists electric current. A low resistivity indicates a material that readily allows electric current. Resistivity is commonly represented by the Greek letter $ρ$ (rho). The SI unit of electrical resistivity is the ohm-meter (Ω⋅m). For example, if a 1 m³ solid cube of material has sheet contacts on two opposite faces, and the resistance between these contacts is 1 Ω, then the resistivity of the material is 1 Ω⋅m.

In mathematics and physics, n-dimensional anti-de Sitter space (AdS_n) is a maximally symmetric Lorentzian manifold with constant negative scalar curvature. Anti-de Sitter space and de Sitter space are named after Willem de Sitter (1872–1934), professor of astronomy at Leiden University and director of the Leiden Observatory. Willem de Sitter and Albert Einstein worked together closely in Leiden in the 1920s on the spacetime structure of the universe.

In quantum physics, a bound state is a quantum state of a particle subject to a potential such that the particle has a tendency to remain localized in one or more regions of space. The potential may be external or it may be the result of the presence of another particle; in the latter case, one can equivalently define a bound state as a state representing two or more particles whose interaction energy exceeds the total energy of each separate particle. One consequence is that, given a potential vanishing at infinity, negative-energy states must be bound. In general, the energy spectrum of the set of bound states is discrete, unlike free particles, which have a continuous spectrum.

In materials science, a dislocation or Taylor's dislocation is a linear crystallographic defect or irregularity within a crystal structure that contains an abrupt change in the arrangement of atoms. The movement of dislocations allow atoms to slide over each other at low stress levels and is known as glide or slip. The crystalline order is restored on either side of a glide dislocation but the atoms on one side have moved by one position. The crystalline order is not fully restored with a partial dislocation. A dislocation defines the boundary between slipped and unslipped regions of material and as a result, must either form a complete loop, intersect other dislocations or defects, or extend to the edges of the crystal. A dislocation can be characterised by the distance and direction of movement it causes to atoms which is defined by the Burgers vector. Plastic deformation of a material occurs by the creation and movement of many dislocations. The number and arrangement of dislocations influences many of the properties of materials.

In materials science, work hardening, also known as strain hardening, is the strengthening of a metal or polymer by plastic deformation. Work hardening may be desirable, undesirable, or inconsequential, depending on the context.

<span class="mw-page-title-main">Rayleigh–Taylor instability</span> Unstable behavior of two contacting fluids of different densities

The Rayleigh–Taylor instability, or RT instability, is an instability of an interface between two fluids of different densities which occurs when the lighter fluid is pushing the heavier fluid. Examples include the behavior of water suspended above oil in the gravity of Earth, mushroom clouds like those from volcanic eruptions and atmospheric nuclear explosions, supernova explosions in which expanding core gas is accelerated into denser shell gas, instabilities in plasma fusion reactors and inertial confinement fusion.

In mathematics, the stability radius of an object at a given nominal point is the radius of the largest ball, centered at the nominal point, all of whose elements satisfy pre-determined stability conditions. The picture of this intuitive notion is this:

Local-density approximations (LDA) are a class of approximations to the exchange–correlation (XC) energy functional in density functional theory (DFT) that depend solely upon the value of the electronic density at each point in space. Many approaches can yield local approximations to the XC energy. However, overwhelmingly successful local approximations are those that have been derived from the homogeneous electron gas (HEG) model. In this regard, LDA is generally synonymous with functionals based on the HEG approximation, which are then applied to realistic systems.

In materials science and engineering, the yield point is the point on a stress-strain curve that indicates the limit of elastic behavior and the beginning of plastic behavior. Below the yield point, a material will deform elastically and will return to its original shape when the applied stress is removed. Once the yield point is passed, some fraction of the deformation will be permanent and non-reversible and is known as plastic deformation.

Channelling is the process that constrains the path of a charged particle in a crystalline solid.

In computational chemistry and computational physics, the embedded atom model, embedded-atom method or EAM, is an approximation describing the energy between atoms and is a type of interatomic potential. The energy is a function of a sum of functions of the separation between an atom and its neighbors. In the original model, by Murray Daw and Mike Baskes, the latter functions represent the electron density. The EAM is related to the second moment approximation to tight binding theory, also known as the Finnis-Sinclair model. These models are particularly appropriate for metallic systems. Embedded-atom methods are widely used in molecular dynamics simulations.

The Navarro–Frenk–White (NFW) profile is a spatial mass distribution of dark matter fitted to dark matter halos identified in N-body simulations by Julio Navarro, Carlos Frenk and Simon White. The NFW profile is one of the most commonly used model profiles for dark matter halos.

A hydrogen-like atom (or hydrogenic atom) is any atom or ion with a single valence electron. These atoms are isoelectronic with hydrogen. Examples of hydrogen-like atoms include, but are not limited to, hydrogen itself, all alkali metals such as Rb and Cs, singly ionized alkaline earth metals such as Ca⁺ and Sr⁺ and other ions such as He⁺, Li²⁺, and Be³⁺ and isotopes of any of the above. A hydrogen-like atom includes a positively charged core consisting of the atomic nucleus and any core electrons as well as a single valence electron. Because helium is common in the universe, the spectroscopy of singly ionized helium is important in EUV astronomy, for example, of DO white dwarf stars.

Methods have been devised to modify the yield strength, ductility, and toughness of both crystalline and amorphous materials. These strengthening mechanisms give engineers the ability to tailor the mechanical properties of materials to suit a variety of different applications. For example, the favorable properties of steel result from interstitial incorporation of carbon into the iron lattice. Brass, a binary alloy of copper and zinc, has superior mechanical properties compared to its constituent metals due to solution strengthening. Work hardening has also been used for centuries by blacksmiths to introduce dislocations into materials, increasing their yield strengths.

The hexatic phase is a state of matter that is between the solid and the isotropic liquid phases in two dimensional systems of particles. It is characterized by two order parameters: a short-range positional and a quasi-long-range orientational (sixfold) order. More generally, a hexatic is any phase that contains sixfold orientational order, in analogy with the nematic phase.

Dislocation creep is a deformation mechanism in crystalline materials. Dislocation creep involves the movement of dislocations through the crystal lattice of the material, in contrast to diffusion creep, in which diffusion is the dominant creep mechanism. It causes plastic deformation of the individual crystals, and thus the material itself.

An electric dipole transition is the dominant effect of an interaction of an electron in an atom with the electromagnetic field.

Geometrically necessary dislocations are like-signed dislocations needed to accommodate for plastic bending in a crystalline material. They are present when a material's plastic deformation is accompanied by internal plastic strain gradients. They are in contrast to statistically stored dislocations, with statistics of equal positive and negative signs, which arise during plastic flow from multiplication processes like the Frank-Read source.

References

1 2 Bragg, Lawrance; Nye, J. F. (1947). "A Dynamical Model of a Crystal Structure" (PDF). Proc. R. Soc. Lond. A. 190 (1023): 474–481. doi: 10.1098/rspa.1947.0089 .
1 2 3 Laboratory Handout in MIT's 3.032: Mechanical Behavior of Materials
↑ Ritacco H, Kiefer F, Langevin D (June 2007). "Lifetime of bubble rafts: cooperativity and avalanches". Phys Rev Lett. 98 (24): 244501. doi:10.1103/PhysRevLett.98.244501. PMID 17677967.
↑ Ritacco HA (November 2020). "Complexity and self-organized criticality in liquid foams. A short review". Adv Colloid Interface Sci. 285: 102282. doi:10.1016/j.cis.2020.102282. PMC 7537653 . PMID 33059304.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[Bragg-1] 1 2 Bragg, Lawrance; Nye, J. F. (1947). "A Dynamical Model of a Crystal Structure" (PDF). Proc. R. Soc. Lond. A. 190 (1023): 474–481. doi: 10.1098/rspa.1947.0089 .

[mit-2] 1 2 3 Laboratory Handout in MIT's 3.032: Mechanical Behavior of Materials

[3] Ritacco H, Kiefer F, Langevin D (June 2007). "Lifetime of bubble rafts: cooperativity and avalanches". Phys Rev Lett. 98 (24): 244501. doi:10.1103/PhysRevLett.98.244501. PMID 17677967.

[4] Ritacco HA (November 2020). "Complexity and self-organized criticality in liquid foams. A short review". Adv Colloid Interface Sci. 285: 102282. doi:10.1016/j.cis.2020.102282. PMC 7537653 . PMID 33059304.

[1]

[2]

[3]

[4]