Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Concurrent effect”

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 289 records · Page 16

Tidal Circularization of Binaries by Resonance Locking. I. The Importance of the Pre-main Sequence

Although tidal dissipation in binary stars has been studied for over a century, theoretical predictions have yet to match the observed properties of binary populations. This work quantitatively examines the recent proposal of tidal circularization by resonance locking, where tidal dissipation arises from resonances between the star’s natural oscillation frequencies and harmonics of the orbital frequency, and where resonances are “locked” for an extended period of time due to concurrent stellar evolution. We focus on tidal resonances with axi-symmetric gravity modes, and examine binaries with primary masses from one to two solar masses. We find that orbital evolution via resonance locking occurs primarily during the star’s pre-main-sequence phase, with the main-sequence phase contributing negligibly. Resonance locking, ignoring nonlinearity, can circularize binaries with peri-center distances out to ∼10 stellar radii, corresponding to circular periods of ∼4–6 days. However, we find resonantly excited gravity modes will become nonlinear in stellar cores, which prevents them from reaching their full, linear amplitudes. We estimate that such a “saturated resonance lock” reduces the circularization period by about a third, but resonance locking remains much more effective than the cumulative actions of equilibrium tides. In a companion paper, we examine recent binary data to compare against theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Pathways for decarbonization of the buildings sector in Ukraine

The paper focuses on Ukraine’s intention to achieve a two-thirds reduction in buildings’ energy consumption for heating and cooling by 2050, concurrently aiming for net zero greenhouse gas emissions and heightened energy security. Here, the study examines the outcomes of retrofitting existing residential, commercial, and public buildings with highly efficient materials, improving construction standards, and transitioning to advanced heating systems. However, Russia’s invasion in 2022 inflicted substantial damage, prompting a shift from retrofit and decarbonization to reconstruction. The Ukrainian government’s Reconstruction Plan emphasizes clean, sustainable, and resilient energy systems. The study employs energy system and integrated assessment models (TIMES-Ukraine and GCAM-Ukraine) to explore scenarios taking into consideration the war, reconstruction, and a net zero CO 2 pathway. Using two models allowed the inter-model comparison. The analysis addresses vital questions on energy resiliency measures and the compounding effects of decarbonization. Findings indicate that Ukraine’s energy goals can be met through strategic retrofitting and economy-wide decarbonization, emphasizing the importance of low-carbon alternatives like district heating with renewable sources. Electrification with renewables and fuel-switching emerges as crucial for achieving building decarbonization. The study offers valuable insights into navigating energy challenges amidst the war and outlines a pathway for Ukraine’s sustainable energy future.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Soil carbon stocks in temperate grasslands differ strongly across sites but are insensitive to decade-long fertilization

We report enhancing soil carbon (C) storage has the potential to offset human-caused increases in atmospheric CO 2 . Rising CO 2 has occurred concurrently with increasing supply rates of biologically limiting nutrients such as nitrogen (N) and phosphorus (P). However, it is unclear how increased supplies of N and P will alter soil C sequestration, particularly in grasslands, which make up nearly a third of non-agricultural land world-wide. Here, we leverage a globally distributed nutrient addition experiment (the Nutrient Network) to examine how a decade of N and P fertilization (alone and in combination) influenced soil C and N stocks at nine grassland sites spanning the continental United States. We measured changes in bulk soil C and N stocks and in three soil C fractions (light and heavy particulate organic matter, and mineral-associated organic matter fractions). Nutrient amendment had variable effects on soil C and N pools that ranged from strongly positive to strongly negative, while soil C and N pool sizes varied by more than an order of magnitude across sites. Piecewise SEM clarified that small increases in plant C inputs with fertilization did not translate to greater soil C storage. Nevertheless, peak season aboveground plant biomass (but not root biomass or production) was strongly positively related to soil C storage at seven of the nine sites, and across all nine sites, soil C covaried with moisture index and soil mineralogy, regardless of fertilization. Overall, we show that site factors such as moisture index, plant productivity, soil texture, and mineralogy were key predictors of cross-site soil C, while nutrient amendment had weaker and site-specific effects on C sequestration. This suggests that prioritizing the protection of highly productive temperate grasslands is critical for reducing future greenhouse gas losses arising from land use change.

59 BASIC BIOLOGICAL SCIENCES↗

Single nucleotide variants drive evolutionary phage-host arms race in anaerobic carbon dioxide-converting microbiome

