Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “thermal response functions”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Widely Tunable Optical and Thermal Properties of Dirac Semimetal Cd 3 As 2

In this paper, a detailed analysis of the temperature-dependent optical properties of epitaxially grown cadmium arsenide (Cd 3 As 2 ), a newly discovered 3D Dirac semimetal is reported. Fermi level tuning—instigated from Pauli-blocking in the linear Dirac cone—and varying Drude response, generate large variations in the mid- and far-infrared optical properties. Thermo-optic shifts larger than those of traditional III–V semiconductors are demonstrated. Electron scattering rate, plasma frequency edge, Fermi level shift, optical conductivity, and electron effective mass analysis of Cd 3 As 2 thin-films are quantified and discussed in detail. The ab initio density functional study and experimental analysis of epitaxially grown Cd 3 As 2 promise applications for nanophotonic and nanoelectronic devices, such as reconfigurable metamaterials and metasurfaces, nanoscale thermal emitters, and on-chip directional antennas.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Time Evolution of Correlation Functions in Quantum Many-Body Systems

We give rigorous analytical results on the temporal behavior of two-point correlation functions—also known as dynamical response functions or Green’s functions—in closed many-body quantum systems. We show that in a large class of translation-invariant models the correlation functions factorize at late times $\langle$A ( t ) B$\rangle$ β → $\langle$ A$\rangle$ β $\langle$ B$\rangle$ β , thus proving that dissipation emerges out of the unitary dynamics of the system. We also show that for systems with a generic spectrum the fluctuations around this late-time value are bounded by the purity of the thermal ensemble, which generally decays exponentially with system size. For autocorrelation functions we provide an upper bound on the timescale at which they reach the factorized late time value. Remarkably, this bound is only a function of local expectation values and does not increase with system size. We give numerical examples that show that this bound is a good estimate in nonintegrable models, and argue that the timescale that appears can be understood in terms of an emergent fluctuation-dissipation theorem. Our study extends to further classes of two point functions such as the symmetrized ones and the Kubo function that appears in linear response theory, for which we give analogous results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Revealing Fast Cu-Ion Transport and Enhanced Conductivity at the CuInP 2 S 6 –In 4/3 P 2 S 6 Heterointerface

Van der Waals layered ferroelectrics, such as CuInP 2 S 6 (CIPS), offer a versatile platform for miniaturization of ferroelectric device technologies. Control of the targeted composition and kinetics of CIPS synthesis enables the formation of stable self-assembled heterostructures of ferroelectric CIPS and nonferroelectric In 4/3 P 2 S 6 (IPS). Here, we use quantitative scanning probe microscopy methods combined with density functional theory (DFT) to explore in detail the nanoscale variability in dynamic functional properties of the CIPS-IPS heterostructure. Within, we report evidence of fast ionic transport which mediates an appreciable out-of-plane electromechanical response of the CIPS surface in the paraelectric phase. Further, we map the nanoscale dielectric and ionic conductivity properties as we thermally stimulate the ferroelectric-paraelectric phase transition, recovering the local dielectric behavior during this phase transition. Finally, aided by DFT, we reveal a substantial and tunable conductivity enhancement at the CIPS/IPS interface, indicating the possibility of engineering its interfacial properties for next generation device applications.

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

Strong anharmonicity dictates ultralow thermal conductivities of type-I clathrates

Type-I clathrate solids have attracted significant interest due to their ultralow thermal conductivities and sub- sequent promise for thermoelectric applications, yet the mechanisms underlying these properties are not well understood. Here, we extend the framework of vibrational dynamical mean-field theory (VDMFT) to calculate temperature-dependent thermal transport properties of solids using a many-body Green’s function approach. When applied to a coarse-grained description of 𝑋 8 Ga 16 Ge 30 , where 𝑋= Ba, Sr, we find that nonresonant scattering between cage acoustic modes and rattling modes leads to a reduction of acoustic phonon lifetimes and thus thermal conductivities. Moreover, we find that the moderate temperature dependence of conductivities above 300 K, which is consistent with experimental measurements, cannot be reproduced by textbook perturbation theory calculations, which predict a 𝑇 −1 dependence. Therefore, we suggest that nonperturbative anharmonic effects, including four- and higher-phonon scattering processes, are responsible for the ultralow thermal conductivities of type-I clathrates.

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

Proton discrimination in CLYC for fast neutron spectroscopy

