Results for “parallel tempering” — Page 2

DOE OSTI2020

BASS: An R Package for Fitting and Performing Sensitivity Analysis of Bayesian Adaptive Spline Surfaces

In this work, we present the R package BASS as a tool for nonparametric regression. The primary focus of the package is fitting fully Bayesian adaptive spline surface (BASS) models and performing global sensitivity analyses of these models. The BASS framework is similar to that of Bayesian multivariate adaptive regression splines (BMARS) from Denison, Mallick, and Smith (1998), but with many added features. The software is built to efficiently handle significant amounts of data with many continuous or categorical predictors and with functional response. Under our Bayesian framework, most priors are automatic but these can be modified by the user to focus on parsimony and the avoidance of overfitting. If directed to do so, the software uses parallel tempering to improve the reversible jump Markov chain Monte Carlo (RJMCMC) methods used to perform inference. We discuss the implementation of these features and present the performance of BASS in a number of analyses of simulated and real data.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Practical applications of machine-learned flows on gauge fields

Normalizing flows are machine-learned maps between different lattice theories which can be used as components in exact sampling and inference schemes. Ongoing work yields increasingly expressive flows on gauge fields, but it remains an open question how flows can improve lattice QCD at state-of-the-art scales. We discuss and demonstrate two applications of flows in replica exchange (parallel tempering) sampling, aimed at improving topological mixing, which are viable with iterative improvements upon presently available flows.

Abbott, Ryan↗

DOE OSTI2018

Sensitivity Analysis and Emulation for Functional Data using Bayesian Adaptive Splines

When a computer code is used to simulate a complex system, one of the fundamental tasks is to assess the sensitivity of the simulator to the different input parameters. In the case of computationally expensive simulators, this is often accomplished via a surrogate statistical model, a statistical output emulator. An effective emulator is one that provides good approximations to the computer code output for wide ranges of input values. In addition, an emulator should be able to handle large dimensional simulation output for a relevant number of inputs; it should flexibly capture heterogeneities in the variability of the response surface; it should be fast to evaluate for arbitrary combinations of input parameters, and it should provide an accurate quantification of the emulation uncertainty. In this paper we discuss the Bayesian approach to multivariate adaptive regression splines (BMARS) as an emulator for a computer model that outputs curves. We introduce modifications to traditional BMARS approaches that allow for fitting large amounts of data and allow for more efficient MCMC sampling. We emphasize the ease with which sensitivity analysis can be performed in this situation. We present a sensitivity analysis of a computer model of the deformation of a protective plate used in pressure-driven experiments. Our example serves as an illustration of the ability of BMARS emulators to fulfill all the necessities of computability, flexibility and reliable calculation on relevant measures of sensitivity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Rapid in situ diversification rates in Rhamnaceae explain the parallel evolution of high diversity in temperate biomes from global to local scales

Summary The macroevolutionary processes that have shaped biodiversity across the temperate realm remain poorly understood and may have resulted from evolutionary dynamics related to diversification rates, dispersal rates, and colonization times, closely coupled with Cenozoic climate change. We integrated phylogenomic, environmental ordination, and macroevolutionary analyses for the cosmopolitan angiosperm family Rhamnaceae to disentangle the evolutionary processes that have contributed to high species diversity within and across temperate biomes. Our results show independent colonization of environmentally similar but geographically separated temperate regions mainly during the Oligocene, consistent with the global expansion of temperate biomes. High global, regional, and local temperate diversity was the result of high in situ diversification rates, rather than high immigration rates or accumulation time, except for Southern China, which was colonized much earlier than the other regions. The relatively common lineage dispersals out of temperate hotspots highlight strong source‐sink dynamics across the cosmopolitan distribution of Rhamnaceae. The proliferation of temperate environments since the Oligocene may have provided the ecological opportunity for rapid in situ diversification of Rhamnaceae across the temperate realm. Our study illustrates the importance of high in situ diversification rates for the establishment of modern temperate biomes and biodiversity hotspots across spatial scales.

