Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “AUTOCORRELATION”

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

Microscopic signatures of yielding in concentrated nanoemulsions under large-amplitude oscillatory shear

Here, we report x-ray photon correlation spectroscopy (XPCS) experiments on a series of concentrated oil-in-water nanoemulsions with varying droplet volume fraction subjected to in situ steady-state large-amplitude oscillatory shear (LAOS). The shear strain causes periodic echoes in the x-ray speckle patterns that lead to peaks in the intensity autocorrelation function. Above an onset strain amplitude that depends on nanoemulsion concentration, the peaks become attenuated, signaling spatially heterogeneous, shear-induced droplet dynamics. These dynamics include irreversible re- arrangements among the droplets that occur in some regions of the nanoemulsions during a given shear cycle and residual strain-like displacements in those regions that do not re-arrange. The wave-vector dependence of the peak attenuation indicates a power-law distribution in the size of regions undergoing shear-induced rearrangement that is similar to that observed previously in LAOS-XPCS measurements on concentrated nanocolloidal gels. The values of the onset strains for re-arrangement correlate with the concentration-dependent macroscopic yielding behavior of the nanoemulsions. Specically, they occur below the strains at which the nanoemulsions become effectively uidized and, except for the lowest concentration nanoemulsion in the study, signicantly above the threshold strain for nonlinear rheological response.

36 MATERIALS SCIENCE↗

Effect of transport mechanism on the coarsening of bicontinuous structures: A comparison between bulk and surface diffusion

Here, coarsening of bicontinuous microstructures is observed in a variety of systems, such as nanoporous metals and mixtures that have undergone spinodal decomposition. To better understand the morphological evolution of these structures during coarsening, we compare the morphologies resulting from two different coarsening mechanisms, surface and bulk diffusion. We perform phase-field simulations of coarsening via each mechanism in a two-phase mixture at nominal volume fractions of 50%-50% and 36%-64%, and the simulated structures are characterized in terms of topology (genus density), the interfacial shape distribution, structure factor, and autocorrelations of phase and mean curvature. We observe self-similar evolution of morphology and topology and agreement with the expected power laws for dynamic scaling, in which the characteristic length scale increases over time proportionally to t 1/4 for surface diffusion and t 1/3 for bulk diffusion. While we observe the expected difference in the coarsening kinetics, we find that differences in self-similar morphology due to coarsening mechanism are relatively small, although typically they are larger at 36% volume fraction than at 50% volume fraction. In particular, we find that bicontinuous structures coarsened via surface diffusion have lower scaled genus densities than structures coarsened via bulk diffusion. We also compare the self-similar morphologies to those in literature and to two model bicontinuous structures, namely, constant-mean-curvature surfaces based on the Schoen G minimal surface and random leveled-wave structures. The average scaled mean curvatures of these model structures agree reasonably with those of the coarsened structures at both 36% and 50%, but we find substantial disagreements in the scaled genus densities and the standard deviations of mean curvature.

36 MATERIALS SCIENCE↗

Cone beam neutron interferometry: From modeling to applications

Phase-grating moiré interferometers (PGMIs) have emerged as promising candidates for the next generation of neutron interferometry, enabling the use of a polychromatic beam and manifesting interference patterns that can be directly imaged by existing neutron cameras. However, the modeling of the various PGMI configurations is limited to cumbersome numerical calculations and backward propagation models which often do not enable one to explore the setup parameters. Here we generalize the Fresnel scaling theorem to introduce a k -space model for PGMI setups illuminated by a cone beam, thus enabling an intuitive forward propagation model for a wide range of parameters and experimental setups. The interference manifested by a PGMI is shown to be a special case of the Talbot effect, and the optimal fringe visibility is shown to occur at the moiré location of the Talbot distances. We derive analytical expressions for the contrast and the propagating intensity profiles in various conditions and provide the first analysis of the PGMI dark-field imaging signal when considering sample characterization. The model's predictions are compared to experimental measurements and good agreement is found between them. Last, we propose and experimentally verify a method to recover contrast at typically inaccessible PGMI autocorrelation lengths. The presented work provides a toolbox for analyzing and understanding existing PGMI setups and their future applications, for example extensions to two-dimensional PGMIs and characterization of samples with nontrivial structures. Published by the American Physical Society 2024

