Search NASA⌕ Search

SEARCH · Search NASA

Results for “ISOTROPISM”

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.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗
DOE OSTI2021

Strong Magnetocrystalline Anisotropy Arising from Metal–Ligand Covalency in a Metal–Organic Candidate for 2D Magnetic Order

Layered metal–organic frameworks are promising candidates for new two-dimensional magnets, as their synthetic programmability of these materials can provide a route to diverse structural and electronic properties. However, such framework materials typically lack the heavy elements that engender magnetocrystalline anisotropy in the monolayer ferromagnets reported to date. Alternative sources of magnetic anisotropy are therefore needed in these materials. Here, we report the synthesis of single crystals of the framework material (NMe 4 ) 2 [Fe 2 L 3 ] (H 2 L = 3,6-dichloro-2,5-dihydroxybenzoquinone) and evaluate the angular dependence of its magnetic properties. Oriented-crystal magnetization measurements reveal strong uniaxial anisotropy, where the easy axis is aligned with the crystallographic c axis. While the spin carriers of this structure are isotropic $S = ^5/_2$ Fe III metal centers and $S = ^1/_2$ organic linkers, the anisotropy energy of the framework material is comparable to that of reported 2D ferromagnets. Density functional theory calculations indicate that the observed magnetocrystalline anisotropy arises from ligand-to-metal charge transfer that enhances the magnetic anisotropy of the otherwise isotropic Fe centers, suggesting that metal–ligand covalency can be utilized as a general additive for the development of 2D magnets. These results in this study show the possibility for (NMe 4 ) 2 [Fe 2 L 3 ] to retain magnetic order down to the 2D monolayer limit. In addition, the combination of large magnetic anisotropy and semiconducting character in (NMe 4 ) 2 [Fe 2 L 3 ] highlights its potential as a new 2D magnetic semiconductor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Probing the Electronic Structure of a Thorium Nitride Complex by Solid-State 15N NMR Spectroscopy

