Search NASA⌕ Search

SEARCH · Search NASA

Results for “density functional theory modeling”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Conserving asphalt resources: Rethinking rejuvenator performance evaluation through peptizing efficiency

Timely rejuvenation and restoration of asphalt are essential conservation practices that help preserve and extend the service life of roads, bridges, driveways, and parking lots. The performance of asphalt rejuvenators is often assessed based on their diffusion rates and softening power, yet these metrics alone fail to capture true rejuvenation potential. Here, this study integrates density functional theory (DFT) modeling with experimental analysis to show that rejuvenation effectiveness is primarily governed by molecular interactions with oxidized asphaltene nanoaggregates. DFT results revealed that amide- and unsaturated-chain compounds, such as hexadecanamide, oleic acid, and 9,17-octadecadienal, act cooperatively to disrupt π–π stacking and exfoliate asphaltene layers, reducing binding strength and enhancing dispersion. Unlike bulky, rigid molecules that remain trapped, these components remain mobile and repeatedly interact with multiple aggregation sites. Experimental validation using rheometry and FTIR confirmed that such molecularly compatible rejuvenators restore the binder flexibility and polydispersity, even when diffusion is relatively slow. Building on this mechanistic foundation, six rejuvenators (A2, A5, A7, A8, A9, and A10) were evaluated using a three-metric performance framework encompassing cracking resistance (Glover–Rowe parameter), UV stability, and surface hydrophobicity. Among the six rejuvenators evaluated, A9 exhibited the highest overall performance, reducing the Glover–Rowe cracking parameter by 85% (from 351 kPa to 53 kPa), demonstrating the greatest resistance to UV-induced aging with a stability index of 4.17 h·kPa⁻¹, and increasing surface hydrophobicity to a contact angle of 103.6°. Its superior performance is primarily attributed to its amide- and unsaturated-chain components, which act cooperatively to disrupt π–π stacking interactions and exfoliate asphaltene layers, thereby promoting molecular deagglomeration and enhancing long-term durability. These results shift the criteria for selecting rejuvenators: effective candidates must pair the electronic capability to unlock aged asphaltenes with sufficient structural stability to resist secondary aging and restore hydrophobicity. Restoring hydrophobicity is critical, as aging reduces the asphalt’s water repellency and increases water diffusion, which in turn accelerates moisture-related damage; consequently, a high-performing rejuvenator must effectively restore the binder’s hydrophobic characteristic. Collectively, these findings provide a framework for the rational design of next-generation of bio-based rejuvenators that enhance pavement longevity and promote long-term sustainability.

Aging↗
DOE OSTI2024

Experimental and Computational Evaluation of Lipidomic In-Source Fragmentation as a Result of Postionization with Matrix-Assisted Laser Desorption/Ionization

Matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI) can provide spatially resolved molecular information about a sample. Recently, a postionization approach (MALDI-2) has been commercially integrated with MALDI-MSI, allowing for bettered sensitivity and consequent improved spatial resolution. While advantages of MALDI-2 have previously been established, we demonstrate here statistically increased in-source fragmentation (ISF) results from postionization with a commercial instrument. Via lipid standard analyses, known MALDI ISF pathways (e.g., loss of trimethylamine) were statistically increased in MALDI-2 compared to MALDI-1 (65–172% increase in fragmentation). Gas phase molecular modeling with density functional theory estimated that the most-weighted virtual orbitals to excite within lipids involve ester and phosphate bonds. Protonated lipid excitation energies are furthermore red-shifted compared to those of other adduct types [e.g., 254 nm for protonated PC(16:0/18:1)] and approach the MALDI-2 laser energy (266 nm). Analysis of rat brain homogenate detected statistically more positive-ion mode peaks with MALDI-2 (1090) than that with MALDI-1 (719), where Kernel density estimations showed that the majority of this enhancement occurs with low m/z ions (i.e., m/z 75–500). Taken together with the lipid standard data, these observations may indicate ISF due to postionization. Finally, while artifact contributions from matrix blanks were also noted, both experimental and computational data sets suggest that the overall extent of ISF is statistically increased in MALDI-2 compared to MALDI-1.

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