Sarenac, D. (ORCID:0000000185753367)↗

Theory of x-ray photon correlation spectroscopy for multiscale flows

Complex multiscale flows associated with instabilities and turbulence are commonly induced under high-energy density (HED) conditions, but accurate measurement of their transport properties has been challenging. X-ray photon correlation spectroscopy (XPCS) with coherent x-ray sources can, in principle, probe material dynamics to infer transport properties using time autocorrelation of density fluctuations. Here we develop a theoretical framework for utilizing XPCS to study material diffusivity in multiscale flows. We extend single-scale shear flow theories to broadband flows using a multiscale analysis that captures shear and diffusion dynamics. Our theory is validated with simulated XPCS for Brownian particles advected in multiscale flows. We demonstrate the versatility of the method over several orders of magnitude in timescale using sequential-pulse XPCS, single-pulse x-ray speckle visibility spectroscopy (XSVS), and double-pulse XSVS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Generalized approach for rapid entropy calculation of liquids and solids

We build a comprehensive methodology for the fast computation of entropy across both solid and liquid phases. The proposed method utilizes a single trajectory of molecular dynamics (MD) to facilitate the calculation of entropy, which is composed of three components. The electronic entropy is determined through the temporal average acquired from density functional theory MD simulations. The vibrational entropy, typically the predominant contributor to the total entropy, even within the liquid state, is evaluated by computing the phonon density of states via the velocity autocorrelation function. The most arduous component to quantify, the configurational entropy, is assessed by probability analysis of the local structural arrangement and atomic distribution. We illustrate, through a variety of examples, that this method is both a versatile and valid technique for characterizing the thermodynamic states of both solids and liquids. Furthermore, this method is employed to expedite the calculation of melting temperatures, demonstrating its practical utility in computational thermodynamics.

36 MATERIALS SCIENCE↗

Superspin renormalization and slow relaxation in random spin systems

We develop an excited-state real-space renormalization group (RSRG-X) formalism to describe the dynamics of conserved densities in randomly interacting spin-12 systems. Our formalism is suitable for systems with U(1) and Z2 symmetries, and we apply it to chains of randomly positioned spins with dipolar XX+YY interactions, as arise in Rydberg quantum simulators and other platforms. The formalism generates a sequence of effective Hamiltonians that provide approximate descriptions for dynamics on successively smaller energy scales. These effective Hamiltonians involve “superspins”: two-level collective degrees of freedom constructed from (anti)aligned microscopic spins. Conserved densities can then be understood as relaxing via coherent collective spin flips. For the well-studied simpler case of randomly interacting nearest-neighbor XX+YY chains, the superspins reduce to single spins. Our formalism also leads to a numerical method capable of simulating the dynamics up to an otherwise inaccessible combination of large system size and late time. Focusing on disorder-averaged infinite-temperature autocorrelation functions, in particular the spin survival probability Sp¯(t), we demonstrate quantitative agreement between our algorithm and exact diagonalization (ED) at low but nonzero frequencies. Such agreement holds for chains with nearest-neighbor, next-nearest-neighbor, and long-range dipolar interactions. Our results indicate decay of Sp¯(t) slower than any power law and feature no significant deviation from the ∼1/ln2(t) asymptote expected from the infinite-randomness fixed-point of the nearest-neighbor model. We also apply the RSRG-X formalism to two-dimensional long-range systems of moderate size and find slow late-time decay of Sp¯(t).

Zhao, Yi J↗

First experimental feasibility study of VIPIC: a custom-made detector for X-ray speckle measurements

The Vertically Integrated Photon Imaging Chip (VIPIC) was custom-designed for X-ray photon correlation spectroscopy, an application in which occupancy per pixel is low but high time resolution is needed. VIPIC operates in a sparsified streaming mode in which each detected photon is immediately read out as a time- and position-stamped event. This event stream can be fed directly to an autocorrelation engine or accumulated to form a conventional image. The detector only delivers non-zero data (sparsified readout), greatly reducing the communications overhead typical of conventional frame-oriented detectors such as charge-coupled devices or conventional hybrid pixel detectors. This feature allows continuous acquisition of data with timescales from microseconds to hours. In this work VIPIC has been used to measure X-ray photon correlation spectroscopy data on polystyrene latex nano-colliodal suspensions in glycerol and on colloidal suspensions of silica spheres in water. Relaxation times of the nano-colloids have been measured for different temperatures. These results demonstrate that VIPIC can operate continuously in the microsecond time frame, while at the same time probing longer timescales.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Correlations and Cascades in Magnetized Turbulence

