Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Graphics”

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 847 records · Page 47

Oleanolic acid activates daf-16 to increase lifespan in Caenorhabditis elegans

Oleanolic acid (OA) is an active ingredient in natural plants. It has been reported to possess a variety of pharmacological activities, but very little is known about its effects of anti-aging. We investigate here whether OA has an impact on longevity in vivo, and more specifically, we have examined effects of OA on the lifespan and stress tolerance in Caenorhabditis elegans (C. elegans). Our results showed that OA could extend the lifespan, increase its stress resistance and reduce the intracellular reactive oxygen species (ROS) in wild-type worms. Moreover, we have found that OA-induced longevity may not be associated with the calorie restriction (CR) mechanism. Our mechanistic studies using daf-16 loss-of-function mutant strains (GR1307) indicated that the extension of lifespan by OA requires daf-16. In addition, OA treatment could also modulate the nuclear localization, and the quantitative real-time PCR results revealed that up-regulation of daf-16 target genes such as sod-3, hsp-16.2 and ctl-1 could prolong lifespan and increase stress response in C. elegans. This study overall uncovers the longevity effect of OA and its underpinning mechanisms. - Graphical abstract: Oleanolic acid modulates the activity of DAF-16 to promote longevity and increase stress resistance in Caenorhabditis elegans. - Highlights: • OA extends the lifespan of wild-type Caenorhabditis elegans. • OA improves the stress resistance and reduces the intracellular ROS level in C. elegans. • OA induces lifespan extension may not proceed through the CR mechanism. • OA extends the lifespan in C. elegans is modulated by daf-16.

60 APPLIED LIFE SCIENCES↗

Logical design of anti-prion agents using NAGARA

To accelerate the logical drug design procedure, we created the program “NAGARA,” a plugin for PyMOL, and applied it to the discovery of small compounds called medical chaperones (MCs) that stabilize the cellular form of a prion protein (PrP{sup C}). In NAGARA, we constructed a single platform to unify the docking simulation (DS), free energy calculation by molecular dynamics (MD) simulation, and interfragment interaction energy (IFIE) calculation by quantum chemistry (QC) calculation. NAGARA also enables large-scale parallel computing via a convenient graphical user interface. Here, we demonstrated its performance and its broad applicability from drug discovery to lead optimization with full compatibility with various experimental methods including Western blotting (WB) analysis, surface plasmon resonance (SPR), and nuclear magnetic resonance (NMR) measurements. Combining DS and WB, we discovered anti-prion activities for two compounds and tegobuvir (TGV), a non-nucleoside non-structural protein NS5B polymerase inhibitor showing activity against hepatitis C virus genotype 1. Binding profiles predicted by MD and QC are consistent with those obtained by SPR and NMR. Free energy analyses showed that these compounds stabilize the PrP{sup C} conformation by decreasing the conformational fluctuation of the PrP{sup C}. Because TGV has been already approved as a medicine, its extension to prion diseases is straightforward. Finally, we evaluated the affinities of the fragmented regions of TGV using QC and found a clue for its further optimization. By repeating WB, MD, and QC recursively, we were able to obtain the optimum lead structure. - Highlights: • NAGARA integrates docking simulation, molecular dynamics, and quantum chemistry. • We found many compounds, e.g., tegobuvir (TGV), that exhibit anti-prion activities. • We obtained insights into the action mechanism of TGV as a medical chaperone. • Using QC, we obtained useful information for optimization of the lead compound, TGV. • NAGARA is a convenient platform for drug discovery and lead optimization.

60 APPLIED LIFE SCIENCES↗

An easily regenerable enzyme reactor prepared from polymerized high internal phase emulsions

