Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “band structure methods”

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

Mirror symmetry breaking and lateral stacking shifts in twisted trilayer graphene

Here, we construct a continuum model of twisted trilayer graphene using ab initio density-functional-theory calculations, and we apply it to address twisted trilayer electronic structure. Our model accounts for moiré variation in site energies, hopping between outside layers, and hopping within layers. We focus on the role of a mirror symmetry present in ABA graphene trilayers with a middle layer twist. The mirror symmetry is lost intentionally when a displacement field is applied between layers, and unintentionally when the top layer is shifted laterally relative to the bottom layer. We calculate two band-structure characteristics that are directly relevant to transport measurements, namely the Drude weight and the weak-field Hall conductivity. We relate them via the Hall density to assess the influence of the accidental lateral stacking shifts currently present in all experimental devices on electronic properties. Finally, we comment on the role of the possible importance of accidental lateral stacking shifts for superconductivity in twisted trilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Twisted symmetric trilayer graphene. II. Projected Hartree-Fock study

The Hamiltonian of the magic-angle twisted symmetric trilayer graphene (TSTG) can be decomposed into a twisted-bilayer-graphene- (TBG-) like flat band Hamiltonian and a high-velocity Dirac fermion Hamiltonian. We use Hartree-Fock mean field approach to study the projected Coulomb interacting Hamiltonian of TSTG developed in Calugaru et al. [Phys. Rev. B 103, 195411 (2021)] at integer fillings ν = -3, -2, -1, and 0 measured from charge neutrality. We study the phase diagram with w 0 /w 1 , the ratio of AA and AB interlayer hoppings, and the displacement field, which introduces an interlayer potential U and hybridizes the TBG-like bands with the Dirac bands. At small U, we find the ground states at all fillings ν are in the same phases as the tensor products of a Dirac semimetal with the filling ν TBG insulator ground states, which are spin-valley polarized at ν = -3, and fully (partially) intervalley coherent at ν = -2, 0(ν = -1) in the flat bands. An exception is ν = -3 with w 0 /w 1 ≳ 0.7 , which possibly becomes a metal with competing orders at small U due to charge transfers between the Dirac and flat bands. At strong U where the bandwidths exceed interactions, all the fillings ν enter a metal phase with small or zero valley polarization and intervalley coherence. Lastly, at intermediate U, semimetal or insulator phases with zero intervalley coherence may arise for ν = -2, -1, 0. Finally, our results provide a simple picture for the electron interactions in TSTG systems, and reveal the connection between the TSTG and TBG ground states.

2-dimensional systems↗
DOE OSTI2024

Adjoint Waveform Tomography for Crustal and Upper Mantle Structure of the Middle East and Southwest Asia for Improved Waveform Simulations Using Openly Available Broadband Data

Here we present a new model of radially anisotropic seismic wavespeeds for the crust and upper mantle of a broad region of the Middle East and Southwest Asia (MESWA) derived from adjoint waveform tomography. The new model enables fully 3D simulations of complete three-component waveforms and provides improved fits that were not possible with previous models. We inverted over 32,000 waveforms from 192 earthquakes recorded by over 1000 openly available broadband seismic stations from permanent and temporary networks in the region with highly uneven coverage. Inversion iterations proceeded from the period band 50–100 s in six stages and 54 total iterations reducing the minimum period to 30 s. Our final model, MESWA, improves waveform fits compared to the starting and other models for both the data used in the inversion and an independent validation set of 66 events. Restitution tests indicate that the model resolves features in the central part of the model to depths of about 150 km. The new model reveals tectonic features imaged by other studies and methods but in a new holistic model of anisotropic shear and compressional wavespeeds (V S and V P , respectively) covering a larger domain with smaller scale length and amplified features. Examples include low crustal V S in the Tethyan belt and low mantle VS following divergent (Gulf of Aden, Red Sea) and transform (Dead Sea fault) margins of the Arabian plate. Low V S is imaged below Cenozoic volcanic centers of the Mecca–Madina–Nafud Line, Arabian Peninsula, and the Türkiye–Iran border region. Elevated V S tracks Makran subduction under southeast Iran with near vertical dip. MESWA could be used as a starting model for further improvements, say, using waveforms from in-country seismic networks that are not currently openly available and/or smaller-scale studies targeting a shorter period. The model could be used to improve earthquake hazard studies and nuclear explosion monitoring.

