Catalog footprint

What is connected

26works
25topics
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

26 published item(s)

preprint2026arXiv

Optimized Three-Dimensional Photovoltaic Structures with LLM guided Tree Search

We present a case study for how AI coding systems can be used to generate novel scientific hypotheses. We combine a generic coding agent (Google's AntiGravity) with an LLM-driven tree search algorithm (Empirical Research Assistance / ERA) to autonomously generate high-efficiency three-dimensional photovoltaic (3DPV) structures that overcome losses limiting flat solar panels at mid-latitudes. These structures operate by presenting favorable angles to the sun throughout the day, and for illustrative purposes we focus on optimizing performance for a single solar day. Our workflow begins by using AntiGravity to reproduce calculations \cite{bernardi2012solar} showing that 3DPV can have energy densities much higher than stationary flat PV panels. We use these initial designs as the starting point for large scale tree search, where we seek improved solutions and score them for their diurnal yield. The initial tree search leads to nominally more efficient solutions, yet they are caused by algorithmic reward hacking, arising from non-physical design features such as structurally levitating disconnected tiers and exploitations of the discretizations in the optics solver. To counteract this, we develop a workflow where the coding agent iteratively patches the physics engine with constraints to eliminate reward hacking. With reward-hacking eliminated, ERA discovers a series of designs with various constraints and improved performance, including optimal designs with different fixed collector areas, optimizing zenith tracking and avoiding self shadowing. Combining coding agents with tree search (ERA) provides a powerful platform for scientific discovery, for problems whose solutions can be empirically evaluated with a score function.

preprint2026arXiv

Prospective multi-pathogen disease forecasting using autonomous LLM-guided tree search

Probabilistic forecasting of infectious diseases is crucial for public health but relies on labor-intensive manual model curation by expert modeling teams. This bespoke development bottlenecks scalability to granular geographic resolutions or emerging pathogens. Here, we present an autonomous system using Large Language Model (LLM)-guided tree search to iteratively generate, evaluate, and optimize executable forecasting software. In a fully prospective, real-time evaluation during the 2025-2026 US respiratory season, the system autonomously discovered methodologically diverse models for influenza, COVID-19, and respiratory syncytial virus (RSV). Aggregating these machine-generated models yielded an ensemble that consistently matched or outperformed the gold-standard, human-curated Centers for Disease Control and Prevention (CDC) hub ensembles out-of-sample. The system successfully navigated data-scarce "cold start" scenarios for RSV. Moreover, controlled retrospective ablations revealed that optimizing log-scale distance metrics prevents reward hacking, while an automated judge-in-the-loop ensures structural fidelity to complex scientific theories. By autonomously translating epidemiological theory into accurate, transparent code, this framework overcomes the modeling labor bottleneck, enabling rapid deployment of expert-level disease forecasting at unprecedented scales.

preprint2022arXiv

Context-Aware Abbreviation Expansion Using Large Language Models

Motivated by the need for accelerating text entry in augmentative and alternative communication (AAC) for people with severe motor impairments, we propose a paradigm in which phrases are abbreviated aggressively as primarily word-initial letters. Our approach is to expand the abbreviations into full-phrase options by leveraging conversation context with the power of pretrained large language models (LLMs). Through zero-shot, few-shot, and fine-tuning experiments on four public conversation datasets, we show that for replies to the initial turn of a dialog, an LLM with 64B parameters is able to exactly expand over 70% of phrases with abbreviation length up to 10, leading to an effective keystroke saving rate of up to about 77% on these exact expansions. Including a small amount of context in the form of a single conversation turn more than doubles abbreviation expansion accuracies compared to having no context, an effect that is more pronounced for longer phrases. Additionally, the robustness of models against typo noise can be enhanced through fine-tuning on noisy data.

preprint2022arXiv

Improved tracking of particles with highly correlated motion

Despite significant advances in particle imaging technologies over the past two decades, few advances have been made in particle tracking, i.e. linking individual particle positions across time series data. The state-of-the-art tracking algorithm is highly effective for systems in which the particles behave mostly independently. However, these algorithms become inaccurate when particle motion is highly correlated, such as in dense or strongly interacting systems. Accurate particle tracking is essential in the study of the physics of dense colloids, such as the study of dislocation formation, nucleation, and shear transformations. Here, we present a new method for particle tracking that incorporates information about the correlated motion of the particles. We demonstrate significant improvement over the state-of-the-art tracking algorithm in simulated data on highly correlated systems.

preprint2022arXiv

Optimal control of nonequilibrium systems through automatic differentiation

Controlling the evolution of nonequilibrium systems to minimize dissipated heat or work is a key goal for designing nanodevices, both in nanotechnology and biology. Progress in computing optimal protocols has thus far been limited to either simple systems or near-equilibrium evolution. Here, we present an approach for computing optimal protocols based on automatic differentiation. Our methodology is applicable to complex systems and multidimensional protocols and is valid arbitrarily far from equilibrium. We validate our method by reproducing theoretical optimal protocols for a Brownian particle in a time-varying harmonic trap. We also compute departures from near-equilibrium behaviour for magnetization reversal on an Ising lattice and for barrier crossing driven by a harmonic trap, which has been used to represent a range of biological processes including biomolecular unfolding reactions. Algorithms based on automatic differentiation outperform the near-equilibrium theory for far-from-equilibrium magnetization reversal and driven barrier crossing. The optimal protocol for crossing an energy landscape barrier of 10kT is found to hasten the approach to, and slow the departure from, the barrier region compared to the near-equilibrium theoretical protocol.

preprint2021arXiv

Machine learning accelerated computational fluid dynamics

Numerical simulation of fluids plays an essential role in modeling many physical phenomena, such as weather, climate, aerodynamics and plasma physics. Fluids are well described by the Navier-Stokes equations, but solving these equations at scale remains daunting, limited by the computational cost of resolving the smallest spatiotemporal features. This leads to unfavorable trade-offs between accuracy and tractability. Here we use end-to-end deep learning to improve approximations inside computational fluid dynamics for modeling two-dimensional turbulent flows. For both direct numerical simulation of turbulence and large eddy simulation, our results are as accurate as baseline solvers with 8-10x finer resolution in each spatial dimension, resulting in 40-80x fold computational speedups. Our method remains stable during long simulations, and generalizes to forcing functions and Reynolds numbers outside of the flows where it is trained, in contrast to black box machine learning approaches. Our approach exemplifies how scientific computing can leverage machine learning and hardware accelerators to improve simulations without sacrificing accuracy or generalization.

preprint2021arXiv

Topological Design of Heterogeneous Self-Assembly

Controlling the topology of structures self-assembled from a set of heterogeneous building blocks is highly desirable for many applications, but is poorly understood theoretically. Here we show that the thermodynamic theory of self-assembly involves an inevitable divergence in chemical potential. The divergence and its detailed structure are controlled by the spectrum of the transfer matrix, which summarizes all of self-assembly design degrees of freedom. By analyzing the transfer matrix, we map out the phase boundary between the designable structures and the unstructured aggregates, driven by the level of cross-talk.

preprint2020arXiv

Predicting ground state configuration of energy landscape ensemble using graph neural network

Many scientific problems seek to find the ground state in a rugged energy landscape, a task that becomes prohibitively difficult for large systems. Within a particular class of problems, however, the short-range correlations within energy minima might be independent of system size. Can these correlations be inferred from small problems with known ground states to accelerate the search for the ground states of larger problems? Here, we demonstrate the strategy on Ising spin glasses, where the interaction matrices are drawn from protein contact maps. We use graph neural network to learn the mapping from an interaction matrix $J$ to a ground state configuration, yielding guesses for the set of most probable configurations. Given these guesses, we show that ground state configurations can be searched much faster than in vanilla simulated annealing. For large problems, a model trained on small $J$ matrices predicts a configurations whose energy is much lower than those obtained by simulated annealing, indicating the size generalizability of the strategy.

preprint2019arXiv

Using Attribution to Decode Dataset Bias in Neural Network Models for Chemistry

Deep neural networks have achieved state of the art accuracy at classifying molecules with respect to whether they bind to specific protein targets. A key breakthrough would occur if these models could reveal the fragment pharmacophores that are causally involved in binding. Extracting chemical details of binding from the networks could potentially lead to scientific discoveries about the mechanisms of drug actions. But doing so requires shining light into the black box that is the trained neural network model, a task that has proved difficult across many domains. Here we show how the binding mechanism learned by deep neural network models can be interrogated, using a recently described attribution method. We first work with carefully constructed synthetic datasets, in which the 'fragment logic' of binding is fully known. We find that networks that achieve perfect accuracy on held out test datasets still learn spurious correlations due to biases in the datasets, and we are able to exploit this non-robustness to construct adversarial examples that fool the model. The dataset bias makes these models unreliable for accurately revealing information about the mechanisms of protein-ligand binding. In light of our findings, we prescribe a test that checks for dataset bias given a hypothesis. If the test fails, it indicates that either the model must be simplified or regularized and/or that the training dataset requires augmentation.

preprint2016arXiv

Pruning to Increase Taylor Dispersion in Physarum polycephalum Networks

How do the topology and geometry of a tubular network affect the spread of particles within fluid flows? We investigate patterns of effective dispersion in the hierarchical, biological transport network formed by Physarum polycephalum. We demonstrate that a change in topology - pruning in the foraging state - causes a large increase in effective dispersion throughout the network. By comparison, changes in the hierarchy of tube radii result in smaller and more localized differences. Pruned networks capitalize on Taylor dispersion to increase the dispersion capability.

preprint2016arXiv

Receptor arrays optimized for natural odor statistics

Natural odors typically consist of many molecules at different concentrations. It is unclear how the numerous odorant molecules and their possible mixtures are discriminated by relatively few olfactory receptors. Using an information-theoretic model, we show that a receptor array is optimal for this task if it achieves two possibly conflicting goals: (i) each receptor should respond to half of all odors and (ii) the response of different receptors should be uncorrelated when averaged over odors presented with natural statistics. We use these design principles to predict statistics of the affinities between receptors and odorant molecules for a broad class of odor statistics. We also show that optimal receptor arrays can be tuned to either resolve concentrations well or distinguish mixtures reliably. Finally, we use our results to predict properties of experimentally measured receptor arrays. Our work can thus be used to better understand natural olfaction and it also suggests ways to improve artificial sensor arrays.

preprint2016arXiv

The Information Capacity of Specific Interactions

Specific interactions are a hallmark feature of self-assembly and signal-processing systems in both synthetic and biological settings. Specificity between components may arise from a wide variety of physical and chemical mechanisms in diverse contexts, from DNA hybridization to shape-sensitive depletion interactions. Despite this diversity, all systems that rely on interaction specificity operate under the constraint that increasing the number of distinct components inevitably increases off-target binding. Here we introduce `capacity', the maximal information encodable using specific interactions, to compare specificity across diverse experimental systems, and to compute how specificity changes with physical parameters. Using this framework, we find that `shape'-coding of interactions has higher capacity than chemical (`color') coding because the strength of off-target binding is strongly sublinear in binding site size for shapes while being linear for colors. We also find that different specificity mechanisms, such as shape and color, can be combined in a synergistic manner, giving a capacity greater than the sum of the parts.