Tailoring Interfaces in Solid-State Batteries Using Interfacial Thermochemistry and Band Alignment

Solid-state lithium ion batteries have enhanced thermal stability compared to batteries using liquid electrolytes. Although their practical implementation is limited by undesired side reactions between the electrode and electrolyte, different modification strategies have been proposed to stabilize these reactive interfaces. These approaches have been primarily based upon bulk materials properties, however, which may not necessarily be representative of the chemistry at the solid/solid interface. In this work, first principles calculations and X-ray scattering experiments are used to elucidate molecular-level reactivity and develop design principles for tailored interfaces based on surface and interfacial properties. Because of its well-known instability, the interface between the Li metal anode and the lithium lanthanum titanate (LLTO) solid electrolyte is applied as a model system. Density functional theory calculations are used to describe bulk, surface, and interfacial thermochemistry of the Li/LLTO system, and interfacial reconstruction is probed using ab initio molecular dynamics. These simulations of the Li/LLTO interface demonstrate facile decomposition concomitant with reduction of Ti 4+ ions. Based on further insights from computational analysis of the surface band edge positions, La 2 O 3 is proposed as an interlayer coating and is shown to provide an energetic barrier for interfacial charge transfer and decomposition reactions from X-ray reflectivity measurements and theoretical calculations. The findings in this work suggest that design strategies based on surface and interfacial properties can be used to kinetically stabilize electrode/electrolyte interfaces in solid-state batteries.

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

Stabilization of Cu 2 O through Site-Selective Formation of a Co 1 Cu Hybrid Single-Atom Catalyst

Single-atom catalysts (SACs) consist of a low coverage of isolated metal atoms dispersed on a metal substrate, called single-atom alloys (SAAs), or alternatively single metal atoms coordinated to oxygen atoms on an oxide support. We present the synthesis of a new type of Co 1 Cu SAC centers on a Cu 2 O(111) support by means of a site-selective atomic layer deposition technique. Isolated metallic Co atoms selectively coordinate to the native oxygen vacancy sites (Cu sites) of the reconstructed Cu 2 O(111) surface, forming a Co 1 Cu SAA with no direct Co–O x bonds. The centers, here referred to as Co 1 Cu hybrid SACs, are found to stabilize the active Cu + sites of the low-cost Cu 2 O catalyst that otherwise is prone to deactivation under reaction conditions. The stability of the Cu 2 O(111) surface was investigated by synchrotron radiation-based ambient-pressure X-ray photoelectron spectroscopy under reducing CO environment. The structure and reduction reaction are modeled by density functional theory calculations, in good agreement with experimental results.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemical Oxidation in Garnet-Type Solid Electrolyte by Formation of Point Defects

All-solid-state batteries hold greater promise for improving safety and energy density over conventional battery technology employing organic liquid electrolytes. One of the required features of a Li + conducting solid electrolyte is electrochemical stability, attained thermodynamically or kinetically, within the targeted operating voltage and temperature ranges. Therefore, understanding of the oxidative or reductive degradation mechanism is important to allow the design of stable solid electrolyte materials. This work contributes to building an understanding of the oxidative degradation mechanism in lithium solid electrolytes at cell operating conditions. Here, we have focused on resolving the oxidative decomposition mechanism of Al-doped lithium garnet Li 6.28 Al 0.24 La 3 Zr 2 O 12 (LLZO) as a state-of-the-art inorganic ceramic electrolyte. By combining experimental and computational analyses, we show that oxidation of LLZO occurs by simultaneous loss of oxygen and lithium from the structure, resulting in substoichiometric LLZO, at a moderate temperature (80 °C) and a high electrode potential (4.3 V vs Li/Li + ). Based on X-ray absorption and diffraction analyses, we find that the zirconium coordination shells in LLZO contract while the crystal structure experiences positive chemical strain upon electrochemical oxidation. The results from ex situ structural characterization of both the local structure and crystal symmetry are supported by a substoichiometric LLZO with lithium and oxygen vacancies, modeled by density functional theory (DFT) calculations. These chemical and structural changes in LLZO suppress effective lithium-ion conductivity by an order of magnitude. Formation of lithium and oxygen vacancies in LLZO upon electrochemical oxidation is different from prior thermodynamic predictions of phase decomposition of LLZO. The difference here is that the experiments were conducted at near-room temperature, which can hinder the kinetics of phase separation, and thus, the resultant LLZO solid electrolyte is still single-phase but substoichiometric in Li and O. In conclusion, these findings contribute an important degradation mechanism of the electrolyte, relevant for practical operational conditions of solid-state batteries.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Luminescent Metal–Organic Framework for the Selective Detection of Aldehydes