A large-scale high-efficient enzyme reactor based on polymerized high internal phase emulsion monolith (polyHIPE) was prepared. First, a porous cross-linked polyHIPE monolith was prepared by in-situ thermal polymerization of a high internal phase emulsion containing styrene, divinylbenzene and polyglutaraldehyde. The enzyme of TPCK-Trypsin was then immobilized on the monolithic polyHIPE. The performance of the resultant enzyme reactor was assessed according to the conversion ability of N{sub α}-benzoyl-L-arginine ethyl ester to N{sub α}-benzoyl-L-arginine, and the protein digestibility of bovine serum albumin (BSA) and cytochrome (Cyt-C). The results showed that the prepared enzyme reactor exhibited high enzyme immobilization efficiency and fast and easy-control protein digestibility. BSA and Cyt-C could be digested in 10 min with sequence coverage of 59% and 78%, respectively. The peptides and residual protein could be easily rinsed out from reactor and the reactor could be regenerated easily with 4 M HCl without any structure destruction. Properties of multiple interconnected chambers with good permeability, fast digestion facility and easily reproducibility indicated that the polyHIPE enzyme reactor was a good selector potentially applied in proteomics and catalysis areas. - Graphical abstract: Schematic illustration of preparation of hypercrosslinking polyHIPE immobilized enzyme reactor for on-column protein digestion. - Highlights: • A reactor was prepared and used for enzyme immobilization and continuous on-column protein digestion. • The new polyHIPE IMER was quite suit for protein digestion with good properties. • On-column digestion revealed that the IMER was easy regenerated by HCl without any structure destruction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Core-6 fucose and the oligomerization of the 1918 pandemic influenza viral neuraminidase

The 1918 H1N1 influenza virus was responsible for one of the most deadly pandemics in human history. Yet to date, the structure component responsible for its virulence is still a mystery. In order to search for such a component, the neuraminidase (NA) antigen of the virus was expressed, which led to the discovery of an active form (tetramer) and an inactive form (dimer and monomer) of the protein due to different glycosylation. In this report, the N-glycans from both forms were released and characterized by mass spectrometry. It was found that the glycans from the active form had 26% core-6 fucosylated, while the glycans from the inactive form had 82% core-6 fucosylated. Even more surprisingly, the stalk region of the active form was almost completely devoid of core-6-linked fucose. These findings were further supported by the results obtained from in vitro incorporation of azido fucose and {sup 3}H-labeled fucose using core-6 fucosyltransferase, FUT8. In addition, the incorporation of fucose did not change the enzymatic activity of the active form, implying that core-6 fucose is not directly involved in the enzymatic activity. It is postulated that core-6 fucose prohibits the oligomerization and subsequent activation of the enzyme. - Graphical abstract: Proposed mechanism for how core-fucose prohibits the tetramerization of the 1918 pandemic viral neuraminidase. Only the cross section of the stalk region with two N-linked glycans are depicted for clarity. (A) Carbohydrate–carbohydrate interaction on non-fucosylated monomer allows tetramerization. (B) Core-fucosylation disrupts the interaction and prevents the tetramerization. - Highlights: • Expressed 1918 pandemic influenza viral neuraminidase has inactive and active forms. • The inactive form contains high level of core-6 fucose, while the active form lacks such modification. • Core fucose could interfere the oligomerization of the neuraminidase and thus its activation. • This discovery may explain why 1918 pandemic influenza caused higher death rate among young population.

60 APPLIED LIFE SCIENCES↗

Ghrelin ameliorates the human alveolar epithelial A549 cell apoptosis induced by lipopolysaccharide

Ghrelin is a gastric acyl-peptide that plays an inhibitory role in cell apoptosis. Herein we investigate the protective effects of ghrelin in LPS-induced apoptosis of human alveolar epithelial A549 cells, along with the possible molecular mechanisms. LPS exposure impaired cell viability and increased apoptosis of A549 cells significantly in concentration- and time-dependent manners embodied in increased Bax and cleaved caspase-3 production, coupled with decreased Bcl-2 levels. Simultaneously, LPS remarkably decreased the expression of phosphatidylinositol 3 kinase/protein kinase B (PI3K/Akt) and extracellular signal-regulated kinas (ERK) in A549 cells. However, ghrelin'pretreatment ameliorated LPS-caused alterations in the ratio of Bax/Bcl-2 and cleaved caspase-3 expression, whereas activated the PI3K/Akt and ERK signaling. These results demonstrate that ghrelin lightens LPS-induced apoptosis of human alveolar epithelial cells partly through activating the PI3K/Akt and ERK pathway and thereby might benefit alleviating septic ALI. -- Graphical abstract: Ghrelin ameliorates the human alveolar epithelial A549 cells apoptosis induced by lipopolysaccharide partly through activating the PI3K/Akt and ERK pathway. Display Omitted -- Highlights: •It has been observed that LPS insult significantly increased apoptosis in A549 cells. •Both Akt and ERK signaling are critical adapter molecules to mediate the ghrelin-mediated proliferative effect. •Ghrelin may have a therapeutic effect in the prevention of LPS-induced apoptosis.

60 APPLIED LIFE SCIENCES↗

