Search NASA⌕ Search

SEARCH · Search NASA

Results for “Physics - Plasma physics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

528 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Osmotic and phoretic competition explains chemotaxic assembly and sorting

Microscale objects responding to chemical gradients by migrating toward or away from a preferred species is a simple yet constitutive mechanism by which transport occurs in biological organisms. Synthetic chemotaxis provides key physical descriptions of simplified systems that can be used in biological models, or in the creation of advanced responsive material systems. In this article, we provide a quantitative framework for understanding synthetic chemotaxis of microparticles which involves a competition between phoresis and osmosis. We present separate quantitative measurements of phoresis and osmosis acting on individual taxing particles, finding that phoresis follows the long-predicted v ∼ 1 / r 2 scaling while the osmotic contribution depends on the geometry and details of the system, and must be solved on a case-by-case basis. Through this, we are able to develop a more accurate picture of particle transport at the single particle level. Equipped with this approach, we go on to describe how high concentrations of particles in a symmetric chemical gradient grow close-packed hives that reach a steady-state size tunable through light intensity or particle size. Last, we demonstrate that mixed particles experiencing the same chemical gradient will selectively migrate toward or away depending on the nature of the particle surface, thereby locally sorting out a particular species. We anticipate these results will be important in describing both biological and synthetic chemotaxis in phoretic systems and should bring a wealth of studies that take advantage of competing osmotic flows to illicit unexpected dynamic active behavior.

Science & Technology - Other Topics↗
DOE OSTI2026

Impacts of Multidimensional Progenitor Perturbations on Core-collapse Supernova Explosions

Numerical studies of core-collapse supernovae have demonstrated the importance of nonradial motions in precollapse progenitors on the explosion outcome. We use the Chimera neutrino radiation hydrodynamics code running seven two-dimensional simulations of 15 M⊙ progenitors with different progenitor structures introduced by different one- and two-dimensional precollapse stellar evolution environments to examine the impacts of stellar structure and nonspherical motion in the precollapse progenitor on the development of explosions. We compare the explosion evolution of these models in terms of shock dynamics, diagnostic energy, neutrino heating, accretion, explosion geometry, nuclear abundances, and turbulent convection. We also analyze how stochastic variation impacts our simulations. Contrary to results reported in prior studies examining the impacts of multidimensional progenitors, we observe similar shock revival times and explosion development in our simulations despite differences in initial compositions and structures. We find no discernible impact from the accretion of nonradial perturbations from a multi-D progenitor onto the stalled shock in the revival and strength of explosion, as fully developed neutrino-driven convection behind the stalled shock is similar for all our models. For models with physically sourced noise in the iron core, a strong oscillation of the shock occurs after bounce and deflects infall laterally, and accelerates the saturation of the lateral turbulent kinetic energy. An examination of model stochasticity shows that any prior expected impacts on explosive outcome due to convection-related perturbations lie below the detectable threshold of numerical variation.

Chen, Chien-Hui [North Carolina State University]↗
DOE OSTI2025

Semiconducting Electrides Derived from Sodalite: A First-Principles Study

Electrides are ionic crystals, with electrons acting as anions occupying well-defined lattice sites. These exotic materials have attracted considerable attention in recent years for potential applications in catalysis, rechargeable batteries, and display technology. Among this class of materials, electride semiconductors can further expand the horizon of potential applications due to the presence of a band gap. However, there are only limited reports on semiconducting electrides, hindering the understanding of their physical and chemical properties. In recent work, we initiated an approach to derive potential electrides via selective removal of symmetric Wyckoff sites of anions from existing complex minerals. Herein, we present a follow-up effort to design semiconducting electrides from parental complex sodalites. Among four candidate compounds, we found that a cubic Ca 4 Al 6 O 12 structure with the I-43m space group symmetry exhibits perfect electron localization at the sodalite cages, with a narrow electronic band gap of 1.8 eV, making it suitable for use in photocatalysis. Analysis of the electronic structures reveals that a lower electronegativity of the surrounding cations drives greater electron localization and promotes the formation of an electride band near the Fermi level. Our work proposes an alternative approach for designing new semiconducting electrides under ambient conditions and offers guidelines for further experimental exploration.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Monolayer Films of Binary Metal Oxide Nanoparticles for Carbon Nanotube Synthesis