58 GEOSCIENCES↗
DOE OSTI2021

Spin Reorientation in Antiferromagnetic Layered FePt 5 P

FePt 5 P, a substitutional variant of the anti-CeCoIn 5 structure type in the space group P4/mmm, was synthesized by a high-temperature solid-state method and structurally characterized by X-ray diffraction. FePt 5 P contains layers of FePt 12 clusters formed by magnetically active Fe and heavy Pt with strong spin-orbit coupling (SOC); the layers are separated by P atoms. The various Fe–Pt distances in FePt 12 clusters generate complex magnetic orders in FePt 5 P. According to temperature-dependent magnetic and specific heat measurements, FePt 5 P shows a stripe-type antiferromagnetic order at T N ≈ 90 K, which is also confirmed by resistivity measurements. Furthermore, a spin reorientation occurs at ~74 and ~68 K in and out of the ab plane based on the specific heat measurements. The temperature-dependent neutron powder diffraction patterns demonstrate the antiferromagnetic order in FePt 5 P, and the spins orientate up to 58.4° with respect to the c axis at 10 K. First-principles calculations of FePt 5 P show the band splitting at the Fermi level by strong SOC and the s–d hybridization between P and Fe/Pt electrons enhances the structural stability and affects the magnetic ordering.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ab initio simulations of metal contacts for graphene-based devices

The precise atomic structure of a metal contact significantly affects the performance of nanoscale electronic devices. We use an accurate, DFT-based non-equilibrium Green’s function method to evaluate various metal contacts with graphene or graphene nanoribbons. Furthermore, for surface metal contacts not chemically bound to graphene, Ti contacts have lower resistance than those of Au, Ca, Ir, Pt, and Sr. However, as an edge contact, Ti has larger resistance than Au. Bridging O atoms at Ti and Au edge contacts lowers the transmission by over 30%.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Physical properties and electronic structure of single-crystal KCo 2 As 2

We present a method for producing high quality KCo 2 As 2 crystals, stable in air and suitable for a variety of measurements. X-ray diffraction, magnetic susceptibility, electrical transport and heat capacity measurements confirm the high quality and an absence of long range magnetic order down to at least 2 K. Residual resistivity values approaching 0.25 µΩ cm are representative of the high quality and low impurity scattering, and a Sommerfeld coefficient γ = 7.3 mJ/mol K 2 signifies weaker correlations than the Fe-based counterparts. Together with Hall effect measurements, angle-resolved photoemission experiments reveal a Fermi surface consisting of electron pockets at the center and corner of the Brillouin zone, in line with theoretical predictions and in contrast to the mixed carrier types of other pnictides with the ThCr 2 Si 2 structure. A large, linear magnetoresistance of 200% at 14 T, together with an observed linear and hyperbolic, rather than parabolic, band dispersions are unusual characteristics of this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fully self-consistent finite-temperature $GW$ in Gaussian Bloch orbitals for solids

In this work, we present algorithmic and implementation details for the fully self-consistent finite-temperature $GW$ method in Gaussian Bloch orbitals for solids. Our implementation is based on the finite-temperature Green's function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of the self-energy are explicitly evaluated. The method is tested by evaluating the band gaps of selected semiconductors and insulators. We show agreement with other, differently formulated, finite-temperature sc ⁢$GW$ implementations when finite-size corrections and basis-set errors are taken into account. By migrating computationally intensive calculations to graphics processing units, we obtain scalable results on large supercomputers with nearly optimal performance. Our work demonstrates the applicability of Gaussian orbital based sc⁢ $GW$ for ab initio correlated material simulations and provides a sound starting point for embedding methods built on top of $GW$.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Transuranium Sulfide via the Boron Chalcogen Mixture Method and Reversible Water Uptake in the NaCu T S 3 Family