The Cs 2 LiYCl 6 :Ce (CLYC) elpasolite scintillator is known for its response to fast and thermal neutrons along with good γ-ray energy resolution. While the 35 Cl(n,p) reaction has been identified as a potential means for CLYC-based fast neutron spectroscopy in the absence of time-of-flight (TOF), previous efforts to functionalize CLYC as a fast neutron spectrometer have been thwarted by the inability to isolate proton interactions from 6 Li(n,α) and 35 Cl(n,α) signals. This work introduces a new approach to particle discrimination in CLYC for fission spectrum neutrons using a multi-gate charge integration algorithm that provides excellent separation between protons and heavier charged particles. Neutron TOF data were collected using a 252 Cf source, an array of EJ-309 organic liquid scintillators, and a 6 Li-enriched CLYC scintillator outfitted with fast electronics. Modal waveforms were constructed corresponding to the different reaction channels, revealing significant differences in the pulse characteristics of protons and heavier charged particles at ultrafast, fast, and intermediate time scales. These findings informed the design of a pulse shape discrimination algorithm, which was validated using the TOF data. This study also proposes an iterative subtraction method to mitigate contributions from confounding reaction channels in proton and heavier charged particle pulse height spectra, opening the door for CLYC-based fast neutron and γ-ray spectroscopy while preserving sensitivity to thermal neutron capture signals.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

A statistical mechanics framework for polymer chain scission, based on the concepts of distorted bond potential and asymptotic matching

To design increasingly tough, resilient, and fatigue-resistant elastomers and hydrogels, the relationship between controllable network parameters at the molecular level (bond type, non-uniform chain length, entanglement density, etc.) to macroscopic quantities that govern damage and failure must be established. Many of the most successful constitutive models for elastomers have been rooted in statistical mechanical treatments of polymer chains. Typically, such constitutive models have used variants of the freely jointed chain model with rigid links. However, since the free energy state of a polymer chain is dominated by enthalpic bond distortion effects as the chain approaches its rupture point, bond extensibility ought to be accounted for if the model is intended to capture chain rupture. To that end, a new bond potential is supplemented to the freely jointed chain model (as derived in the u FJC framework of Buche and Silberstein (2021) and Buche et al. (2022)), which we have extended to yield a tractable, closed-form model of single chain behavior that should be amenable to continuum-level constitutive model development. Inspired by the asymptotically matched u FJC model response in both the low/intermediate chain force and high chain force regimes, a simple, quasi-polynomial bond potential energy function is derived. This bond potential exhibits harmonic behavior near the equilibrium state and anharmonic behavior for large bond stretches tending to a characteristic energy plateau (akin to the Lennard-Jones and Morse bond potentials). Using this bond potential, approximate yet highly-accurate analytical functions for bond stretch and chain force dependent upon chain stretch are established. Then, using this polymer chain model, a stochastic thermal fluctuation-driven chain rupture framework is developed. This framework is based upon a force-modified tilted bond potential that accounts for distortional bond potential energy, allowing for the derivation and subsequent calculation of the dissipated chain scission energy. Here, the cases of rate-dependent and rate-independent scission are accounted for throughout the rupture framework. The impact of Kuhn segment number on chain rupture behavior is also investigated. The model is fit to single chain mechanical response data collected from atomic force microscopy tensile tests for validation and to glean deeper insight into the molecular physics taking place. Due to their analytical nature, this polymer chain model and the associated rupture framework can in the future be implemented in finite element models accounting for fracture and fatigue in polydisperse elastomer networks.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Structural and physical properties of two distinct 2D lead halides with intercalated Cu( II )

Transition metal cation intercalation between the layers of two-dimensional (2D) metal halides is an underexplored research area. In this work we focus on the synthesis and physical property characterizations of two layered hybrid lead halides: a new compound [Cu(O 2 C–CH 2 –NH 2 ) 2 ]Pb 2 Br 4 and the previously reported [Cu(O 2 C–(CH 2 ) 3 –NH 3 ) 2 ]PbBr 4 . These compounds exhibit 2D layered crystal structures with incorporated Cu 2+ between the metal halide layers, which is achieved by combining Cu(II) and lead bromide with suitable amino acid precursors. The resultant [Cu(O 2 C–(CH 2 ) 3 –NH 3 ) 2 ]PbBr 4 adopts a 2D layered perovskite structure, whereas the new compound [Cu(O 2 C–CH 2 –NH 2 ) 2 ]Pb 2 Br 4 crystallizes with a new structure type based on edge-sharing dodecahedral PbBr5O3 building blocks. [Cu(O 2 C–CH 2 –NH 2 ) 2 ]Pb 2 Br 4 is a semiconductor with a bandgap of 3.25 eV. It shows anisotropic charge transport properties with a semiconductor resistivity of 1.44 × 10 10 Ω cm (measured along the a-axis) and 2.17 × 10 10 Ω cm (along the bc-plane), respectively. The fabricated prototype detector based on this material showed response to soft low-energy X-rays at 8 keV with a detector sensitivity of 1462.7 μCGy –1 cm –2 , indicating its potential application for ionizing radiation detection. Finally, these encouraging results are discussed together with the results from density functional theory calculations, optical, magnetic, and thermal property characterization experiments.

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