Binary metal oxide nanoparticles (biMO-NPs) combining transition and main-group metals with well-organized atomic architectures have unique potential as robust and economical heterogeneous catalysts. Herein, metal oxide nanoparticles (CoAlxOy) using cobalt and aluminum were synthesized, and their structure, morphology, and chemical composition were investigated using various characterization techniques. The biMO-NPs, which had a Gaussian-type size distribution (∼6 nm), were assembled on plain silicon substrates modified chemically with linker molecules bearing suitable end groups. Surface analysis revealed an ordered, uniform, close-packed arrangement of biMO-NPs forming a monolayer architecture with nanoscale thickness (∼12 nm), high surface uniformity, and negligible defects. The nanoparticle monolayer film (NP-MF) was employed as a catalyst to grow carbon nanotube (CNT) forests via the thermal chemical vapor deposition (CVD) method, similar to those grown from catalyst thin films deposited by physical vapor deposition. Structural characterization confirmed the growth of a dense, uniform, high-quality vertically aligned carbon nanotube (VA-CNT) forest with a length of ∼125 µm, which is comparable to the VA-CNTs grown from expensive thin films. Thus, CNT growth was successfully catalyzed by the biMO-NPs, highlighting a simple and inexpensive alternative for catalyst deposition. In this innovative wet-chemistry approach, the catalyst and catalyst support are combined within a single nanoparticle before NP-MF assembly. Interestingly, this approach provides a scalable and cost-effective pathway for CNT growth, eliminating the need for expensive thin deposition techniques.

Regmi, Bishow [University of Cincinnati]↗
DOE OSTI2026

Insights into elevated temperature tensile deformation mechanisms and kink banding in additively manufactured tungsten

Additive manufacturing (AM) potentially enables fabrication and repair of plasma facing components (PFCs) of tungsten. However, deformation mechanisms in AM-fabricated tungsten under service-relevant thermomechanical conditions remain unknown majorly due to challenges in achieving crack-free W. This study reveals the deformation mechanisms operative under tension at 800 °C and 1200 °C in the electron beam melting powder bed fusion (EBM-PBF) fabricated W. Textured columnar grains with mixed <001 >/<111 >|| build-direction were observed. Strong anisotropy in tensile properties persisted across test temperatures, and all specimens exhibited extensive deformation-induced banding phenomena. Grain boundary character analysis of the deformed specimens indicated that these deformation bands were kink bands with tilt grain boundaries prominently present at the band-matrix interface. Interestingly, the material within the kink bands rotated toward higher resolved shear stress, providing insights into the origins of kink band formation in body-centered cubic refractory metals. Intragranular misorientation axis analysis revealed that plasticity was dominated by {110} <111 > slip at 800 °C, whereas additional (213)[11⁢̄1] and (112)[11⁢̄1] slip systems activated at 1200 °C. A dense low-angle boundary networks observed near fractured surface indicated the onset of dynamic recrystallization. Results reveal plasticity governing mechanisms in AM-fabricated W under power plant-relevant conditions, and provide insights into kink band attributes in refractory metals.

Mayes, Riley [ORNL] (ORCID:000900089307010X)↗
DOE OSTI2026

Improving the Performance of NEML2 with Modern Graph Compilation Backends