Plant Sciences↗

DOE OSTI2018

Evolutionary and demographic correlates of Pleistocene coastline changes in the Sicilian wall lizard Podarcis wagleriana

Abstract Aim Emergence of coastal lowlands during Pleistocene ice ages might have provided conditions for glacial expansions (demographic and spatial), rather than contraction, of coastal populations of temperate species. Here, we tested these predictions in the insular endemic Sicilian wall lizard Podarcis wagleriana . Location Sicily and neighbouring islands. Methods We sampled 179 individuals from 45 localities across the whole range of P. wagleriana . We investigated demographic and spatial variations through time using Bayesian coalescent models (Bayesian phylogeographic reconstruction, Extended Bayesian Skyline plots, Isolation‐with‐migration models) based on multilocus DNA sequence data. We used species distribution modelling to reconstruct present and past habitat suitability. Results We found two main lineages distributed in the east and west portions of the current species range and a third lineage restricted to a small area in the north of Sicily. Multiple lines of evidence from palaeogeographic (shorelines), palaeoclimatic (species distribution models), and multilocus genetic data (demographic and spatial Bayesian reconstructions) indicate that these lineages originated in distinct refugia, located in the north‐western and south‐eastern coastal lowlands, during Middle Pleistocene interglacial phases, and came into secondary contact following demographic and spatial expansions during the last glacial phase. Main conclusions This scenario of interglacial contraction and glacial expansion is in sharp contrast with patterns commonly observed in temperate species on the continent but parallels recent findings on other Mediterranean island endemics. Such a reverse expansion–contraction ( EC ) dynamic has been likely associated with glacial increases of climatically suitable coastal lowlands, suggesting this might be a general pattern in Mediterranean island species and also in other coastal regions strongly affected by glacial marine regressions during glacial episodes. This study provides explicit predictions and some methodological recommendations for testing the reverse EC model in other region and taxa.

Senczuk, Gabriele↗

DOE OSTI2018

Reduced Free Energy Perturbation/Hamiltonian Replica Exchange Molecular Dynamics Method with Unbiased Alchemical Thermodynamic Axis

We report that replica-exchange molecular dynamics (REMD) has been proven to efficiently improve the convergence of free energy perturbation (FEP) calculations involving considerable reorganization of their surrounding. We previously introduced the FEP/(λ,H)-REMD algorithm for ligand binding, in which replicas along the alchemical thermodynamic coupling axis λ were expanded as a series of Hamiltonian boosted replicas along a second axis to form a two-dimensional (2D) replica-exchange exchange map [Jiang, W.; Roux, B., J. Chem. Theory Comput. 2010, 6 (9), 2559-2565]. Aiming to achieve a similar performance at a lower computational cost, we propose here a modified version of this algorithm in which only the end-states along the alchemical axis are augmented by boosted replicas. The reduced FEP/(λ,H)-REMD method with one-dimensional (1D) unbiased alchemical thermodynamic coupling axis λ is implemented on the basis of generic multiple copy algorithm (MCA) module of the biomolecular simulation program NAMD. The flexible MCA framework of NAMD enables a user to design customized replica-exchange patterns through Tcl scripting in the context of a highly parallelized simulation program without touching the source code. Two Hamiltonian tempering boosting scheme were examined with the new algorithm: a first one based on potential energy rescaling of a pre-identified “solute”, and a second one via the introduction of flattening torsional free energy barriers. As two illustrative examples with reliable experiment data, the absolute binding free energies of pxylene and n-butylbenzene to the nonpolar cavity of the L99A mutant of T4 lysozyme were calculated. Lastly, the tests demonstrate that the new protocol efficiently enhances the sampling of torsional motions for backbone and side chains around the binding pocket and accelerates the convergence of the free energy computations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2018

Austenite and precipitation in secondary-hardening ultra-high-strength stainless steel

