Oak Ridge National Laboratory

GROMACS is one of the most widely used open-source and free software codes in chemistry, used primarily for dynamical simulations of biomolecules. It provides a rich set of calculation types, preparation and analysis tools. Several advanced techniques for free-energy calculations are supported. In version 5, it reaches new performance heights, through several new and enhanced parallelization algorithms. These work on every level; SIMD registers inside cores, multithreading, heterogeneous CPU-GPU acceleration, state-of-the-art 3D domain decomposition, and ensemble-level parallelization through built-in replica exchange and the separate Copernicus framework. The latest best-in-class compressed trajectory storage format is supported.

BACKGROUND: The quality of automated gene prediction in microbial organisms has improved steadily over the past decade, but there is still room for improvement. Increasing the number of correct identifications, both of genes and of the translation initiation sites for each gene, and reducing the overall number of false positives, are all desirable goals. RESULTS: With our years of experience in manually curating genomes for the Joint Genome Institute, we developed a new gene prediction algorithm called Prodigal (PROkaryotic DYnamic programming Gene-finding ALgorithm). With Prodigal, we focused specifically on the three goals of improved gene structure prediction, improved translation initiation site recognition, and reduced false positives. We compared the results of Prodigal to existing gene-finding methods to demonstrate that it met each of these objectives. CONCLUSION: We built a fast, lightweight, open source gene prediction program called Prodigal http://compbio.ornl.gov/prodigal/. Prodigal achieved good results compared to existing methods, and we believe it will be a valuable asset to automated microbial annotation pipelines.

Studies of the human microbiome have revealed that even healthy individuals differ remarkably in the microbes that occupy habitats such as the gut, skin and vagina. Much of this diversity remains unexplained, although diet, environment, host genetics and early microbial exposure have all been implicated. Accordingly, to characterize the ecology of human-associated microbial communities, the Human Microbiome Project has analysed the largest cohort and set of distinct, clinically relevant body habitats so far. We found the diversity and abundance of each habitat’s signature microbes to vary widely even among healthy subjects, with strong niche specialization both within and among individuals. The project encountered an estimated 81–99% of the genera, enzyme families and community configurations occupied by the healthy Western microbiome. Metagenomic carriage of metabolic pathways was stable among individuals despite variation in community structure, and ethnic/racial background proved to be one of the strongest associations of both pathways and microbes with clinical metadata. These results thus delineate the range of structural and functional configurations normal in the microbial communities of a healthy population, enabling future characterization of the epidemiology, ecology and translational applications of the human microbiome. The Human Microbiome Project Consortium reports the first results of their analysis of microbial communities from distinct, clinically relevant body habitats in a human cohort; the insights into the microbial communities of a healthy population lay foundations for future exploration of the epidemiology, ecology and translational applications of the human microbiome. The Human Microbiome Project (HMP), supported by the National Institutes of Health Common Fund, has the goal of characterizing the microbial communities that inhabit and interact with the human body in sickness and in health. In two Articles in this issue of Nature, the HMP Consortium presents the first population-scale details of the organismal and functional composition of the microbiota across five areas of the body. An associated News & Views discusses the initial results — which, along with those of a series of co-publications, already constitute the most extensive catalogue of organisms and genes related to the human microbiome yet published — and highlights some of the major questions that the project will tackle in the next few years.