Microbial bioconversions are shaped by environmental perturbations and the adaptation of resident microbiomes. Prokaryotes coexist with bacteriophages, yet their coevolutionary trajectories remain underexplored. Here, we investigate the effects of a cultivation vessel leak on an anaerobic consortium performing carbon dioxide reduction. Using time-series shotgun metagenomic sequencing, we reconstruct microbial and viral genomes to track community shifts. We further apply single-nucleotide variant profiling and CRISPR array analysis to monitor viral microdiversity and host defense mechanisms. After bioaugmentation restores bioconversion efficiency, the consortium undergoes pronounced restructuring, with new dominant taxa emerging from the rare biosphere. We identify patterns consistent with phage predation selectively removing certain species, while others exhibit resilience to infection. This shift aligns with a widespread viral outbreak and a transient increased frequency of single nucleotide variants in bacterial CRISPR–Cas defense genes. Expansion of CRISPR spacers further supports that CRISPR-mediated processes influence microbial resilience. Concurrently, phages infecting resilient hosts exhibited adaptive evolution, marked by high genetic heterogeneity. Selective pressure varies across their genomes, targeting infectivity genes and protospacer-adjacent motifs. These findings highlight a dynamic evolutionary arms race driven by the selection of beneficial genetic variants, providing a mechanistic framework for multi-omics investigations, and informing biotechnological applications, including phage-based microbiome manipulation.

Ghiotto, G↗

Concurrent probing of electron-lattice dephasing induced by photoexcitation in 1 T -TaSeTe using ultrafast electron diffraction

It has been technically challenging to concurrently probe the electrons and the lattices in materials during nonequilibrium processes, allowing their correlations to be determined. Here, in a single set of ultrafast electron diffraction patterns taken on the charge-density-wave (CDW) material 1T-TaSeTe, we discover a temporal shift in the diffraction intensity measurements as a function of scattering angle. With the help of dynamic models and theoretical calculations, we show that the ultrafast electrons probe both the valence-electron and lattice dynamic processes, resulting in the temporal shift measurements. Our results demonstrate unambiguously that the CDW is not merely a result of the periodic lattice deformation ever present in 1T-TaSeTe but has significant electronic origin. This method demonstrates an approach for studying many quantum effects that arise from electron-lattice dephasing in molecules and crystals for next-generation devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Correlating Crystallography, Magnetism, and Electronic Structure Across Anhysteretic First-Order Phase Transition in Pr 2 In

Here, temperature-dependent powder X-ray diffraction and magnetization measurements of Pr 2 In conclusively prove that the unusual anhysteretic first-order paramagnetic-ferromagnetic phase transition in the compound is related to concurrent changes in both the magnetic and crystallographic lattices. At the same time, the hexagonal Ni 2 In-type structure is stable at least between 6 and 298 K, including at T C = ~57 K. From the density functional theory calculations, the electronic structure of the compound is extraordinarily sensitive to minor changes in lattice parameters that occur across the phase transition, revealing the origin of strong magnetoelastic coupling. In the vicinity of T C , the maximum entropy change, ΔS Max = -16 J Kg -1 K -1 induced by a moderate magnetic field change of 20 kOe (ΔS Max = -20 J Kg -1 K -1 for 50 kOe magnetic field change) is comparable to other known potentially functional materials demonstrating large cryogenic magnetocaloric effect.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

ECP ST Capability Assessment Report (CAR) for VTK-m (FY20)

The ECP/VTK-m project is providing the core capabilities to perform scientific visualization on Exascale architectures. The ECP/VTK-m project fills the critical feature gap of performing visualization and analysis on processors like graphics-based processors. The results of this project will be delivered in tools like ParaView, Vislt, and Ascent as well as in stand-alone form. Moreover, these projects are depending on this ECP effort to be able to make effective use of ECP architectures. One of the biggest recent changes in high-performance computing is the increasing use of accelerators. Accelerators contain processing cores that independently are inferior to a core in a typical CPU, but these cores are replicated and grouped such that their aggregate execution provides a very high computation rate at a much lower power. Current and future CPU processors also require much more explicit parallelism. Each successive version of the hardware packs more cores into each processor, and technologies like hyper threading and vector operations require even more parallel processing to leverage each core's full potential. VTK-m is a toolkit of scientific visualization algorithms for emerging processor architectures. VTK-m supports the fine-grained concurrency for data analysis and visualization algorithms required to drive extreme scale computing by providing abstract models for data and execution that can be applied to a variety of algorithms across many different processor architectures. The ECP/VTK-m project is building up the VTK-m codebase with the necessary visualization algorithm implementations that run across the varied hardware platforms to be leveraged at the Exascale. We will be working with other ECP projects, such as ALPINE, to integrate the new VTK-m code into production software to enable visualization on our HPC systems.