Highlights: • Austenite and precipitation in Ferrium S53 steel were analyzed. • Austenite decreased to less than 5% after cryogenic treatment and tempering. • M{sub 2}C carbide could segregate inside martensite lath or in the larger parallel microstructure. - Abstract: Ferrium S53 is a new secondary-hardening ultra-high-strength stainless steel, with more Cr than other Co-Ni ultra-high-strength steels. Austenite has a great effect on the fracture toughness of this metal, and the carbides that are precipitated are the main strengthening phase of Ferrium S53 steel; in this study, its microstructural changes during heat treatment were observed and analyzed by X-ray diffraction (XRD), Optical microscope (OM), SEM, and transmission electron microscopy (TEM). About 40.85% retained austenite in the quenched sample is decreased to less than 5% after cryogenic treatment and tempering. In the final sample, block austenite is retained austenite, and film austenite around 5 nm in width is reverted austenite. Inside the martensite lath, parallel twin-like structures are actually segregated M{sub 2}C precipitates, and most of them are smaller than 20 nm. These precipitates can be aggregated into a several-micron-long parallel zonal microstructure.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Genomic dissection of anthracnose resistance response in sorghum [Sorghum bicolor (L.) Moench]

Sorghum [Sorghum bicolor (L.) Moench] is the fifth most important grain crop behind maize, wheat, rice, and barley. Today, it is of interest as a source of fermentable sugars for the production of renewable fuels and chemicals, and as a source of biomass for co-firing. The productivity and profitability of sorghum are limited by several biotic constraints, most notably anthracnose caused by the fungal pathogen Colletotrichum sublineolum. The most cost-effective and environmentally benign strategy to control anthracnose is through the incorporation of resistance genes. Over the last three years, our research efforts have been directed to identify new sources of resistance in temperate adapted and tropical germplasm, and to delimited genomic regions associated with the observe anthracnose resistant response. Three biparental mapping populations derived from the resistant lines SC112-14, QL3 and IS18760 were evaluated for anthracnose resistance response in Texas, Georgia, Florida and Puerto Rico. In parallel, three high density recombination maps were constructed and used to identify resistant loci. Anthracnose resistant response in line SC112-14 is controlled by a major locus on chromosome 5. Segregation analysis of 1,500 progenies delimited the resistance locus on chromosome 5 to a 23-kb region harboring three candidate genes, including Sobic.005G17230 identified by GWAS of the sorghum association panel (SAP). The latter gene belongs to a family of genes encoding F-box proteins indicating that this resistance response involved in signaling cascades and transcriptional reprograming, rather than recognition of pathotype-associated molecular patterns. In contrast, anthracnose resistant response in lines QL3 and IS18760 is controlled by multiple small-effect genes. Greenhouse evaluation of a representative subset of the three mapping populations against nine pathotypes found that lines susceptible in the field could be resistant to a single pathotype in the greenhouse. Thus, the activation of a resistance response system by a single pathotype could not provide a broader resistance response against multiple pathotypes. The screening of 1,801 sweet sorghum accessions from the National Plant Germplasm System identified 654 accessions with Brix value larger than 10, which in turn was used to select a subset of 233 accessions for evaluation of anthracnose resistant response. Even though most of the accessions were not completely infected by anthracnose, 28 accessions were completely resistant against pathotypes from Texas, Georgia, Florida and Puerto Rico. Genotyping-by-sequencing analysis of this subset identified 157,843 single nucleotide polymorphisms. Population structure analysis of the subset based on a subset of 2,345 unlinked SNPs found that the genetic diversity could be divided into four populations. The genetic relatedness among accessions within populations suggests most of the resistant germplasm may contain few different resistance sources. These resistance sources present in sweet sorghum germplasm could expedite the development of new resistant sweet sorghum cultivars and hybrids by avoiding time-consuming introgression breeding approaches with non-sweet sorghums serving as donor of the resistance alleles.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2018

