Search NASA⌕ Search

SEARCH · Search NASA

Results for “small bodies”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

We report a spin-free formulation of the multireference (MR) driven similarity renormalization group (DSRG) based on the ensemble normal ordering of Mukherjee and Kutzelnigg (1997). This ensemble averages over all microstates of a given total spin quantum number, and therefore, it is invariant with respect to SU(2) transformations. As such, all equations may be reformulated in terms of spin-free quantities and they closely resemble those of spin-adapted closed-shell coupled cluster (CC) theory. The current implementation is used to assess the accuracy of various truncated MR-DSRG methods (perturbation theory up to third order and iterative methods with single and double excitations) in computing the constants of 33 first-row diatomic molecules. The accuracy trends for these first-row diatomics are consistent with our previous benchmark on a small subset of closed-shell diatomic molecules. We then present the first MR-DSRG application on transition-metal complexes by computing the spin splittings of the [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ molecules. A focal point analysis (FPA) shows that third-order perturbative corrections are essential to achieve reasonably converged energetics. The FPA based on the linearized MR-DSRG theory with one- and two-body operators and up to a quintuple- ζ basis set predicts the spin splittings of [Fe(H 2 O) 6 ] 2+ and [Fe(NH 3 ) 6 ] 2+ to be –35.7 and –17.1 kcal mol –1 , respectively, showing good agreement with the results of local CC theory with singles, doubles, and perturbative triples.

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

A time-dependent momentum-resolved scattering approach to core-level spectroscopies

While new light sources allow for unprecedented resolution in experiments with X-rays, a theoretical understanding of the scattering cross-section is lacking. In the particular case of strongly correlated electron systems, numerical techniques are quite limited, since conventional approaches rely on calculating a response function (Kramers-Heisenberg formula) that is obtained from a perturbative analysis of scattering processes in the frequency domain. This requires a knowledge of a full set of eigenstates in order to account for all intermediate processes away from equilibrium, limiting the applicability to small tractable systems. In this work, we present an alternative paradigm, recasting the problem in the time domain and explicitly solving the time-dependent Schrödinger equation without the limitations of perturbation theory: a faithful simulation of the scattering processes taking place in actual experiments, including photons and core electrons. We show how this approach can yield the full time and momentum resolved Resonant Inelastic X-Ray Scattering (RIXS) spectrum of strongly interacting many-body systems. We demonstrate the formalism with an application to Mott insulating Hubbard chains using the time-dependent density matrix renormalization group method, which does not require a priory knowledge of the eigenstates and can solve very large systems with dozens of orbitals. This approach can readily be applied to systems out of equilibrium without modification and generalized to other spectroscopies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Planar Systems for Quantum Information

This project aims to develop two‐dimensional (2D) moiré materials as a quantum simulator to implement model Hamiltonians and their phase diagrams. Progress in quantum information science (QIS) requires the development of advanced quantum materials systems. The rich family of layered van der Waals materials and their heterostructures present opportunities to create previously unrealized types of applications for QIS. Specifically, when two layers of van der Waals materials are overlaid with a small twist angle or/and lattice mismatch, a moiré superlattice with a period of about ten nanometers is formed. This provides a periodic trapping potential for electrons. Electrons can tunnel between the traps and repel each other by their mutual Coulomb interactions. The platform of 2D moiré materials provides many attractive features, including tunability of length and energy scales, charge density, and even lattice symmetry. It presents new possibilities for realizing quantum simulation of the many-body physics in a solid-state platform. This integrated team of six investigators seeks to develop relevant theoretical treatments to link ab-initio studies of 2D moiré materials to model Hamiltonians and to evaluate correlated phases predicted by these model Hamiltonians in the relevant regimes. On the experimental side, the team aims to develop methods to realize a homogeneous and highly controlled potential landscape for the electrons and to initiate, protect, and measure their quantum many-body states.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Formation, cooling history and age of impact events on the IIE iron parent body: Evidence from the Miles meteorite

Most iron meteorites formed in planetary cores during differentiation, but the IIE iron meteorites have chemical and physical features that are inconsistent with this origin. By combining mineral chemistry, mineral modes and three-dimensional petrography, we reconstruct the bulk chemistry of the felsic silicate-bearing Miles IIE iron meteorite and demonstrate that the silicate inclusion compositions are similar to partial melts produced experimentally from an H chondrite composition. We use the reconstructed bulk composition, mineralogy and thermodynamic modelling to show that melting above ~ 1200°C under reducing conditions formed metal (Fe-Ni alloy) and felsic silicate partial melts. Upon cooling, the melts crystallized Mg-rich pyroxenes, Na- and K-rich feldspars, and tridymite. Importantly, this mechanism enriches cosmochemically volatile elements (i.e., those with a 50% condensation temperature of ~430–830°C, like Na and K) to the level found in the felsic silicate inclusions. The presence of crystallographically disordered srilankite (only stable above 1160°C) and an absence of Widmanstätten texture require both high peak temperatures and rapid cooling, which cannot be explained by core formation. Instead, they point to small melt volumes, a transient heat pulse, and small thermal mass, and imply efficient physical segregation of silicate and metallic melts through buoyancy separation followed by rapid cooling that arrested the separation of metal and silicate liquid phases. In situ 207 Pb/ 206 Pb age of 4542.3 ± 4.0 Ma in Zr-oxide and phosphate minerals dates the melting event that formed the silicate inclusions. This age aligns with the earliest ages found in other IIE iron meteorite silicates and requires a heating event ~25 million years after the solar system formed. We found 39 Ar/ 40 Ar ages of 3495 ± 52 Ma (low-T) and 4303 ± 7 Ma (high-T) in a K-feldspar grain, with the 3495 Ma age aligning with later thermal events recorded in other IIE iron meteorites. Dating reveals the complex petrogenetic and thermal history of Miles and the IIE iron meteorites. This is the first IIE iron meteorite found to record evidence of heating at 4.5 and 3.5 Ga likely from impact events. We propose that high-velocity impact(s) into an iron-rich, porous chondritic parent body at ~4.54 Ga produced immiscible metal and silicate melts that cooled rapidly and trapped low density silicate inclusions within high density metal. Other IIE irons that formed at lower peak temperatures (900–1000°C) contain chondritic silicate inclusions and relict chondrules, supporting this conceptual model. In conclusion, this petrogenesis is consistent with thermodynamic modelling, experimental data and the wide range of peak temperatures and cooling rates observed in the IIE iron meteorites.

58 GEOSCIENCES↗
DOE OSTI2024

Liquid state theory of the structure of model polymerized ionic liquids

We employ polymer integral equation theory to study a simplified model of semiflexible polymerized ionic liquids (PolyILs) that interact via hard core repulsions and short range screened Coulomb interactions. The multi-scale structure in real and Fourier space of PolyILs (ions chosen to mimic Li, Na, K, Br, PF 6 , and TFSI) are determined as a function of melt density, Coulomb interaction strength, and ion size. Comparisons with a homopolymer melt, a neutral polymer–solvent-like athermal mixture, and an atomic ionic liquid are carried out to elucidate the distinct manner that ions mediate changes of polymer packing, the role of excluded volume effects, and the influence of chain connectivity, respectively. The effect of Coulomb strength depends in a rich manner on ion size and density, reflecting the interplay of steric packing, ion adsorption, and charge layering. Ion-mediated bridging of monomers is found, which intensifies for larger ions. Intermediate range charge layering correlations are characterized by a many-body screening length that grows with PolyIL density, cooling, and Coulomb strength, in disagreement with Debye–Hückel theory, but in accord with experiments. Qualitative differences in the collective structure, including an ion-size-dependent bifurcation of the polymer structure factor peak and pair correlation function, are predicted. The monomer cage order parameter increases significantly, but its collective ion counterpart decreases, as ions become smaller. Such behaviors allow one to categorize PolyILs into two broad classes of small and large ions. Furthermore, dynamical implications of the predicted structural results are qualitatively discussed.

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

Non-Hermitian Quantum Mechanics Approach for Extracting and Emulating Continuum Physics Based on Bound-State-like Calculations

Here, this Letter introduces a unified emulation framework for studying continuum physics in finite quantum systems. Using a reduced basis method, we construct powerful emulators for the inhomogeneous Schrödinger equation that operate in a combined parameter space of complex energy (𝐸) and other inputs (𝜽). Within the space, the emulators simultaneously perform analytical continuation in 𝐸—extracting continuum physics from numerically simpler bound-state-like calculations—and interpolate this entire process across 𝜽. This yields a small, non-Hermitian system whose properties (e.g., resonances and scattering observables) can be rapidly predicted for any 𝜽. Crucially, the complex-𝐸 emulation provides a pathway to compute continuum observables for complex systems where advanced bound-state methods exist but direct continuum calculations are yet to be developed, while the 𝜽 emulation enables rapid parameter-space exploration and can be adapted to accelerate other existing continuum calculations. Demonstrations with two- and three-body systems highlight the method’s effectiveness and suggest its connection to (near-)optimal rational approximation. This Letter presents the key results, with further details reserved for a companion paper.

ab initio calculations↗
DOE OSTI2022

Beware of fake ν ’s: The effect of massive neutrinos on the nonlinear evolution of cosmic structure

Massive neutrinos suppress the growth of cosmic structure on small, nonlinear scales. It is thus often proposed that using statistics beyond the power spectrum can tighten constraints on the neutrino mass by extracting additional information from these nonlinear scales. We study the information content regarding neutrino mass at the field level, quantifying how much of this information arises from the difference in nonlinear evolution between a cosmology with one fluid [cold dark matter (CDM)] and two fluids ( CDM + neutrinos ). We do so by running two N -body simulations, one with and one without massive neutrinos, both with the same phases, and matching their linear power spectrum at a given low redshift. This effectively isolates the information encoded in the linear initial conditions from the nonlinear cosmic evolution. We demonstrate that, for k ≲ 1 h / Mpc , and for a single redshift, there is negligible difference in the real-space CDM field between the two simulations. This suggests that all the information regarding neutrino mass is in the linear power spectrum set by the initial conditions. Thus, any probe based on the CDM field alone will have negligible constraining power beyond that which exists at the linear level over the same range of scales. Consequently, any probe based on the halo field will contain little information beyond the linear power. We find similar results for the matter field responsible for weak lensing. We also demonstrate that there may be much information beyond the power spectrum in the 3D matter field; however, this is not observable in modern surveys via dark matter halos or weak lensing. Finally, we show that there is additional information to be found in redshift space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Minimal two-body quantum absorption refrigerator

Here, we study the phenomenon of absorption refrigeration, where refrigeration is achieved by heating instead of work, in two different setups: a minimal set up based on coupled qubits and two nonlinearly coupled resonators. Considering ZZ interaction between the two qubits, we outline the basic ingredients required to achieve cooling. Using local as well as global master equations, we observe that inclusion of an XX type term in the qubit-qubit coupling is detrimental to cooling. We compare the cooling effect obtained in the qubit case with that of nonlinearly coupled resonators (multilevel system) where the ZZ interaction translates to a Kerr-type nonlinearity. For small to intermediate strengths of nonlinearity, we observe that multilevel quantum systems, for example qutrits, give better cooling effect compared to the qubits. Using Keldysh nonequilibrium Green's function formalism, we go beyond first order sequential tunneling processes and study the effect of higher order processes on refrigeration. We find a reduced cooling effect compared to the master equation calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Hamiltonian simulation in the low-energy subspace

Abstract We study the problem of simulating the dynamics of spin systems when the initial state is supported on a subspace of low energy of a Hamiltonian H . This is a central problem in physics with vast applications in many-body systems and beyond, where the interesting physics takes place in the low-energy sector. We analyze error bounds induced by product formulas that approximate the evolution operator and show that these bounds depend on an effective low-energy norm of H . We find improvements over the best previous complexities of product formulas that apply to the general case, and these improvements are more significant for long evolution times that scale with the system size and/or small approximation errors. To obtain these improvements, we prove exponentially decaying upper bounds on the leakage to high-energy subspaces due to the product formula. Our results provide a path to a systematic study of Hamiltonian simulation at low energies, which will be required to push quantum simulation closer to reality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Galaxy-multiplet clustering from DESI DR2

We present an efficient estimator for higher-order galaxy clustering using small groups of nearby galaxies, or multiplets. Using the Luminous Red Galaxy (LRG) sample from the Dark Energy Spectroscopic Instrument (DESI) Data Release 2, we identify galaxy multiplets as discrete objects and measure their cross-correlations with the general galaxy field. Our results show that the multiplets exhibit stronger clustering bias as they trace more massive dark matter halos than individual galaxies. When comparing the observed clustering statistics with the mock catalogs generated from the N-body simulation AbacusSummit, we find that the mocks underpredict multiplet clustering despite reproducing the galaxy two-point auto-correlation reasonably well. This discrepancy indicates that the standard Halo Occupation Distribution (HOD) model is insufficient to describe the properties of galaxy multiplets, revealing the greater constraining power of this higher-order statistic on galaxy-halo connection and the possibility that multiplets are specific to additional assembly bias. We demonstrate that incorporating secondary biases into the HOD model improves agreement with the observed multiplet statistics, specifically by allowing galaxies to preferentially occupy halos in denser environments. Our results highlight the potential of utilizing multiplet clustering, beyond traditional two-point correlation measurements, to break degeneracies in models describing the galaxy-dark matter connection.

cosmology↗
DOE OSTI2021

Exploring the Effects of Emplacement Conditions on Explosion P/S Ratios across Local to Regional Distances

High-frequency (~> 2 Hz) seismic P/S amplitude ratios are well-established as a discriminant to distinguish between natural earthquakes and underground explosions at regional distances (~200–1500 km). As research shifts toward identifying lower-yield events, work has begun to investigate the potential of this discriminant for use at local distances (<200 km), in which initial results raise questions about its effectiveness. Here, we utilize data from several chemical explosion experiment series at the Nevada National Security Site in southern Nevada in the United States to study explosion Pg/Lg ratios across the range of local to regional distances. Additionall, the experiments are conducted over differing emplacement conditions, with contrasting geologies and a variety of yields and depths of burial, including surface explosions. We first establish the similarities of Pg/Lg ratios from chemical explosions to those from historic nuclear tests and conclude that, as previous data have suggested, chemical explosion ratios are good proxies for nuclear tests. We then examine Pg/Lg ratios from the new experiment series as functions of distance, yield, depth of burial, and scaled depth of burial (SDOB). At far-local and regional distances, we observe consistently higher ratios from hard-rock explosions compared to ones in a weaker dry alluvium medium, consistent with prior regional distance results. No other trends with yield, depth of burial, or SDOB are strongly evident. Scatter in the observed ratios is very high, particularly at the shortest event-to-station distances, suggesting that small-scale path effects play a significant role. On average, the local distance explosion Pg/Lg ratios show remarkable consistency across all the variations in emplacement. Explosion source models will need to reproduce these results.

58 GEOSCIENCES↗
DOE OSTI2020

Lattice QCD Inputs for nuclear double beta decay

Second order β -decay processes with and without neutrinos in the final state are key probes of nuclear physics and of the nature of neutrinos. Neutrinoful double- β decay is the rarest Standard Model process that has been observed and provides a unique test of the understanding of weak nuclear interactions. Observation of neutrinoless double- β decay would reveal that neutrinos are Majorana fermions and that lepton number conservation is violated in nature. While significant progress has been made in phenomenological approaches to understanding these processes, establishing a connection between these processes and the physics of the Standard Model and beyond is a critical task as it will provide input into the design and interpretation of future experiments. The strong-interaction contributions to double- β decay processes are non-perturbative and can only be addressed systematically through a combination of lattice Quantum Chromoodynamics (LQCD) and nuclear many-body calculations. In this review, current efforts to establish the LQCD connection are discussed for both neutrinoful and neutrinoless double- β decay. LQCD calculations of the hadronic contributions to the neutrinoful process $nn &rarr; ppe^-e^-\bar{v}_e\bar{v}_e$ and to various neutrinoless pionic transitions are reviewed, and the connections of these calculations to the phenomenology of double- β decay through the use of effective field theory (EFTs) is highlighted. At present, LQCD calculations are limited to small nuclear systems, and to pionic subsystems, and require matching to appropriate EFTs to have direct phenomenological impact. However, these calculations have already revealed qualitatively that there are terms in the EFTs that can only be constrained from double- β decay processes themselves or using inputs from LQCD. Finally, future prospects for direct calculations in larger nuclei are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

The Simons Observatory: A fully remote controlled calibration system with a sparse wire grid for cosmic microwave background telescopes

For cosmic microwave background (CMB) polarization observations, calibration of detector polarization angles is essential. We have developed a fully remote controlled calibration system with a sparse wire grid that reflects linearly polarized light along the wire direction. The new feature is a remote-controlled system for regular calibration, which has not been possible in sparse wire grid calibrators in past experiments. The remote control can be achieved by two electric linear actuators that load or unload the sparse wire grid into a position centered on the optical axis of a telescope between the calibration time and CMB observation. Furthermore, the sparse wire grid can be rotated by using a motor. A rotary encoder and a gravity sensor are installed on the sparse wire grid to monitor the wire direction. They allow us to achieve detector polarization angle calibration with an expected systematic error of 0.08°. The calibration system will be installed in small-aperture telescopes at Simons Observatory.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Structure, Hydration, and Interactions of Native and Hydrophobically Modified Phytoglycogen Nanoparticles

Phytoglycogen is a highly branched polymer of glucose produced as soft, compact nanoparticles by sweet corn. Properties such as softness, porosity, and mechanical integrity, combined with nontoxicity and biodegradability, make phytoglycogen nanoparticles ideal for applications involving the human body, ranging from skin moisturizing and rejuvenation agents in personal care formulations to functional therapeutics in biomedicine. To further broaden the range of applications, phytoglycogen nanoparticles can be chemically modified with hydrophobic species such as octenyl succinic anhydride (OSA). In this study, we present a self-consistent model of the particle structure, water content, and degree of chemical modification of the particles, as well as the emergence of well-defined interparticle spacings in concentrated dispersions, based on small-angle neutron scattering (SANS) measurements of aqueous dispersions of native phytoglycogen nanoparticles and particles that were hydrophobically modified using octenyl succinic anhydride (OSA) in both its protiated (pOSA) and deuterated (dOSA) forms. Measurements on native particles with reduced polydispersity have allowed us to refine the particle morphology, which is well described by a hairy particle (core-chain) geometry with short chains decorating the surface of the particles. The isotopic variants of OSA-modified particles enhanced the scattering contrast for neutrons, revealing lightly modified hairy chains for small degrees of substitution (DS) of OSA, and a raspberry particle geometry for the largest DS value, where the OSA-modified hairy chains collapse to form small seeds on the surface of the particles. This refined model of native and OSA-modified phytoglycogen nanoparticles establishes a quantitative basis for the development of new applications of this promising sustainable nanotechnology.

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

Harnessing Neural Networks for Elucidating X-ray Absorption Structure–Spectrum Relationships in Amorphous Carbon

Improved understanding of structural and chemical properties through local experimental probes, such as X-ray absorption near-edge structure (XANES) spectroscopy, is crucial for the understanding and design of functional materials. In recent years, significant advancements have been made in the development of data science approaches for the automated interpretation of XANES structure–spectrum relationships. However, existing studies have primarily focused on crystalline solids and small molecules, while fewer efforts have been devoted to disordered systems. Thus, in this work, we demonstrate the development of neural network models for predicting and interpreting XANES spectra of amorphous carbon (a-C) from local structural descriptors. Comparison between different structural descriptors expectedly shows that the inclusion of both bond length and bond angle information is necessary for an accurate prediction of the spectra. Among the descriptors considered in this work, we find that the local many-body tensor representation yields the highest accuracy and greatest interpretability so that it can be leveraged to understand the importance of structural motifs in determining XANES spectra. Furthermore, we also discuss performance of neural network models for predicting both local structure features, such as bond lengths and bond angles, and global chemical composition, such as the sp:sp 2 :sp 3 ratio.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Many-Body Effects in FeN4 Center Embedded in Graphene

We introduce a computational approach to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce generalized gradient approximation (GGA)+U results when ignoring inter-orbital Coulomb repulsion U ′ and Hund exchange J. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra- and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN 4 centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation, yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, U ′ changes the relative position of the peaks in the density of states, particularly on the iron d z 2 orbital, which is expected to affect the binding of ligands greatly.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Plasticity of Behavioral and Growth Responses to Two Feeding Regimes and Implications for the Domestication of Steelhead Trout

Abstract The effects of domestication selection on fitness are of serious concern in hatchery-reared anadromous salmonids, which often experience selection for rapid growth and associated deviations from natural agonistic and feeding behaviors. The high rations that are used to raise age-1 smolts and the reduced rations that are used to produce age-2 smolts potentially result in very different behavioral tactics and social dynamics of competition in hatchery tanks. Consequently, steelhead Oncorhynchus mykiss hatchery programs with recovery or conservation goals may need to change their aquaculture practices to minimize domestication selection. This study investigated how the feeding rate (proportion of presented food items that were successfully captured) and relative competitive ability (ability to monopolize food items in small groups) of age-1 and age-2 steelhead corresponded to feeding location (surface or bottom) and growth rate in hatchery rearing tanks. We raised steelhead from a common broodstock to smolt at age 1 or age 2 by manipulating their rations (age-1 fish received the high ration, age-2 the low ration). The fish in the age-1 treatment had a lower feeding rate than did those in the age-2 treatment. Feeding rate was positively correlated with both body size and with growth rate in the rearing tanks at both rations. The age-2 steelhead were more likely to be detected at a single feeding location than were the age-1 fish, which were more frequently detected at both locations. The more competitive age-1 fish were more frequently detected at both feeding locations, and the less competitive fish were detected at only one location. Relative competitive ability was not associated with feeding location in the age-2 treatment, nor was it associated with growth rate in either treatment. Age-2 steelhead that were detected only at the bottom or at both feeding locations had higher growth rates than those that were detected only at the surface, suggesting greater scramble competition near the point of food delivery under the restricted ration. Feeding location was not associated with growth rate in the age-1 treatment. The results indicate that ration directly affects feeding rate and that ration and relative competitive ability interact to influence feeding rate. Because feeding rate was positively associated with growth rate at both ration levels in the hatchery tanks, behaviors that are associated with feeding motivation may be an important trait under selection in salmonid hatcheries.

Tatara, Christopher P.↗
DOE OSTI2021

Moiréless correlations in ABCA graphene

Atomically thin van der Waals materials stacked with an interlayer twist have proven to be an excellent platform toward achieving gate-tunable correlated phenomena linked to the formation of flat electronic bands. In this work we demonstrate the formation of emergent correlated phases in multilayer rhombohedral graphene––a simple material that also exhibits a flat electronic band edge but without the need of having a moiré superlattice induced by twisted van der Waals layers. We show that two layers of bilayer graphene that are twisted by an arbitrary tiny angle host large (micrometer-scale) regions of uniform rhombohedral four-layer (ABCA) graphene that can be independently studied. Scanning tunneling spectroscopy reveals that ABCA graphene hosts an unprecedentedly sharp van Hove singularity of 3–5-meV half-width. We demonstrate that when this van Hove singularity straddles the Fermi level, a correlated many-body gap emerges with peak-to-peak value of 9.5 meV at charge neutrality. Mean-field theoretical calculations for model with short-ranged interactions indicate that two primary candidates for the appearance of this broken symmetry state are a charge-transfer excitonic insulator and a ferrimagnet. Finally, we show that ABCA graphene hosts surface topological helical edge states at natural interfaces with ABAB graphene which can be turned on and off with gate voltage, implying that small-angle twisted double-bilayer graphene is an ideal programmable topological quantum material.

42 ENGINEERING↗