SEARCH · Engineering Papers

Results for “Derived Matrices”

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.

SourceResearch connectionFrom yearThrough year

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

138 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Immobilization of cesium and iodine into Cs 3 Bi 2 I 9 perovskite-silica composites and core-shell waste forms with high waste loadings and chemical durability

Cs$_3$Bi$_2$I$_9$, a defect perovskite derivative, is a potential host phase to immobilize iodine and cesium with high waste loadings. In this work, two strategies were explored to form Cs$_3$Bi$_2$I$_9$-silica composites and a core-shell structure in order to improve chemical durability of waste form materials meanwhile maintaining high waste loadings. Cs$_3$Bi$_2$I$_9$ loadings as high as 70 wt.% were incorporated into a silica matrix to form silica-ceramic composites, and 20 wt.% Cs$_3$Bi$_2$I$_9$ was encapsulated into silica to form a core–shell structure by low temperature spark plasma sintering. Chemical durability of the composite and core-shell waste forms was evaluated by semi-dynamic leaching experiments, and Cs and I were incongruently released from waste form matrices. A BiOI alteration layer formed, acting as a passivation layer to reduce the release of radionuclides. The long-term iodine release rate was low (30 mg m$^{-2}$ day$^{-1}$) for the 70 wt.%Cs$_3$Bi$_2$I$_9$–silica composite leached in deionized water at 90 °C, which can be further reduced to 5 × 10$^{-3}$ mg m$^{-2}$ day$^{-1}$ for the 20 wt.% core-shell structure. This work highlights a robust way to immobilize the highly mobile radionuclides with high waste loadings through encapsulation into durable matrices and a surface passivating mechanism that can greatly reduce the elemental transport from waste form materials and significantly enhance their chemical durability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Selective Removal of Barium and Hardness Ions from Brackish Water with Chemically Enhanced Electrodialysis

The use of a nonconventional water resource for energy and industrial applications often requires extraction of low-level undesirable ions from matrices of benign dominant ions. In this study, the selective extraction of Ba 2+ and Mg 2+ from synthetic brine solutions was evaluated using strong acid cation-exchange membranes modified with the surface deposition of macrocyclic molecules including crown ethers and calixarenes. Compared to the bare membrane, vinylbenzo-18-crown-6 (VB18C6) and calix[4]arene-amended membranes showed increased selectivity for Ba 2+ and Mg 2+ with respect to the dominant ion (Na + ) by up to fourfolds, with the calix[4]arene-modified membrane achieving more selective separation than VB18C6. Optimal selectivity was achieved at a moderate-to-high current density (3.1–6.3 mA/cm 2 ), which was attributed to the alleviation of transport limitation in the boundary layer by the surface modification. Amendment of a calix[4]arene derivative with a crown-6 “boot strap” and two carboxylic groups resulted in reduced selectivity due to the formation of strong complexes with divalent ions. Furthermore, these results show that surface deposition of ion sequestrants can be a versatile approach to improve the membrane’s selectivity; however, the performance is sensitive to feed water salinity, current loading, and the macrocycle-ion chemistry, where there is a trade-off between ion affinity and mobility.

36 MATERIALS SCIENCE↗

DOE OSTI2023

The engineered single guide RNA structure as a biomarker for gene-editing reagent exposure

Abstract CRISPR arrays and CRISPR-associated (Cas) proteins comprise a prevalent adaptive immune system in bacteria and archaea. These systems defend against exogenous parasitic mobile genetic elements. The adaption of single effector CRISPR-Cas systems has massively facilitated gene-editing due to the reprogrammable guide RNA. The guide RNA affords little priming space for conventional PCR-based nucleic acid tests without foreknowledge of the spacer sequence. Further impeding detection of gene-editor exposure, these systems are derived from human microflora and pathogens ( Staphylococcus pyogenes , Streptococcus aureus , etc.) that contaminate human patient samples. The single guide RNA—formed from the CRISPR RNA (crRNA) and transactivating RNA (tracrRNA)—harbors a variable tetraloop sequence between the two RNA segments, complicating PCR assays. Identical single effector Cas proteins are used for gene-editing and naturally by bacteria. Antibodies raised against these Cas proteins are unable to distinguish CRISPR-Cas gene-editors from bacterial contaminant. To overcome the high potential for false positives, we have developed a DNA displacement assay to specifically detect gene-editors. We leveraged the single guide RNA structure as an engineered moiety for gene-editor exposure that does not cross-react with bacterial CRISPRs. Our assay has been validated for five common CRISPR systems and functions in complex sample matrices.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Extracting the asymptotic behavior of S-matrix elements from their phases