The detection of toxic, hazardous chemical species is an important task because they pose serious risks to either the environment or human health. Luminescent metal–organic frameworks (LMOFs) as alternative sensors offer rapid and sensitive detection of chemical species. Interactions between chemical species and LMOFs result in changes in the photoluminescence (PL) profile of the LMOFs which can be readily detected using a simple fluorometer. Herein, we report the use of a robust, Zn-based LMOF, [Zn 5 (μ 3 -OH) 2 (adtb) 2 (H 2 O) 5 ·5 DMA] (Zn-adtb, LMOF-341), for the selective detection of benzaldehyde. Upon exposure to benzaldehyde, Zn-adtb experiences significant luminescent quenching, as characterized through PL experiments. Photoluminescent titration experiments reveal that LMOF-341 has a detection limit of 64 ppm and a Ksv value of 179 M –1 for benzaldehyde. Furthermore, we study the guest–host interactions that occur between LMOF-341 and benzaldehyde through in situ Fourier transform infrared and computational modeling employing density functional theory. The results show that benzaldehyde interacts more strongly with LMOF-341 compared to formaldehyde and propionaldehyde. Furthermore, our combined studies also reveal that the mechanism of luminescence quenching originates from an electron-transfer process.

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

Static and Dynamic Correlations in Water: Comparison of Classical Ab Initio Molecular Dynamics at Elevated Temperature with Path Integral Simulations at Ambient Temperature

It is a common practice in ab initio molecular dynamics (AIMD) simulations of water to use an elevated temperature to overcome the overstructuring and slow diffusion predicted by most current density functional theory (DFT) models. The simulation results obtained in this distinct thermodynamic state are then compared with experimental data at ambient temperature based on the rationale that a higher temperature effectively recovers nuclear quantum effects (NQEs) that are missing in the classical AIMD simulations. In this work, we systematically examine the foundation of this assumption for several DFT models as well as for the many-body MB-pol model. Here, we find for the cases studied that a higher temperature does not correctly mimic NQEs at room temperature, which is especially manifest in significantly different three-molecule correlations as well as hydrogen bond dynamics. In many of these cases, the effects of NQEs are the opposite of the effects of carrying out the simulations at an elevated temperature.

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

Complex Fluorine Chemical Potential Effects on the Shape and Compositional Heterogeneity of KTa 1– x Nb x O 3 Nanoparticles

When kinetic conditions dominate and dictate the growth in a nanoparticle synthesis, properties of the synthesis environment can have considerable effects on the properties of the products. Such effects were studied here, where the solution environment was changed via the addition of KF to the hydrothermal syntheses of KTa 1– x Nb x O 3 and KTaO 3 . One result demonstrated the straightforward cause-and-effect relationship between the solution and reaction kinetics: KF directly increased the solution stability of the Ta species and therefore decreased its reaction rate, resulting in a change in composition heterogeneity of Ta and Nb in the KTa 1– x Nb x O 3 particles. However, not all effects are so simple; changing the chemical potential of the solution with KF can also promote the formation of particles with anisotropic defect enhanced kinetic Wulff shapes instead of cuboidal shapes. The increased F chemical potential in the solution enabled the formation of planar defects in the bulk, which accelerated growth in-plane to form particles characterized by flat rectangular flake geometries. In this work, thermodynamic modeling with density functional theory calculations confirmed that sufficient KF concentrations can drive the formation of a defect phase K n +1 Ta n O 3 n F.

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