Many terrestrial and astrophysical plasmas encompass very large dynamical ranges in space and time, which are not accessible by direct numerical simulations. Thus, idealized subvolumes are often used to study small-scale effects including the dynamics of magnetized turbulence. A significant aspect of magnetized turbulence is the transfer of energy from large to small scales, in part through the operation of a turbulent cascade. In this work, we present a new shell-to-shell energy transfer analysis framework for understanding energy transfer within magnetized turbulence and in particular, through the cascade. We demonstrate the viability of this framework through application to a series of isothermal subsonic and supersonic simulations of compressible magnetized turbulence and utilize results from this analysis to establish a non-linear benchmark for compressible magnetized turbulence in the subsonic regime. We further study how the autocorrelation time of the driving and its normalization systematically change properties of compressible magnetized turbulence. For example, we find that 6-in-time forcing with a constant energy injection leads to a steeper slope in kinetic energy spectrum and less efficient small-scale dynamo action. We examine how these results can impact a range of diagnostics relevant for a range of terrestrial and astrophysical applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fine root morphology is phylogenetically structured, but nitrogen is related to the plant economics spectrum in temperate trees

Summary Plant functional traits have revealed trade‐offs related to life‐history adaptations, geographical distributions, and ecosystem processes. Fine roots are essential in plant resource acquisition and play an important role in soil carbon cycling. Nonetheless, root trait variation is still poorly quantified and rarely related to the rest of the plant. We examined chemical and morphological traits of 34 temperate arbuscular mycorrhizal tree species, representing three main angiosperm clades (super‐orders asterid, magnoliid and rosid). We tested to what extent fine root chemical and morphological traits were correlated similarly to the leaf economical spectrum ( LES ) or were structured by ancestral affiliations among species. Root traits did not display the same trade‐offs as leaves (e.g. specific root length was not correlated with root N, whereas specific leaf area was correlated with leaf N). Moreover, 75% of below‐ground traits were phylogenetically structured according to Pagel's λ and Abouheif's C mean autocorrelation tests, as opposed to 28% of above‐ground traits. Magnoliids showed thicker, less branched roots than asterids or rosids, but rosid roots exhibited lower N and higher non‐acid‐hydrolysable (e.g. lignin) content than other species. In contrast, leaf traits did not differ significantly among super‐orders. At the whole‐tree level, chemical traits such as nitrogen tissue content and lignin content were correlated between above and below‐ground organs. The distribution of root traits in woody temperate trees was better explained by shared ancestry than by the nutrient content and structural trade‐offs expected by the LES hypothesis. Root chemistry and morphology differed substantially among species belonging to different super‐orders, suggesting deep divergences in resource acquisition strategies among major angiosperm groups. Although we found partial support for the idea of whole‐plant integration based on corresponding nitrogen content across all organs (i.e. a plant economics spectrum), our study stresses phylogenetic affiliation as the primary driver of root trait distributions among angiosperms, a pattern that could be easily overlooked based solely on above‐ground observations.

Valverde‐Barrantes, Oscar J.↗

Time series methods for the analysis of soundscapes and other cyclical ecological data

Biodiversity monitoring has entered an era of ‘big data’, exemplified by a near-continuous collection of sounds, images, chemical and other signals from organisms in diverse ecosystems. Such data streams have the potential to help identify new threats, assess the effectiveness of conservation interventions, as well as generate new ecological insights. However, appropriate analytical methods are often still missing, particularly with respect to characterizing cyclical temporal patterns. Here, we present a framework for characterizing and analysing ecological responses that represent nonstationary, complex temporal patterns and demonstrate the value of using Fourier transforms to decorrelate continuous data points. In our example, we use a framework based on three approaches (spectral analysis, magnitude squared coherence, and principal component analysis) to characterize differences in tropical forest soundscapes within and across sites and seasons in Gabon. By reconstructing the underlying, cyclic behaviour of the soundscape for each site, we show how one can identify circadian patterns in acoustic activity. Soundscapes in the dry season had a complex diel cycle, requiring multiple harmonics to represent daily variation, while in the wet season there was less variance attributable to the daily cyclic patterns. Our framework can be applied to most continuous, or near-continuous ecological data collected at a fine temporal resolution, allowing ecologists to explore patterns of temporal autocorrelation at multiple levels for biologically meaningful trends. Such methods will become indispensable as biological big data are used to understand the impact of anthropogenic pressures on biodiversity and to inform efforts to mitigate them.