Genomic dissection of anthracnose resistance response in sorghum [ Sorghum bicolor (L.) Moench]

Sorghum [Sorghum bicolor (L.) Moench] is the fifth most important grain crop behind maize, wheat, rice, and barley. Today, it is of interest as a source of fermentable sugars for the production of renewable fuels and chemicals, and as a source of biomass for co-firing. The productivity and profitability of sorghum are limited by several biotic constraints, most notably anthracnose caused by the fungal pathogen Colletotrichum sublineolum. The most cost-effective and environmentally benign strategy to control anthracnose is through the incorporation of resistance genes. Over the last three years, our research efforts have been directed to identify new sources of resistance in temperate adapted and tropical germplasm, and to delimited genomic regions associated with the observe anthracnose resistant response. Three biparental mapping populations derived from the resistant lines SC112-14, QL3 and IS18760 were evaluated for anthracnose resistance response in Texas, Georgia, Florida and Puerto Rico. In parallel, three high density recombination maps were constructed and used to identify resistant loci. Anthracnose resistant response in line SC112-14 is controlled by a major locus on chromosome 5. Segregation analysis of 1,500 progenies delimited the resistance locus on chromosome 5 to a 23-kb region harboring three candidate genes, including Sobic.005G17230 identified by GWAS of the sorghum association panel (SAP). The latter gene belongs to a family of genes encoding F-box proteins indicating that this resistance response involved in signaling cascades and transcriptional reprograming, rather than recognition of pathotype-associated molecular patterns. In contrast, anthracnose resistant response in lines QL3 and IS18760 is controlled by multiple small-effect genes. Greenhouse evaluation of a representative subset of the three mapping populations against nine pathotypes found that lines susceptible in the field could be resistant to a single pathotype in the greenhouse. Thus, the activation of a resistance response system by a single pathotype could not provide a broader resistance response against multiple pathotypes. The screening of 1,801 sweet sorghum accessions from the National Plant Germplasm System identified 654 accessions with Brix value larger than 10, which in turn was used to select a subset of 233 accessions for evaluation of anthracnose resistant response. Even though most of the accessions were not completely infected by anthracnose, 28 accessions were completely resistant against pathotypes from Texas, Georgia, Florida and Puerto Rico. Genotyping-by-sequencing analysis of this subset identified 157,843 single nucleotide polymorphisms. Population structure analysis of the subset based on a subset of 2,345 unlinked SNPs found that the genetic diversity could be divided into four populations. The genetic relatedness among accessions within populations suggests most of the resistant germplasm may contain few different resistance sources. These resistance sources present in sweet sorghum germplasm could expedite the development of new resistant sweet sorghum cultivars and hybrids by avoiding time-consuming introgression breeding approaches with non-sweet sorghums serving as donor of the resistance alleles.

09 BIOMASS FUELS↗

DOE OSTI2018

Calibrated simulation analysis for integration of evaporative cooling and radiant cooling system for different Indian climatic zones

Radiant cooling system (RCS) has proven to be an energy efficient system for meeting a building's cooling requirements. RCS is energy efficient and provides better thermal comfort compared to conventional all-air heating ventilation air conditioning (HVAC) system. To further improve the efficiency of RCS, a parallel evaporative cooling system (cooling tower) is coupled with RCS and analyzed for different climatic condition using calibrated model. A statistical analysis of weather files (based on wet bulb temperature), was used to identify the availability of useful water for cooling tower integration with RCS. A comprehensive simulation feasibility study of the application of cooling tower in RCS was performed for different cities to cover every climatic zones of India. It was found that in summer, the wet bulb temperature (WBT) of different climatic zones, except warm-humid, is suitable for the integration of cooling tower with RCS. An experimental setup was designed and developed for integration of parallel chiller and cooling tower with panel-based RCS for cooling and a dedicated outdoor air system (DOAS) for dehumidification and ventilation. Experiments were conducted for chiller and cooling tower operated RCS in Composite climate of Jaipur, Rajasthan, India. Building Energy Models (BEM) were developed for both the chiller-operated RCS and cooling tower-operated RCS in EnergyPlus and calibrated with the measured data. Using the calibrated models, performance of the system was analyzed for different climatic zones of India. A chiller-operated RCS was considered as a baseline to compare the annual energy saving potential and monthly performance of the cooling tower integrated RCS. Here, in the cooling tower–operated RCS, a total annual savings of 7% in hot and dry climates, 11% in composite climates and 20% in temperate climates was achieved compared to the chiller-operated RCS.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Linking Dissolved Organic Matter Composition to Landscape Properties in Wetlands Across the United States of America