Inelastic Neutron Scattering Study of Magnetic Exchange Pathways in MnS

We report an investigation of the magnetic structure and magnetic exchange pathways in MnS via neutron scattering methods, aided by density functional theory (DFT) modeling. The material has been confirmed to undergo antiferromagnetic (AFM) ordering at 152 K, with the magnetic structure representing AFM stacking of ferromagnetic (FM) (111) planes of Mn magnetic moments. Correspondingly, the magnetic structure is described by a propagation vector k = (1/2, 1/2, 1/2), with the volume of the magnetic unit cell being 8 times larger than the volume of the nuclear unit cell. Analysis of inelastic neutron scattering (INS) data collected on a powder sample of MnS revealed that the next-nearest-neighbor magnetic exchange constant (J 2 ) exceeds the nearest-neighbor exchange constant (J 1 ) by more than 3 times, while in the case of MnO, which exhibits the same nuclear and magnetic structures as MnS, the J 2 /J 1 ratio was reported to be below 1.5. Although for MnO the signs of both J 1 and J 2 indicated AFM exchange interactions, machine-learning INS data analysis in combination with DFT calculations suggests that the INS data collected on MnS are best described with J 1 < 0 and J 2 > 0, corresponding to FM and AFM exchange couplings, respectively. To achieve a satisfactory fit to the experimentally observed data, the Hamiltonian used to model the INS spectra also included the next-next-nearest-neighbor magnetic exchange constant (J 3 ). The best-fit model has been obtained with the values of the exchange constants J 1 = -0.27 meV, J 2 = 1.05 meV, and J 3 = -0.19 meV.

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

Electronic Structure of Molecular Cobalt Catalysts for H 2 Production Revealed by Multifrequency EPR

Improving the design of hydrogen evolution reaction (HER) catalysts requires a detailed understanding of the catalyst's structure and environment. Here wwe have used multifrequency electron paramagnetic resonance (EPR) spectroscopy (X-band, Q-band, and D-band) in combination with density functional theory (DFT) modeling to elucidate the electronic structure of three Co(II) molecular catalysts for H 2 evolution to determine correlations of electronic structure parameters with catalytic activity and stability. The Co(II) coordination complexes investigated here are representative of two types of H 2 -evolving molecular catalysts: cobaloximes, formed from cobalt coordination to two glyoxime ligands, and cobalt poly(pyridyl) complexes, which commonly have either four or five pyridyl ligands coordinating a single Co(II) ion. These Co(II) complexes were studied in a variety of solvents with a range of polarities and stoichiometric amounts of potential ligands to the cobalt ion to assess the strength and nature of ligand- cobalt interactions. Using magnetic resonance parameters obtained by EPR spectroscopy as reference data, we reconstructed the electronic structure of the Co(II) complexes with the help of DFT modeling. For all three catalysts, we uncovered a linear correlation between the Co(II) A z value and the g x value that is dependent on the strength of the axial ligands. The correlations that we have established between the electronic structure and metal coordination environment provide insight into the structural and electronic factors underpinning the observed trends in HER catalysis activity.

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

Revealing Intermetallic Active Sites of PtNi Nanocatalysts for Reverse Water Gas Shift Reaction

Bimetallic nanoparticles have gained great interest due to their potential applications in catalysis. By combining Lattice Boltzmann Modeling (LBM), Density Functional Theory (DFT) calculations, and kinetic Monte Carlo (KMC) simulation, here we rationalize the superior performance of supported Pt-Ni alloy nanostructures over Pt and Ni alone during the reverse Water Gas Shift (rWGS) reaction, which was observed previously. The LBM simulation including the diffusion and flow of the reacting gas molecules within the supporting nanoporous channels in SBA-15 indicates that the simple combination of contributions from Ni and Pt alone is not enough to result in the observed promotion in rWGS activity and selectivity of the Pt-Ni alloy. Instead, the synergy between Pt and Ni via a unique intermetallic motif is responsible according to the DFT and KMC simulation. In conclusion, compared to the monometallic counterparts, the intermetallic-introduced ensemble effect at the low-coordinated Pt-Ni hybrid site of Pt-Ni alloy nanocatalysts enables the enhanced activation of CO 2 and promoted removal of CO/H 2 O as gas phase product.

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