The solid-state 15N NMR powder spectra of the thorium nitride complex, [K(18-crown-6)(THF)2][(R2N)3Th(m-15N)(Th(NR2)3] ([K][1-15N], R = SiMe3), and the thorium amide complex, [Th(NR2)3(15NH2)] (2-15N) were recorded. The spectrum for [K][1-15N] represents the first reported solid-state 15N NMR data for an actinide nitride complex. The experimentally measured tensor spans are 847 ppm for [K][1-15N] and 237 ppm for 2-15N. Both shielding tensors exhibit a near-zero asymmetry parameter, which for [K][1-15N] is consistent with a local rotational symmetry of its 15N-labelled nitride ligand. For 2-15N, the lack of asymmetry can be rationalized by a quasi-free Th-NH2 bond rotation in the solid-state. DFT calculations overestimate the tensor span somewhat for [K][1-15N], but provide isotropic shifts in good agreement with both the solid-state and solution values for both complexes. Natural localized molecular orbital (NLMO) analyses of the nuclear shielding reveal that the larger tensor span in [K][1-15N] vs. 2-15N is primarily a consequence of more pronounced covalency of the N-Th bonds, and large spin-orbit coupling due to significant Th 5f orbital contribution to those bonds, impacting the principal components of the shielding tensor perpendicular to the Th-N-Th axis. Overall, our analysis confirms the involvement of the 5f orbitals in Th-N multiple bonds, and further demonstrates the value of solid-state NMR spectroscopy for interrogating actinide-ligand bonding.

actinides, covalent bonding, Nuclear Magnetic Reso↗
DOE OSTI2025

Precursor Identity and Surfactant Concentration Influence Shape of UO 2 Nanoparticles

A fundamental step toward studying the unique properties of actinide nanomaterials is control over the shape of actinide nanoparticles. Toward this goal, this work demonstrates the effects of precursor identity and surfactant concentration on the shape of uranium dioxide (UO 2 ) nanoparticles. UO 2 nanoparticles were synthesized by thermal decomposition of different precursors in the presence of oleylamine and oleic acid as surfactants. The size, shape, phase, and chemical composition of the nanoparticles was evaluated using transmission electron microscopy (TEM), energy dispersive X-ray spectroscopy (EDS), electron diffraction (ED), and U L 3 -edge X-ray absorption fine-structure (XAFS) spectroscopy. Anisotropic UO 2 nanocubes of ∼ 4 nm were obtained only with low surfactant concentrations, while increasing the surfactant concentration resulted in formation of nanoparticles with an isotropic, sphere morphology. The importance of precursor identity was also investigated by employing U(hfa) 4 , U(acac) 4 , UO 2 (acac) 2 , and UO 2 (hfa) 2 ·H 2 O (where hfa = hexafluoracetylacetone and acac = acetylacetone) as precursors. The nanocube morphology was only observed when U(hfa) 4 was used as a precursor. XAS allowed for comparison of the disorder in anisotropic vs isotropic UO 2 nanoparticles.

Nanocubes↗
DOE OSTI2021

Cylindrical Confinement of Nanocolloidal Cholesteric Liquid Crystal

The organization of nanocolloidal liquid crystals in constrained geometries has fundamental and practical importance, since under confinement, liquid crystals contain stable topological defects that can serve as templates for nanoparticle organization. Three-dimensional confinement of cholesteric (Ch) liquid crystals formed by cellulose nanocrystals (CNCs) have been extensively studied; however, their two-dimensional confinement remains underinvestigated. In this work, we report the results of systematic experimental studies of two-dimensional confinement of Ch–CNC liquid crystal in cylindrical capillaries with varying inner diameters. Confinement resulted in phase separation of the Ch–CNC liquid crystal into a Ch shell formed by concentric CNC pseudolayers with the helicoidal axis perpendicular to the inner surface of the capillary walls, and a micrometer-diameter isotropic core thread running parallel to the long axis of the capillary. The morphology of the confined Ch–CNC liquid crystal varied when progressively increasing the degree confinement. Finally, we show that phase separation of the Ch–CNC liquid crystal into a Ch shell and an isotropic core is preserved in flexible capillary tubing, suggesting the applicability of this system for the fabrication of flexible optical waveguides.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Detecting Reactive Products in Carbon Capture Polymers with Chemical Shift Anisotropy and Machine Learning

Aminopolymers are attractive sorbents for CO 2 direct air capture applications due to their high density of amine groups, which can readily react with atmospheric levels of CO 2 to form chemisorbed species. The identity of these chemisorbed species and the functional groups that form upon oxidative degradation depends on both material properties and processing conditions, forming a variety of carbonyl-type sites such as ammonium carbamates, bicarbonates, carbonates, carbamic acids, ureas, and amides. 13 C solid-state nuclear magnetic resonance (NMR) is often used to help elucidate the identity of these reacted species, but it is challenging due to the narrow chemical shift range of carbonyl sites. Herein, we demonstrate the application of a two-dimensional (2D) chemical shift anisotropy (CSA) recoupling pulse sequence (ROCSA) to obtain CSA tensor values at each isotropic chemical shift, overcoming limitations of isotropic peak resolution. CSA tensor values describe the local chemical environment and can readily differentiate between the chemisorbed and degradation products. To aid identification, we also developed a k-nearest neighbor (kNN) classification model to distinguish the functional groups via their CSA tensor parameters. This methodology was demonstrated on poly(ethylenimine) in γ-Al 2 O 3 exposed to CO 2 and showed that the chemisorbed products are ammonium carbamate and a mixed carbamate–carbamic acid species. The sample was analyzed again after desorption at 100 °C inducing mild degradation, and the remaining products were strongly bound carbamate and urea species. In conclusion, the combination of 2D CSA measurements coupled with a kNN classification model enhances the ability to accurately identify chemisorbed or degradation products in complex carbon capture materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Modulating the Extent of Anisotropic Cuprophilicity via High Pressure with Piezochromic Luminescence Sensitization

We report metallophilicity has been widely studied as a fundamental supramolecular interaction. However, the extent and directionality thereof remain controversial. A major obstacle lies in the difficulty to separately control the geometry and chemical composition. Herein, we address this challenge by modulating metallophilicity with mechanical pressure. Using a multinuclear Cu(I) complex as model system, we report anomalous anisotropies of (supra)molecular structures, vibrations, and interaction energies upon isotropic compression as well as concomitant (essentially turn-on) piezochromic luminescence enhancement with ~10 3 modulation. The in situ characterizations indicate opposite behaviors of contact distances and cuprophilic interactions for intermolecular vs intramolecular Cu-Cu pairs under pressure. Theoretical calculations break down the attractive and repulsive forces associated with cuprophilicity, its spontaneous 4p-3d hybridization origin, and direction-dependent interaction strength. The use of isotropic mechanical force reveals the intrinsic anisotropy of metallophilicity in multinuclear systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Extraordinary Magnetic Response of an Anisotropic 2D Antiferromagnet via Site Dilution

A prominent characteristic of 2D magnetic systems is the enhanced spin fluctuations, which reduce the ordering temperature. We report that a magnetic field of only 1000th of the Heisenberg superexchange interaction can induce a crossover, which for practical purposes is the effective ordering transition, at temperatures about 6 times the Néel transition in a site-diluted two-dimensional anisotropic quantum antiferromagnet. Such a strong magnetic response is enabled because the system directly enters the antiferromagnetically ordered state from the isotropic disordered state, skipping the intermediate anisotropic stage. The underlying mechanism is achieved on a pseudospin-half square lattice realized in the [(SrIrO3)1/(SrTiO3)2] superlattice thin film that is designed to linearly couple the staggered magnetization to external magnetic fields by virtue of the rotational symmetry-preserving Dzyaloshinskii–Moriya interaction. Furthermore, our model analysis shows that the skipping of the anisotropic regime despite finite anisotropy is due to the enhanced isotropic fluctuations under moderate dilution.

2D antiferromagnet↗
DOE OSTI2021

Record-Low Thermal Boundary Resistance between Diamond and GaN-on-SiC for Enabling Radiofrequency Device Cooling

The implementation of 5G-and-beyond networks requires faster, high-performance, and power-efficient semiconductor devices, which are only possible with materials that can support higher frequencies. Gallium nitride (GaN) power amplifiers are essential for 5G-and-beyond technologies since they provide the desired combination of high frequency and high power. These applications along with terrestrial hub and backhaul communications at high power output can present severe heat removal challenges. The cooling of GaN devices with diamond as the heat spreader has gained significant momentum since device self-heating limits GaN’s performance. However, one of the significant challenges in integrating polycrystalline diamond on GaN devices is maintaining the device performance while achieving a low diamond/GaN channel thermal boundary resistance. In this study, we achieved a record-low thermal boundary resistance of around 3.1 ± 0.7 m 2 K/GW at the diamond/Si 3 N 4 /GaN interface, which is the closest to theoretical prediction to date. The diamond was integrated within ~1 nm of the GaN channel layer without degrading the channel’s electrical behavior. Furthermore, we successfully minimized the residual stress in the diamond layer, enabling more isotropic polycrystalline diamond growth on GaN with thicknesses >2 μm and a ~1.9 μm lateral grain size. More isotropic grains can spread the heat in both vertical and lateral directions efficiently. Using transient thermoreflectance, the thermal conductivity of the grains was measured to be 638 ± 48 W/m K, which when combined with the record-low thermal boundary resistance makes it a leading-edge achievement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Exploring the Role of the Central Carbide of the Nitrogenase Active-Site FeMo-cofactor through Targeted 13 C Labeling and ENDOR Spectroscopy

Mo-dependent nitrogenase is a major contributor to global biological N 2 reduction, which sustains life on Earth. Its multi-metallic active-site FeMo-cofactor (Fe 7 MoS 9 C-homocitrate) contains a carbide (C 4– ) centered within a trigonal prismatic CFe 6 core resembling the structural motif of the iron carbide, cementite. The role of the carbide in FeMo-cofactor binding and activation of substrates and inhibitors is unknown. Here, to explore this role, the carbide has been in effect selectively enriched with 13 C, which enables its detailed examination by ENDOR/ESEEM spectroscopies. 13 C-carbide ENDOR of the S = 3/2 resting state (E 0 ) is remarkable, with an extremely small isotropic hyperfine coupling constant, C a = +0.86 MHz. Turnover under high CO partial pressure generates the S = 1/2 hi-CO state, with two CO molecules bound to FeMo-cofactor. This conversion surprisingly leaves the small magnitude of the 13 C carbide isotropic hyperfine-coupling constant essentially unchanged, C a = –1.30 MHz. This indicates that both the E 0 and hi-CO states exhibit an exchange-coupling scheme with nearly cancelling contributions to C a from three spin-up and three spin-down carbide-bound Fe ions. In contrast, the anisotropichyperfine coupling constant undergoes a symmetry change upon conversion of E 0 to hi-CO that may be associated with bonding and coordination changes at Fe ions. In combination with the negligible difference between CFe 6 core structures of E 0 and hi-CO, these results suggest that in CO-inhibited hi-CO the dominant role of the FeMo-cofactor carbide is to maintain the core structure, rather than to facilitate inhibitor binding through changes in Fe-carbide covalency or tretching/breaking of carbide–Fe bonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Pitch Angle Dependence of Electron and Ion Flux Changes During Local Magnetic Dipolarization Inside Geosynchronous Orbit

In this study we have statistically examined flux changes of 1–1,000-keV electrons, and hydrogen, helium, and oxygen ions during local magnetic dipolarization inside geosynchronous orbit, with a focus on their pitch angle dependence. Here, using 144 dipolarization events that were selected in a previous study with Van Allen Probes observations in 2012–2016, we have performed a superposed epoch analysis of differential flux changes after the dipolarization onset for each species. On average the electron flux increases primarily around pitch angle (α) = 90° at >80 keV and almost isotropically at 10–50 keV. The electron flux at <5 keV increases at α = 0° and 180° but slightly decreases (or remains unchanged) at α = 90°. On the other hand, the ion flux at >80 keV increases around α = 90°, while at <30 keV it decreases nearly independent of pitch angle. Only the low-energy (<5 keV) helium and oxygen ion flux changes indicate a strong field-aligned enhancement, which could be attributed to their outflows from the topside ionosphere. After the dipolarization onset, >5-keV electron and >30-keV ion fluxes exhibit a perpendicular anisotropy, which is most pronounced for 50–200-keV electrons. Both distributions become more isotropic with decreasing energy. A noticeable field-aligned anisotropy is seen for <5-keV ion fluxes. These statistical pitch angle-dependent electron and ion properties during dipolarizations may be explained by combining various processes, including adiabatic acceleration and/or transport, wave-particle interactions, and ion outflow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer

Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The Viscoplastic Behavior of Natural Hydrate‐Bearing Sandy‐Silts Under Uniaxial Strain Compression ( K 0 Loading)

Abstract The in‐situ stress state and geomechanical properties of hydrate‐bearing sediments impact hydrate formation and gas production strategies. We explore the uniaxial strain compression and stress evolution of natural hydrate‐bearing sandy‐silts from Green Canyon Block 955 in the deep‐water Gulf of Mexico. We performed constant rate of strain uniaxial strain experiments, interrupted by periods where we held the axial stress constant, to explore the vertical deformation and the evolution of the ratio of lateral to axial effective stress ( K 0 ) with time. The hydrate‐bearing sandy‐silt is stiffer and has a larger K 0 than the equivalent hydrate‐free sediment upon loading. During stress holds, the void ratio decreases sigmoidally with the log of time, and K 0 converges to isotropic conditions. We interpret that during loading, the hydrate bears the load and deforms. With time, the hydrate redistributes the load and K 0 increases. We used a viscoelastic model to describe the behavior. The model accurately captures deformation and K 0 trends but does not reproduce all the complex interactions of the hydrate with the porous skeleton. We anticipate that viscous effects within hydrate sediments will impact reservoir compression and stresses during production (hours to days), result in isotropic stress state over geological timescales, and explain the creeping movement in submarine landslides.

58 GEOSCIENCES↗
DOE OSTI2023

Orbital character of the spin-reorientation transition in TbMn 6 Sn 6

Ferromagnetic (FM) order in a two-dimensional kagome layer is predicted to generate a topological Chern insulator without an applied magnetic field. The Chern gap is largest when spin moments point perpendicular to the kagome layer, enabling the capability to switch topological transport properties, such as the quantum anomalous Hall effect, by controlling the spin orientation. In TbMn 6 Sn 6 , the uniaxial magnetic anisotropy of the Tb 3+ ion is effective at generating the Chern state within the FM Mn kagome layers while a spin-reorientation (SR) transition to easy-plane order above T SR = 310 K provides a mechanism for switching. Here, we use inelastic neutron scattering to provide key insights into the fundamental nature of the SR transition. The observation of two Tb excitations, which are split by the magnetic anisotropy energy, indicates an effective two-state orbital character for the Tb ion, with a uniaxial ground state and an isotropic excited state. The simultaneous observation of both modes below T SR confirms that orbital fluctuations are slow on magnetic and electronic time scales < ps and act as a spatially-random orbital alloy. A thermally-driven critical concentration of isotropic Tb ions triggers the SR transition.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Superior zero thermal expansion dual-phase alloy via boron-migration mediated solid-state reaction

Rapid progress in modern technologies demands zero thermal expansion (ZTE) materials with multi-property profiles to withstand harsh service conditions. Thus far, the majority of documented ZTE materials have shortcomings in different aspects that limit their practical utilization. Here, we report on a superior isotropic ZTE alloy with collective properties regarding wide operating temperature windows, high strength-stiffness, and cyclic thermal stability. A boron-migration-mediated solid-state reaction (BMSR) constructs a salient “plum pudding” structure in a dual-phase Er-Fe-B alloy, where the precursor ErFe 10 phase reacts with the migrated boron and transforms into the target Er 2 Fe 14 B (pudding) and α-Fe phases (plum). The formation of such microstructure helps to eliminate apparent crystallographic texture, tailor and form isotropic ZTE, and simultaneously enhance the strength and toughness of the alloy. These findings suggest a promising design paradigm for comprehensive performance ZTE alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Two-component nematic superconductivity in 4Hb-TaS 2

Most superconductors have an isotropic, single component order parameter and are well described by the standard (BCS) theory for superconductivity. Unconventional, multiple-component superconductors are exceptionally rare and are much less understood. Here, we combine scanning tunneling microscopy and angle-resolved macroscopic transport for studying the candidate chiral superconductor, 4Hb-TaS 2 . We reveal quasi-periodic one-dimensional modulations in the tunneling conductance accompanied by two-fold symmetric superconducting critical field. The strong modulation of the in-plane critical field, H c2 , points to a nematic, unconventional order parameter. However, the imaged vortex core is isotropic at low temperatures. We suggest a model that reconciles this apparent discrepancy and takes into account previously observed spontaneous time-reversal symmetry breaking at low temperatures. The model describes a competition between a dominating chiral superconducting order parameter and a nematic one. The latter emerges close to the normal phase. Our results strongly support the existence of two-component superconductivity in 4Hb-TaS 2 and can provide valuable insights into other systems with coexistent charge order and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Evidence of a distinct collective mode in Kagome superconductors

The collective modes of the superconducting order parameter fluctuation can provide key insights into the nature of the superconductor. Recently, a family of superconductors has emerged in non-magnetic kagome materials AV 3 Sb 5 (A = K, Rb, Cs), exhibiting fertile emergent phenomenology. However, the collective behaviors of Cooper pairs have not been studied. Here, we report a distinct collective mode in CsV 3-x Ta x Sb 5 using scanning tunneling microscope/spectroscopy. The spectral line-shape is well-described by one isotropic and one anisotropic superconducting gap, and a bosonic mode due to electron-mode coupling. With increasing x, the two gaps move closer in energy, merge into two isotropic gaps of equal amplitude, and then increase synchronously. The mode energy decreases monotonically to well below 2Δ and survives even after the charge density wave order is suppressed. We propose the interpretation of this collective mode as Leggett mode between different superconducting components or the Bardasis-Schrieffer mode due to a subleading superconducting component.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Many-Body Potential for Simulating Self-Assembly of Polymer-Grafted Nanoparticles in a Polymer Matrix

Many-body interactions between polymer-grafted nanoparticles (NPs) play a key role in promoting their assembly into low-dimensional structures within polymer melts, even when the particles are spherical and isotropically grafted. However, capturing such interactions in simulations of NP assembly is very challenging because explicit modeling of the polymer grafts and melt chains is highly computationally expensive, even using coarse-grained models. Here, we develop a many-body potential for describing the effective interactions between spherical polymer-grafted NPs in a polymer matrix through a machine-learning approach. The approach involves using permutationally invariant polynomials to fit two- and three-body interactions derived from the potential of mean force calculations. The potential developed here reduces the computational cost by several orders of magnitude, thereby, allowing us to explore assembly behavior over large length and time scales. We show that the potential not only reproduces previously known assembled phases such as 1D strings and 2D hexagonal sheets, which generally cannot be achieved using isotropic two-body potentials, but can also help discover interesting phases such as networks, clusters, and gels. We demonstrate how each of these assembly morphologies intrinsically arises from a competition between two- and three-body interactions. Our approach for deriving many-body effective potentials can be readily extended to other colloidal systems, enabling researchers to make accurate predictions of their behavior and dissect the role of individual interaction energy terms of the overall potential in the observed behavior.

36 MATERIALS SCIENCE↗