Abstract Wetlands are integral to the global carbon cycle, serving as both a source and a sink for organic carbon. Their potential for carbon storage will likely change in the coming decades in response to higher temperatures and variable precipitation patterns. We characterized the dissolved organic carbon (DOC) and dissolved organic matter (DOM) composition from 12 different wetland sites across the USA spanning gradients in climate, landcover, sampling depth, and hydroperiod for comparison to DOM in other inland waters. Using absorption spectroscopy, parallel factor analysis modeling, and ultra‐high resolution mass spectroscopy, we identified differences in DOM sourcing and processing by geographic site. Wetland DOM composition was driven primarily by differences in landcover where forested sites contained greater aromatic and oxygenated DOM content compared to grassland/herbaceous sites which were more aliphatic and enriched in N and S molecular formulae. Furthermore, surface and porewater DOM was also influenced by properties such as soil type, organic matter content, and precipitation. Surface water DOM was relatively enriched in oxygenated higher molecular weight formulae representing HUP High O/C compounds than porewaters, whose DOM composition suggests abiotic sulfurization from dissolved inorganic sulfide. Finally, we identified a group of persistent molecular formulae (3,489) present across all sites and sampling depths (i.e., the signature of wetland DOM) that are likely important for riverine‐to‐coastal DOM transport. As anthropogenic disturbances continue to impact temperate wetlands, this study highlights drivers of DOM composition fundamental for understanding how wetland organic carbon will change, and thus its role in biogeochemical cycling.

Environmental Sciences & Ecology↗

DOE OSTI2019

Winter hardiness of Miscanthus (III): Genome‐wide association and genomic prediction for overwintering ability in Miscanthus sinensis

Abstract Overwintering ability is an important selection criterion for Miscanthus breeding in temperate regions. Insufficient overwintering ability of the currently leading Miscanthus biomass cultivar, M . × giganteus (M×g) ‘1993–1780', in regions where average annual minimum temperatures are −26.1°C (USDA hardiness zone 5) or lower poses a pressing need to develop new cultivars with superior cold tolerance. To facilitate breeding of Miscanthus , this study characterized phenotypic and genetic variation of overwintering ability in an M. sinensis germplasm panel consisting of 564 accessions, evaluated in field trials at three locations in North America and two in Asia. Genome‐wide association (GWA) and genomic prediction analyses were performed. The Korea/N China M. sinensis genetic group is a valuable gene pool for cold tolerance. The Yangtze‐Qinling, Southern Japan, and Northern Japan genetic groups were also potential sources of cold tolerance. A total of 73 marker–trait associations were detected for overwintering ability. Estimated breeding value for overwintering ability based on these 73 markers could explain 55% of the variation for first winter overwintering ability among M. sinensis . Average genomic prediction ability for overwintering ability across 50 fivefold cross‐validations was high (~0.73) after accounting for population structure. Common genomic regions for overwintering ability were detected by GWA analyses and a previous parallel QTL mapping study using three interconnected biparental F 1 populations. One QTL on Miscanthus LG 8 encompassed five GWA hits and a known cold‐responsive gene, COR47 . The other overwintering ability QTL on Miscanthus LG 11 contained two GWA hits and three known cold stress‐related genes, carboxylesterase 13 ( CEX13 ), WRKY2 transcription factor, and cold shock domain ( CSDP1 ). Miscanthus accessions collected from high latitude locations with cold winters had higher rates of overwintering, and more alleles for overwintering, than accessions collected from southern locations with mild winters.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2019