The behavior of 5f electrons in soft ligand environments makes actinides, and especially transuranium chalcogenides, an intriguing class of materials for fundamental studies. Due to the affinity of actinides for oxygen, however, it is a challenge to synthesize actinide chalcogenides using non-metallic reagents. Using the Boron Chalcogen Mixture (BCM) method, we achieved the synthesis of the transuranium sulfide NaCuNpS 3 starting from the oxide reagent, NpO 2 . Via the same synthetic route, the isostructural composition of NaCuUS 3 was synthesized and the material contrasted with NaCuNpS 3 . Single crystals of the U-analog, NaCuUS 3 , were found to undergo an unexpected reversible hydration process to form NaCuUS 3 ·xH 2 O (x ≈ 1.5). Here, a large combination of techniques was used to fully characterize the structure, hydration process, and electronic structures, specifically a combination of single crystal, powder, high temperature powder X-ray diffraction, extended X-ray ab-sorption fine structure, infrared, and inductively coupled plasma spectroscopies, thermogravimetric analysis, and density functional theory calculations. The outcome of these analyses enabled us to determine the composition of NaCuUS 3 ·xH 2 O and obtain a structural model that demonstrated the retention of the local structure within the [CuUS 3 ] – layers throughout the hydration-dehydration process. Band structure, density of states, and Bader charge calculations for NaCuUS 3 , NaCu-US 3 ·xH 2 O, and NaCuNpS 3 along with X-ray absorption near edge structure, UV-vis-NIR and work function measurements on ACuUS 3 (A = Na, K, and Rb) and NaCuUS 3 ·xH 2 O samples were carried out to demonstrate that electronic properties arise from the [CuTS 3 ] – layers and show surprisingly little dependence on the interlayer distance.

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

A far‐red‐emitting (Gd,Y) 3 (Ga,Al) 5 O 12 :Mn 2+ ceramic phosphor with enhanced thermal stability for plant cultivation

Abstract As for plants, far‐red (FR) light with wavelength from 700 nm to 740 nm is critical for processes of photosynthesis and photomorphogenesis. Light‐controlled development depends on light to control cell differentiation, structural and functional changes, and finally converge into the formation of tissues and organs. Phosphor converted FR emission under LED excitation is a cost‐effective and high‐efficient way to provide artificial FR light source. With the aim to develop an efficient FR phosphor that can promote the plant growth, a series of gadolinium yttrium gallium garnet (GYGAG) transparent ceramic phosphors co‐doped with Mn 2+ and Si 4+ have been fabricated via chemical co‐precipitation method, followed sintered in O 2 and hot isostatic pressing in this work. Under UV excitation, the phosphor exhibited two bright and broadband red emission spectra due to Mn 2+ : 4 T 1 → 6 A 1 spin‐forbidden transition, and one of which located in the right FR region. And then, Ce 3+ ions were co‐doped as the activator to enhance the absorption at blue light region and the emission of Mn 2+ . It turns out that the emission band of GYGAG transparent ceramic phosphors matches well with the absorption band of phytochrome P FR , which means they are promising to be applied in plant cultivation light‐emitting diodes (LEDs) for modulating plant growth. Besides, the thermal stability of this material was investigated systematically, and an energy transferring model involves defects was also proposed to explain the phenomenon of abnormal temperature quenching.

Qiu, Xuting↗
DOE OSTI2020

Importance of van der Waals interactions for ab initio studies of topological insulators

We investigate the lattice and electronic structures of the bulk and surface of the prototypical layered topological insulators Bi 2 Se 3 and Bi 2 Te 3 using ab initio density functional methods, and systematically compare the results of different methods of including van der Waals (vdW) interactions. We show that the methods utilizing semi-empirical energy corrections yield accurate descriptions of these materials, with the most precise results obtained by properly accounting for the long-range tail of the vdW interactions. The bulk lattice constants, distances between quintuple layers and the Dirac velocity of the topological surface states (TSS) are all in excellent agreement with experiment. In Bi 2 Te 3 , hexagonal warping of the energy dispersion leads to complex spin textures of the TSS at moderate energies, while in Bi 2 Se 3 these states remain almost perfectly helical away from the Dirac point, showing appreciable signs of hexagonal warping at much higher energies, above the minimum of the bulk conduction band. Furthermore, our results establish a framework for unified and systematic self-consistent first principles calculations of topological insulators in bulk, slab and interface geometries, and provides the necessary first step toward ab initio modeling of topological heterostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