97 MATHEMATICS AND COMPUTING↗

Layer-Parallel Training of Deep Residual Neural Networks

Residual neural networks (ResNets) are a promising class of deep neural networks that have shown excellent performance for a number of learning tasks, e.g., image classification and recognition. Mathematically, ResNet architectures can be interpreted as forward Euler discretizations of a nonlinear initial value problem whose time-dependent control variables represent the weights of the neural network. Hence, training a ResNet can be cast as an optimal control problem of the associated dynamical system. For similar time-dependent optimal control problems arising in engineering applications, parallel-in-time methods have shown notable improvements in scalability. This paper demonstrates the use of those techniques for efficient and effective training of ResNets. The proposed algorithms replace the classical (sequential) forward and backward propagation through the network layers with a parallel nonlinear multigrid iteration applied to the layer domain. This adds a new dimension of parallelism across layers that is attractive when training very deep networks. From this basic idea, we derive multiple layer-parallel methods. The most efficient version employs a simultaneous optimization approach where updates to the network parameters are based on inexact gradient information in order to speed up the training process. Finally, using numerical examples from supervised classification, we demonstrate that the new approach achieves a training performance similar to that of traditional methods, but enables layer-parallelism and thus provides speedup over layer-serial methods through greater concurrency.

97 MATHEMATICS AND COMPUTING↗

Micro-beam bending of FCC bicrystals: A comparison between defect dynamics simulations and experiments

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. In conclusion, simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

36 MATERIALS SCIENCE↗

Fast redox switches lead to rapid transformation of goethite in humid tropical soils: A Mössbauer spectroscopy study

Abstract Humid tropical forest soils experience frequent rainfall, which limits oxygen diffusion and creates redox heterogeneity in upland soils. In this study we gauged the effect of short‐term anoxic conditions on Fe mineralogy of relatively Fe‐ and C‐rich surface soils (C/Fe mole ratio ∼5) from a humid tropical forest in Puerto Rico. Soils subjected to 4‐d oxic/anoxic oscillation were characterized by selective chemical extractions, Mössbauer spectroscopy (MBS), and X‐ray diffraction. Chemical extraction data suggested that rapidly switching redox conditions had subtle effects on bulk Fe mineralogy. Mössbauer, on the other hand, indicated that (a) the soil Fe is a mixture of goethites of varying characteristics with minor contributions from ferrihydrite (<5%) and Fe(III)‐organic matter (OM) phases (∼10%), and (b) anoxic conditions rapidly transformed all forms of goethite to relatively stable Fe oxides. Such fast changes in goethite features could be due to rapid depletion and sorption of bio‐released structural Al or sorbed OM onto residual soil components. The rapid temporal changes in MBS parameters and corresponding pore water nominal oxidation state of C values suggest that Fe‐C transformations in these upland tropical soils are complex and intricately coupled. A comprehensive understanding of the fate of Fe, Al, and OM (Fe‐organic moieties) during redox switches and concurrent changes in pore water chemistry is critical for development of robust transport models in humid tropical soils, which are subject to episodic low‐redox events.

54 ENVIRONMENTAL SCIENCES↗

Climate Challenges and Nonproliferation: Addressing the Issue through Technical Cooperation

Climate change is an urgent global challenge that is already severely impacting many regions of the world. The 2015 Paris Agreement was a landmark achievement, while the 2023 UN COP 28 Climate Conference in Dubai recognized nuclear energy and its applications as a proven and sustainable means to help societies adapt to climate change and take measures to mitigate and, possibly, reverse its effects. Twenty-two countries, supported by over 120 companies, pledged to triple the share of global nuclear energy generation by 2050. The global interest in nuclear applications for peaceful uses has never been higher. The deployment of advanced reactors around the globe – focused primarily on lowering the carbon footprint and allowing for socio-economic development – must be addressed responsibly through setting priorities for its implementation. Tackling the challenges that new technologies, such as advanced reactors, bring to the nonproliferation regime is at the top of the list. The regime stands on the three pillars of the Nuclear Nonproliferation Treaty (NPT): nonproliferation, peaceful uses of nuclear energy, and disarmament. Recognizing the inherent risks of expanding applications of nuclear materials and technologies for peaceful uses, all such efforts must concurrently strengthen the nonproliferation norm enshrined in the NPT. The challenges of adapting to and mitigating climate change with the use of advanced reactors is already impacting the discussions and expectations about the future of the NPT. The dialog with non-traditional domestic, and international partners on a strategic technical cooperation approach is key to proposing and implementing scientific and technical solutions for urgent climate issues, while framing the discussion within nonproliferation requirements and commitments, and demonstrating the underlying value of the NPT in support of peaceful nuclear applications. This paper addresses adaptation to climate challenges and mitigation of them through advanced reactors and establishes nonproliferation linkages that derive from the deployment of this technology. This paper presents an assessment of the benefits of mechanisms for technical cooperation and peaceful uses of nuclear technology in the framework of Article IV of the NPT.