Data from Winter Hardiness of Miscanthus (III): Genome‐Wide Association and Genomic Prediction for Overwintering Ability in Miscanthus sinensis

Overwintering ability is an important selection criterion for Miscanthus breeding in temperate regions. Insufficient overwintering ability of the currently leading Miscanthus biomass cultivar, M. ×giganteus (M×g) ‘1993–1780′, in regions where average annual minimum temperatures are −26.1°C (USDA hardiness zone 5) or lower poses a pressing need to develop new cultivars with superior cold tolerance. To facilitate breeding of Miscanthus , this study characterized phenotypic and genetic variation of overwintering ability in an M. sinensis germplasm panel consisting of 564 accessions, evaluated in field trials at three locations in North America and two in Asia. Genome‐wide association (GWA) and genomic prediction analyses were performed. The Korea/N China M. sinensis genetic group is a valuable gene pool for cold tolerance. The Yangtze‐Qinling, Southern Japan, and Northern Japan genetic groups were also potential sources of cold tolerance. A total of 73 marker–trait associations were detected for overwintering ability. Estimated breeding value for overwintering ability based on these 73 markers could explain 55% of the variation for first winter overwintering ability among M. sinensis . Average genomic prediction ability for overwintering ability across 50 fivefold cross‐validations was high (~0.73) after accounting for population structure. Common genomic regions for overwintering ability were detected by GWA analyses and a previous parallel QTL mapping study using three interconnected biparental F1 populations. One QTL on Miscanthus LG 8 encompassed five GWA hits and a known cold‐responsive gene, COR47 . The other overwintering ability QTL on Miscanthus LG 11 contained two GWA hits and three known cold stress‐related genes, carboxylesterase 13 (CEX13), WRKY2 transcription factor, and cold shock domain (CSDP1). Miscanthus accessions collected from high latitude locations with cold winters had higher rates of overwintering, and more alleles for overwintering, than accessions collected from southern locations with mild winters.

Biomass Analytics↗

DOE OSTI2020

Decreased Soil Organic Matter in a Long-Term Soil Warming Experiment Lowers Soil Water Holding Capacity and Affects Soil Thermal and Hydrological Buffering

Long-term soil warming can decrease soil organic matter (SOM), resulting in self-reinforcing feedback to the global climate system. We investigated additional consequences of SOM reduction for soil water holding capacity (WHC) and soil thermal and hydrological buffering. At a long-term soil warming experiment in a temperate forest in the northeastern United States, we suspended the warming treatment for 104 days during the summer of 2017. The formerly heated plot remained warmer (+0.39 °C) and drier (-0.024 cm 3 H 2 O cm -3 soil) than the control plot throughout the suspension. We measured decreased SOM content (-0.184 g SOM g -1 for O horizon soil, -0.010 g SOM g -1 for A horizon soil) and WHC (-0.82 g H 2 O g -1 for O horizon soil, -0.18 g H2O g -1 for A horizon soil) in the formerly heated plot relative to the control plot. Reduced SOM content accounted for 62% of the WHC reduction in the O horizon and 22% in the A horizon. We investigated differences in SOM composition as a possible explanation for the remaining reductions with Fourier transform infrared (FTIR) spectra. We found FTIR spectra that correlated more strongly with WHC than SOM, but those particular spectra did not differ between the heated and control plots, suggesting that SOM composition affects WHC but does not explain treatment differences in this study. We conclude that SOM reductions due to soil warming can reduce WHC and hydrological and thermal buffering, further warming soil and decreasing SOM. This feedback may operate in parallel, and perhaps synergistically, with carbon cycle feedbacks to climate change.

54 ENVIRONMENTAL SCIENCES↗