DFT Plus DMFT Study of Spin-Charge-Lattice Coupling in Covalnt LaCoO{siub 3}.

We study energetics and the nature of both homogeneous and mixed spin (MS) states in LaCoO3 incorporating structural changes of the crystal volume expansion and the Co-O bond disproportionation (BD) during the spin-state transition using the density functional theory plus dynamical mean field theory (DFT + DMFT) method. DFT + DMFT predicts that energetics of both excited spin states are almost the same while DFT + U calculations of the same structures energetically favor the MS states and produce various metastable solutions whose energetics depend sensitively on final spin states. Within DFT + DMFT, the homogeneous spin state in the expanded crystal volume shows the multiconfigurational nature with non-negligible occupancy probabilities of both high spin (HS) and low spin (LS) states along with d(6) and d(7) charge configurations indicating the dynamically fluctuating nature of spin and charge states due to the Co-O covalency. The nature of the MS state under the BD structure reveals that Co sites with the long Co-O bonds develop a Mott insulating state and favor HS with a d(6) configuration, while more covalent Co sites with the short Co-O bonds occupy more LS states with a d(7) configuration and behave as a band insulator, as a result, charge ordering is induced in the BD structure from the spin-state ordering. We also find that both energetics and electronic structure sensitively depend on the Co-O covalency effect, which can be tuned by changing the double counting potential and the resulting d occupancy (N-d), and N-d close to 6.7 is consistent with the nature of the spin-state transition. Our results show that structural changes during the spin-state transition can play an important role in understanding energetics and electronic structure of LaCoO3.

Park, Hyowon↗
DOE OSTI2024

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Cadmium Selective Etching in CdTe Solar Cells Produces Detrimental Narrow-Gap Te in Grain Boundaries

Recent advances in design and processing technology have made possible commercialization of polycrystalline (px)-CdTe as a photovoltaic absorber. Grain boundaries (GBs) are the most prominent structural defects in these devices and undergo significant changes during device fabrication. However, the effects of device fabrication processes on these GBs are not entirely understood. Prevailing models of GBs in thin-film photovoltaics consider individual GBs to have homogeneous properties in their area. Here, using an aberration-corrected scanning transmission electron microscope (STEM)-based low-loss and core-loss electron energy-loss spectroscopy (EELS), we show that back-surface etching of CdTe leads to inhomogeneity within individual grain boundaries. We observe that etching the back surface leads to the conversion of a region of GBs from CdTe to an elemental Te, which has an only 0.33 eV band gap, as deep as 1 μm from the back surface. The presence of elemental Te in GBs this deep into the absorber layer will increase recombination in the absorber layer and limit the extractable open-circuit voltage, thus reducing device efficiency. However, additive methods for back contact formation such as deposition of Te, ZnTe, or other materials preserve the CdTe stoichiometry of the GBs. Thus, especially for the next generations of CdTe-based cells having longer minority carrier diffusion length and/or thinner absorber layers, additive back contacting methods are superior.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Metastable piezoelectric group-IV monochalcogenide monolayers with a buckled honeycomb structure