Vibrational study of CO, O 2 , and H 2 Adsorbed on the CoCrFeNi (110) High Entropy Alloy Surface

The vibrational properties of CO, O 2 , and H 2 molecularly or dissociatively adsorbed on a CoCrFeNi(110) surface have been probed using highresolution energy loss spectroscopy (HREELS) and modeled using density functional theory (DFT) calculations. Large (~20 mm3) single-crystal, quaternary face-centered cubic CoCrFeNi was synthesized via a modified Czochralski technique. We show strong evidence that CO adsorbs primarily on bridge and on-top sites in compositionally varied local environments, which reflect the random, multielemental surface composition inherent in a high entropy alloy. A variation of adsorption sites is also found with oxygen, which exhibits two broad groups of modes. Comparison to previous photoemission and theoretical studies suggests that the higher energy modes consist primarily of local CrO x species, while the lower energy modes are due to oxygen atoms adsorbed on other metal sites. Unlike CO and O 2 , HREELS upon H 2 adsorption shows only two much narrower modes and is consistent with atomic adsorption on 3-fold hollow sites. The hypothesized adsorption sites for all three species are directly corroborated by our DFT calculations.

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

Ship-in-a-Bottle Synthesis of High Concentration of N2 Molecules in a Cage-Structured Electride

We report the formation of neutral nitrogen molecules in the cages of [Ca12Al14O32]2+ (C12A7) framework compensated by extra-framework anions. NH3 treatment of C12A7 electride (C12A7:e–) at 800 °C leads to the formation of N2 and NH2– species in the C12A7 cages. N2 and NHx species in the cages are identified using the Raman spectroscopy of 14NH3 and 15NH3-treated C12A7:e–. The concentration of H and N in the C12A7 cages after NH3 treatment is ~1021 cm–3. We propose a two-step mechanism, supported by density functional theory (DFT) modeling, of N2 incorporation into the C12A7 cages: incorporation of NH2– formed from decomposition of NH3 at C12A7:e– surface followed by the NH2- species reacting to form N2 molecules. Encapsulation of neutral molecules, as opposed to negatively charged species reported in C12A7 previously, offers new opportunities for trapping and storing gaseous substances in nano-porous materials.

Nakao, Takuya↗
DOE OSTI2024

Coupled Experimental–Theoretical Characterization of a Carbon Electrode in Vanadium Redox Flow Batteries using X-ray Absorption Spectroscopy

Vanadium redox flow batteries (VRFBs) have emerged as promising solutions for stationary grid energy storage due to their high efficiency, scalability, safety, near room-temperature operation conditions, and the ability to independently size power and energy capacities. The performance of VRFBs heavily relies on the redox couple reactions of V 2+ /V 3+ and VO 2+ /VO 2 + on carbon electrodes. Therefore, a thorough understanding of the surface functionality of carbon electrodes and their propensity for degradation during electrochemical cycles is crucial for designing VRFBs with extended lifespans. In this study, we present a coupled experimental–theoretical approach based on carbon K edge X-ray absorption spectroscopy (XAS) to characterize carbon electrodes prepared under different conditions and identify relevant functional groups that contribute to unique spectroscopic features. Atomic models were created to represent functional groups, such as hydroxyl, carboxyl, methyl, and aldehyde, bonded to carbon atoms in either sp 2 or sp 3 environments. The interactions between functionalized carbon and various solvated vanadium complexes were modeled using density functional theory. A library of carbon K-edge XAS spectra was generated for distinct carbon atoms in different functional groups, both before and after interacting with solvated vanadium complexes. Here we demonstrate how these simulated spectra can be used to deconvolve ex situ experimental spectra measured from carbon electrodes and to track changes in the electrode composition following immersion in different electrolytes or extended cycling within a functional VRFB. By doing so, we identify the active species present on the carbon electrodes, which play a crucial role in determining their electrochemical performance.