NEML2 vectorizes constitutive-model evaluation for large-scale multiphysics simulation, using PyTorch as its tensor backend so that a batch of material-point updates runs on CPU or GPU through a single implementation. In the two prior reports in this series it was a C++-native library, deployed through TorchScript tracing and just-in-time (JIT) compilation; it has since been rewritten from the ground up into a Python-native library deployed through Ahead-of-Time Inductor (AOTInductor), a modern PyTorch graph-compilation backend. The rewrite is driven by a persistent tension, not a language preference: NEML2 composes constitutive models at runtime from a registry of small, independently-authored pieces, and that flexibility is difficult to reconcile with the compile-time knowledge an efficient GPU kernel needs. This report documents the rewrite and the investment that accompanied it: the AOTInductor export pipeline that turns a Python-authored model into a portable, Python-free compiled artifact loadable from pure C++; the eager and compiled runtimes and the new implicit solver layer built on them; a head-to-head benchmark of legacy JIT against AOTInductor; the physics-model catalog and its worked examples; the developer tooling; and the corresponding overhaul of MOOSE’s NEML2 integration that lets MOOSE consume it. A central objective is to examine whether modern PyTorch graph-compilation backends are effective for MOOSE GPU integration. The benchmark answers directly: AOTInductor outperforms legacy JIT on every GPU scenario measured, by 1.0–4.5×. Modern graph-compilation backends are effective for MOOSE GPU integration, and AOTInductor specifically – not compilation in the abstract – is why.