Electronic cigarette aerosols and copper nanoparticles induce mitochondrial stress and promote DNA fragmentation in lung fibroblasts

Oxidants or nanoparticles have recently been identified as constituents of aerosols released from various styles of electronic cigarettes (E-cigs). Cells in the lung may be directly exposed to these constituents and harbor reactive properties capable of incurring acute cell injury. Our results show mitochondria are sensitive to both E-cig aerosols and aerosol containing copper nanoparticles when exposed to human lung fibroblasts (HFL-1) using an Air-Liquid Interface culture system, evident by elevated levels of mitochondrial ROS (mtROS). Increased mtROS after aerosol exposure is associated with reduced stability of OxPhos electron transport chain (ETC) complex IV subunit and nuclear DNA fragmentation. Increased levels of IL-8 and IL-6 in HFL-1 conditioned media were also observed. These findings reveal both mitochondrial, genotoxic, and inflammatory stresses are features of direct cell exposure to E-cig aerosols which are ensued by inflammatory duress, raising a concern on deleterious effect of vaping. - Graphical abstract: Oxidants and possibly reactive properties of metal particles in E-cig aerosols impart mitochondrial oxidative stress and DNA damage. These biological effects accompany inflammatory response which may raise concern regarding long term E-cig use. Mitochondria may be particularly sensitive to reactive properties of E-cig aerosols in addition to the potential for them to induce genotoxic stress by generating increased ROS. - Highlights: • Mitochondria are sensitive to both E-cig aerosols and metal nanoparticles. • Increased mtROS by E-cig aerosol is associated with disrupted mitochondrial energy. • E-cig causes nuclear DNA fragmentation. • E-cig aerosols induce pro-inflammatory response in human fibroblasts.

60 APPLIED LIFE SCIENCES↗

Functional characterization of archaeal homologs of human nuclear RNase P proteins Rpp21 and Rpp29 provides insights into the molecular basis of their cooperativity in catalysis

Ribonuclease P (RNase P) is a ribonucleoprotein that catalyzes the processing of 5′ leader sequences of precursor tRNAs (pre-tRNA). RNase P proteins PhoRpp21 and PhoRpp29 in the hyperthermophilic archaeon Pyrococcus horikoshii, homologs of human nuclear RNase P proteins Rpp21 and Rpp29 respectively, fold into a heterodimeric structure and synergistically function in the activation of the specificity domain (S-domain) in RNase P RNA (PhopRNA). To elucidate the molecular basis for their cooperativity, we first analyzed binding ability to PhopRNA using a pull-down assay. The result showed that PhoRpp21 is able to bind to PhopRNA in the absence of PhoRpp29, whereas PhoRpp29 alone has reduced affinity to PhopRNA, suggesting that PhoRpp21 primarily functions as a binding element for PhopRNA in the PhoRpp21-PhoRpp29 complex. Mutational analyses suggested that although individual positively charged clusters contribute little to the PhopRNA binding, Lys53, Lys54, and Lys56 at the N-terminal helix (α2) in PhoRpp21 and 10 C-terminal residues in PhoRpp29 are essential for PhopRNA activation. Moreover, deletion of a single stranded loop linking P11 and P12 helices in the PhopRNA S-domain impaired the PhoRpp21-PhoRpp29 complex binding to PhopRNA. Collectively, the present results suggest that PhoRpp21 binds the loop between P11 and P12 helices through overall positively charged clusters on the surface of the complex and serves as a scaffold for PhoRpp29 to optimize structural conformation of its N-terminal helix (α2) in PhoRpp21, as well as C-terminal residues in PhoRpp29, for RNase P activity. - Graphical abstract: A mechanistic model for the synergistic activation of PhopRNA by PhoRpp21 and PhoRpp29. PhoRpp21 binds to the loop between P11 and P12 helices in PhopRNA and serves as a scaffold for PhoRpp29 to optimize a structural conformation of a positively charged edge composed of lysine residues (●) at α2 in PhoRpp21and C-terminal residues (â–‹) in PhoRpp29 on the protein complex which probably stabilize stacking interactions between PhopRNA and pre-tRNA. ▲ indicates the processing site in pre-tRNA. Display Omitted - Highlights: • RNase P proteins PhoRpp21 and PhoRpp29 synergistically function in RNase P catalysis. • PhoRpp21 functions as an RNA-binding element in the PhoRpp21-PhoRpp29 complex. • PhoRpp21 serves as a scaffold for PhoRpp29 in RNase P. • PhoRpp21 binds the loop between P11 and P12 helices in the PhopRNA S-domain. • The molecular basis for the cooperativity of PhoRpp21 and PhoRpp29 is presented.