25 ENERGY STORAGE↗
DOE OSTI2020

Enhancement of van der Waals Interlayer Coupling through Polar Janus MoSSe

Interlayer coupling plays essential roles in the quantum transport, polaritonic, and electrochemical properties of stacked van der Waals (vdW) materials. In this work, we report the unconventional interlayer coupling in vdW heterostructures (HSs) by utilizing an emerging 2D material, Janus transition metal dichalcogenides (TMDs). In contrast to conventional TMDs, monolayer Janus TMDs have two different chalcogen layers sandwiching the transition metal and thus exhibit broken mirror symmetry and an intrinsic vertical dipole moment. Such a broken symmetry is found to strongly enhance the vdW interlayer coupling by as much as 13.2% when forming MoSSe/MoS 2 HS as compared to the pristine MoS 2 counterparts. Our noncontact ultralow-frequency Raman probe, linear chain model, and density functional theory calculations confirm the enhancement and reveal the origins as charge redistribution in Janus MoSSe and reduced interlayer distance. Furthermore, our results uncover the potential of tuning interlayer coupling strength through Janus heterostacking.

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

Efficient Intermolecular Charge Transport in π-Stacked Pyridinium Dimers Using Cucurbit[8]uril Supramolecular Complexes

Intermolecular charge transport through π-conjugated molecules plays an essential role in biochemical redox processes and energy storage applications. In this work, we observe highly efficient intermolecular charge transport upon dimerization of pyridinium molecules in the cavity of a synthetic host (cucurbit[8]uril, CB[8]). Stable, homoternary complexes are formed between pyridinium molecules and CB[8] with high binding affinity, resulting in an offset stacked geometry of two pyridiniums inside the host cavity. The charge transport properties of free and dimerized pyridiniums are characterized using a scanning tunneling microscope-break junction (STM-BJ) technique. Our results show that π-stacked pyridinium dimers exhibit comparable molecular conductance to isolated, single pyridinium molecules, despite a longer transport pathway and a switch from intra- to intermolecular charge transport. Control experiments using a CB[8] homologue (cucurbit[7]uril, CB[7]) show that the synthetic host primarily serves to facilitate dimer formation and plays a minimal role on molecular conductance. Molecular modeling using density functional theory (DFT) reveals that pyridinium molecules are planarized upon dimerization inside the host cavity, which facilitates charge transport. In addition, the π-stacked pyridinium dimers possess large intermolecular LUMO–LUMO couplings, leading to enhanced intermolecular charge transport. Altogether, this work demonstrates that supramolecular assembly can be used to control intermolecular charge transport in π-stacked molecules.

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

Magic in twisted transition metal dichalcogenide bilayers

The long-wavelength moiré superlattices in twisted 2D structures have emerged as a highly tunable platform for strongly correlated electron physics. We study the moiré bands in twisted transition metal dichalcogenide homobilayers, focusing on WSe 2 , at small twist angles using a combination of first principles density functional theory, continuum modeling, and Hartree-Fock approximation. We reveal the rich physics at small twist angles θ < 4°, and identify a particular magic angle at which the top valence moiré band achieves almost perfect flatness. In the vicinity of this magic angle, we predict the realization of a generalized Kane-Mele model with a topological flat band, interaction-driven Haldane insulator, and Mott insulators at the filling of one hole per moiré unit cell. The combination of flat dispersion and uniformity of Berry curvature near the magic angle holds promise for realizing fractional quantum anomalous Hall effect at fractional filling. We also identify twist angles favorable for quantum spin Hall insulators and interaction-induced quantum anomalous Hall insulators at other integer fillings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