54 ENVIRONMENTAL SCIENCES↗

Landscape‐scale lidar analysis of aboveground biomass distribution in secondary Brazilian Atlantic Forest

Abstract Secondary forests account for more than half of tropical forests and represent a growing carbon sink, but rates of biomass accumulation vary by a factor of two or more even among plots in the same landscape. To better understand the drivers of this variability, we used airborne lidar to measure forest canopy height and estimate biomass over 4529 ha at Serra do Conduru Park in Southern Bahia, Brazil. We measured trees in 30 georeferenced field plots (0.25‐ha each) to estimate biomass using allometry. Then we estimated aboveground biomass density ( ABD ) across the lidar study area using a statistical model developed from our field plots. This model related the 95th percentile of the distribution of lidar return heights to ABD . We overlaid this map of ABD on a Landsat‐derived forest age map to determine rates of biomass accumulation. We found rapid initial biomass regeneration (~6 Mg/ha yr), which slowed as forests aged. We also observed high variability in both height and biomass across the landscape within forests of similar age. Nevertheless, a regression model that accounted for spatial autocorrelation and included forest age, slope, and distance to roads or open areas explained 62 and 77 percent of the landscape variation in ABD and canopy height, respectively. Thus, while there is high spatial heterogeneity in forest recovery, and the drivers of this heterogeneity warrant further investigation, we suggest that a relatively simple set of predictor variables is sufficient to explain the majority of variance in both height and ABD in this landscape.

Becknell, Justin M.↗

Quantifying climate–growth relationships at the stand level in a mature mixed‐species conifer forest

Abstract A range of environmental factors regulate tree growth; however, climate is generally thought to most strongly influence year‐to‐year variability in growth. Numerous dendrochronological (tree‐ring) studies have identified climate factors that influence year‐to‐year variability in growth for given tree species and location. However, traditional dendrochronology methods have limitations that prevent them from adequately assessing stand‐level (as opposed to species‐level) growth. We argue that stand‐level growth analyses provide a more meaningful assessment of forest response to climate fluctuations, as well as the management options that may be employed to sustain forest productivity. Working in a mature, mixed‐species stand at the Howland Research Forest of central Maine, USA , we used two alternatives to traditional dendrochronological analyses by (1) selecting trees for coring using a stratified (by size and species), random sampling method that ensures a representative sample of the stand, and (2) converting ring widths to biomass increments, which once summed, produced a representation of stand‐level growth, while maintaining species identities or canopy position if needed. We then tested the relative influence of seasonal climate variables on year‐to‐year variability in the biomass increment using generalized least squares regression, while accounting for temporal autocorrelation. Our results indicate that stand‐level growth responded most strongly to previous summer and current spring climate variables, resulting from a combination of individualistic climate responses occurring at the species‐ and canopy‐position level. Our climate models were better fit to stand‐level biomass increment than to species‐level or canopy‐position summaries. The relative growth responses (i.e., percent change) predicted from the most influential climate variables indicate stand‐level growth varies less from to year‐to‐year than species‐level or canopy‐position growth responses. By assessing stand‐level growth response to climate, we provide an alternative perspective on climate–growth relationships of forests, improving our understanding of forest growth dynamics under a fluctuating climate.

Teets, Aaron↗

Joint effects of climate, tree size, and year on annual tree growth derived from tree–ring records of ten globally distributed forests