preprint2015arXiv

Two-Dimensional Clusters of Colloidal Spheres: Ground States, Excited States, and Structural Rearrangements

We study experimentally what is arguably the simplest yet non-trivial colloidal system: two-dimensional clusters of 6 spherical particles bound by depletion interactions. These clusters have multiple, degenerate ground states whose equilibrium distribution is determined by entropic factors, principally the symmetry. We observe the equilibrium rearrangements between ground states as well as all of the low-lying excited states. In contrast to the ground states, the excited states have soft modes and low symmetry, and their occupation probabilities depend on the size of the configuration space reached through internal degrees of freedom, as well as a single "sticky parameter" encapsulating the depth and curvature of the potential. Using a geometrical model that accounts for the entropy of the soft modes and the diffusion rates along them, we accurately reproduce the measured rearrangement rates. The success of this model, which requires no fitting parameters or measurements of the potential, shows that the free-energy landscape of colloidal systems and the dynamics it governs can be understood geometrically.

preprint2014arXiv

Multifarious Assembly Mixtures: Systems Allowing Retrieval of Diverse Stored Structures

Self-assembly materials are traditionally designed so that molecular or meso-scale components form a single kind of large structure. Here, we propose a scheme to create "multifarious assembly mixtures", which self-assemble many different large structures from a set of shared components. We show that the number of multifarious structures stored in the solution of components increases rapidly with the number of different types of components. Yet, each stored structure can be retrieved by tuning only a few parameters, the number of which is only weakly dependent on the size of the assembled structure. Implications for artificial and biological self-assembly are discussed.

