SEARCH · Search NASA

Results for “analyzable function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Anion Exchange Ionomers: Impact of Chemistry on Thin-Film Properties

Ionomer thin-films (i.e., 20–100 nm) on supports serve as model systems to understand ionomer-catalyst interfacial behavior as well as the confinement-driven deviation in properties from bulk membranes. While ionomer thin-films have been examined for proton exchange ionomers, the thin-film properties of anion exchange ionomers (AEIs) remain largely unexplored. More importantly, delineating the convoluted impact of chemistry and confinement on thin-film morphology and hydration is of interest to advancing the field on functional ionic interfaces. In this work, these aspects are studied by using AEIs of different backbones (perfluorinated, aliphatic, and aromatic) and side chains (various lengths, and single versus dual functional groups). Quartz-crystal microbalance and spectroscopic ellipsometry are used to analyze density and coupled with calculated free volume fraction of thin-films to provide insights on their gas transport properties. AEI side-chain's chemical character plays a key role in how confinement modulates hydration (in thin-film versus bulk). Overall, the results elucidate the effects of backbone, side-chain chemistry versus anion/cation type in the confinement-driven changes in thin-film morphology and swelling. This study also provides new insights for tuning AEI transport functionalities at interfaces via chemistry, which can benefit the design and development of electrode-ionomers for alkaline membrane-based energy systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

How Precisely Can Individual Molecules Be Analyzed? A Case Study on Locally Quantifying Forces and Energies Using Scanning Probe Microscopy

Recent advances in scanning probe microscopy methodology have enabled the measurement of tip-sample interactions with picometer accuracy in all three spatial dimensions, thereby providing a detailed site-specific and distance-dependent picture of the related properties. This paper explores the degree of detail and accuracy that can be achieved in locally quantifying probe-molecule interaction forces and energies for adsorbed molecules. Toward this end, cobalt phthalocyanine (CoPc), a promising CO 2 reduction catalyst, was studied on Ag(111) as a model system using low-temperature, ultrahigh vacuum noncontact atomic force microscopy. Data were recorded as a function of distance from the surface, from which detailed three-dimensional maps of the molecule's interaction with the tip for normal and lateral forces as well as the tip-molecule interaction potential were constructed. The data were collected with a CO molecule at the tip apex, which enabled a detailed visualization of the atomic structure. Determination of the tip-substrate interaction as a function of distance allowed isolation of the molecule-tip interactions; when analyzing these in terms of a Lennard-Jones-type potential, the atomically resolved equilibrium interaction energies between the CO tethered to the tip and the CoPc molecule could be recovered. Interaction energies peaked at less than 160 meV, indicating a physisorption interaction. As expected, the interaction was weakest at the aromatic hydrogens around the periphery of the molecule and strongest surrounding the metal center. In conclusion, the interaction, however, did not peak directly above the Co atom but rather in pockets surrounding it.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Towards functionalized lignin and its derivatives for high-value material applications