Hu, Gary (Tianchen) [Argonne National Laboratory (↗
DOE OSTI2026

Optimal binning of correlated measurements

Experimental measurements are commonly represented on a discrete grid, requiring a balance between granularity and statistical noise. Two strategies have traditionally been used to improve such representations: selecting an appropriate bin width to control discretization error and applying kernel-based smoothing to suppress fluctuations. Despite their shared goal, these approaches have largely developed independently, without a unified statistical description of how discretization and correlation jointly determine measurement precision. Here, we extend the discussion of optimal interval averaging to a correlation-aware setting by Gaussian process regression, which explicitly accounts for correlations among neighboring bins. Starting from first principles, we derive the mean-squared error of discretized measurements and obtain closed-form asymptotic expressions for the optimal bin width and correlation length. When recast in reduced variables, the theory reveals distinct universal scaling laws governing the error in the correlation-free and correlation-controlled regimes. Characterized by intrinsically smooth intensity profiles and counting-based statistics, neutron scattering measurements are well suited for demonstrating the enhanced error contraction enabled by inter-bin correlations. We show that such improvement is achievable over the experimentally accessible Q-range and across multiple instruments and material systems. These results show that explicitly accounting for correlations systematically reshapes the limits of precision in discretized, noise-limited measurements. More broadly, the framework provides a transferable statistical foundation for optimizing data representation, inference, and experimental design across the physical and data sciences.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗
DOE OSTI2026

Automated scanning probe microscopy of combinatorial ferroelectric libraries: Gaussian-process-guided exploration and noise-aware experiment planning

Combinatorial materials libraries provide an efficient route for mapping composition–property relationships, but their broader impact depends on rapid, quantitative, and functionally relevant characterization. Scanning Probe Microscopy (SPM), including piezoresponse force microscopy (PFM), offers significant potential for quantitative, functionally relevant combi-library readouts. Here, we implement a fully automated SPM workflow for ferroelectric combinatorial libraries and benchmark Gaussian-process-based Bayesian optimization strategies for autonomous experiment planning. The workflow integrates automated probe motion, contact optimization, imaging, and dual amplitude resonance tracking-PFM spectroscopy, and uses scalarized spectroscopic observables to guide subsequent measurements. Stage motion, probe engagement, in-contact tuning, imaging, spectroscopy, and the choice of the next measurement location all proceed without human input. We demonstrate the approach on Sm-doped BiFeO 3 and Zn x Mg 1−x O libraries. By comparing vanilla Bayesian optimization with a measured-noise variant, we show that explicit treatment of local reproducibility can improve modeling of composition-dependent response when the measured variance is physically meaningful, but can also reduce robustness when variability is dominated by outliers or topographic artifacts. Furthermore, these results establish automated SPM as a bridge between combinatorial synthesis and quantitative functional characterization.

Liu, Yu [University of Tennessee, Knoxville, TN (U↗
DOE OSTI2026

Spectroscopic Demarcation of Emergent Photons and Spinons in a Dipolar-Octupolar Quantum Spin Liquid

The identification of fractionalized excitations in quantum spin liquids (QSLs) remains a central challenge in condensed matter physics. In dipolar-octupolar (DO) pyrochlores, such as Ce 2⁢ Zr 2 ⁢O 7 , the candidate 𝜋-flux quantum spin ice (QSI) state is predicted to host both gapless emergent photons and a continuum of spinons. However, resolving these modes at zero field is complicated by their spectral overlap and the presence of nonmagnetic scattering near zero energy. Here, in this study, we report neutron scattering experiments on Ce 2 ⁢Zr 2 ⁢O 7 under a magnetic field along the [1, 1, 1] direction. In contrast to previous unpolarized studies at zero field that relied on high-temperature subtraction, we use a same-temperature high-field subtraction protocol to isolate the photon and spinon modes. Leveraging the selective coupling of the magnetic field to the dipolar degrees of freedom, we demonstrate the spectroscopic demarcation of these excitations. We observe that weak fields (≈ 0.15 T) suppress the low-energy photon weight while leaving the high-energy spinon continuum robust, albeit hardened. Our results, supported by gauge mean-field theory and exact diagonalization calculations, provide strong evidence for the 𝜋-flux QSI state and introduce a powerful field-tuning protocol for investigating DO-QSLs.

Gao, Bin [Rice University, Houston, TX (United Sta↗
DOE OSTI2026

Material-dependent photon ionizing radiation effects in Si and GaAs PIN diodes: A numerical investigation

We present a finite-element drift-diffusion-Poisson model in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework to compare the radiation response of silicon (Si) and gallium arsenide (GaAs) PIN diodes under high-energy photon irradiation. The model solves coupled carrier continuity and Poisson’s equations with Shockley-Read-Hall recombination, and is verified against standard analytical J-V behavior. Using a simplified 1D geometry with ideal Ohmic contacts, we quantify device response under forward and reverse bias with a 100 MeV photon flux. Under forward bias, Si exhibits markedly greater radiation sensitivity than GaAs, including larger increases in current density, stronger local field and carrier-product perturbations, and higher recombination. Under reverse bias, GaAs shows larger radiation-induced photocurrent and broader current-density peaks near junctions, indicating an advantage for photodetection. Integrated steady-state recombination is consistently higher in Si across voltages. Under periodic photon pulses, GaAs produces higher-amplitude photoresponse and settles more rapidly than Si. These results highlight material-dependent trade-offs for radiation-tolerant, high-speed optoelectronics and provide guidance for selecting PIN architectures in aerospace, nuclear, and high-energy physics environments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dressed-State Hamiltonian Engineering in a Strongly Interacting Solid-State Spin Ensemble

In quantum science applications, ranging from many-body physics to quantum metrology, dipolar interactions in spin ensembles are often controlled via Floquet engineering. However, this technique typically reduces the interaction strength between spins and effectively weakens the coupling to a target sensing field, limiting the metrological sensitivity. In this Letter, we develop and demonstrate an alternative method that directly tunes the native dipolar interaction in an ensemble of nitrogen-vacancy (NV) centers in diamond, thereby overcoming these limitations inherent to Floquet engineering. Our approach utilizes dressed-state qubit encoding under a bias magnetic field applied perpendicular to the crystal lattice orientation. This method leads to a 3.2× enhancement of the dimensionless coherence parameter JT 2 compared to state-of-the-art Floquet engineering and a 2.6× (8.3 dB) enhanced sensitivity in ac magnetometry. Furthermore, our results provide a powerful Hamiltonian engineering tool for future studies with NV ensembles and other interacting higher-spin (S > $\frac{1}{2}$) systems.

Quantum control↗
DOE OSTI2026

AutoBEM: A scalable framework for nationwide building energy simulation and retrofit evaluation in the United States

This paper presents AutoBEM, an integrated, automated framework for nationwide building energy modeling and retrofit evaluation in the United States. Unlike prior UBEM platforms that either rely primarily on representative stock sampling or operate at city scale, AutoBEM automates the generation of building-resolved, physics-based EnergyPlus/OpenStudio simulation models at national scale using GIS-derived geometry, prototype-based assumptions, and standardized scalable workflows. Leveraging the Model America dataset and high-performance computing, AutoBEM generates and simulates energy models for 122.9 million buildings, representing 97.8% of the U.S. building stock. These models are being made publicly and freely available as the Model America v1.0 (MAv1) dataset. AutoBEM supports detailed, building-level assessments of energy consumption, CO2 emissions, and post-processed anthropogenic heat emissions (AHE), and evaluates 151 energy conservation measures (ECMs) using localized utility pricing and building characteristics. In addition, AutoBEM incorporates both typical and future climate conditions through integration with Typical Meteorological Year (TMY) and Future TMY (fTMY) weather data derived from IPCC scenarios. In a case study of Phoenix, Arizona, AutoBEM identified several high-efficiency HVAC upgrades and selected envelope measures with short modeled payback periods (1.5 years) for certain building types and standards. Simulations under future climate scenarios (SSP5–RCP8.5) project an 11.3% increase in electricity use and a 32% reduction in natural gas demand by 2100, underscoring the need for climate-adaptive retrofit planning. By enabling reproducible, bottom-up, and location-specific analysis at scale, AutoBEM provides a step toward a national digital twin of the built environment and supports data-driven screening and planning for decarbonization, resilience, and energy equity.

Li, Hang [ORNL] (ORCID:0000000306001920)↗
DOE OSTI2026

In situ studies on heavy ion irradiation and partial oxidation induced stacking faults and nanograins in tungsten nanolaminates

Understanding the microstructural evolution of tungsten (W) under extreme irradiation environments is critical for its application as a plasma-facing material in future fusion reactors. Heavy ion irradiation study offers expedited damage accumulation and thus allows irradiation tolerance property prediction in a short period of time. In this study we explore in situ Kr ion irradiation of W nanolaminates at 800 °C and uncover a complex interplay of irradiation-induced transformations, including the emergence of stacking faults, nanograin formation, and generation of thickened grain boundaries defined as GB regions with the thickness of up to 15 nm. High-resolution transmission electron microscopy studies reveal the formation of extended planar faults and a metastable hexagonal close-packed (hcp) phase within the body-centered cubic (bcc) matrix. These structural transitions are facilitated by irradiation-induced shear. They are further stabilized by the presence of oxygen, which promotes stacking fault formation and vacancy trapping. Elevated temperatures enhance defect mobility, enabling dynamic recrystallization into ultrafine grains and the thickening of GBs due to defect absorption and impurity segregation. These findings unravel non-equilibrium phase transformation pathways in heavy ion irradiated W and highlight the critical role of impurity-mediated defect dynamics in governing its radiation tolerance properties.

Wazeer, Adil [Purdue University]↗
DOE OSTI2025

Crystal Growth Scale-Up and Stability of RbSr 2 X 5 :Eu Scintillators

The discovery and development of new scintillation materials support national security needs. In these applications, large volumes of scintillator crystals are needed to achieve efficient screening for contraband. Therefore, one important step in the discovery of new scintillator compositions is testing their feasibility for scale-up and their stability. In this work, high-quality Ø22 mm crystals of two new scintillators RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu were grown via the Vertical Bridgman method, and their scintillation properties were characterized. Here, the Ø22 mm RbSr 2 I 5 :Eu crystals could be grown with fast translation rates up to 3.5 mm/h. Both RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu had high scintillation performance, including light yields of 46,000 and 61,000 ph/MeV, respectively, for Ø22 × 35 mm crystals. Additionally, properties related to physical stability were investigated, including the coefficients of thermal expansion via high-temperature X-ray diffraction (HTXRD) as well as moisture sensitivity. HTXRD confirmed the absence of solid–solid phase transitions and showed that RbSr 2 Br 5 had minimal thermal expansion anisotropy compared to RbSr 2 I 5 and some other inorganic metal halide scintillators, which favors the growth of large-sized crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Charm and bottom hadrons in hot hadronic matter

Heavy quarks, and the hadrons containing them, are excellent probes of the QCD medium formed in high-energy heavy-ion collisions, as they provide essential information on the transport properties of the medium and how quarks color-neutralize into hadrons. Large theoretical and phenomenological efforts have been dedicated thus far to assess the diffusion of charm and bottom quarks in the quark–gluon plasma and their subsequent hadronization into heavy-flavor (HF) hadrons. However, the fireball formed in heavy-ion collisions also features an extended hadronic phase, and therefore any quantitative analysis of experimental observables needs to account for the rescattering of charm and bottom hadrons. This is further reinforced by the presence of a QCD cross-over transition and the notion that the interaction strength is maximal in the vicinity of the pseudo-critical temperature. We review existing approaches for evaluating the interactions of open HF hadrons in a hadronic heat bath and the pertinent results for scattering amplitudes, spectral functions and transport coefficients. While most of the work to date has focused on D -mesons, we also discuss excited states as well as HF baryons and the bottom sector. Both the HF hadro-chemistry and bottom observables will play a key role in future experimental measurements. We also conduct a survey of transport calculations in heavy-ion collisions that have included effects of hadronic HF diffusion and assess its impact on various observables.

effective field theories↗
DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗
DOE OSTI2024

Unravelling ultralow thermal conductivity in perovskite Cs2AgBiBr6: dominant wave-like phonon tunnelling and strong anharmonicity

Abstract Understanding the lattice dynamics and heat transport physics in the lead-free halide double perovskites remains an outstanding challenge due to their lattice dynamical instability and strong anharmonicity. In this work, we investigate the microscopic mechanisms of anharmonic lattice dynamics and thermal transport in lead-free halide double perovskite Cs 2 AgBiBr 6 from first principles. We combine self-consistent phonon calculations with bubble diagram correction and a unified theory of lattice thermal transport that considers both the particle-like phonon propagation and wave-like tunnelling of phonons. An ultra-low thermal conductivity at room temperature (~0.21 Wm −1 K −1 ) is predicted with weak temperature dependence( ~ T −0.34 ), in sharp contrast to the conventional ~T −1 dependence. Particularly, the vibrational properties of Cs 2 AgBiBr 6 are featured by strong anharmonicity and wave-like tunnelling of phonons. Anharmonic phonon renormalization from both the cubic and quartic anharmonicities are found essential in precisely predicting the phase transition temperature in Cs 2 AgBiBr 6 while the negative phonon energy shifts induced by cubic anharmonicity has a significant influence on particle-like phonon propagation. Further, the contribution of the wave-like tunnelling to the total thermal conductivity surpasses that of the particle-like propagation above around 310 K, indicating the breakdown of the phonon gas picture conventionally used in the Peierls-Boltzmann Transport Equation. Importantly, further including four-phonon scatterings is required in achieving the dominance of wave-like tunnelling, as compared to the dominant particle-like propagation channel when considering only three-phonon scatterings. Our work highlights the importance of lattice anharmonicity and wave-like tunnelling of phonons in the thermal transport in lead-free halide double perovskites.

Chemistry↗
DOE OSTI2025

Patchy nanoparticles by atomic stencilling

Stencilling, in which patterns are created by painting over masks, has ubiquitous applications in art, architecture and manufacturing. Modern, top-down microfabrication methods have succeeded in reducing mask sizes to under 10 nm, enabling ever smaller microdevices as today’s fastest computer chips. Meanwhile, bottom-up masking using chemical bonds or physical interactions has remained largely unexplored, despite its advantages of low cost, solution-processability, scalability and high compatibility with complex, curved and three-dimensional (3D) surfaces. Here we report atomic stencilling to make patchy nanoparticles (NPs), using surface-adsorbed iodide submonolayers to create the mask and ligand-mediated grafted polymers onto unmasked regions as ‘paint’. We use this approach to synthesize more than 20 different types of NP coated with polymer patches in high yield. Polymer scaling theory and molecular dynamics (MD) simulation show that stencilling, along with the interplay of enthalpic and entropic effects of polymers, generates patchy particle morphologies not reported previously. These polymer-patched NPs self-assemble into extended crystals owing to highly uniform patches, including different non-closely packed superlattices. We propose that atomic stencilling opens new avenues in patterning NPs and other substrates at the nanometre length scale, leading to precise control of their chemistry, reactivity and interactions for a wide range of applications, such as targeted delivery, catalysis, microelectronics, integrated metamaterials and tissue engineering.

36 MATERIALS SCIENCE↗