Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting methods”

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 181 records · Page 10

A combined experimental and computational analysis of failure mechanisms in open-hole cross-ply laminates under flexural loading

In this work, integrated experimental tests and computational modeling are proposed to investigate the failure mechanisms of open-hole cross-ply carbon fiber reinforced polymer (CFRP) laminated composites. In particular, we propose two effective methods, which include width-tapered double cantilever beam (WTDCB) and fixed-ratio mixed-mode end load split (FRMMELS) tests, to obtain the experimental data more reliably. We then calibrate the traction-separation laws of cohesive zone model (CZM) used among laminas of the composites by leveraging these two methods. The experimental results of fracture energy, i.e. G Ic and G Tc , obtained from WTDCB and FRMMELS tests are generally insensitive to the crack length thus requiring no effort to accurately measure the crack tip. Moreover, FRMMELS sample contains a fixed mixed-mode ratio of G IIc /G Tc depending on the width taper ratio. Examining comparisons between experimental results of FRMMELS tests and failure surface of B–K failure criterion predicted from a curve fitting, good agreement between the predictions and experimental data has been found, indicating that FRMMELS tests are an effective method to determine mixed-mode fracture criterion. In addition, a coupled experimental-computational modeling of WTDCB, edge notched flexure, and FRMMELS tests are adopted to calibrate and validate the interfacial strengths. Finally, failure mechanisms of open-hole cross-ply CFRP laminates under flexural loading have been studied systematically using experimental and multi-scale computational analyses based on the developed CZM model. The initiation and propagation of delamination, the failure of laminated layers as well as load-displacement curves predicted from computational analyses are in good agreement with what we have observed experimentally.

36 MATERIALS SCIENCE↗

Electronic structure, magnetic properties, and exchange splitting of gadolinium intermetallics

Here, we report here the electronic structure calculations of the magnetic moment, magneto-crystalline anisotropy energy (MAE), and density of states (DOS) of Gd-metal, GdAl 2 , GdAl 3 , GdCo 2 , GdCo 3 , and GdCo 5 . From the DOS, the exchange splitting energy is derived to map ferrimagnetic to paramagnetic phase transition (Curie temperature, T C ). The employed theoretical method here is based on advanced density functional theory, specifically the linearized augmented plane wave method (LAPW) within the local spin density approximation (LSDA), including spin–orbit coupling (SOC) and Hubbard model parameter (U). DOS analysis of both Gd-Al and Gd-Co systems reveals a strong correlation between exchange splitting energy and T C . The Gd-Co system exhibits an enormous exchange splitting energy among Co-3d electrons, vastly overwhelming the vanishing Gd 5d exchange splitting. Finally, these findings suggest that exchange splitting of Co-3d electrons plays a key role in the ground state magnetism and magnetic transition temperatures in these materials.

36 MATERIALS SCIENCE↗

scPlantAnnotate: an accurate and robust transformer-based model for plant cell type annotation

Accurate cell type annotation remains a major bottleneck in plant single-cell RNA sequencing (scRNA-seq), where existing tools are often adapted from animal studies and perform sub-optimally on plant data. The lack of plant-specific computational frameworks limits the construction of plant cell atlases and downstream biological discovery. We develop and evaluate scPlantAnnotate, a Transformer-based reference annotation framework tailored for plant scRNA-seq data, and benchmark it against state-of-the-art deep learning and conventional methods across multiple plant species. Species-specific scPlantAnnotate models were trained using curated datasets from Arabidopsis thaliana, Zea mays, Oryza sativa, and Glycine max. We compared scPlantAnnotate with leading baselines under both standard random-split evaluation and a more stringent leave-one-dataset-out setting, which tests robustness to completely unseen datasets and tissue types. scPlantAnnotate consistently outperforms existing approaches across all four species under random-split evaluation. In the leave-one-dataset-out setting for A. thaliana, where performance drops markedly for all methods due to strong batch effects and dataset heterogeneity, scPlantAnnotate nonetheless achieves the highest Accuracy, Macro-F1, Balanced Accuracy, and Macro-AUROC on average and ranks first on most held-out datasets. These results demonstrate improved robustness to dataset shifts, a critical yet underexplored challenge in plant scRNA-seq analysis. A freely accessible web server enables users to annotate their own datasets using pretrained models. scPlantAnnotate provides a plant-specific, Transformer-based framework for single-cell annotation that delivers state-of-the-art performance and enhanced robustness to unseen datasets. By addressing limitations of existing tools and enabling scalable reference-based annotation, scPlantAnnotate supports the development of comprehensive plant cell atlases and facilitates broader use of single-cell genomics in plant biology.