MOTIVATION: Molecular simulation has historically been a low-throughput technique, but faster computers and increasing amounts of genomic and structural data are changing this by enabling large-scale automated simulation of, for instance, many conformers or mutants of biomolecules with or without a range of ligands. At the same time, advances in performance and scaling now make it possible to model complex biomolecular interaction and function in a manner directly testable by experiment. These applications share a need for fast and efficient software that can be deployed on massive scale in clusters, web servers, distributed computing or cloud resources. RESULTS: Here, we present a range of new simulation algorithms and features developed during the past 4 years, leading up to the GROMACS 4.5 software package. The software now automatically handles wide classes of biomolecules, such as proteins, nucleic acids and lipids, and comes with all commonly used force fields for these molecules built-in. GROMACS supports several implicit solvent models, as well as new free-energy algorithms, and the software now uses multithreading for efficient parallelization even on low-end systems, including windows-based workstations. Together with hand-tuned assembly kernels and state-of-the-art parallelization, this provides extremely high performance and cost efficiency for high-throughput as well as massively parallel simulations. AVAILABILITY: GROMACS is an open source and free software available from http://www.gromacs.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The promise of quantum computers is that certain computational tasks might be executed exponentially faster on a quantum processor than on a classical processor1. A fundamental challenge is to build a high-fidelity processor capable of running quantum algorithms in an exponentially large computational space. Here we report the use of a processor with programmable superconducting qubits2–7 to create quantum states on 53 qubits, corresponding to a computational state-space of dimension 253 (about 1016). Measurements from repeated experiments sample the resulting probability distribution, which we verify using classical simulations. Our Sycamore processor takes about 200 seconds to sample one instance of a quantum circuit a million times—our benchmarks currently indicate that the equivalent task for a state-of-the-art classical supercomputer would take approximately 10,000 years. This dramatic increase in speed compared to all known classical algorithms is an experimental realization of quantum supremacy8–14 for this specific computational task, heralding a much-anticipated computing paradigm. Quantum supremacy is demonstrated using a programmable superconducting processor known as Sycamore, taking approximately 200 seconds to sample one instance of a quantum circuit a million times, which would take a state-of-the-art supercomputer around ten thousand years to compute.

High-entropy alloys are equiatomic, multi-element systems that can crystallize as a single phase, despite containing multiple elements with different crystal structures. A rationale for this is that the configurational entropy contribution to the total free energy in alloys with five or more major elements may stabilize the solid-solution state relative to multiphase microstructures. We examined a five-element high-entropy alloy, CrMnFeCoNi, which forms a single-phase face-centered cubic solid solution, and found it to have exceptional damage tolerance with tensile strengths above 1 GPa and fracture toughness values exceeding 200 MPa·m(1/2). Furthermore, its mechanical properties actually improve at cryogenic temperatures; we attribute this to a transition from planar-slip dislocation activity at room temperature to deformation by mechanical nanotwinning with decreasing temperature, which results in continuous steady strain hardening.

Ever since its discovery the notion of Berry phase has permeated through all branches of physics. Over the past three decades it was gradually realized that the Berry phase of the electronic wave function can have a profound effect on material properties and is responsible for a spectrum of phenomena, such as polarization, orbital magnetism, various (quantum, anomalous, or spin) Hall effects, and quantum charge pumping. This progress is summarized in a pedagogical manner in this review. A brief summary of necessary background is given and a detailed discussion of the Berry phase effect in a variety of solid-state applications. A common thread of the review is the semiclassical formulation of electron dynamics, which is a versatile tool in the study of electron dynamics in the presence of electromagnetic fields and more general perturbations. Finally, a requantization method is demonstrated that converts a semiclassical theory to an effective quantum theory. It is clear that the Berry phase should be added as an essential ingredient to our understanding of basic material properties.

We show that inversion symmetry breaking together with spin-orbit coupling leads to coupled spin and valley physics in monolayers of MoS2 and other group-VI dichalcogenides, making possible controls of spin and valley in these 2D materials. The spin-valley coupling at the valence-band edges suppresses spin and valley relaxation, as flip of each index alone is forbidden by the valley-contrasting spin splitting. Valley Hall and spin Hall effects coexist in both electron-doped and hole-doped systems. Optical interband transitions have frequency-dependent polarization selection rules which allow selective photoexcitation of carriers with various combination of valley and spin indices. Photoinduced spin Hall and valley Hall effects can generate long lived spin and valley accumulations on sample boundaries. The physics discussed here provides a route towards the integration of valleytronics and spintronics in multivalley materials with strong spin-orbit coupling and inversion symmetry breaking.