Multiple two-dimensional materials are being naïvely termed stable on the grounds of displaying phonon dispersions with no negative frequencies and of not collapsing on molecular dynamics calculations at fixed volume. But, if these phases do not possess the smallest possible structural energy, how does one understand and establish their actual meta stability? To answer this question, twelve two-dimensional group-IV monochalcogenide monolayers (SiS, SiSe, SiTe, GeS, GeSe, GeTe, SnS, SnSe, SnTe, PbS, PbSe, and PbTe) with a buckled honeycomb atomistic structure—belonging to symmetry group P3m1—displaying an out-of-plane intrinsic electric polarization are shown to be metastable by three independent methods. First, we uncover a coordination-preserving structural transformation from the low-buckled honeycomb structure onto the lower-energy Pnm2 1 (or Pmmn for PbS, PbSe, and PbTe) phase to estimate energy barriers E B that must be overcome during such structural transformation. Using the curvature of the local minima and E B as inputs to Kramers escape formula, large escape times are found, implying the structural metastability of the buckled honeycomb phase (with the exception of PbS and PbSe, these phases display escape times ranging from 700 years to multiple times the age of the universe and can be considered “stable” for practical purposes in that relative sense). The second demonstration is provided by phonon dispersion relations that include the effect of long-range Coulomb forces and display no negative vibrational modes. The third and final demonstration of structural metastability is furnished by room-temperature ab initio molecular dynamics for selected compounds. Here, the magnitude of the electronic band gap evolves with chemical composition. Different from other binary two-dimensional compounds such as transition metal dichalcogenide monolayers and hexagonal boron nitride monolayers which only develop an in-plane piezoelectric response, the twelve group-IV monochalcogenide monolayers with a buckled honeycomb structure also display out-of-plane piezoelectric properties.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Point defects in Ga 2 O 3

In the field of high-power electronics, gallium oxide (Ga 2 O 3 ) is attracting attention due to its wide band gap and ability to be doped n-type. Point defects, including vacancies, impurities, and dopants, play important roles in optimizing device performance. This tutorial discusses the fundamental properties of point defects in monoclinic β-Ga 2 O 3 and the methods employed to study them. Oxygen vacancies are deep donors that do not cause n-type conductivity but may compensate acceptors. Gallium vacancies are deep acceptors that can be partially passivated by hydrogen. Substitutional magnesium is a promising acceptor that produces semi-insulating material and also forms a complex with hydrogen. Calcium and iron also have deep acceptor levels. Iridium deep donors are introduced into crystals grown from a melt in an Ir crucible. Here, other defects are introduced by irradiation with energetic particles such as neutrons or protons. In addition to altering the electronic properties, defects give rise to UV/visible emission bands in photoluminescence and cathodoluminescence spectra.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultradiffuse Galaxies (UDGs) with Hyper Suprime-Cam. I. Revised Catalog of Coma Cluster UDGs*

This is the first in a series of papers on the properties of ultradiffuse galaxies (UDGs) in clusters of galaxies. We present an updated catalog of UDGs in the Coma Cluster using g- and r-band images obtained with Hyper Suprime-Cam (HSC) of the Subaru telescope. We develop a method to find UDGs even in the presence of contaminating objects, such as halos and background galaxies. This study expands upon our previous works that covered about half the area of the Coma Cluster. The HSC observations covered the whole Coma Cluster up to the virial radius and beyond (an area twice as large as the previous studies) and doubled the numbers of UDGs (r eff , r ≥ 1.5 kpc) and sub-UDGs (1.0 ≤ r eff , r < 1.5 kpc) to 774 and 729, respectively. The new UDGs show internal properties consistent with those of previous studies (e.g., a Sérsic index of approximately 1), and are distributed across the cluster, with a concentration around the cluster center. The whole cluster coverage clearly revealed an excess of their distribution toward the east to southwest direction along the cluster center, where Coma connects to other large-scale structure, and where a known substructure exists (the NGC 4839 subgroup). The alignment of the UDG distribution along the large-scale structure around Coma supports the interpretation that most of them lie at the distance of the Coma Cluster and the NGC 4839 subgroup.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

α-Quartz Phase Stabilization, Surface Texturing, and Tunable Optical Properties of Nanocrystalline GeO 2 Films Made by Pulsed-Laser Deposition: Implications for Optical and Optoelectronic Applications