Bioinformatics↗

Toward non-Si electronics: From remote epitaxy to layer splitting of 2D materials for mixed dimensional heterostructures

The current electronics has been mainly dominated by Si-based devices due to their mature processing system and exceptional cost-effectiveness. However, next generation electronics needs novel functionalities that cannot be realized by Si because of intrinsic limitation of Si. Accordingly, demand for non-Si electronics has been getting substantially high. Unfortunately, current methodology requires extremely high cost for non-Si materials, which impedes the progress in developing the non-Si based electronics. Here, I will discuss about our group’s efforts to address this issue. Our team recently conceived a new crystalline growth, termed as “remote epitaxy”, which can copy/paste crystalline information from substrates remotely through graphene, thus generating single-crystalline films on graphene. As interfacial binding energy is attenuated by inserting graphene at interface, the single-crystalline films can be easily exfoliated from the slippery graphene surface. Also, the graphene-coated substrates can be, in principle, reused infinitely to produce single-crystalline films. Thus, the remote epitaxy can produce non-Si semiconductor films with unprecedented cost efficiency while allowing additional flexible device functionality required for current ubiquitous electronics. Next, I will discuss about a layer splitting technique which can be a potential solution to overcome the problem in obtaining large-scale and monolayer 2D materials. A 2D material-based heterostructure has been intensively studied because of its unique device functionalities and novel physics. However, it is extremely challenging to secure large-scale and monolayer 2D materials because of following issues: 1) poor scalability for laboratory fabrication processes of 2D heterostructures and 2) lack of well-defined control parameters for kinetics of 2D materials and predictable number of layers of 2D materials. To resolve this issue, we conceived a new approach called “layer-resolved splitting” which obtains multiple monolayer from multilayer 2D materials by controlling interfacial toughness contrast. As this method is versatile and universal, we can, in principle, apply to all 2D materials. We succeeded in having large-scale, monolayer 2D materials through our approach and, thereby 2D heterostructures were demonstrated for functional devices. Lastly, I would like to discuss opportunities of mixed-dimensional heterostructure demonstrated by remote epitaxy and layer-resolved splitting. As they produce freestanding 3D bulk films and 2D atomic layers, a new type of 3D/2D heterostructures can be realized where a new physics and new device architecture are revealed. Therefore, I believe that a new opportunity will be discovered through the mixed-dimensional heterostructures.

Bae, Sang-Hoon↗

Prediction of Feasibility of Polaronic OER on (110) Surface of Rutile TiO 2

The polaronic effects at the atomic level hold paramount significance for advancing the efficacy of transition metal oxides in applications pertinent to renewable energy. The lattice–distortion mediated localization of photoexcited carriers in the form of polarons plays a pivotal role in the photocatalysis. This investigation focuses on rutile TiO 2 , an important material extensively explored for solar energy conversion in artificial photosynthesis, specifically targeting the generation of green H 2 through photoelectrochemical (PEC) H 2 O splitting. By employing Hubbard-U corrected and hybrid density functional theory (DFT) methods, we systematically probe the polaronic effects in the catalysis of oxygen evolution reaction (OER) on the (110) surface of rutile TiO 2 . Theoretical understanding of polarons within the surface, coupled with simulations of OER at distinct titanium (Ti) and oxygen (O) active sites, reveals diverse polaron formation energies within the lattice sites with strong preference for bulk and surface bridge (O b ) oxygen sites. Moreover, we provide the evidence for the facilitative role of polarons in OER. Here, we find that hole polarons situated at the equatorial oxygen sites near the Ti–active site, along with bridge site hole polarons distal from the O b active site yield a small reduction in OER overpotential by ~0.06 eV and ~0.12 eV, respectively. However, subsurface, equatorial, and bridge site hole polarons significantly reduce the Ti-active site OER overpotential by ~0.4 eV through the peroxo–type oxygen pathway. We also observe that the presence of hole polarons stabilizes the *OH, *O, and *OOH intermediate species compared to the scenario without hole polarons. Overall, this study provides a detailed mechanistic insight into polaron–mediated OER, offering a promising avenue for improving the catalytic activity of transition metal oxide-based photocatalysts catering to renewable energy requisites.

14 SOLAR ENERGY↗

Open Specy 1.0: Automated (Hyper)spectroscopy for Microplastics