The asymptotic kinematic limits of S-matrices are dominated by large logarithms which, roughly speaking, fall into two categories: those which are controlled by a renormalization group (RG) scale, which we may think of as logs involving ratios of invariant mass scales, and those which are functions of ratios of rapidities, so called “rapidity logs”. It has been pointed out by Caron-Huot and Wilhlem [1] that RG anomalous dimension can be extracted from the phase of the S-matrix, which can greatly simplify calculations via unitarity methods. In this paper we generalize the results of [1] to show that the phase can be used to reconstruct rapidity anomalous dimensions, by performing a special type of complex boost. The methodology introduced allows one to calculate without the need for a rapidity regulator which can lead to significant simplifications. We demonstrate the use of this method to derive the rapidity anomalous dimensions in the Sudakov form factor and the two parton soft function at two loops order.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2026

Applications of Decellularized Plant Tissues in Regenerative Medicine and Tissue Engineering

The development of biomaterials capable of supporting complex tissue growth remains a central challenge in regenerative medicine and tissue engineering, particularly in replicating the structural, mechanical, and transport functions of native extracellular matrices. While decellularized animal tissues have demonstrated significant success as scaffolds for tissue engineering, they are still constrained by cost, immunogenicity, and ethical concerns. In recent years, decellularized plant tissues have emerged as a compelling alternative scaffold platform due to their inherent vascular architectures, ethical sourcing, tunable mechanical properties, cytocompatibility, and sustainability. This review summarizes current strategies for the decellularization of plant tissues, including chemical, enzymatic, and physical approaches, and discusses how these methods preserve plant cell wall structure while removing immunogenic components. Advances in surface loading and functionalization, including protein coatings, oxidation, nanoparticle incorporation, peptide conjugation, and bioactive molecule loading, have further enhanced cell adhesion, differentiation, biodegradability, and immunomodulation. Recent applications of decellularized plant scaffolds in cardiac, skeletal muscle, bone, nerve, and wound healing contexts are reviewed, highlighting proof-of-concept successes and remaining challenges. Beyond therapeutic applications, plant-derived scaffolds have also enabled physiologically relevant in vitro models for vascular biology, mechanotransduction, cancer, metabolic tissues, and drug response studies. Collectively, these advances position decellularized plant tissues as versatile, low-cost, and ethically favorable biomaterials with growing relevance for both regenerative medicine and tissue modeling.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Orbital dependent complications for close vs well-separated electrons in diradicals

We investigate two limits in open-shell diradical systems: O3, in which the interesting orbitals are in close proximity to one another, and (C21H13)2, where there is a significant spatial separation between the two orbitals. In accord with earlier calculations, we find that standard density-functional approximations do not predict the open-shell character for the former case but uniformly predict the open-shell character for the latter case. We trace the qualitatively incorrect behavior in O3 predicted by these standard density functional approximations to self-interaction error and use the Fermi–Löwdin-orbital-self-interaction-corrected formalism to determine accurate triplet, closed-shell singlet, and open-shell broken-spin-symmetry electronic configurations. Analysis of the resulting many-electron overlap matrices allows us to unambiguously show that the broken-spin-symmetry configurations do not participate in the representation of the Ms = 0 triplet states and allows us to reliably extract the singlet–triplet splitting in O3 by analyzing the energy as a function of Fermi-orbital-descriptor permutations. The results of these analyses predict the percentage of open-shell character in O3, which agrees well with conventional wavefunction-based methods. While these techniques are expected to be required in cases near the Coulson–Fischer point, we find that they will be less necessary in diradical systems with well-separated electrons, such as (C21H13)2. Results based on energies from self-interaction-corrected generalized gradient, local density, and Hartree–Fock approximations and experimental results are in generally good agreement for O3. These results help form the basis for deriving extended Heisenberg-like Hamiltonians that are needed for descriptions of molecular magnets when there are competing low-energy electronic configurations.