preprint2012arXiv

A geometrical approach to computing free energy landscapes from short-ranged potentials

Particles interacting with short-ranged potentials have attracted increasing interest, partly for their ability to model mesoscale systems such as colloids interacting via DNA or depletion. We consider the free energy landscape of such systems as the range of the potential goes to zero. In this limit, the landscape is entirely defined by geometrical manifolds, plus a single control parameter. These manifolds are fundamental objects that do not depend on the details of the interaction potential, and provide the starting point from which any quantity characterizing the system -- equilibrium or non-equilibrium -- can be computed for arbitrary potentials. To consider dynamical quantities we compute the asymptotic limit of the Fokker-Planck equation, and show that it becomes restricted to the low-dimensional manifolds connected by "sticky" boundary conditions. To illustrate our theory, we compute the low-dimensional manifolds for n<=8 identical particles, providing a complete description of the lowest-energy parts of the landscape including floppy modes with up to 2 internal degrees of freedom. The results can be directly tested on colloidal clusters. This limit is a novel approach for understanding energy landscapes, and our hope is that it can also provide insight into finite-range potentials.

preprint2012arXiv

Development of Knife-Edge Ridges on Ion-Bombarded Surfaces

We demonstrate in both laboratory and numerical experiments that ion bombardment of a modestly sloped surface can create knife-edge like ridges with extremely high slopes. Small pre-fabricated pits expand under ion bombardment, and the collision of two such pits creates knife-edge ridges. Both laboratory and numerical experiments show that the pit propagation speed and the precise shape of the knife edge ridges are universal, independent of initial conditions, as has been predicted theoretically. These observations suggest a novel method of fabrication in which a surface is pre-patterned so that it dynamically evolves to a desired target pattern made of knife-edge ridges.