Microplastic spectral analysis is one of the most time-consuming processes in studying microplastic pollution, often requiring days per sample. Researchers are transitioning to automated batch and hyperspectral image analysis techniques to enhance efficiency. Open Specy, initially aimed at manual single-spectrum analysis, has now integrated automated methods. This updated version, Open Specy 1.0, introduces several new features, including two algorithms for automated processing (smoothing and particle compression), an extensive library containing over 40,000 open-source Raman and FTIR spectra, and two machine learning classifiers (logistic regression and k medoids) developed from this library. Furthermore, it includes a revamped user interface, an R package, and a benchmark data set for testing future advancements in automated techniques. Researchers evaluated various configurations for hyperspectral smoothing, particle identification, compression, and splitting, to achieve combined recovery rates between 50 and 150% particle counts, identities, and sizes with a coefficient of variation (CV) of less than 40% (the accredited standard). Mean absorbance times the standard deviation provided a consistent particle identification. Hyperspectral smoothing led to a 96% combined recovery rate and reduced variability (CV = 38%) compared to the 86% recovery (CV = 83%) of nonsmoothed controls. Additionally, compressing spectra for particles was significantly faster (>3x) and showed similar accuracy but with reduced variability than processing each pixel individually. Key challenges persist in automating spectral analysis, particularly in refining particle splitting algorithms, and improving identification routines to minimize false positives and negatives. In conclusion, new methods in sample preparation for better stabilization and dispersion of particles could overcome some of these issues.

13 HYDRO ENERGY↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Novel Chalcopyrites for Advanced Photoelectrochemical Water Splitting

With the support of DoE’s EERE office, our team has established a unique tool-chest of capabilities, including theoretical modeling (Lawrence Livermore National Laboratory: LLNL), state-of-the-art synthesis (Hawaii Natural Energy Institute: HNEI, Stanford, and the National Renewable Energy Laboratory: NREL) and advanced materials and interfaces characterization (University of Nevada, Las Vegas: UNLV, and Lawrence Berkeley National Laboratory: LBNL), to accelerate the development of high efficiency and durable chalcopyrite materials for advanced photoelectrochemical (PEC) water splitting. Using this synergistic approach, we have successfully created new wide bandgap chalcopyrite photocathodes generating over 10 mA/cm 2 , developed innovative strategies to protect them from corrosion, and engineered novel integration methods to circumvent thin film materials mechanical, chemical and thermal incompatibility. In Task 1 “Modeling and synthesis of chalcopyrite photocathodes”, we expanded our library of wide bandgap chalcopyrites for PEC water splitting. With support from LLNL’s “Computational Materials Diagnostics and Optimization of PEC Devices”, LBNL’s “photophysical” and NREL’s “I-III-VI Compound Semiconductors for Water-Splitting” nodes, we investigated two new chalcopyrite candidates for PEC water splitting: Cu(In,Al)Se 2 and Cu(In,B)Se 2 . We also further developed ordered vacancy compounds, such as CuGa 3 Se 5 , with unprecedented durability during PEC waters splitting in acidic solutions. In Task 2 “Interfaces engineering for enhanced efficiency and durability”, we addressed both the non-ideal band-edge positions of chalcopyrites with respect to water redox potentials, as well as their chemical instability under PEC water splitting, with a buried-junctions approach. With help from NREL’s “High-Throughput Experimental Thin Film Combinatorial Capabilities” and “Corrosion Analysis of Materials” nodes, we engineered environmentally friendly n-type buffers, including Mn x Zn 1-x O, to adjust the chalcopyrite band-edge positions and achieved photovoltages as high as 925 mV. Also, we integrated non-precious catalytic-protecting layers, such as WO 3 , to enhance the water splitting long-term stability of chalcopyrite absorbers. Finally, in Task 3 “Hybrid photoelectrode device integration”, we proposed an innovative method to bond wide bandgap photocathodes onto narrow bandgap PV drivers at room temperature using conductive polymers. Our semi-monolithic approach addressed fundamental processing incompatibility issues, as both the photocathode and the PV driver are processed separately. Proof-of-concept whole-chalcopyrite tandems were obtained by consecutive exfoliation and transfer of fully integrated 1.85 eV CuGa 3 Se 5 and 1.13 eV CuInGaSe 2 stacks from their Mo/SLG substrates onto a new single FTO host substrate.

08 HYDROGEN↗

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗

Effect of Ligand Polarity on the Internal Dipoles and Ferroelectric Distortion in BaTiO 3 Nanocubes