Prah, Christina↗

Effects of secondary void-initiating particles on the steady-state crack growth resistance of high-strength steel

In this investigation, we focus on the dependence of stable crack growth resistance, as measured by the short-rod fracture toughness test, on secondary void-initiating particles such as AlN and Ti-based grain-refining precipitates in a variety of high-strength steels with tempered martensitic microstructures. The model developed by Ritchie and Thompson is modified to illustrate the significant amounts of toughening that can result from the refinement of secondary particles. Analysis of the data suggests that material strength is a predominant factor in increasing the short-rod fracture toughness relative to linear-elastic measures of initiation fracture toughness, but the extent of toughening is limited by the size and number density of secondary particles in the microstructure. The variation in estimates of secondary microvoid initiation and growth strains with precipitate size reinforce the notion that primary fracture at non-metallic inclusions and secondary fracture at smaller particles occur as sequential processes with a degree of concurrence that is dependent on the state of precipitation in both particle dispersions. Toughening in this connection is maximized by increases in microvoid growth strain that result from decreases in the size and areal number density of secondary void-initiating particles. Finally, the occurrence of transient instabilities during crack extension in short-rod specimens is explained with a phenomenological model that relates crack growth stability to natural variations in the dispersion of secondary void-initiating particles in the microstructure.

42 ENGINEERING↗

Bacterial community dynamics as a result of growth-yield trade-off and multispecies metabolic interactions toward understanding the gut biofilm niche

Abstract Bacterial communities are ubiquitous, found in natural ecosystems, such as soil, and within living organisms, like the human microbiome. The dynamics of these communities in diverse environments depend on factors such as spatial features of the microbial niche, biochemical kinetics, and interactions among bacteria. Moreover, in many systems, bacterial communities are influenced by multiple physical mechanisms, such as mass transport and detachment forces. One example is gut mucosal communities, where dense, closely packed communities develop under the concurrent influence of nutrient transport from the lumen and fluid-mediated detachment of bacteria. In this study, we model a mucosal niche through a coupled agent-based and finite-volume modeling approach. This methodology enables us to model bacterial interactions affected by nutrient release from various sources while adjusting individual bacterial kinetics. We explored how the dispersion and abundance of bacteria are influenced by biochemical kinetics in different types of metabolic interactions, with a particular focus on the trade-off between growth rate and yield. Our findings demonstrate that in competitive scenarios, higher growth rates result in a larger share of the niche space. In contrast, growth yield plays a critical role in neutralism, commensalism, and mutualism interactions. When bacteria are introduced sequentially, they cause distinct spatiotemporal effects, such as deeper niche colonization in commensalism and mutualism scenarios driven by species intermixing effects, which are enhanced by high growth yields. Moreover, sub-ecosystem interactions dictate the dynamics of three-species communities, sometimes yielding unexpected outcomes. Competitive, fast-growing bacteria demonstrate robust colonization abilities, yet they face challenges in displacing established mutualistic systems. Bacteria that develop a cooperative relationship with existing species typically obtain niche residence, regardless of their growth rates, although higher growth yields significantly enhance their abundance. Our results underscore the importance of bacterial niche dynamics in shaping community properties and succession, highlighting a new approach to manipulating microbial systems.

Microbiology↗

funcX: Federated Function as a Service for Science

Here, funcX is a distributed function as a service (FaaS) platform that enables flexible, scalable, and high performance remote function execution. Unlike centralized FaaS systems, funcX decouples the cloud-hosted management functionality from the edge-hosted execution functionality. funcX's endpoint software can be deployed, by users or administrators, on arbitrary laptops, clouds, clusters, and supercomputers, in effect turning them into function serving systems. funcX's cloud-hosted service provides a single location for registering, sharing, and managing both functions and endpoints. It allows for transparent, secure, and reliable function execution across the federated ecosystem of endpoints-enabling users to route functions to endpoints based on specific needs. funcX uses containers (e.g., Docker, Singularity, and Shifter) to provide common execution environments across endpoints. funcX implements various container management strategies to execute functions with high performance and efficiency on diverse funcX endpoints. funcX also integrates with an in-memory data store and Globus for managing data that may span endpoints. We motivate the need for funcX, present our prototype design and implementation, and demonstrate, via experiments on two supercomputers, that funcX can scale to more than 130000 concurrent workers. We show that funcX's container warming-aware routing algorithm can reduce the completion time for 3,000 functions by up to 61% compared to a randomized algorithm and the in-memory data store can speed up data transfers by up to 3x compared to a shared file system.