Similar to its popular older cousins the fullerene, the carbon nanotube, and graphene, the latest form of nanocarbon, the carbon nanodot, is inspiring intensive research efforts in its own right. These surface-passivated carbonaceous quantum dots, so-called C-dots, combine several favorable attributes of traditional semiconductor-based quantum dots (namely, size- and wavelength-dependent luminescence emission, resistance to photobleaching, ease of bioconjugation) without incurring the burden of intrinsic toxicity or elemental scarcity and without the need for stringent, intricate, tedious, costly, or inefficient preparation steps. C-dots can be produced inexpensively and on a large scale (frequently using a one-step pathway and potentially from biomass waste-derived sources) by many approaches, ranging from simple candle burning to in situ dehydration reactions to laser ablation methods. In this Review, we summarize recent advances in the synthesis and characterization of C-dots. We also speculate on their future and discuss potential developments for their use in energy conversion/storage, bioimaging, drug delivery, sensors, diagnostics, and composites.

The ecological consequences of biodiversity loss have aroused considerable interest and controversy during the past decade. Major advances have been made in describing the relationship between species diversity and ecosystem processes, in identifying functionally important species, and in revealing underlying mechanisms. There is, however, uncertainty as to how results obtained in recent experiments scale up to landscape and regional levels and generalize across ecosystem types and processes. Larger numbers of species are probably needed to reduce temporal variability in ecosystem processes in changing environments. A major future challenge is to determine how biodiversity dynamics, ecosystem processes, and abiotic factors interact.

We report the draft genome of the black cottonwood tree, Populus trichocarpa. Integration of shotgun sequence assembly with genetic mapping enabled chromosome-scale reconstruction of the genome. More than 45,000 putative protein-coding genes were identified. Analysis of the assembled genome revealed a whole-genome duplication event; about 8000 pairs of duplicated genes from that event survived in the Populus genome. A second, older duplication event is indistinguishably coincident with the divergence of the Populus and Arabidopsis lineages. Nucleotide substitution, tandem gene duplication, and gross chromosomal rearrangement appear to proceed substantially more slowly in Populus than in Arabidopsis. Populus has more protein-coding genes than Arabidopsis, ranging on average from 1.4 to 1.6 putative Populus homologs for each Arabidopsis gene. However, the relative frequency of protein domains in the two genomes is similar. Overrepresented exceptions in Populus include genes associated with lignocellulosic wall biosynthesis, meristem development, disease resistance, and metabolite transport.

Background Lignin, nature’s dominant aromatic polymer, is found in most terrestrial plants in the approximate range of 15 to 40% dry weight and provides structural integrity. Traditionally, most large-scale industrial processes that use plant polysaccharides have burned lignin to generate the power needed to productively transform biomass. The advent of biorefineries that convert cellulosic biomass into liquid transportation fuels will generate substantially more lignin than necessary to power the operation, and therefore efforts are underway to transform it to value-added products. Advances Bioengineering to modify lignin structure and/or incorporate atypical components has shown promise toward facilitating recovery and chemical transformation of lignin under biorefinery conditions. The flexibility in lignin monomer composition has proven useful for enhancing extraction efficiency. Both the mining of genetic variants in native populations of bioenergy crops and direct genetic manipulation of biosynthesis pathways have produced lignin feedstocks with unique properties for coproduct development. Advances in analytical chemistry and computational modeling detail the structure of the modified lignin and direct bioengineering strategies for targeted properties. Refinement of biomass pretreatment technologies has further facilitated lignin recovery and enables catalytic modifications for desired chemical and physical properties. Outlook Potential high-value products from isolated lignin include low-cost carbon fiber, engineering plastics and thermoplastic elastomers, polymeric foams and membranes, and a variety of fuels and chemicals all currently sourced from petroleum. These lignin coproducts must be low cost and perform as well as petroleum-derived counterparts. Each product stream has its own distinct challenges. Development of renewable lignin-based polymers requires improved processing technologies coupled to tailored bioenergy crops incorporating lignin with the desired chemical and physical properties. For fuels and chemicals, multiple strategies have emerged for lignin depolymerization and upgrading, including thermochemical treatments and homogeneous and heterogeneous catalysis. The multifunctional nature of lignin has historically yielded multiple product streams, which require extensive separation and purification procedures, but engineering plant feedstocks for greater structural homogeneity and tailored functionality reduces this challenge.

