Geometry and dynamics link form, function, and evolution of finch beaks

Darwin’s finches are a classic example of adaptive radiation, exemplified by their adaptive and functional beak morphologies. To quantify their form, we carry out a morphometric analysis of the three-dimensional beak shapes of all of Darwin’s finches and find that they can be fit by a transverse parabolic shape with a curvature that increases linearly from the base toward the tip of the beak. The morphological variation of beak orientation, aspect ratios, and curvatures allows us to quantify beak function in terms of the elementary theory of machines, consistent with the dietary variations across finches. Finally, to explain the origin of the evolutionary morphometry and the developmental morphogenesis of the finch beak, we propose an experimentally motivated growth law at the cellular level that simplifies to a variant of curvature-driven flow at the tissue level and captures the range of observed beak shapes in terms of a simple morphospace. Altogether, our study illuminates how a minimal combination of geometry and dynamics allows for functional form to develop and evolve.

Nonadiabatic Dynamics Algorithms with Only Potential Energies and Gradients: Curvature-Driven Coherent Switching with Decay of Mixing and Curvature-Driven Trajectory Surface Hopping

Direct dynamics by mixed quantum–classical nonadiabatic methods is an important tool for understanding processes involving multiple electronic states. Very often, the computational bottleneck of such direct simulation comes from electronic structure theory. For example, at every time step of a trajectory, nonadiabatic dynamics requires potential energy surfaces, their gradients, and the matrix elements coupling the surfaces. The need for the couplings can be alleviated by employing the time derivatives of the wave functions, which can be evaluated from overlaps of electronic wave functions at successive timesteps. However, evaluation of overlap integrals is still expensive for large systems. In addition, for electronic structure methods for which the wave functions or the coupling matrix elements are not available, nonadiabatic dynamics algorithms become inapplicable. In this work, building on recent work by Baeck and An, we propose new nonadiabatic dynamics algorithms that only require adiabatic potential energies and their gradients. The new methods are named curvature- driven coherent switching with decay of mixing (κCSDM) and curvature-driven trajectory surface hopping (κTSH). We show how powerful these new methods are in terms of computer time and good agreement with methods employing nonadiabatic coupling vectors computed in conventional ways. The lowering of the computational cost will allow longer nonadiabatic trajectories and greater ensemble averaging to be affordable, and the ability to calculate the dynamics without electronic structure coupling matrix elements extends the dynamics capability to new classes of electronic structure methods.

Curvature flows, scaling laws and the geometry of attrition under impacts

Impact induced attrition processes are, beyond being essential models of industrial ore processing, broadly regarded as the key to decipher the provenance of sedimentary particles. Here we establish the first link between microscopic, particle-based models and the mean field theory for these processes. Based on realistic computer simulations of particle-wall collision sequences we first identify the well-known damage and fragmentation energy phases, then we show that the former is split into the abrasion phase with infinite sample lifetime (analogous to Sternberg’s Law) at finite asymptotic mass and the cleavage phase with finite sample lifetime, decreasing as a power law of the impact velocity (analogous to Basquin’s Law). This splitting establishes the link between mean field models (curvature-driven partial differential equations) and particle-based models: only in the abrasion phase does shape evolution emerging in the latter reproduce with startling accuracy the spatio-temporal patterns (two geometric phases) predicted by the former.

Curvature effects and radial homoclinic snaking

In this work, we revisit some general results on the dynamics of circular fronts between homogeneous states and the formation of localized structures in two dimensions (2D). We show how the bifurcation diagram of axisymmetric structures localized in radius fits within the framework of collapsed homoclinic snaking. In 2D, owing to curvature effects, the collapse of the snaking structure follows a different scaling that is determined by the so-called nucleation radius. Moreover, in the case of fronts between two symmetry-related states, the precise point in parameter space to which radial snaking collapses is not a ‘Maxwell’ point but is determined by the curvature-driven dynamics only. In this case, the snaking collapses to a ‘zero surface tension’ point. Near this point, the breaking of symmetry between the homogeneous states tilts the snaking diagram. A different scaling law is found for the collapse of the snaking curve in each case. Curvature effects on axisymmetric localized states with internal structure are also discussed, as are cellular structures separated from a homogeneous state by a circular front. While some of these results are well understood in terms of curvature-driven dynamics and front interactions, a proper mathematical description in terms of homoclinic trajectories in a radial spatial dynamics description is lacking.