Enhancing the efficiency of time-dependent density functional theory calculations of dynamic response properties

X-ray Thomson scattering (XRTS) constitutes an essential technique for diagnosing material properties under extreme conditions, such as high pressures and intense laser heating. Time-dependent density functional theory (TDDFT) is one of the most accurate available ab initio methods for modeling XRTS spectra, as well as a host of other dynamic material properties. However, strong thermal excitations, along with the need to account for variations in temperature and density as well as the finite size of the detector significantly increase the computational cost of TDDFT simulations compared to ambient conditions. In this work, we present a broadly applicable method for optimizing and enhancing the efficiency of TDDFT calculations. Our approach is based on a one-to-one mapping between the dynamic structure factor and the imaginary time density–density correlation function, which naturally emerges in Feynman’s path integral formulation of quantum many-body theory. Specifically, we combine rigorous convergence tests in the imaginary time domain with a constraints-based attenuation of narrow-band fluctuations to improve the efficiency of TDDFT modeling without the introduction of any significant bias. As a result, we can report a speed-up by up to an order of magnitude, thus substantially reducing the burden of computational cost required for XRTS analysis.

Moldabekov, Zhandos A. [Helmholtz-Zentrum Dresden-↗
DOE OSTI2024

Phase Selection During Solidification and Solid-State Phase Transformations in an Al-10Ce-8Mn (wt pct) Alloy

In multicomponent Al-Ce alloys, and especially after additive manufacturing (AM), complex and metastable solidification microstructures are frequently observed. Here, in this research, the relationship between solidification conditions and phase selection are explored for an Al-10Ce-8Mn (wt pct) alloy using a systematic study of laser melting conditions. Three solidification modes were observed: primary Al 10 Mn 2 Ce; primary Al 20 Mn 2 Ce; and eutectic FCC Al + Al 20 Mn 2 Ce. These solidification modes were correlated to specific liquid-solid interface velocities using a simple thermal model, showing the change in primary solidification phase for low (< 6.8 × 10 −4 m/s), moderate (between 8.2 × 10 −4 and 5.9 × 10 −2 m/s) and high solidification velocities (> 6.2 x 10 −2 m/s) for the above three solidification microstructures, respectively. These results were rationalized by using interface response function (IRF) theory to describe the solidification undercooling for the possible primary intermetallic phases. The implication of the local phase selection from differing solidification conditions is summarized by a comparison of hardness which demonstrates the potential variance of Vickers hardness from 101 to 242 (VHV) by changing the laser velocity from 1 to 83 mm/s. Interestingly, on heat treatment at 400°C, the decomposition pathways of the solidification microstructure and hardness were also found to be different, thereby opening multiple pathways for spatial microstructure and property control within AM components.

Sisco, Kevin↗
DOE OSTI2024

Thermal aging effects on crosslinked polyethylene cable insulation with decabromodiphenyl ether flame retardant alternative

Decabromodiphenyl ether (decaBDE) has been extensively used as a flame retardant in several applications, including nuclear electrical cable insulation. However, decaBDE has been identified as a persistent, bioaccumulative and toxic (PBT) substance, leading to regulatory scrutiny. The Environmental Protection Agency (EPA) published a regulation on January 6, 2021, aimed at phasing out the manufacturing, processing, and distribution of decaBDE. This rule set a compliance deadline of March 8, 2021, for the manufacture and processing of decaBDE, and an extended deadline of January 6, 2023, for specific applications including wire and cable insulation in nuclear power generation facilities. In response to such regulations, RSCC, a major supplier of safety-related electrical cables and associated products to the U.S. nuclear industry updated the formula of their crosslinked polyethylene (XLPE) insulation to replace the historically used decaBDE flame retardant with an acceptable alternative. This change from the previous decaBDE-containing XLPE prompted interest in comparative performance of the two material formulations, especially with respect to characteristics relevant to safety-related function such as thermal and radiation resistance. RSCC graciously provided samples of wire insulated with the decaBDE-containing XLPE formulation and corresponding wire insulated with XLPE of the new formulation, containing a decaBDE alternative. In this work we compare characteristics of the two formulations and a previously produced commercial version of the RSCC decaBDE-containing XLPE insulation subjected to thermal aging at 150 °C and 165 °C. The comparison was focused on mechanical durability, thermal stability in the oxidative environment, and chemical structures. Briefly, • Tensile elongation at break (EAB) results showed loss of mechanical elasticity with longer aging time, as expected. Aging time dependence of EAB did not differ between the decaBDE-containing and decaBDE-alternative samples. • Subtle differences between the two materials can be detected from Fourier-transform infrared spectroscopy (FTIR) absorbance spectra in the range below 1700 cm -1 , are assumed to be related to decomposition of flame retardant additives during thermal aging. • The oxidation induction time (OIT) data seemed to show that the unaged decaBDE-containing XLPE material is more thermally stable than the unaged decaBDE-alternative material, but the discrepancy in OIT decreased with aging time and the OIT values of the two materials became similar starting with the 4 th day of aging at 165 °C. This thermal aging investigation confirmed that the mechanical durability, a key property monitored for cable qualification, was not significantly affected by the modification of the formulation with a decaBDE alternative flame-retardant system in the investigated thermal aging conditions. Further studies on the same sets of materials exposed to thermal and gamma radiation aging would further inform comparison of the materials safety-related function.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular origin of viscoelasticity and influence of methylation in mesophase pitch

The viscoelastic and thermomechanical properties of pitches are responsible for their melt-spinning behavior, which is a critical step for manufacturing high-performance pitch-based carbon fibers. Here, we systematically explore the impact of methyl group modifications on the viscoelastic and thermal properties of mesophase pitches. We employ a range of atomistic modeling approaches, including Density Functional Theory (DFT), Density Functional Tight Binding (DFTB), and Classical Molecular Mechanics (MM), to provide detailed insights into the molecular interactions and structural changes. Further, our results revealed the molecular mechanisms that promote layered structures leading to the anisotropic nature of the viscoelastic behavior of mesophase pitch. Furthermore, we propose a modified molecular representation of naphthalene-based mesophase pitch based on the analysis of x-ray diffraction measurements. This study provides fundamental insights into the molecular structures of mesophase pitch and the role of methyl groups controlling its viscosity, which offer valuable insights into mesophase-based carbon fiber production.

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

The developments in modifying functionality of lignin and its application in biocomposites

With the advancement of sustainable material innovations, renewable natural biopolymers are gradually replacing traditional metal and petroleum-based synthetic materials due to their environmental friendliness, biodegradability, and economic advantages. Lignin, the second most abundant natural aromatic polymer in the plant kingdom, has emerged as a key candidate raw material for the development of green polymer systems because of its unique phenylpropane unit network structure, high carbon content, and rich functional group characteristics. However, challenges such as the inherent structural complexity, chemical inertness, and uneven molecular weight distribution of lignin limit its direct application. By employing modification strategies such as chemical functionalization and physical regulation, researchers can precisely control its reactivity, thermal stability, and interfacial compatibility, enabling the preparation of high-performance lignin-based functional composites. Here, this paper systematically reviews the principles and methodological advancements in lignin's multi-dimensional modification technology. It analyzes the mechanisms by which various chemical and physical modification techniques enhance the mechanical properties, functional responsiveness, and environmental adaptability of materials, and discusses innovative applications in fields such as intelligent packaging, biomedical materials, and energy storage devices. Furthermore, this review addresses the key challenges encountered in the high-value transformation of lignin, with the aim of offering a theoretical framework and technical pathway for the transformative development of lignin from agricultural and forestry by-products to functional material platforms.

Functional composites↗
DOE OSTI2020

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗
DOE OSTI2021

Stability and molecular pathways to the formation of spin defects in silicon carbide

Abstract Spin defects in wide-bandgap semiconductors provide a promising platform to create qubits for quantum technologies. Their synthesis, however, presents considerable challenges, and the mechanisms responsible for their generation or annihilation are poorly understood. Here, we elucidate spin defect formation processes in a binary crystal for a key qubit candidate—the divacancy complex (VV) in silicon carbide (SiC). Using atomistic models, enhanced sampling simulations, and density functional theory calculations, we find that VV formation is a thermally activated process that competes with the conversion of silicon (V Si ) to carbon monovacancies (V C ), and that VV reorientation can occur without dissociation. We also find that increasing the concentration of V Si relative to V C favors the formation of divacancies. Moreover, we identify pathways to create spin defects consisting of antisite-double vacancy complexes and determine their electronic properties. The detailed view of the mechanisms that underpin the formation and dynamics of spin defects presented here may facilitate the realization of qubits in an industrially relevant material.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Super-expansive thermo-reversible interstitial solid solution of nanocrystal superlattices with mesogens

Designing superlattices of nanocrystals to mimic and extend the properties of atomic crystals has been a long-standing motivation in materials chemistry. Interstitial solid solutions, such as steel, are well-studied atomic lattices in which mobile components move among the interstices. These materials exhibit unique properties, including reversible structural changes and phase transitions. Interstitial solid solutions possess unique dynamic structures and reversible responses, which motivate the creation of their colloidal equivalents. Here, in this study, we report a fully thermo-reversible colloidal interstitial solid solution by combining liquid crystals and nanocrystals functionalized with promesogenic ligands. Mesogen molecules fill and diffuse among the interstices of a superlattice, resulting in a super-large thermal expansivity. The approach uses a modular design of interparticle interactions, allowing control of interparticle distance, microstructure and transition between crystallographic forms.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Applying Dynamic Strain on Thin Oxide Films Immobilized on a Pseudoelastic Nickel-Titanium Alloy

Direct alteration of materials structure/function through strain is a burgeoning area of research that has allowed for novel properties of materials to emerge. Tuning materials structure can be achieved by controlling an external force imposed on materials and inducing stress-strain responses (i.e., applying dynamic strain). Electroactive thin films are typically deposited on shape or volume tunable elastic substrates, where mechanical loading (i.e. compression or tension) can affect films structure and function through imposed strain. Here we summarize methods for straining n-type doped titanium dioxide (TiO2) films prepared by a thermal treatment of a pseudo-elastic nickel-titanium alloy. The main purpose of the described methods is to study how strain affects electrocatalytic activities of metal oxide, specifically hydrogen evolution and oxygen evolution reactions. The same system can be adapted to study the effect of strain more broadly. Strain engineering can be applied for optimization of a material function, as well as for design of adjustable, multifunctional (photo)electrocatalytic materials under external stress control.

92 SC - Basic Energy Sciences (BES) (SC-22)↗
DOE OSTI2024

Enhancing actuation control in monolithic liquid crystal elastomer films with dual heat and light responsiveness

This study presents a dual heat and light-responsive monolithic liquid crystal elastomers (LCEs) film designed to enhance actuation control. The LCE films were capable of various motion types, including curvature reversal, oscillation, and sustainable flipping. By incorporating an azobenzene compound, the films exhibited additional functionalities such as enhanced curvature reversal, accelerated oscillation, and the ability to lock continuous flipping motion upon light stimulation. Furthermore, by adjusting the aspect ratio (AR) of a single film, different thermal motion states were realized, demonstrating the flexibility of this approach. Observations using an NIR camera confirmed temperature fluctuations corresponding to thermal exchange, and the influence of UV light on motion regulation was quantitatively analyzed. In conclusion, the ability to regulate motion through dual stimuli paves the way for improved dynamic motion control and thermal management in advanced robotic applications.

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

Responses of Boreal Forest Ecosystems and Permafrost to Climate Change and Disturbances: A Modeling Perspective

The north circumpolar region contains a large amount of carbon. This carbon storage is vulnerable due to permafrost degradation and wildfire disturbances under ongoing and projected climate change. Climate warming and wildfires change soil organic horizons gradually or abruptly, and modify permafrost thermal-hydrology and biogeochemistry, ecosystem structures, functions, and capability of sequestrating rising atmospheric CO2. Land models do not fully take accounts of these interactions and its complexity in the high latitude. This chapter describes a terrestrial ecosystem model with dynamic organic soil module (DOS-TEM) and its unique freezing-thawing algorithm, and presents key results of its applications mainly in boreal forests of Alaska. The DOS-TEM explicitly considers interactions of soil thermal and hydrological processes, permafrost degradation and the direct and indirect effects of wildfire disturbances, in addition to soil–plant C and N cycles. We first introduce four modules of DOS-TEM, focusing on its disturbance module and coupling with a dynamic organic soil module. Then we describe and validate DOS-TEM’s freezing-thawing algorithm and development based on two-directional Stefan algorithm (TDSA). Finally, we apply the DOS-TEM at site and region scales, with a focus on model ability to dynamically simulate soil organic thickness under warming and wildfires, and consequent impacts on permafrost in the Yukon River Basin. We conclude that land surface model development is urgently needed to include other critical landscape processes, such as thermalkarst and other disturbances, to synchronize thermal-hydrological-biogeochemical processes, and to incorporate an advanced understanding of biospheric feedbacks to atmosphere and ecosystems. Such a complexity of modeling scope is plausible with advancement of high performance computing.

Yi, Shuhua↗