Tree rings provide an invaluable long-term record for understanding how climate and other drivers shape tree growth and forest productivity. However, conventional tree-ring analysis methods were not designed to simultaneously test effects of climate, tree size, and other drivers on individual growth. This has limited the potential to test ecologically relevant hypotheses on tree growth sensitivity to environmental drivers and their interactions with tree size. Here, we develop and apply a new method to simultaneously model nonlinear effects of primary climate drivers, reconstructed tree diameter at breast height (DBH), and calendar year in generalized least squares models that account for the temporal autocorrelation inherent to each individual tree's growth. We analyze data from 3811 trees representing 40 species at 10 globally distributed sites, showing that precipitation, temperature, DBH, and calendar year have additively, and often interactively, influenced annual growth over the past 120 years. Growth responses were predominantly positive to precipitation (usually over ≥3-month seasonal windows) and negative to temperature (usually maximum temperature, over ≤3-month seasonal windows), with concave-down responses in 63% of relationships. Climate sensitivity commonly varied with DBH (45% of cases tested), with larger trees usually more sensitive. Trends in ring width at small DBH were linked to the light environment under which trees established, but basal area or biomass increments consistently reached maxima at intermediate DBH. Accounting for climate and DBH, growth rate declined over time for 92% of species in secondary or disturbed stands, whereas growth trends were mixed in older forests. These trends were largely attributable to stand dynamics as cohorts and stands age, which remain challenging to disentangle from global change drivers. By providing a parsimonious approach for characterizing multiple interacting drivers of tree growth, our method reveals a more complete picture of the factors influencing growth than has previously been possible.

generalized least squares (GLS)↗

Multiple Changepoint Detection for Non‐Gaussian Time Series

ABSTRACT This article combines methods from existing techniques to identify multiple changepoints in non‐Gaussian autocorrelated time series. A transformation is used to convert a Gaussian series into a non‐Gaussian series, enabling penalized likelihood methods to handle non‐Gaussian scenarios. When the marginal distribution of the data is continuous, the methods essentially reduce to the change of variables formula for probability densities. When the marginal distribution is count‐oriented, Hermite expansions and particle filtering techniques are used to quantify the scenario. Simulations demonstrating the efficacy of the methods are given and two data sets are analyzed: 1) the proportion of home runs hit by Major League Baseball batters from 1920 to 2023 and 2) a six‐dimensional series of tropical cyclone counts from the Earth's basins of generation from 1980 to 2023. In the first series, beta marginal distributions are used to describe the proportions; in the second, Poisson marginal distributions seem appropriate.