60 APPLIED LIFE SCIENCES↗

Substituted (E)-2-(2-benzylidenehydrazinyl)-4-methylthiazole-5-carboxylates as dual inhibitors of 15-lipoxygenase & carbonic anhydrase II: synthesis, biochemical evaluation and docking studies

15-Lipoxygenase (15-LOX) plays a major role in many inflammatory lung diseases including chronic obstructive pulmonary disease (COPD), asthma and chronic bronchitis. Over-expression of 15-LOX is related with some specific carcinomas including pancreatic, gastric and brain tumor. Similarly among different isozymes of carbonic anhydrase (CA), CA II is expressed in pancreatic, gastric carcinomas as well as in brain tumors. Therefore, novel potent inhibitors of both 15-LOX and CA II are required to explore the role of these enzymes further and to enable the drug discovery efforts. For this purpose, a series of benzyledinyl-hydrazinyl substituted thiazole derivatives were designed, synthesized and characterized by FTIR, {sup 1}H, &{sup 13}C NMR spectroscopy. The derivatives were then evaluated for their potential to inhibit 15-LOX and bovine carbonic anhydrase II (bCA II). Most of these compounds showed excellent inhibitory potential for 15-LOX with an IC{sub 50} of 0.12 ± 0.002 to 0.69 ± 0.5 μM and showed moderate inhibition potency for bCA II with compound 5h (IC{sub 50} = 1.26 ± 0.24 μM) being the most active. The most potent compound 5a that emerged as a dual inhibitor of both enzymes, exhibiting 24 times greater selectivity for 15-LOX over bCA II. Compound 5a exhibited dual potent inhibitory activity against both 15-LOX and bCA II enzymes having IC{sub 50} values of 0.12 ± 0.002 and 2.93 ± 0.22 μM, respectively. Molecular docking studies of potent as well as dual inhibitors were also carried out to provide an insight into the binding site interactions. - Graphical abstract: Putative binding mode of most potent compound in the active pocket of 15-lipoxygenase. Display Omitted.

60 APPLIED LIFE SCIENCES↗

Biomarker identification for statin sensitivity of cancer cell lines

Highlights: • Graphical modeling is used to find a biomarker signature for statin sensitivity. • Predictions from the model were experimentally validated in four human cell lines. • Enrichment analysis and predicted combination therapies show relevance of signature. Statins are potent cholesterol reducing drugs that have been shown to reduce tumor cell proliferation in vitro and tumor growth in animal models. Moreover, retrospective human cohort studies demonstrated decreased cancer-specific mortality in patients taking statins. We previously implicated membrane E-cadherin expression as both a marker and mechanism for resistance to atorvastatin-mediated growth suppression of cancer cells; however, a transcriptome-profile-based biomarker signature for statin sensitivity has not yet been reported. Here, we utilized transcriptome data from fourteen NCI-60 cancer cell lines and their statin dose-response data to produce gene expression signatures that identify statin sensitive and resistant cell lines. We experimentally confirmed the validity of the identified biomarker signature in an independent set of cell lines and extended this signature to generate a proposed statin-sensitive subset of tumors listed in the TCGA database. Finally, we predicted drugs that would synergize with statins and found several predicted combination therapies to be experimentally confirmed. The combined bioinformatics-experimental approach described here can be used to generate an initial biomarker signature for anticancer drug therapy.

60 APPLIED LIFE SCIENCES↗

Multichannel Analysis of Surface Waves Accelerated (MASWAccelerated): Software for efficient surface wave inversion using MPI and GPUs