preprint2011arXiv

Algorithm for a microfluidic assembly line

Microfluidic technology has revolutionized the control of flows at small scales giving rise to new possibilities for assembling complex structures on the microscale. We analyze different possible algorithms for assembling arbitrary structures, and demonstrate that a sequential assembly algorithm can manufacture arbitrary 3D structures from identical constituents. We illustrate the algorithm by showing that a modified Hele-Shaw cell with 7 controlled flowrates can be designed to construct the entire English alphabet from particles that irreversibly stick to each other.

preprint2011arXiv

Deriving Finite Sphere Packings

Sphere packing problems have a rich history in both mathematics and physics; yet, relatively few analytical analyses of sphere packings exist, and answers to seemingly simple questions are unknown. Here, we present an analytical method for deriving all packings of n spheres in R^3 satisfying minimal rigidity constraints (>= 3 contacts per sphere and >= 3n-6 total contacts). We derive such packings for n <= 10, and provide a preliminary set of maximal contact packings for 10 < n <= 20. The resultant set of packings has some striking features, among them: (i) all minimally rigid packings for n <= 9 have 3n-6 contacts, (ii) non-rigid packings satisfying minimal rigidity constraints arise for n >= 9, (iii) the number of ground states (i.e. packings with the maximal number of contacts) oscillates with respect to n, (iv) for 10 <= n <= 20 there are only a small number of packings with the maximal number of contacts, and for 10 <= n < 13 these are all commensurate with the HCP lattice. The general method presented here may have applications to other related problems in mathematics, such as the Erdos repeated distance problem and Euclidean distance matrix completion problems.