In pursuit of low-carbon development, the production of diverse functional materials from renewable resources, especially lignocellulosic biomass, is of vital importance. Lignin is a major component of lignocellulose, and it is nature's only true high-volume aromatic polymer, which has the advantages of biodegradability, biocompatibility, and low acquisition cost. After modification by a variety of physicochemical strategies, it can be applied as an alternative to basic industrial materials, including polyurethanes, phenolic, and epoxy resins, which stands comparably to conventional synthetic polymers in terms of both performance and reduced cost. To further unlock the potential of lignin, an alluring opportunity is the exploitation of lignin for high-value materials, especially nanomaterials, that find extensive applications in energy and environment areas. This review provides a systematic summary and perspective of research that has been devoted to lignin transformation to high-value materials, ranging from industrially well-established engineering materials to the recently emerged nanomaterials that were used for energy and environment applications. Latest cutting-edge innovations on lignin modification, controlled depolymerization, and assembly during the last five years are summarized. Structure-function relationships of lignin materials in terms of their specific applications are analyzed. Furthermore, challenges and future opportunities for lignin conversion to high-value materials are also provided. We wish this review will stimulate further advances in lignin based high-value materials, and promote “waste” into “wealth”.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Study of Self-Interaction Errors in Density Functional Calculations of Magnetic Exchange Coupling Constants Using Three Self-Interaction Correction Methods

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three non-empirical density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. Phys. Rev. B, 1981, 23, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. et al., J. Chem. Phys. 2006, 124, 094108], and the recent local scaling method [Zope, R. R. et al., J. Chem. Phys. 2019, 151, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H-He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Furthermore, analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Microstructure of amide-functionalized polyethylenes determined by NMR relaxometry

Amidation of polyethylenes creates a range of amide-containing materials with enhanced properties, but the effect of these functional groups on the microstructure of these new materials is not known. Here we employ solid-state nuclear magnetic resonance (NMR) techniques to analyze the microstructure of amide-modified polyethylenes. While a decrease in crystallinity was observed with increasing amounts of functionalization, we found by measuring the chain mobility of the crystalline, amorphous, and interphasial regions of the polyethylenes with NMR relaxation techniques that the grafted amidyl groups partition into the rigid amorphous fraction (RAF) between the crystalline and amorphous regions. The chemical specificity of these NMR experiments creates precise assessments of the location of functional groups within the materials. Together, these insights into the microstructure and morphology of amide-containing polyethylenes lay a foundation for a deeper understanding of the structure and properties of functional polyolefins.

Haber, Shira [Lawrence Berkeley National Laborator↗

DOE OSTI2022

Multilevel Analysis, Design, and Modeling of Coupling Advanced Nuclear Reactors and Thermal Energy Storage in an Integrated Energy System

This report discusses the different options for coupling thermal energy storage (TES) systems to advanced nuclear power plants (A-NPPs) in order to enable flexible and hybrid plant operation. An advanced light-water reactor (A LWR), a high-temperature gas-cooled reactor (HTGR) and a liquid-metal fast reactor (LMFR) were selected as the initial use cases for demonstrating a thermally balanced energy storage coupling design for thermal power extraction. The models presented herein showcase several design considerations, focusing on optimal deployment methodologies for achieving steady-state and transient-state operation with minimum disruption to the nuclear power cycle. This first part of the study presents steady-state models developed using Aspen HYSYS®, with the thermal energy bypass for NPP-TES coupling being varied at up to 50%. The various components were sized using the Aspen Process Economic Analyzer (APEA) and Aspen Exchanger Design and Rating (EDR), when applicable. Cost functions from these models were developed using the latest publicly available data obtained from APEA V11. The TES-coupled A-NPP steady-state models and cost functions then provided a baseline for dynamic operation and process optimization by using Idaho National Laboratory (INL)’s Framework for Optimization of Resources and Economics (FORCE) tools. A stochastic optimization of the various energy storage systems coupled to the A-NPPs was then performed using the Risk Analysis Virtual Environment (RAVEN) and its dispatch optimization plugin, the Holistic Energy Resource Optimization Network (HERON). The signal processing and synthetic history capabilities of RAVEN were used to account for the unpredictable behavior of electricity markets. An autoregressive moving average (ARMA) model was used to analyze price signals from the Pennsylvania-New Jersey-Maryland (PJM) market and were applied to the HERON analysis in order to optimize a system with the best economics. Transient modeling evaluation was then performed using Modelica models within the HYBRID repository, which was developed at INL for the Department of Energy Integrated Energy Systems program for the characterization of dynamic integrated system behavior and feedback. This includes evaluation of the TES-coupled A-LWR systems’ impact on physical and thermal system response during imposed system demands. Additional TES-coupled reactor types, coupling approaches, markets, and TES technologies will be evaluated in future work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

Computationally efficient and error aware surrogate construction for numerical solutions of subsurface flow through porous media

Limiting the injection rate to restrict the pressure below a threshold at a critical location can be an important goal of simulations that model the subsurface pressure between injection and extraction wells. The pressure is approximated by the solution of Darcy’s partial differential equation for a given permeability field. The subsurface permeability is modeled as a random field since it is known only up to statistical properties. This induces uncertainty in the computed pressure. Solving the partial differential equation for an ensemble of random permeability simulations enables estimating a probability distribution for the pressure at the critical location. These simulations are computationally expensive, and practitioners often need rapid online guidance for real-time pressure management. An ensemble of numerical partial differential equation solutions is used to construct a Gaussian process regression model that can quickly predict the pressure at the critical location as a function of the extraction rate and permeability realization. The Gaussian process surrogate analyzes the ensemble of numerical pressure solutions at the critical location as noisy observations of the true pressure solution, enabling robust inference using the conditional Gaussian process distribution. Our first novel contribution is to identify a sampling methodology for the random environment and matching kernel technology for which fitting the Gaussian process regression model scales as O ( n log n ) instead of the typical O ( n 3 ) rate in the number of samples n used to fit the surrogate. The surrogate model allows almost instantaneous predictions for the pressure at the critical location as a function of the extraction rate and permeability realization. Our second contribution is a novel algorithm to calibrate the uncertainty in the surrogate model to the discrepancy between the true pressure solution of Darcy’s equation and the numerical solution. Finally, although our method is derived for building a surrogate for the solution of Darcy’s equation with a random permeability field, the framework broadly applies to solutions of other partial differential equations with random coefficients.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Double-sided electron energy analyzer for measurement of non-Maxwellian electron energy distributions

Here, a double-sided electron energy analyzer is developed for studies of magnetic reconnection. It can measure electron energy distribution functions along two directions opposite to each other at the same time. Each side is composed of a floating reference grid, an energy selector grid, and a collector plate. The voltage of the selector grid is swept from –40 to 0 V with respect to the reference grid with a frequency of 1 MHz. This fast sweeping is required to resolve sub-Alfvénic changes in plasma quantities of the Magnetic Reconnection Experiment, where the typical Alfvénic time is a few microseconds. The reliability of the energy analyzer is checked in Maxwellian plasmas away from the reconnection region. In this case, the electron temperature computed from the electron energy distribution function agrees with measurements of a reference triple Langmuir probe. When it is located near the reconnection region, the temperatures of the tail electron population from both sides, facing into the electron flow and facing away from it, exceed the bulk electron temperature measured by the Langmuir probe by a factor of about 2.

47 OTHER INSTRUMENTATION↗

DOE OSTI2018

Elemental mapping and geochemical characterization of gas producing shales by laser induced breakdown spectroscopy

Laser induced breakdown spectroscopy (LIBS) has been used for the analysis of hydrocarbon bearing shale samples. Shale samples taken from a Marcellus gas well at depths ranging from 7498 to 7551 ft (2285.4–2301.5 m) were analyzed by LIBS using an 81 × 81 grid pattern covering an 8 × 8 mm area. The data collected from these experiments were used to construct 2D elemental maps for each sample including the hydrocarbon forming elements C and H. Results show that the spatial elemental composition of the shale varies due to the matrix of the rock, and as a function of the sample depth. The accuracy of the analysis was confirmed by analyzing a shale sample of known elemental concentrations and obtaining a good agreement between analyzed and reference values. It has been shown that LIBS can be used to determine elemental composition variations in hydrocarbon bearing shales in a laboratory setting. In conclusion, extending this capability into wellbores will enable producers to rapidly target resources with greater accuracy.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2022

Impact of acute temperature and air pollution exposures on adult lung function: A panel study of asthmatics

Background Individuals with respiratory conditions, such as asthma, are particularly susceptible to adverse health effects associated with higher levels of ambient air pollution and temperature. This study evaluates whether hourly levels of fine particulate matter (PM 2.5 ) and dry bulb globe temperature (DBGT) are associated with the lung function of adult participants with asthma. Methods and findings Global positioning system (GPS) location, respiratory function (measured as forced expiratory volume at 1 second (FEV 1 )), and self-reports of asthma medication usage and symptoms were collected as part of the Exposure, Location, and Lung Function (ELF) study. Hourly ambient PM 2.5 and DBGT exposures were estimated by integrating air quality and temperature public records with time-activity patterns using GPS coordinates for each participant (n = 35). The relationships between acute PM 2.5 , DBGT, rescue bronchodilator use, and lung function collected in one week periods and over two seasons (summer/winter) were analyzed by multivariate regression, using different exposure time frames. In separate models, increasing levels in PM 2.5 , but not DBGT, were associated with rescue bronchodilator use. Conversely DBGT, but not PM 2.5 , had a significant association with FEV 1 . When DBGT and PM 2.5 exposures were placed in the same model, the strongest association between cumulative PM 2.5 exposures and the use of rescue bronchodilator was identified at the 0–24 hours (OR = 1.030; 95% CI = 1.012–1.049; p-value = 0.001) and 0–48 hours (OR = 1.030; 95% CI = 1.013–1.057; p-value = 0.001) prior to lung function measure. Conversely, DBGT exposure at 0 hours (β = 3.257; SE = 0.879; p-value>0.001) and 0–6 hours (β = 2.885; SE = 0.903; p-value = 0.001) hours before a reading were associated with FEV 1 . No significant interactions between DBGT and PM 2.5 were observed for rescue bronchodilator use or FEV 1 . Conclusions Short-term increases in PM 2.5 were associated with increased rescue bronchodilator use, while DBGT was associated with higher lung function (i.e. FEV 1 ). Further studies are needed to continue to elucidate the mechanisms of acute exposure to PM 2.5 and DBGT on lung function in asthmatics.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2019

Inverse regression for ridge recovery: a data-driven approach for parameter reduction in computer experiments

Parameter reduction can enable otherwise infeasible design and uncertainty studies with modern computational science models that contain several input parameters. In statistical regression, techniques for sufficient dimension reduction (SDR) use data to reduce the predictor dimension of a regression problem. A computational scientist hoping to use SDR for parameter reduction encounters a problem: a computer prediction is best represented by a deterministic function of the inputs, so data comprised of computer simulation queries fail to satisfy the SDR assumptions. To address this problem, here we interpret SDR methods sliced inverse regression (SIR) and sliced average variance estimation (SAVE) as estimating the directions of a ridge function, which is a composition of a low-dimensional linear transformation with a nonlinear function. Within this interpretation, SIR and SAVE estimate matrices of integrals whose column spaces are contained in the ridge directions’ span; we analyze and numerically verify convergence of these column spaces as the number of computer model queries increases. Moreover, we show example functions that are not ridge functions but whose inverse conditional moment matrices are low-rank. Consequently, the computational scientist should beware when using SIR and SAVE for parameter reduction, since SIR and SAVE may mistakenly suggest that truly important directions are unimportant.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Online peak-aware energy scheduling with untrusted advice

This paper studies the online energy scheduling problem in a hybrid model where the cost of energy is proportional to both the volume and peak usage, and where energy can be either locally generated or drawn from the grid. Inspired by recent advances in online algorithms with Machine Learned (ML) advice, we develop parameterized deterministic and randomized algorithms for this problem such that the level of reliance on the advice can be adjusted by a trust parameter. We then analyze the performance of the proposed algorithms using two performance metrics: robustness that measures the competitive ratio as a function of the trust parameter when the advice is inaccurate, and consistency for competitive ratio when the advice is accurate. Since the competitive ratio is analyzed in two different regimes, we further investigate the Pareto optimality of the proposed algorithms. Our results show that the proposed deterministic algorithm is Pareto-optimal, in the sense that no other online deterministic algorithms can dominate the robustness and consistency of our algorithm. Furthermore, we show that the proposed randomized algorithm dominates the Pareto-optimal deterministic algorithm. Our large-scale empirical evaluations using real traces of energy demand, energy prices, and renewable energy generations highlight that the proposed algorithms outperform worst-case optimized algorithms and fully data-driven algorithms.

Lee, Russell↗

DOE OSTI2021

Determination of Carbon Functional Groups in Pyrolysis Bio-Oils using 13 C NMR: Laboratory Analytical Procedure (LAP)

Pyrolysis is a process that can be used to convert biomass to solid, liquid and gaseous products for use as renewable chemicals and fuels. The liquid fraction, known as “bio-oil” is complex and challenging to characterize, particularly by means of GC/MS, GPC, LC and FT-IR. NMR is capable of analyzing whole bio-oil samples and can provide quantitative results to characterize different functional groups or types of carbon present in bio-oil. This Laboratory Analytical Procedure (LAP) quantifies different carbon functional groups in whole bio-oil samples. This information can be used for comparisons between different pyrolysis experimental conditions or different upgrading processes and catalysts, and also allow for comparisons between bio-oils produced at different facilities.

09 BIOMASS FUELS↗

DOE OSTI2025

Bottomonium properties in the quark-gluon plasma from a lattice-QCD informed 𝑇-matrix approach

Recent computations of bottomonium correlation functions with extended sources in lattice-discretized quantum chromodynamics (lQCD) provide new insights into heavy-quark dynamics at distance scales which are of the order of the inverse temperature. We analyze these results employing the thermodynamic 𝑇-matrix framework, in a continued effort to interpret lQCD data for quarkonium correlation functions in a nonperturbative and self-consistently solved quantum-many-body approach to a strongly coupled quark-gluon plasma (QGP). Its key inputs are the in-medium driving kernel (potential) of the scattering equation and an interference function which implements 3-body effects in the quarkonium coupling to the thermal medium. A simultaneous description of lQCD results for the bottomonium correlators with extended operators and the previously analyzed Wilson line correlators only requires minor refinements of the potential but calls for stronger interference effects at larger separation of the bottom quark and antiquark. We then analyze the poles of the self-consistent 𝑇 matrices on the real axis to assess the survival of the various bound states. Here, we estimate the pertinent temperatures where the poles disappear for the various bottomonium states and discuss the relation to the corresponding peaks in the bottomonium spectral functions. We also recalculate the spatial diffusion coefficient of the QGP and find it to be similar to that in our previous study.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Comparative Density Functional Theory Study of Magnetic Exchange Couplings in Dinuclear Transition-Metal Complexes

Multicenter transition-metal complexes (MCTMs) with magnetically interacting ions have been proposed as components for information-processing devices and storage units. For any practical application of MCTMs as magnetic units, it is crucial to characterize their magnetic behavior, and in particular, the isotropic magnetic exchange coupling, J, between its magnetic centers. Due to the large size of typical MCTMs, density functional theory is the only practical electronic structure method for evaluating the J coupling. Here, we assess the accuracy of different density functional approximations for predicting the magnetic couplings of eight dinuclear transition-metal complexes, including five dimanganese, two dicopper, and one divanadium with known reliable experimental J couplings spanning from ferromagnetic to strong antiferromagnetic. The density functionals considered include global hybrid functionals which mix semilocal density functional approximations and exact exchange with a fixed admixing parameter, six local hybrid functionals where the admixing parameters are extended to be spatially dependent, the SCAN and r 2 SCAN meta-generalized gradient approximations (GGAs), and two widely used GGAs. We found that global hybrids tested in this work have a tendency to over-correct the error in magnetic coupling parameters from the Perdew–Burke–Ernzerhof (PBE) GGA as seen for manganese complexes. The performance of local hybrid density functionals shows no improvement in terms of bias and is scattered without a clear trend, suggesting that more efforts are needed for the extension from global to local hybrid density functionals for this particular property. The SCAN and r 2 SCAN meta-GGAs are found to perform as well as benchmark global hybrids on most tested complexes. We further analyze the charge density redistribution of meta-GGAs as well as global and local hybrid density functionals with respect to that of PBE, in connection to the self-interaction error or delocalization error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Characterization of nuclear breakup as a function of hard-scattering kinematics using dijets measured by ATLAS in p+Pb collisions

This Letter analyzes the sensitivity of event geometry estimators to the initial-state kinematics of hard scattering in proton–lead collisions. This analysis uses dijets as a proxy for the parton–parton scattering configuration, correlating it with event geometry estimators, namely the energy deposited in the Zero-Degree Calorimeter and the transverse energy recorded in the Forward Calorimeter in the Pb-going direction. The analysis uses data recorded by the ATLAS detector at the Large Hadron Collider with a nucleon–nucleon center-of-mass energy of 8.16 TeV, corresponding to an integrated luminosity of 56 nb − 1 . The jets are measured within the pseudorapidity interval −2.8 < η < 4.5, where positive η values correspond to the direction of the proton beam. Results are presented as a function of the Bjorken-x of the parton originating from the proton, xp . Both event geometry estimators are found to be dependent on xp , with the energy deposited in the Zero-Degree Calorimeter about six times less sensitive to xp compared with the transverse energy deposited in the Forward Calorimeter.

Aad, G↗

DOE OSTI2023

Multielectrode electrochemical cell for in situ structural characterization of amorphous thin-film catalysts using high-energy X-ray scattering

A multielectrode-based electrochemical cell allows the structural characterization of an amorphous thin-film water oxidation catalyst under various electrochemical potentials using high-energy X-ray scattering and atomic pair distribution function (PDF) techniques. A multielectrode with five electrodes provides a sufficiently low background signal to enable high-energy X-ray scattering (HEXS) measurements and amplifies the extremely low HEXS signals from samples for high-resolution PDF analysis of in situ data from thin-film catalysts. Glassy carbon (GC) creates a relatively low intensity HEXS pattern and is used as a working electrode. Instead of a three-dimensional (3D) porous electrode architecture, the flat geometry of the electrode enables various deposition techniques to be used for the preparation of a highly conductive metal oxide layer. PDF analysis demonstrates high spatial resolution for a 230 nm thick amorphous iridium oxide film deposited on two roughened 60 µm thick GC electrodes. The PDF analysis resolves the domain size and distinguishes changes in fine structure which are directly correlated with the structure and function of the catalysts. In conclusion, the results bring the opportunity to analyze the structure of nanometre-scale amorphous thin-film catalysts in an electrolyte-compatible and compact 3D-printed electrochemical cell in a three-electrode configuration.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Genome-Wide Analysis of Terpene Synthase Gene Family in Mentha longifolia and Catalytic Activity Analysis of a Single Terpene Synthase

Terpenoids are a wide variety of natural products and terpene synthase (TPS) plays a key role in the biosynthesis of terpenoids. Mentha plants are rich in essential oils, whose main components are terpenoids, and their biosynthetic pathways have been basically elucidated. However, there is a lack of systematic identification and study of TPS in Mentha plants. In this work, we genome-widely identified and analyzed the TPS gene family in Mentha longifolia, a model plant for functional genomic research in the genus Mentha. A total of 63 TPS genes were identified in the M. longifolia genome sequence assembly, which could be divided into six subfamilies. The TPS-b subfamily had the largest number of genes, which might be related to the abundant monoterpenoids in Mentha plants. The TPS-e subfamily had 18 members and showed a significant species-specific expansion compared with other sequenced Lamiaceae plant species. The 63 TPS genes could be mapped to nine scaffolds of the M. longifolia genome sequence assembly and the distribution of these genes is uneven. Tandem duplicates and fragment duplicates contributed greatly to the increase in the number of TPS genes in M. longifolia. The conserved motifs (RR(X)8W, NSE/DTE, RXR, and DDXXD) were analyzed in M. longifolia TPSs, and significant differentiation was found between different subfamilies. Adaptive evolution analysis showed that M. longifolia TPSs were subjected to purifying selection after the species-specific expansion, and some amino acid residues under positive selection were identified. Furthermore, we also cloned and analyzed the catalytic activity of a single terpene synthase, MlongTPS29, which belongs to the TPS-b subfamily. MlongTPS29 could encode a limonene synthase and catalyze the biosynthesis of limonene, an important precursor of essential oils from the genus Mentha. This study provides useful information for the biosynthesis of terpenoids in the genus Mentha.

59 BASIC BIOLOGICAL SCIENCES↗