Multichannel Analysis of Surface Waves (MASW) is a technique frequently used in geotechnical engineering and engineering geophysics to infer 1D layered models of seismic shear wave velocities in the top tens to hundreds of meters of the subsurface. We aim to accelerate MASW calculations by capitalizing on modern computer hardware available in the workstations of most engineers: multiple cores and graphics processing units (GPUs). We propose new parallel and GPU accelerated algorithms for computing 1D MASW inversion, and provide software implementations in C using Message Passing Interface (MPI) and CUDA. These algorithms take advantage of sparsity that arises in the problem, and the work balance between processes considers typical data trends. We compare our methods to an existing open source Matlab MASW tool. Our serial C implementation achieves a 2x speedup over the Matlab software, and we continue to see improvements by parallelizing the problem with MPI. Here we see nearly perfect strong and weak scaling for uniform data, and improve strong scaling for realistic data by repartitioning the problem to process mapping. By utilizing GPUs available on most modern workstations, we observe an additional 1.3x speedup over the serial C implementation on the first use of the method. We typically repeatedly evaluate theoretical dispersion curves as part of an optimization procedure, and on the GPU the kernel can be cached for faster reuse on later runs. We observe a 3.2x speedup on the cached GPU runs compared to the serial C runs. This work is the first open-source parallel or GPU-accelerated software tool for MASW imaging, and should enable geotechnical engineers to fully utilize all computer hardware at their disposal.

58 GEOSCIENCES↗

An integrated modeling framework with open architecture for phase field simulation of multi-component alloys

An integrated modeling framework (PanPhaseField) has been developed, which enables a direct and fast coupling between CALPHAD calculations and large-scale phase field simulations for multi-component alloys. Further, it adopts an open architecture allowing for integration of user-defined phase field models in a plug-and-play manner by taking full advantage of the user-friendly graphical interface of Pandat software. The developed modeling platform becomes an enabling tool that can be used to simulate the evolution of spatially varying microstructures of industrial complex alloys for various engineering applications.

36 MATERIALS SCIENCE↗

Historical experiment to measure irradiation-induced creep of graphite

This paper presents historical results of graphite irradiation-induced creep experiments that were performed at Oak Ridge National Laboratory from the 1950's to the 1970's. These experiments were performed at temperatures from 150°C to 1000 °C, and bend stresses ranging from 500 to 5000 psi (~3.3–34.5 MPa). The experimental setup utilized in-situ measurement of specimen displacement, on-line applied stress control, and the ability to change stress during the experiment. The different stress conditions showed that the primary creep strain and the steady-state creep rates both have a linear stress dependence. The temperature range used in this work resulted in trends that have not be previously presented in the literature: 1) a linear dependence of primary creep strain on temperature, and 2) the shape of steady state creep rate versus temperature (see graphical abstract). The maximum dose in the specimens was 0.9 dpa, which is sufficient to achieve steady-state creep without the structural changes that alter the observed creep behavior. The results from this experiment provide evidence that dispels that the pinning-unpinning model describes the mechanism of irradiation creep in graphite. Instead these results suggest a dislocation climb mechanism is the probable mechanism for creep within the crystalline regions.

36 MATERIALS SCIENCE↗

The crystal and coordination chemistry of neptunium in all its oxidation states: An expanded structural hierarchy of neptunium compounds

We report structures of 388 neptunium-bearing crystal structures are sorted into an expanded structural hierarchy and examined. The majority of structures contain Np(V) (149), Np(VI) (99), or Np(IV) (91) whereas structures containing Np(II), Np(III), Np(VII), or more than one oxidation state are rare. Coordination geometry of Np(V) and Np(VI) compounds is dominated by the approximately linear neptunyl cation. Actinyl–actinyl interactions (AAIs) are most prevalent in Np(V) compounds and are observed in 48% of these extended structures. Prismatic coordination geometries are commonly observed in Np(IV) compounds. Np(III) organometallic complexes display a variety of coordination modes. Coordination chemistry of Np(VII) compounds is dominated by the tetra-oxido core. Updated average bond length values are provided for Np(IV), Np(V), Np(VI), and Np(VII) compounds and updated bond valence parameters for Np(IV) are reported. A graphical notation system for neptunyl structures coordinated by organic ligands was developed to facilitate discussion and identify similarities between neptunyl compounds coordinated by oxyanions and neptunyl compounds coordinated by organic ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Production of H{sub 2} by water radiolysis in cement paste under electron irradiation: A joint experimental and theoretical study

Long-term confinement of nuclear waste is one of the main challenges faced by the nuclear industry. Fission products such as {sup 90}Sr and {sup 137}Cs, both β{sup −} emitters known to induce serious health hazards, represent the largest fraction of nuclear waste. Cement is a good candidate to store them, provided it can resist the effects of irradiation over time. Here, we have investigated the effects of β{sup −} decay on cement by performing electron irradiation experiments on different samples. We show that H{sub 2} production in cement, the main effect of water radiolysis, depends strongly on composition and relative humidity. First-principles calculations indicate that the water-rich interlayer regions with Ca{sup 2+} ions act as electron traps that promote the formation of H{sub 2}. They also show that holes localize in water-rich regions in low Ca content samples and are then able to participate in H{sub 2} production. This work provides new understanding of radiolysis effects in cements. - Graphical abstract: H{sub 2} produced by radiolysis in cement paste is investigated by experimental and computational studies.