preprint2011arXiv

Skating on a Film of Air: Drops Impacting on a Surface

Drops impacting on a surface are ubiquitous in our everyday experience. This impact is understood within a commonly accepted hydrodynamic picture: it is initiated by a rapid shock and a subsequent ejection of a sheet leading to beautiful splashing patterns. However, this picture ignores the essential role of the air that is trapped between the impacting drop and the surface. Here we describe a new imaging modality that is sensitive to the behavior right at the surface. We show that a very thin film of air, only a few tens of nanometers thick, remains trapped between the falling drop and the surface as the drop spreads. The thin film of air serves to lubricate the drop enabling the fluid to skate on the air film laterally outward at surprisingly high velocities, consistent with theoretical predictions. Eventually this thin film of air must break down as the fluid wets the surface. We suggest that this occurs in a spinodal-like fashion, and causes a very rapid spreading of a wetting front outwards; simultaneously the wetting fluid spreads inward much more slowly, trapping a bubble of air within the drop. Our results show that the dynamics of impacting drops are much more complex than previously thought and exhibit a rich array of unexpected phenomena that require rethinking classical paradigms.

preprint2010arXiv

Breakup of diminutive Rayleigh jets

Discharging a liquid from a nozzle at sufficient large velocity leads to a continuous jet that due to capillary forces breaks up into droplets. Here we investigate the formation of microdroplets from the breakup of micron-sized jets with ultra high-speed imaging. The diminutive size of the jet implies a fast breakup time scale $τ_\mathrm{c} = \sqrt{ρr^3 / γ}$ of the order of 100\,ns{}, and requires imaging at 14 million frames per second. We directly compare these experiments with a numerical lubrication approximation model that incorporates inertia, surface tension, and viscosity [Eggers and Dupont, J. Fluid Mech. 262, 205 (1994); Shi, Brenner, and Nagel, Science 265, 219 (1994)]. The lubrication model allows to efficiently explore the parameter space to investigate the effect of jet velocity and liquid viscosity on the formation of satellite droplets. In the phase diagram we identify regions where the formation of satellite droplets is suppressed. We compare the shape of the droplet at pinch-off between the lubrication approximation model and a boundary integral (BI) calculation, showing deviations at the final moment of the pinch-off. Inspite of this discrepancy, the results on pinch-off times and droplet and satellite droplet velocity obtained from the lubrication approximation agree with the high-speed imaging results.

preprint2010arXiv

MD-Predicted Phase diagrams for Pattern Formation due to Ion Irradiation