Here, surface adsorbates and surrounding matrix species have been demonstrated to affect the properties of nanoscale ferroelectrics and nanoscale ferroelectric composites; potentially counteracting performance losses that can occur in small particle sizes. In this work, the effects of nonpolar oleic acid (OA) and polar tetrafluoroborate (BF 4 - ) ligand capping on the surface of various sizes of BaTiO 3 nanocubes have been investigated with combined neutron diffraction and neutron pair distribution function (PDF), density functional theory (DFT), and ab initio molecular dynamics (AIMD) methods. The low real space PDF region provides an unobstructed view of rhombohedral (split short and long) Ti-O distances in BaTiO 3 nanocubes, mimicking the well-established order-disorder local structure found in bulk BaTiO 3 . Interestingly, the intermediate-range order in nanocubes is found to be orthorhombic, rather than tetragonal. It is concluded that polar ligands adsorbed at BaTiO 3 surfaces stabilize the correlation length scale of local rhombohedral distortions in ferroelectric nanoparticles relative to nonpolar ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Kβ X-ray Emission Spectra Analysis Using Bayesian Optimization

The Kβ X-ray emission spectrum of 3 d transition metals is rich with electronic and structural information due to strong exchange interactions with the valence shell of the metal, and has become crucial for understanding their spin and oxidation states. The spectrum is commonly treated using crystal-field multiplet theory, a semi-empirical theory that uses tunable parameters to control the strength of the effects present in X-ray emission spectroscopy (XES). However, determining the experimental values of these parameters remains a challenge. We present a methodology that applies Bayesian optimization to crystal-field multiplet theory to determine parameter values. The algorithm is tested on the X-ray emission spectra of a collection of Mn, Co, and Ni oxides. We are able to find optimal values for the four most impactful parameters: Slater−Condon reduction factors F dd , F pd , and G pd , and crystal field splitting 10 Dq . The algorithm produces significantly improved accuracy compared to current analysis methods, and probes interparameter dependencies by modeling the error landscape. This advancement enhances XES analysis by offering an approach of obtaining quantitative electronic structural information on 3 d transition metal valence shells, facilitating applications across various scientific fields.

Bayesian optimization↗

CATMoS: Collaborative Acute Toxicity Modeling Suite

Background: Humans are exposed to tens of thousands of chemical substances that need to be assessed for their potential toxicity. Acute systemic toxicity testing serves as the basis for regulatory hazard classification, labeling, and risk management. However, it is cost- and time-prohibitive to evaluate all new and existing chemicals using traditional rodent acute toxicity tests. In silico models built using existing data facilitate rapid acute toxicity predictions without using animals. Objectives: The U.S. Interagency Coordinating Committee on the Validation of Alternative Methods Acute Toxicity Workgroup organized an international collaboration to develop in silico models for predicting acute oral toxicity based on five different endpoints: LD50 value, U.S. Environmental Protection Agency hazard categories, Globally Harmonized System for Classification and Labelling hazard categories, very toxic chemicals (LD50 =50 mg/kg), and non-toxic chemicals (LD50 >2000 mg/kg). Methods: An acute oral toxicity data inventory for 11,992 chemicals was compiled, split into training and evaluation sets, and made available to 35 participating international research groups that submitted a total of 139 predictive models. Predictions that fell within the applicability domains of the submitted models were evaluated using external validation sets. These were then combined into consensus models to leverage strengths of individual approaches. Results: The resulting consensus predictions, which leverage the collective strengths of each individual model, form the Collaborative Acute Toxicity Modeling Suite (CATMoS). CATMoS demonstrated high performance in terms of accuracy and robustness when compared to in vivo results. Discussion: CATMoS is being evaluated by regulatory agencies for its utility and applicability as a potential replacement for in vivo rat acute oral toxicity studies. CATMoS predictions for over 800,000 chemicals have been made available via the NTP’s Integrated Chemical Environment. The models are also implemented in a free, standalone open-source tool, OPERA, which allows predictions of new and untested chemicals to be made.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

Integration of Concentrating Solar Power with High Temperature Electrolysis for Hydrogen Production: Preprint