Chemistry↗

DOE OSTI2021

Measurement report: Molecular composition, optical properties, and radiative effects of water-soluble organic carbon in snowpack samples from northern Xinjiang, China

Abstract. Water-soluble organic carbon (WSOC) in the cryosphere has an important impact on the biogeochemistry cycling and snow–ice surface energy balance through changes in the surface albedo. This work reports on the chemical characterization of WSOC in 28 representative snowpack samples collected across a regional area of northern Xinjiang, northwestern China. We employed multimodal analytical chemistry techniques to investigate both bulk and molecular-level composition of WSOC and its optical properties, informing the follow-up radiative forcing (RF) modeling estimates. Based on the geographic differences and proximity of emission sources, the snowpack collection sites were grouped as urban/industrial (U), rural/remote (R), and soil-influenced (S) sites, for which average WSOC total mass loadings were measured as 1968 ± 953 ng g−1 (U), 885 ± 328 ng g−1 (R), and 2082 ± 1438 ng g−1 (S), respectively. The S sites showed the higher mass absorption coefficients at 365 nm (MAC365) of 0.94 ± 0.31 m2 g−1 compared to those of U and R sites (0.39 ± 0.11 m2 g−1 and 0.38 ± 0.12 m2 g−1, respectively). Bulk composition of WSOC in the snowpack samples and its basic source apportionment was inferred from the excitation–emission matrices and the parallel factor analysis featuring relative contributions of one protein-like (PRLIS) and two humic-like (HULIS-1 and HULIS-2) components with ratios specific to each of the S, U, and R sites. Additionally, a sample from site 120 showed unique pollutant concentrations and spectroscopic features remarkably different from all other U, R, and S samples. Molecular-level characterization of WSOC using high-resolution mass spectrometry (HRMS) provided further insights into chemical differences among four types of samples (U, R, S, and 120). Specifically, many reduced-sulfur-containing species with high degrees of unsaturation and aromaticity were uniquely identified in U samples, suggesting an anthropogenic source. Aliphatic/protein-like species showed the highest contribution in R samples, indicating their biogenic origin. The WSOC components from S samples showed high oxygenation and saturation levels. A few unique CHON and CHONS compounds with high unsaturation degree and molecular weight were detected in the 120 sample, which might be anthraquinone derivatives from plant debris. Modeling of the WSOC-induced RF values showed warming effects of 0.04 to 0.59 W m−2 among different groups of sites, which contribute up to 16 % of that caused by black carbon (BC), demonstrating the important influences of WSOC on the snow energy budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

A road map to cosmological parameter analysis with third-order shear statistics: III. Efficient estimation of third-order shear correlation functions and an application to the KiDS-1000 data

Context. Third-order lensing statistics contain a wealth of cosmological information that is not captured by second-order statistics. However, the computational effort it takes to estimate such statistics in forthcoming stage IV surveys is prohibitively expensive. Aims. We derive and validate an efficient estimation procedure for the three-point correlation function (3PCF) of polar fields such as weak lensing shear. We then use our approach to measure the shear 3PCF and the third-order aperture mass statistics on the KiDS-1000 survey. Methods We constructed an efficient estimator for third-order shear statistics that builds on the multipole decomposition of the 3PCF. We then validated our estimator on mock ellipticity catalogs obtained from N -body simulations. Finally, we applied our estimator to the KiDS-1000 data and presented a measurement of the third-order aperture statistics in a tomographic setup. Results. Our estimator provides a speedup of a factor of ∼100–1000 compared to the state-of-the-art estimation procedures. It is also able to provide accurate measurements for squeezed and folded triangle configurations without additional computational effort. We report a significant detection of tomographic third-order aperture mass statistics in the KiDS-1000 data (S/N = 6.69). Conclusions. Our estimator will make it computationally feasible to measure third-order shear statistics in forthcoming stage IV surveys. Furthermore, it can be used to construct empirical covariance matrices for such statistics.

Astronomy & Astrophysics↗

DOE OSTI2021

Physics of the Inverted Harmonic Oscillator: From the lowest Landau level to event horizons