Energetic particle irradiation of solids can cause surface ultra-smoothening, self-organized nanoscale pattern formation, or degradation of the structural integrity of nuclear reactor components. Periodic patterns including high-aspect ratio quantum dots, with occasional long-range order and characteristic spacing as small as 7 nm, have stimulated interest in this method as a means of sub-lithographic nanofabrication. Despite intensive research there is little fundamental understanding of the mechanisms governing the selection of smooth or patterned surfaces, and precisely which physical effects cause observed transitions between different regimes has remained a matter of speculation. Here we report the first prediction of the mechanism governing the transition from corrugated surfaces to flatness, using only parameter-free molecular dynamics simulations of single-ion impact induced crater formation as input into a multi-scale analysis, and showing good agreement with experiment. Our results overturn the paradigm attributing these phenomena to the removal of target atoms via sputter erosion. Instead, the mechanism dominating both stability and instability is shown to be the impact-induced redistribution of target atoms that are not sputtered away, with erosive effects being essentially irrelevant. The predictions are relevant in the context of tungsten plasma-facing fusion reactor walls which, despite a sputter erosion rate that is essentially zero, develop, under some conditions, a mysterious nanoscale topography leading to surface degradation. Our results suggest that degradation processes originating in impact-induced target atom redistribution effects may be important, and hence that an extremely low sputter erosion rate is an insufficient design criterion for morphologically stable solid surfaces under energetic particle irradiation.

preprint2004arXiv

Shear-dependent apparent slip on hydrophobic surfaces: The Mattress Model

Recent experiments (Zhu & Granick (2001) Phys. Rev. Lett. 87 096105) have measured a large shear dependent fluid slip at partially wetting fluid-solid surfaces. We present a simple model for such slip, motivated by the recent observations of nanobubbles on hydrophobic surfaces. The model considers the dynamic response of bubbles to change in hydrodynamic pressure due to the oscillation of a solid surface. Both the compression and diffusion of gas in the bubbles decrease the force on the oscillating surface by a ``leaking mattress'' effect, thereby creating an apparent shear-dependent slip. With bubbles similar to those observed by atomic force microscopy to date, the model is found to lead to force decreases consistent with the experimental measurements of Zhu & Granick.

preprint1999arXiv

Collapsing Bacterial Cylinders

Under special conditions bacteria excrete an attractant and aggregate. The high density regions initially collapse into cylindrical structures, which subsequently destabilize and break up into spherical aggregates. This paper presents a theoretical description of the process, from the structure of the collapsing cylinder to the spacing of the final aggregates. We show that cylindrical collapse involves a delicate balance in which bacterial attraction and diffusion nearly cancel, leading to corrections to the collapse laws expected from dimensional analysis. The instability of a collapsing cylinder is composed of two distinct stages: Initially, slow modulations to the cylinder develop, which correspond to a variation of the collapse time along the cylinder axis. Ultimately, one point on the cylinder pinches off. At this final stage of the instability, a front propagates from the pinch into the remainder of the cylinder. The spacing of the resulting spherical aggregates is determined by the front propagation.

preprint1998arXiv

Electrostatic Edge Instability of Lipid Membranes

This paper considers the effect of electrostatics on the stability of a charged membrane. We show that at low ionic strength and high surface charge density, repulsion between charges on the membrane renders it unstable to the formation of holes. A free straight edge is unstable to modulations with wavelength longer than the Debye screening length. Hence at low ionic strength, membranes will disintegrate into vesicles. We use these results to interpret stable holes in red blood cell ghosts (Steck et. al., Science 168:255(1970)).

preprint1996arXiv

Acoustic Energy Storage in Single Bubble Sonoluminescence

Single bubble sonoluminescence is understood in terms of a shock focusing towards the bubble center. We present a mechanism for significantly enhancing the effect of shock focusing, arising from the storage of energy in the acoustic modes of the gas. The modes with strongest coupling are not spherically symmetric. The storage of acoustic energy gives a framework for understanding how light intensities depend so strongly on ambient gases and liquids and suggests that the light intensities of successive flashes are highly correlated.

preprint1995arXiv

Bubble Shape Oscillations and the Onset of Sonoluminescence

An air bubble trapped in water by an oscillating acoustic field undergoes either radial or nonspherical pulsations depending on the strength of the forcing pressure. Two different instability mechanisms (the Rayleigh--Taylor instability and parametric instability) cause deviations from sphericity. Distinguishing these mechanisms allows explanation of many features of recent experiments on sonoluminescence, and suggests methods for finding sonoluminescence in different parameter regimes.