FLUXNET is a global network of micrometeorological flux measurement sites that measure the exchanges of carbon dioxide, water vapor, and energy between the biosphere and atmosphere. At present over 140 sites are operating on a long-term and continuous basis. Vegetation under study includes temperate conifer and broadleaved (deciduous and evergreen) forests, tropical and boreal forests, crops, grasslands, chaparral, wetlands, and tundra. Sites exist on five continents and their latitudinal distribution ranges from 70°N to 30°S. FLUXNET has several primary functions. First, it provides infrastructure for compiling, archiving, and distributing carbon, water, and energy flux measurement, and meteorological, plant, and soil data to the science community. (Data and site information are available online at the FLUXNET Web site, http://www-eosdis.ornl.gov/FLUXNET/.) Second, the project supports calibration and flux intercomparison activities. This activity ensures that data from the regional networks are intercomparable. And third, FLUXNET supports the synthesis, discussion, and communication of ideas and data by supporting project scientists, workshops, and visiting scientists. The overarching goal is to provide information for validating computations of net primary productivity, evaporation, and energy absorption that are being generated by sensors mounted on the NASA Terra satellite. Data being compiled by FLUXNET are being used to quantify and compare magnitudes and dynamics of annual ecosystem carbon and water balances, to quantify the response of stand-scale carbon dioxide and water vapor flux densities to controlling biotic and abiotic factors, and to validate a hierarchy of soil–plant–atmosphere trace gas exchange models. Findings so far include 1) net CO2 exchange of temperate broadleaved forests increases by about 5.7 g C m−2 day−1 for each additional day that the growing season is extended; 2) the sensitivity of net ecosystem CO2 exchange to sunlight doubles if the sky is cloudy rather than clear; 3) the spectrum of CO2 flux density exhibits peaks at timescales of days, weeks, and years, and a spectral gap exists at the month timescale; 4) the optimal temperature of net CO2 exchange varies with mean summer temperature; and 5) stand age affects carbon dioxide and water vapor flux densities.

The prohibitive cost of platinum for catalyzing the cathodic oxygen reduction reaction (ORR) has hampered the widespread use of polymer electrolyte fuel cells. We describe a family of non-precious metal catalysts that approach the performance of platinum-based systems at a cost sustainable for high-power fuel cell applications, possibly including automotive power. The approach uses polyaniline as a precursor to a carbon-nitrogen template for high-temperature synthesis of catalysts incorporating iron and cobalt. The most active materials in the group catalyze the ORR at potentials within ~60 millivolts of that delivered by state-of-the-art carbon-supported platinum, combining their high activity with remarkable performance stability for non-precious metal catalysts (700 hours at a fuel cell voltage of 0.4 volts) as well as excellent four-electron selectivity (hydrogen peroxide yield <1.0%).

A computational model is a representation of some physical or other system of interest, first expressed mathematically and then implemented in the form of a computer program; it may be viewed as a function of inputs that, when evaluated, produces outputs. Motivation for this article comes from computational models that are deterministic, complicated enough to make classical mathematical analysis impractical and that have a moderate-to-large number of inputs. The problem of designing computational experiments to determine which inputs have important effects on an output is considered. The proposed experimental plans are composed of individually randomized one-factor-at-a-time designs, and data analysis is based on the resulting random sample of observed elementary effects, those changes in an output due solely to changes in a particular input. Advantages of this approach include a lack of reliance on assumptions of relative sparsity of important inputs, monotonicity of outputs with respect to inputs, or adequacy of a low-order polynomial as an approximation to the computational model.

Consider the ridge estimate (λ) for β in the model unknown, (λ) = (X T X + nλI)−1 X T y. We study the method of generalized cross-validation (GCV) for choosing a good value for λ from the data. The estimate is the minimizer of V(λ) given by where A(λ) = X(X T X + nλI)−1 X T . This estimate is a rotation-invariant version of Allen's PRESS, or ordinary cross-validation. This estimate behaves like a risk improvement estimator, but does not require an estimate of σ2, so can be used when n − p is small, or even if p ≥ 2 n in certain cases. The GCV method can also be used in subset selection and singular value truncation methods for regression, and even to choose from among mixtures of these methods.

1. Crystal lattices. General theory 2. . Crystal lattices. Applications 3. Interaction of light with non-conducting crystals 4. Electrons in a perfect lattice 5. Cohesive forces in metals 6. Transport phenomena 7. Magnetic properties of metals 8. Ferromagnetism 9. Interaction of light with electrons in solids 10. Semi-conductors and luminescence 11. Superconductivity

The fourth version of the Community Climate System Model (CCSM4) was recently completed and released to the climate community. This paper describes developments to all CCSM components, and documents fully coupled preindustrial control runs compared to the previous version, CCSM3. Using the standard atmosphere and land resolution of 1° results in the sea surface temperature biases in the major upwelling regions being comparable to the 1.4°-resolution CCSM3. Two changes to the deep convection scheme in the atmosphere component result in CCSM4 producing El Niño–Southern Oscillation variability with a much more realistic frequency distribution than in CCSM3, although the amplitude is too large compared to observations. These changes also improve the Madden–Julian oscillation and the frequency distribution of tropical precipitation. A new overflow parameterization in the ocean component leads to an improved simulation of the Gulf Stream path and the North Atlantic Ocean meridional overturning circulation. Changes to the CCSM4 land component lead to a much improved annual cycle of water storage, especially in the tropics. The CCSM4 sea ice component uses much more realistic albedos than CCSM3, and for several reasons the Arctic sea ice concentration is improved in CCSM4. An ensemble of twentieth-century simulations produces a good match to the observed September Arctic sea ice extent from 1979 to 2005. The CCSM4 ensemble mean increase in globally averaged surface temperature between 1850 and 2005 is larger than the observed increase by about 0.4°C. This is consistent with the fact that CCSM4 does not include a representation of the indirect effects of aerosols, although other factors may come into play. The CCSM4 still has significant biases, such as the mean precipitation distribution in the tropical Pacific Ocean, too much low cloud in the Arctic, and the latitudinal distributions of shortwave and longwave cloud forcings.

A Weyl semimetal is a new state of matter that hosts Weyl fermions as emergent quasiparticles and admits a topological classification that protects Fermi arc surface states on the boundary of a bulk sample. This unusual electronic structure has deep analogies with particle physics and leads to unique topological properties. We report the experimental discovery of a Weyl semimetal, tantalum arsenide (TaAs). Using photoemission spectroscopy, we directly observe Fermi arcs on the surface, as well as the Weyl fermion cones and Weyl nodes in the bulk of TaAs single crystals. We find that Fermi arcs terminate on the Weyl fermion nodes, consistent with their topological character. Our work opens the field for the experimental study of Weyl fermions in physics and materials science.

This report identifies twelve building block chemicals that can be produced from sugars via biological or chemical conversions. The twelve building blocks can be subsequently converted to a number of high-value bio-based chemicals or materials. Building block chemicals, as considered for this analysis, are molecules with multiple functional groups that possess the potential to be transformed into new families of useful molecules. The twelve sugar-based building blocks are 1,4-diacids (succinic, fumaric and malic), 2,5-furan dicarboxylic acid, 3-hydroxy propionic acid, aspartic acid, glucaric acid, glutamic acid, itaconic acid, levulinic acid, 3-hydroxybutyrolactone, glycerol, sorbitol, and xylitol/arabinitol.