In this work, we present the inverted harmonic oscillator (IHO) Hamiltonian as a paradigm to understand the quantum mechanics of scattering and time-decay in a diverse set of physical systems. As one of the generators of area preserving transformations, the IHO Hamiltonian can be studied as a dilatation generator, squeeze generator, a Lorentz boost generator, or a scattering potential. In establishing these different forms, we demonstrate the physics of the IHO that underlies phenomena as disparate as the Hawking–Unruh effect and scattering in the lowest Landau level (LLL) in quantum Hall systems. We derive the emergence of the IHO Hamiltonian in the LLL in a gauge invariant way and show its exact parallels with the Rindler Hamiltonian that describes quantum mechanics near event horizons. This approach of studying distinct physical systems with symmetries described by isomorphic Lie algebras through the emergent IHO Hamiltonian enables us to reinterpret geometric response in the lowest Landau level in terms of relativistic effects such as Wigner rotation. Further, the analytic scattering matrix of the IHO points to the existence of quasinormal modes (QNMs) in the spectrum, which have quantized time-decay rates. We present a way to access these QNMs through wave packet scattering, thus proposing a novel effect in quantum Hall point contact geometries that parallels those found in black holes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Structural and Dynamical Roles of Bound Polymer Chains in Rubber Reinforcement

The addition of nanofillers to rubber matrices is a powerful route to improve the mechanical properties. Here, we focus on a molecular understanding of basic mechanisms that are important for the reinforcement in rubbers. The key role in this process is ascribed to bound rubber (BR) that engages with the matrix as well as with adjacent nanofillers. To date, this understanding has been impeded by the lack of experimental tools to directly probe the BR chains buried in a polymer matrix composed of the same polymer. To tackle this challenge, we combine neutron scattering/spectroscopy techniques with isotope-labeling and molecular dynamics simulations. The system is a simplified carbon-black-filled polybutadiene. The combined experimental and computational results provide new insights into the local structural and dynamical heterogeneities of BR chains and their interactions with the matrix polymer, highlighting (i) the structural partition of the bound chains into three components (i.e., trains, loops, and tails) and their fractions; (ii) their dynamical hierarchies, i.e., the trains that remain immobile on the filler surface, the loops that are fairly large and hence allow the interdigitation of matrix chains, and the tails with their unique characteristics to reach far out into the matrix and entangle with matrix chains. These multiple roles of the constituent components of the BR chains promote the formation of a well-developed adhesive polymer–filler interface, enhancing the elastic property of a filled rubber. Finally, the comprehensive understanding derived and validated by the model rubber will be translatable to many other polymer nanocomposites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

BEYONDPLANCK X. Planck Low Frequency Instrument frequency maps with sample-based error propagation

We present Planck Low Frequency Instrument (LFI) frequency sky maps derived within the BEYONDPLANCK framework. This framework draws samples from a global posterior distribution that includes instrumental, astrophysical, and cosmological parameters, and the main product is an entire ensemble of frequency sky map samples, each of which corresponds to one possible realization of the various modeled instrumental systematic corrections, including correlated noise, time-variable gain, as well as far sidelobe and bandpass corrections. This ensemble allows for computationally convenient end-to-end propagation of low-level instrumental uncertainties into higher-level science products, including astrophysical component maps, angular power spectra, and cosmological parameters. We show that the two dominant sources of LFI instrumental systematic uncertainties are correlated noise and gain fluctuations, and the products presented here support – for the first time – full Bayesian error propagation for these effects at full angular resolution. We compared our posterior mean maps with traditional frequency maps delivered by the Planck Collaboration, and find generally good agreement. The most important quality improvement is due to significantly lower calibration uncertainties in the new processing, as we find a fractional absolute calibration uncertainty at 70 GHz of Δg 0 /g 0 = 5 x 10 -5 , which is nominally 40 times smaller than that reported by Planck 2018. However, we also note that the original Planck 2018 estimate has a nontrivial statistical interpretation, and this further illustrates the advantage of the new framework in terms of producing self-consistent and well-defined error estimates of all involved quantities without the need of ad hoc uncertainty contributions. We describe how low-resolution data products, including dense pixel-pixel covariance matrices, may be produced from the posterior samples directly, without the need for computationally expensive analytic calculations or simulations. We conclude that posterior-based frequency map sampling provides unique capabilities in terms of low-level systematics modeling and error propagation, and may play an important role for future Cosmic Microwave Background (CMB) B-mode experiments aiming at nanokelvin precision.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