36 MATERIALS SCIENCE↗

edxia: Microstructure characterisation from quantified SEM-EDS hypermaps

The characterisation of cement paste microstructure is an important step towards understanding durability mechanisms in cementitious materials. Scanning electron microscopy (SEM) coupled with energy dispersive spectroscopy (EDS) is a widely used technique to analyse the microstructure at the micron-scale. However, it is challenging, notably because the characteristic size of many phases is found on a scale smaller than the EDS interaction volume. This work presents a new image analysis framework to identify phases and quantify the microstructure of cementitious materials from SEM-EDS hypermaps. By leveraging domain knowledge, representative points are attributed to phases and mixtures of phases based on ratio plots. Then, quantitative analysis of the microstructure can be carried out (chemical composition, particle size distributions, volume fractions, …). We demonstrate the abilities of the framework, and we present possible applications and extensions of the method. The framework is available as both a graphical interface and a Python code.

36 MATERIALS SCIENCE↗

Nonlinear elasticity with the Shifted Boundary Method

Here, we propose a new unfitted/immersed computational framework for nonlinear solid mechanics, which bypasses the complexities associated with the generation of CAD representations and subsequent body-fitted meshing. This approach allows to speed up the cycle of design and analysis in complex geometry and requires relatively simple computer graphics representations of the surface geometries to be simulated, such as the Standard Tessellation Language (STL format). Complex data structures and integration on cut elements are avoided by means of an approximate boundary representation and a modification (shifting) of the boundary conditions to maintain optimal accuracy. An extensive set of computational experiments in two and three dimensions is included.

97 MATHEMATICS AND COMPUTING↗

On the application of betweenness centrality in chemical network analysis: Computational diagnostics and model reduction

Here, the concept of shortest path and betweenness centrality (BC) is introduced in combustion systems, with its application demonstrated in chemical network analysis and skeletal mechanism generation. After establishing the graphical representation of a chemical network for given pressure, temperature, and species concentrations, a metric BC is defined to rank the shortest paths passed by other nodes between the pair of source and target node, and as such captures the dominant indirect kinetic pathways between any pair of nodes in the network for computational diagnostics. As a result, the controlling pathways are retained by collecting the species with larger BC values. Such a notion then indicates that the index of importance of species could then be assigned based on their BC values to further guide skeletal mechanism generation. Unlike existing methods, the betweenness centrality approach takes into account of both the fluxes between species and their relative positioning within the chemical network. To demonstrate its potential utility to combustion studies, the approach was applied to the GRI-3.0, LLNL and USC-Mech II mechanisms to identify the important pathways in the chemical network at each local reaction state, and develop skeletal mechanisms from all reaction state samples in auto-ignition and perfectly stirred reactor (PSR) simulations. The performance of the BC ranking is compared to the methods of directed relation graph (DRG), DRG with error propagations (DRGEP) and sensitivity analysis (SA), and is shown to possess sufficient utility in producing skeletal mechanisms with good accuracy and flexibility for the cases studied.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating dissipative particle dynamics simulations for soft matter systems

Dissipative particle dynamics (DPD) is a coarse-grained particle-based simulation method that offers microscopic-scale insights into soft matter systems. We present an efficient implementation of a DPD model for graphical processing units (GPUs). As implemented in the LAMMPS molecular dynamics package, it can run effectively on current-generation supercomputers which often have hybrid nodes containing multi-core CPUs and (one or more) GPUs. Using efficient communication of information between the CPUs and GPUs, DPD interactions can be computed on the GPU while other portions of a full simulation model (boundary conditions, constraints, bonded interactions, diagnostic calculations, etc.) can be performed on the CPU. Our GPU-enhanced runs show a speedup of up to 9.5× versus many-core CPU simulations, and can run scalably across thousands of compute nodes. We briefly discuss how the new GPU implementation was validated against the CPU version for thermodynamics, diffusion, and hydrodynamic behavior. We also highlight large-scale models which the faster DPD implementation has enabled, for studies of monolayer self-assembly and thin-film instabilities.

36 MATERIALS SCIENCE↗