Hydrogen (H2) has been identified as a leading sustainable contender to replace fossil fuels in transportation and electricity generation. H2 production can be achieved by concentrating solar thermal power (CSP) systems collecting thermal energy from the sun to various chemical processes for fuel production. Fuel production via solar thermal chemical processes integrated with CSP uses the full spectrum of sunlight compared with photovoltaic power conversion and stores solar energy directly and efficiently [1]. The solar fuel production can be realized by thermochemical processes (e.g., water splitting for H2 production, carbon dioxide reduction, or methane reforming) or thermal electrochemical methods (e.g., integration with solid oxide electrolysis cell). Technology development for CSP-integrated solar fuel production requires broad technological bases from solar energy collection to chemical energy conversion. H2 generated from renewable sources can be an energy carrier for a carbon-free economy. Integrating CSP with high temperature electrolysis (HTE) using solid oxide electrolysis cells (SOEC) provides a renewable path for H2 generation. The CSP-HTE integration approach provides the benefit of thermal energy storage (TES) for continuous operation, improved capacity, and SOEC life. H2 gas has low energy density for transportation, pipeline networks are expensive, and H2 liquefaction is energy intensive. However, an alternative method for H2 distribution is to use carbon dioxide (CO2) capture and liquid hydrocarbon synthesis to convert solar energy into liquid fuels that are compatible with the existing fossil fuel infrastructure.

concentrating solar thermal power↗

Baryon fraction from the BAO amplitude: a consistent approach to parameterizing perturbation growth

Galaxy clustering constrains the baryon fraction Omega_b/Omega_m through the amplitude of baryon acoustic oscillations and the suppression of perturbations entering the horizon before recombination. This produces a different pre-recombination distribution of baryons and dark matter. After recombination, the gravitational potential responds to both components in proportion to their mass, allowing robust measurement of the baryon fraction. This is independent of new-physics scenarios altering the recombination background (e.g. Early Dark Energy). The accuracy of such measurements does, however, depend on how baryons and CDM are modeled in the power spectrum. Previous template-based splitting relied on approximate transfer functions that neglected part of information. We present a new method that embeds an extra parameter controlling the balance between baryons and dark matter in the growth terms of the perturbation equations in the CAMB Boltzmann solver. This approach captures the baryonic suppression of CDM prior to recombination, avoids inconsistencies, and yields a clean parametrization of the baryon fraction in the linear power spectrum, separating out the simple physics of growth due to the combined matter potential. We implement this framework in an analysis pipeline using Effective Field Theory of Large-Scale Structure with HOD-informed priors and validate it against noiseless LCDM and EDE cosmologies with DESI-like errors. The new scheme achieves comparable precision to previous splitting while reducing systematic biases, providing a more robust way to baryon-fraction measurements. In combination with BBN constraints on the baryon density and Alcock-Paczynski estimates of the matter density, these results strengthen the use of baryon fraction measurements to derive a Hubble constant from energy densities, with future DESI and Euclid data expected to deliver competitive constraints.

Crespi, Andrea [U. Waterloo (main); Waterloo U., I↗

Pooled PPIseq: Screening the SARS-CoV-2 and human interface with a scalable multiplexed protein-protein interaction assay platform

Protein-Protein Interactions (PPIs) are a key interface between virus and host, and these interactions are important to both viral reprogramming of the host and to host restriction of viral infection. In particular, viral-host PPI networks can be used to further our understanding of the molecular mechanisms of tissue specificity, host range, and virulence. At higher scales, viral-host PPI screening could also be used to screen for small-molecule antivirals that interfere with essential viral-host interactions, or to explore how the PPI networks between interacting viral and host genomes co-evolve. Current high-throughput PPI assays have screened entire viral-host PPI networks. However, these studies are time consuming, often require specialized equipment, and are difficult to further scale. Here, we develop methods that make larger-scale viral-host PPI screening more accessible. This approach combines the mDHFR split-tag reporter with the iSeq2 interaction-barcoding system to permit massively-multiplexed PPI quantification by simple pooled engineering of barcoded constructs, integration of these constructs into budding yeast, and fitness measurements by pooled cell competitions and barcode-sequencing. We applied this method to screen for PPIs between SARS-CoV-2 proteins and human proteins, screening in triplicate >180,000 ORF-ORF combinations represented by >1,000,000 barcoded lineages. Our results complement previous screens by identifying 74 putative PPIs, including interactions between ORF7A with the taste receptors TAS2R41 and TAS2R7, and between NSP4 with the transmembrane KDELR2 and KDELR3. We show that this PPI screening method is highly scalable, enabling larger studies aimed at generating a broad understanding of how viral effector proteins converge on cellular targets to effect replication.

60 APPLIED LIFE SCIENCES↗

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