97 MATHEMATICS AND COMPUTING↗

Online task-space motion control for positioner-coordinated multi-robot manufacturing systems

Incorporating multiple robotic manipulators into large-scale manufacturing systems enhances production efficiency and expands manufacturing capabilities beyond those of single-robot systems. Workpiece positioners in robotic manufacturing have demonstrated significant benefits for process optimization, but coordination strategies for multi-robot systems with shared positioners have received limited attention. This work presents a task-space coordinated trajectory-tracking control framework for multi-robot manufacturing systems, in which robots coordinate their motions within a shared, dynamic workpiece positioning frame. A workpiece positioner actively adjusts the pose of the manufactured component to enable greater operational concurrency and improve overall production efficiency. The proposed motion-coordination scheme employs a distributed and scalable architecture, supporting coordination across heterogeneous multi-robot systems. Two optimization methodologies are introduced to manage kinematic redundancies and maintain continuous, near-optimal operation throughout the manufacturing process. The first strategy exploits a task-space dimensionality reduction to achieve locally optimal configurations by leveraging symmetry-axis rotations of the tool. The second strategy utilizes the workpiece positioner to drive the coordinated robots toward stable and kinematically favorable configurations. For both optimization strategies, multiple objectives are defined to improve key performance metrics, including manipulability, configuration consistency, proximity to mechanical limits, and motion efficiency. Addressing a key limitation of existing coordination approaches, the framework is designed around online setpoint modification, allowing coordinated robots to respond effectively to in-situ process feedback. The proposed control framework is validated using the Robot Operating System (ROS) middleware on a combination of physical and simulated multi-robot system hardware.

Arbogast, Alex [ORNL] (ORCID:0000000154740723)↗

Sheath transitions in a cylindrical filament discharge: Axisymmetric 1D3V PIC-MCC simulations

We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.

Electrical hysteresis↗

Distributed optimization for multi-commodity urban traffic control

A distributed method for concurrent traffic signal and routing control of traffic networks is proposed. The method is based on the multi-commodity store-and-forward model, in which the destinations are the commodities. The system benefits from the communication between vehicles and infrastructure, providing optimal signal timings to intersections and routes to vehicles on a link-by-link basis. Using the augmented Lagrangian to model the constraints into the objective, the baseline centralized problem is decomposed into a set of objective-coupled subproblems, one for each intersection, enabling the solution to be computed by a distributed- gradient projection algorithm. Further, the intersection agents only need to communicate and coordinate with neighboring intersections to ensure convergence to the optimal solution while tolerating suboptimal iterations that offer more flexibility, unlike other distributed approaches. Through microsimulation, we demonstrate the effectiveness of the proposed algorithm in traffic networks with time-varying demand. Computational analysis shows that the distributed problem is suitable for real-time applications. A robustness analysis show that the distributed formulation enables a graceful degradation of the system in case of failure.

Augmented Lagrangian↗

Multiplex knockout of trichome-regulating MYB duplicates in hybrid poplar using a single gRNA

As the focus for CRISPR/Cas-edited plants moves from proof-of-concept to real-world applications, precise gene manipulation will increasingly require concurrent multiplex editing for polygenic traits. A common approach for editing across multiple sites is to design one guide RNA (gRNA) per target; however, this complicates construct assembly and increases the possibility of off-target mutations. In this study, we utilized one gRNA to target MYB186, a known positive trichome regulator, as well as its paralogs MYB138 and MYB38 at a consensus site for mutagenesis in hybrid poplar (Populus tremula × P. alba INRA 717-1B4). Unexpected duplications of MYB186 and MYB138 resulted in eight alleles for the three targeted genes in the hybrid poplar. Deep sequencing and polymerase chain reaction analyses confirmed editing across all eight targets in nearly all of the resultant glabrous mutants, ranging from small indels to large genomic dropouts, with no off-target activity detected at four potential sites. This highlights the effectiveness of a single gRNA targeting conserved exonic regions for multiplex editing. Additionally, cuticular wax and whole-leaf analyses showed a complete absence of triterpenes in the trichomeless mutants, hinting at a previously undescribed role for the nonglandular trichomes of poplar.

59 BASIC BIOLOGICAL SCIENCES↗