Lund, Robert [Department of Statistics University ↗

Imaging of electrical response of NiO x under controlled environment with sub-25-nm resolution

The spatially resolved electrical response of rf-sputtered polycrystalline NiO x films composed of 40 nm crystallites was investigated under different relative humidity levels (RH). The topological and electrical properties (surface potential and resistance) were characterized using Kelvin probe force microscopy (KPFM) and conductive scanning probe microscopy at 0%, 50%, and 80% relative humidity with sub 25nm resolution. The surface potential of NiO x decreased by about 180 mV and resistance decreased in a nonlinear fashion by about 2 G when relative humidity was increased from 0% to 80%. The dimensionality of surface features obtained through autocorrelation analysis of topological, surface potential and resistance maps increased linearly with increased relative humidity as water was adsorbed onto the film surface. Spatially resolved surface potential and resistance of the NiO x films were found to be heterogeneous, with distinct features that grew in size from about 60 nm to 175 nm between 0% and 80% RH levels, respectively. Here, we find that the changes in the heterogeneous character of the NiO films are consistent through the topological, surface potential, and resistance measurements, suggesting that the nanoscale surface potential and resistance properties converge with the mesoscale properties as water is adsorbed onto the NiO x film.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Quantum noise properties of CT images with anatomical textured backgrounds across reconstruction algorithms: FBP and SAFIRE

Purpose: Quantum noise properties of CT images are generally assessed using simple geometric phantoms with uniform backgrounds. Such phantoms may be inadequate when assessing nonlinear reconstruction or postprocessing algorithms. The purpose of this study was to design anatomically informed textured phantoms and use the phantoms to assess quantum noise properties across two clinically available reconstruction algorithms, filtered back projection (FBP) and sinogram affirmed iterative reconstruction (SAFIRE). Methods: Two phantoms were designed to represent lung and soft-tissue textures. The lung phantom included intricate vessel-like structures along with embedded nodules (spherical, lobulated, and spiculated). The soft tissue phantom was designed based on a three-dimensional clustered lumpy background with included low-contrast lesions (spherical and anthropomorphic). The phantoms were built using rapid prototyping (3D printing) technology and, along with a uniform phantom of similar size, were imaged on a Siemens SOMATOM Definition Flash CT scanner and reconstructed with FBP and SAFIRE. Fifty repeated acquisitions were acquired for each background type and noise was assessed by estimating pixel-value statistics, such as standard deviation (i.e., noise magnitude), autocorrelation, and noise power spectrum. Noise stationarity was also assessed by examining the spatial distribution of noise magnitude. The noise properties were compared across background types and between the two reconstruction algorithms. Results: In FBP and SAFIRE images, noise was globally nonstationary for all phantoms. In FBP images of all phantoms, and in SAFIRE images of the uniform phantom, noise appeared to be locally stationary (within a reasonably small region of interest). Noise was locally nonstationary in SAFIRE images of the textured phantoms with edge pixels showing higher noise magnitude compared to pixels in more homogenous regions. For pixels in uniform regions, noise magnitude was reduced by an average of 60% in SAFIRE images compared to FBP. However, for edge pixels, noise magnitude ranged from 20% higher to 40% lower in SAFIRE images compared to FBP. SAFIRE images of the lung phantom exhibited distinct regions with varying noise texture (i.e., noise autocorrelation/power spectra). Conclusions: Quantum noise properties observed in uniform phantoms may not be representative of those in actual patients for nonlinear reconstruction algorithms. Anatomical texture should be considered when evaluating the performance of CT systems that use such nonlinear algorithms.

60 APPLIED LIFE SCIENCES↗

TU-CD-BRB-09: Prediction of Chemo-Radiation Outcome for Rectal Cancer Based On Radiomics of Tumor Clinical Characteristics and Multi-Parametric MRI

Purpose: To evaluate the tumor clinical characteristics and quantitative multi-parametric MR imaging features for prediction of response to chemo-radiation treatment (CRT) in locally advanced rectal cancer (LARC). Methods: Forty-three consecutive patients (59.7±6.9 years, from 09/2013 – 06/2014) receiving neoadjuvant CRT followed by surgery were enrolled. All underwent MRI including anatomical T1/T2, Dynamic Contrast Enhanced (DCE)-MRI and Diffusion-Weighted MRI (DWI) prior to the treatment. A total of 151 quantitative features, including morphology/Gray Level Co-occurrence Matrix (GLCM) texture from T1/T2, enhancement kinetics and the voxelized distribution from DCE-MRI, apparent diffusion coefficient (ADC) from DWI, along with clinical information (carcinoembryonic antigen CEA level, TNM staging etc.), were extracted for each patient. Response groups were separated based on down-staging, good response and pathological complete response (pCR) status. Logistic regression analysis (LRA) was used to select the best predictors to classify different groups and the predictive performance were calculated using receiver operating characteristic (ROC) analysis. Results: Individual imaging category or clinical charateristics might yield certain level of power in assessing the response. However, the combined model outperformed than any category alone in prediction. With selected features as Volume, GLCM AutoCorrelation (T2), MaxEnhancementProbability (DCE-MRI), and MeanADC (DWI), the down-staging prediciton accuracy (area under the ROC curve, AUC) could be 0.95, better than individual tumor metrics with AUC from 0.53–0.85. While for the pCR prediction, the best set included CEA (clinical charateristics), Homogeneity (DCE-MRI) and MeanADC (DWI) with an AUC of 0.89, more favorable compared to conventional tumor metrics with an AUC ranging from 0.511–0.79. Conclusion: Through a systematic analysis of multi-parametric MR imaging features, we are able to build models with improved predictive value over conventional imaging or clinical metrics. This is encouraging, suggesting the wealth of imaging radiomics should be further explored to help tailor the treatment into the era of personalized medicine. This work is supported by the National Science Foundation of China (NSFC Grant No. 81201091), National High Technology Research and Development Program of China (863 program, Grant No. 2015AA020917), and Fund Project for Excellent Abroad Scholar Personnel in Science and Technology.

60 APPLIED LIFE SCIENCES↗