FUN-BioCROP model with litter decomposition parameters derived from the LIDET dataset

This repository contains the code and data necessary to run the FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms, and Protection) model with litter decomposition parameters derived from a modified Monte Carlo simulation that used the Long-term Intersite Decomposition Experiment Team dataset. Related publication:Juice, S.M., Ridgeway, J.R., Hartman, M.D., Parton, W.J., Berardi, D.M., Sulman, B.N., Allen, K.E., & Brzostek, E.R. Reparameterizing litter decomposition using a simplified Monte Carlo method improves litter decay simulated by a microbial model and alters bioenergy soil carbon estimates. Description of Files: FUNBioCROP_LIDET Study.Rmd R code with FUN-BioCROP model that can be run with 10 different sets of parameters for litter decomposition (Baseline, LIDET, or eight other best parameter sets identified in the modified Monte Carlo simulation. CORPSE Functions_Bioenergy_V2.R Code with CORPSE model functions, called by FUNBioCROP_LIDET Study.Rmd Model Input Data: bulk.csv, bulk_till.csv, rhizo.csv, rhizo_till.csv, litter.csv Initial C and N (kg C or N/m2) pool values for each soil compartment, final values from spin up. All five files have the same columns: (Column - Description - Units) uFastC - Unprotected fast decomposing carbon - kg carbon/m2 uSlowC - Unprotected slow decomposing carbon - kg carbon/m2 uNecroC - Unprotected necromass carbon - kg carbon/m2 pFastC - Protected fast decomposing carbon - kg carbon/m2 pSlowC - Protected slow decomposing carbon - kg carbon/m2 pNecroC - Protected necromass carbon - kg carbon/m2 livingMicrobeC - Carbon in living microbial biomass - kg carbon/m2 uFastN - Unprotected fast decomposing nitrogen - kg nitrogen/m2 uSlowN - Unprotected slow decomposing nitrogen - kg nitrogen/m2 uNecroN - Unprotected necromass nitrogen - kg nitrogen/m2 pFastN - Protected fast decomposing nitrogen - kg nitrogen/m2 pSlowN - Protected slow decomposing nitrogen - kg nitrogen/m2 pNecroN - Protected necromass nitrogen - kg nitrogen/m2 inorganicN - Inorganic nitrogen - kg nitrogen/m2 CO2 - Carbon in carbon dioxide - kg carbon/m2 livingMicrobeN - Nitrogen in living microbial biomass - kg nitrogen/m2 Model Input Data: FluxTower_AvgSoilT.csv: Average daily soil temperature (oC) at 10 cm depth at University of Illinois Urbana-Champaign (UIUC) Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: FluxTower_AvgSoilVWC.csv: Average daily soil volumetric water content (VWC) at 10 cm depth at UIUC Energy Farm flux tower from 7/2008-3/2016. (One year of averaged data) Model Input Data: input_CCS_LIDET Study.csv: This file has daily data to run FUN-BioCROP (Column - Description - Units): yr - calendar year - year doy - day of year (1 to 365) (no leap year) - day anpp - aboveground NPP (DayCent) - kg C/m2/day bnpp - belowground NPP (DayCent) - kg C/m2/day aglivc - live aboveground biomass carbon (DayCent) - kg C/m2 bglivcj - live juvenile fine root biomass carbon (DayCent) - kg C/m2 bglivcm - live mature fine root biomass carbon (DayCent) - kg C/m2 aglivn - live aboveground biomass nitrogen (DayCent) - kg N/m2 bglivnj - live juvenile fine root biomass nitrogen (DayCent) - kg N/m2 bglivnm - live mature fine root biomass nitrogen (DayCent) - kg N/m2 nyr - simulation year - year cult - indicates a cultivation event (0 or 1) crop - indicates a new crop (0 or 1) fert - indicates a fertilizer event (0 or 1) frst - indicates the first day of the growing season (0 or 1) harv - indicates a harvest event (0 or 1) last - indicates the end of the growing season (0 or 1) croptype - crop type (0=none; 1=alfalfa; 2=corn; 3=grass clover pasture; 4=soybean; 5=wheat) cropsrl - crop specific root length - mm/g root cultrhizmix - fraction of rhizosphere mixed with bulk soil during cultivation (0.0-1.0) - fraction cultlitmix - fraction of litter mixed with bulk soil during cultivation (0.0-1.0) - fraction harvremov - fraction of above ground biomass removed during harvest (0.0-1.0) - fraction fertamt - fertilization amount - g N/m2 lifehist - plant life history (0 = annual, 1 = perennial) froot_turnover_c - amount of C in fine root turnover - kg C/m2 froot_turnover_n - amount of N in fine root turnover - kg N/m2 agrd_turnover_c - amount of C in aboveground biomass turnover - kg C/m2 agrd_turnover_n - amount of N in aboveground biomass turnover - kg N/m2 leaf_litter_fastfrac - Fast decomposing fraction of leaf litter (0.0-1.0) - fraction root_litter_fastfrac - Fast decomposing fraction of root litter (0.0-1.0) - fraction root_diameter - root diameter - mm root_length - root length - mm root/m2 rhizo_frac - fraction of total soil volume that is rhizosphere (0.0 - 1.0) - fraction date - date in format YYYY-MM-DD Instructions: Save the model code ("FUN-BioCROP_LIDET Study.Rmd") and accompanything files (data streams and CORPSE function code) in the same folder. In model code "Chunk 3: Load CORPSE Data Streams" set the working directory (setwd) to the folder with the files saved in step #1. In "Chunk 5: Define LIDET parameter sets" select the litter decomposition parameter set to be used in the run, and comment out all other sets. If changing any parameter values, edit them in "Chunk 6: Load parameters." Run all chunks up to and including "Chunk 10: Prepare Data for Export." In "Chunk 11: Export Output Data" edit data frames for export and filenames, as necessary. "Chunk 12: Graph Total Soil C" makes a figure of C remaining over the model run period. Description of each model chunk (in file FUN-BioCROP_LIDET Study.Rmd): Chunk 1: Remove all functions, clear memory. Removes all functions from R environment, clears the memory. Chunk 2: Load Packages. Loads packages necessary to run the code. Chunk 3: Load CORPSE Data Streams. Sets the working directory and loads the data files necessary to run CORPSE. Chunk 4: Load CORPSE Functions. Loads the R script with CORPSE functions from the working directory, "CORPSE Functions_Bioenergy_V2.R". Chunk 5: Define LIDET parameter sets. Has ten different parameter sets for litter decomposition tested in this study: Baseline parameters, LIDET parameters, and the other 8 best performing parameter sets identified in the modified Monte Carlo. To run the model, all but one parameter set must be commented out. Chunk 6: Load Parameters. Loads all fixed parameters to run the model. Data frame with definitions of parameters is in the CORPSE function script "CORPSE Functions_Bioenergy_V2.R" Chunk 7: Prepare Data Streams. Takes data streams loaded in Chunk 3 and puts them in the format necessary to run the model. The model is coded to run at least two sites at a time, so if only one site is being run it must be run in duplicate. Individual data tables of daily values are created in this chunk from the input data file. Chunk 8: Set Initial Conditions. Creates data tables of soil C and N pools for each soil compartment (rhizo_till, rhizo, bulk_till, bulk, litter) and loads initial values into the data tables. Creates lists for each soil compartment to hold model output. Chunk 9: Load FUN Data and Set Up Matrices. Uses DayCent data to calculate FUN input data: root and leaf N demand, total N demand, plant CN, leaf N available for retranslocation, and litter production. Creates matrices for FUN model outputs. Chunk 10: Run Model. Runs the model. Chunk 11: Prepare Data for Export. Combines data from each day saved as lists into data frames for each soil compartment. Adds values from all soil compartments together to calculate total soil values, creates separate data frames for each soil C and N pool (e.g., protected slow C) for the total soil value. Adds different C and N pools together to calculate total soil C and N for all layers. Creates data frame of ratio of protected to unprotected SOC. Organizes FUN data for export. Chunk 12: Export Results. Exports CSV files of model results to the working directory. Chunk 13: Graph Total Soil C. Makes figure of C remaining over time. Related Links: Original FUN-BioCROP model: https://github.com/BrzostekEcologyLab/FUN-BioCROP LIDET dataset: https://andlter.forestry.oregonstate.edu/data/abstract.aspx?dbcode=TD023

Juice, Stephanie↗