Germanium oxide (GeO 2 ) has great potential in multifunctional devices and next-generation power electronics due to its high thermal conductivity and ambipolar doping capability. However, the complexity of synthesizing the desirable polymorph with a controlled phase, surface/interface quality, microstructure, and functional properties is the main barrier to GeO 2 utilization in advanced applications. Here, in this regard, we present a method to realize the hexagonal (h) or a-quartz type GeO 2 with nano-textured surface morphology on sapphire substrates using a hybrid synthesis strategy that comprises pulsed laser deposition (PLD) and post-deposition thermal annealing. We performed a comprehensive study to investigate the effect of annealing temperature, which was varied in a wide range (600-1100 °C), on the crystal structure, phase, surface morphology, chemical stoichiometry, defect states, and optical properties of PLD-grown GeO 2 films. As-deposited GeO 2 films at 500 °C were amorphous. Upon annealing, the GeO 2 films induced an amorphous-to-crystalline phase transformation; GeO 2 films annealed at higher annealing temperatures (≥900 °C) stabilized in the hexagonal phase and demonstrated excellent crystal quality and chemical stability. Thermally activated growth process showed increased average crystallite size, which was varied in the range of 20-130 (±2) nm, while the surface roughness followed a similar trend. The spectral transmittance and band gap also increased with increasing annealing temperature. The resulting h-GeO 2 films, particularly those obtained at annealing temperatures in the 900-1100 °C range, had a higher band gap of 6.2-6.3 eV and displayed excellent optical transmittance in the visible region. Moreover, the absence of extended valence band maxima and reduced optical defect density support the quality improvement upon annealing. When considering phase-pure bulk and nanostructured GeO 2 as a possible candidate for ultra wide band gap semiconductors in cutting-edge technological applications, the results of the current work can be beneficial to realize high structural and optical quality a-quartz structured GeO 2 films.

GeO2↗
DOE OSTI2023

Energy Exascale Earth System Model v2.1.0

First release of version 2.1 of the Energy Exascale Earth System Model. [ATM] The atmosphere component remains EAM. There are no major changes in the default configuration since 2.0. New features include: A semi-lagrangian tracer transport for theta-l dycore, a new algorithm for finding the tropopause, new RRM mesh configurations. Add and update SSP370 and SSP585 cases. Restore the FIDEAL case. [LAND] The land component is ELM. There are no major changes in the default configuration since 2.0. Several option features have been added including: implementation of topography-based subgrid structure (topounits) and accompanying parameterizations and atmospheric forcing downscaling methods; a new plant hydraulics scheme; two-way land-river hydrological coupling through the infiltration of floodplain water; an implementation of perennial crops; updates to the SNICAR-AD snow radiative transfer model; and implementation of soil erosion and sediment yield in ELM-Erosion. Each of these new changes is modular in design and can be turned on or off as the user specifies; they are currently being tested in different “BGC” configurations. [OCEAN] The ocean component remains MPAS-Ocean. Major change since version 2.0 include the addition of the Fox-Kemper et al. 2011 parameterization for submesocale eddies, a correction for barotropic thickness consistency that reduces divergence noise, and the addition of an ocean carbon conservation analysis member. [SEAICE] The sea-ice component remains MPAS-Seaice. Major changes since version 2.0 include: A correction to how shortwave parameters are interpolated in the snicar-ad 5-band radiation scheme, the addition of a sea ice carbon conservation analysis member, updates to the default sea ice biogeochemistry namelist parameters to be consistent with version 2.0 improvements to nitrogen cycling and a correction in the ice-ocean dissolved organic nitrogen coupling. [LAND ICE] The land-ice component remains MPAS-Albany-landIce (MALI). Major changes since 2.0 include an update to the MALI version and the Greenland mesh.[RIVER] The river model is MOSART. There are no major changes in the default configuration since 2.0. A major new optional feature is two-way river-ocean hydrological coupling between MOSART and MPAS-O. This change can be turned on or off as the user specifies, and is being tested in different configurations. [COUPLER] The coupler remains cpl7/MCT. Major changes since version 2.0 include: Carbon budget calculated when heat/water budgets active. Fix a bug in land-atm fluxes for tri-grid configurations. [OTHER] a small bug in the zenith angle calculation was fixed in the data models.

ECP↗