Daniel O.B. Jones

Reproducibility in density functional theory calculations of solids

A comparison of DFT methods Density functional theory (DFT) is now routinely used for simulating material properties. Many software packages are available, which makes it challenging to know which are the best to use for a specific calculation. Lejaeghere et al. compared the calculated values for the equation of states for 71 elemental crystals from 15 different widely used DFT codes employing 40 different potentials (see the Perspective by Skylaris). Although there were variations in the calculated values, most recent codes and methods converged toward a single value, with errors comparable to those of experiment. Science, this issue p. 10.1126/science.aad3000; see also p. 1394 A survey of recent density functional theory methods shows a convergence to more accurate property calculations. [Also see Perspective by Skylaris] INTRODUCTION The reproducibility of results is one of the underlying principles of science. An observation can only be accepted by the scientific community when it can be confirmed by independent studies. However, reproducibility does not come easily. Recent works have painfully exposed cases where previous conclusions were not upheld. The scrutiny of the scientific community has also turned to research involving computer programs, finding that reproducibility depends more strongly on implementation than commonly thought. These problems are especially relevant for property predictions of crystals and molecules, which hinge on precise computer implementations of the governing equation of quantum physics. RATIONALE This work focuses on density functional theory (DFT), a particularly popular quantum method for both academic and industrial applications. More than 15,000 DFT papers are published each year, and DFT is now increasingly used in an automated fashion to build large databases or apply multiscale techniques with limited human supervision. Therefore, the reproducibility of DFT results underlies the scientific credibility of a substantial fraction of current work in the natural and engineering sciences. A plethora of DFT computer codes are available, many of them differing considerably in their details of implementation, and each yielding a certain “precision” relative to other codes. How is one to decide for more than a few simple cases which code predicts the correct result, and which does not? We devised a procedure to assess the precision of DFT methods and used this to demonstrate reproducibility among many of the most widely used DFT codes. The essential part of this assessment is a pairwise comparison of a wide range of methods with respect to their predictions of the equations of state of the elemental crystals. This effort required the combined expertise of a large group of code developers and expert users. RESULTS We calculated equation-of-state data for four classes of DFT implementations, totaling 40 methods. Most codes agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Even in the case of pseudization approaches, which largely depend on the atomic potentials used, a similar precision can be obtained as when using the full potential. The remaining deviations are due to subtle effects, such as specific numerical implementations or the treatment of relativistic terms. CONCLUSION Our work demonstrates that the precision of DFT implementations can be determined, even in the absence of one absolute reference code. Although this was not the case 5 to 10 years ago, most of the commonly used codes and methods are now found to predict essentially identical results. The established precision of DFT codes not only ensures the reproducibility of DFT predictions but also puts several past and future developments on a firmer footing. Any newly developed methodology can now be tested against the benchmark to verify whether it reaches the same level of precision. New DFT applications can be shown to have used a sufficiently precise method. Moreover, high-precision DFT calculations are essential for developing improvements to DFT methodology, such as new density functionals, which may further increase the predictive power of the simulations. Recent DFT methods yield reproducible results. Whereas older DFT implementations predict different values (red darts), codes have now evolved to mutual agreement (green darts). The scoreboard illustrates the good pairwise agreement of four classes of DFT implementations (horizontal direction) with all-electron results (vertical direction). Each number reflects the average difference between the equations of state for a given pair of methods, with the green-to-red color scheme showing the range from the best to the poorest agreement. The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Does Presence of a Mid-Ocean Ridge Enhance Biomass and Biodiversity?

In contrast to generally sparse biological communities in open-ocean settings, seamounts and ridges are perceived as areas of elevated productivity and biodiversity capable of supporting commercial fisheries. We investigated the origin of this apparent biological enhancement over a segment of the North Mid-Atlantic Ridge (MAR) using sonar, corers, trawls, traps, and a remotely operated vehicle to survey habitat, biomass, and biodiversity. Satellite remote sensing provided information on flow patterns, thermal fronts, and primary production, while sediment traps measured export flux during 2007–2010. The MAR, 3,704,404 km2 in area, accounts for 44.7% lower bathyal habitat (800–3500 m depth) in the North Atlantic and is dominated by fine soft sediment substrate (95% of area) on a series of flat terraces with intervening slopes either side of the ridge axis contributing to habitat heterogeneity. The MAR fauna comprises mainly species known from continental margins with no evidence of greater biodiversity. Primary production and export flux over the MAR were not enhanced compared with a nearby reference station over the Porcupine Abyssal Plain. Biomasses of benthic macrofauna and megafauna were similar to global averages at the same depths totalling an estimated 258.9 kt C over the entire lower bathyal north MAR. A hypothetical flat plain at 3500 m depth in place of the MAR would contain 85.6 kt C, implying an increase of 173.3 kt C attributable to the presence of the Ridge. This is approximately equal to 167 kt C of estimated pelagic biomass displaced by the volume of the MAR. There is no enhancement of biological productivity over the MAR; oceanic bathypelagic species are replaced by benthic fauna otherwise unable to survive in the mid ocean. We propose that globally sea floor elevation has no effect on deep sea biomass; pelagic plus benthic biomass is constant within a given surface productivity regime.

Marine litter distribution and density in European seas, from the shelves to deep basins

Anthropogenic litter is present in all marine habitats, from beaches to the most remote points in the oceans. On the seafloor, marine litter, particularly plastic, can accumulate in high densities with deleterious consequences for its inhabitants. Yet, because of the high cost involved with sampling the seafloor, no large-scale assessment of distribution patterns was available to date. Here, we present data on litter distribution and density collected during 588 video and trawl surveys across 32 sites in European waters. We found litter to be present in the deepest areas and at locations as remote from land as the Charlie-Gibbs Fracture Zone across the Mid-Atlantic Ridge. The highest litter density occurs in submarine canyons, whilst the lowest density can be found on continental shelves and on ocean ridges. Plastic was the most prevalent litter item found on the seafloor. Litter from fishing activities (derelict fishing lines and nets) was particularly common on seamounts, banks, mounds and ocean ridges. Our results highlight the extent of the problem and the need for action to prevent increasing accumulation of litter in marine environments.

Biotic and Human Vulnerability to Projected Changes in Ocean Biogeochemistry over the 21st Century

Mora and colleagues show that ongoing greenhouse gas emissions are likely to have a considerable effect on several biogeochemical properties of the worlds oceans, with potentially serious consequences for biodiversity and human welfare.

Global reductions in seafloor biomass in response to climate change

Seafloor organisms are vital for healthy marine ecosystems, contributing to elemental cycling, benthic remineralization, and ultimately sequestration of carbon. Deep-sea life is primarily reliant on the export flux of particulate organic carbon from the surface ocean for food, but most ocean biogeochemistry models predict global decreases in export flux resulting from 21st century anthropogenically induced warming. Here we show that decadal-to-century scale changes in carbon export associated with climate change lead to an estimated 5.2% decrease in future (2091–2100) global open ocean benthic biomass under RCP8.5 (reduction of 5.2 Mt C) compared with contemporary conditions (2006–2015). Our projections use multi-model mean export flux estimates from eight fully coupled earth system models, which contributed to the Coupled Model Intercomparison Project Phase 5, that have been forced by high and low representative concentration pathways (RCP8.5 and 4.5, respectively). These export flux estimates are used in conjunction with published empirical relationships to predict changes in benthic biomass. The polar oceans and some upwelling areas may experience increases in benthic biomass, but most other regions show decreases, with up to 38% reductions in parts of the northeast Atlantic. Our analysis projects a future ocean with smaller sized infaunal benthos, potentially reducing energy transfer rates though benthic multicellular food webs. More than 80% of potential deep-water biodiversity hotspots known around the world, including canyons, seamounts, and cold-water coral reefs, are projected to experience negative changes in biomass. These major reductions in biomass may lead to widespread change in benthic ecosystems and the functions and services they provide.

Jelly-falls historic and recent observations: a review to drive future research directions

The biological pump describes the transport of particulate matter from the sea surface to the ocean’s interior including the seabed. The contribution by gelatinous zooplankton bodies as particulate organic matter (POM) vectors (“jelly-falls”) has been neglected owing to technical and spatiotemporal sampling limitations. Here, we assess the existing evidence on jelly-falls from early ocean observations to present times. The seasonality of jelly-falls indicates that they mostly occur after periods of strong upwelling and/or spring blooms in temperate/subpolar zones and during late spring/early summer. A conceptual model helps to define a jelly-fall based on empirical and field observations of biogeochemical and ecological processes. We then compile and discuss existing strategic and observational oceanographic techniques that could be implemented to further jelly-falls research. Seabed video- and photography-based studies deliver the best results, and the correct use of fishing techniques, such as trawling, could provide comprehensive regional datasets. We conclude by considering the possibility of increased gelatinous biomasses in the future ocean induced by upper ocean processes favouring their populations, thus increasing jelly-POM downward transport. We suggest that this could provide a “natural compensation” for predicted losses in pelagic POM with respect to fuelling benthic ecosystems.

Single-molecule decoding of combinatorially modified nucleosomes

Deciphering the histone code Regulatory information is stored both in DNA sequences and in the chromatin proteins that package the genome. The covalent modification of histones plays a critical role in signaling whether a gene or genomic region should be active or inactive. Shema et al. used high-throughput single-molecule imaging to reveal the different combinations of modifications on millions of individual nucleosomes. Single-molecule DNA sequencing determined the exact location of these modified nucleosomes in the genome. Science, this issue p. 717 Determining the location and epigenetic modification state of individual nucleosomes will help unpick the “histone code.” Different combinations of histone modifications have been proposed to signal distinct gene regulatory functions, but this area is poorly addressed by existing technologies. We applied high-throughput single-molecule imaging to decode combinatorial modifications on millions of individual nucleosomes from pluripotent stem cells and lineage-committed cells. We identified definitively bivalent nucleosomes with concomitant repressive and activating marks, as well as other combinatorial modification states whose prevalence varies with developmental potency. We showed that genetic and chemical perturbations of chromatin enzymes preferentially affect nucleosomes harboring specific modification states. Last, we combined this proteomic platform with single-molecule DNA sequencing technology to simultaneously determine the modification states and genomic positions of individual nucleosomes. This single-molecule technology has the potential to address fundamental questions in chromatin biology and epigenetic regulation.

Environmental Impacts of the Deep-Water Oil and Gas Industry: A Review to Guide Management Strategies

The industrialization of the deep sea is expanding worldwide. Expanding oil and gas exploration activities in the absence of sufficient baseline data in these ecosystems has made environmental management challenging. Here, we review the types of activities that are associated with global offshore oil and gas development in water depths over 200 m, the typical impacts of these activities, some of the more extreme impacts of accidental oil and gas releases, and the current state of management in the major regions of offshore industrial activity including 18 exclusive economic zones. Direct impacts of infrastructure installation, including sediment resuspension and burial by seafloor anchors and pipelines, are typically restricted to a radius of approximately 100 m on from the installation on the seafloor. Discharges of water-based and low-toxicity oil-based drilling muds and produced water can extend over 2 km, while the ecological impacts at the population and community levels on the seafloor are most commonly on the order of 200-300 m from their source. These impacts may persist in the deep sea for many years and likely longer for its more fragile ecosystems, such as cold-water corals. This synthesis of information provides the basis for a series of recommendations for the management of offshore oil and gas development. An effective management strategy, aimed at minimizing risk of significant environmental harm, will typically encompass regulations of the activity itself (e.g. discharge practices, materials used), combined with spatial (e.g. avoidance rules and marine protected areas) and temporal measures (e.g. restricted activities during peak reproductive periods). Spatial management measures that encompass representatives of all of the regional deep-sea community types is important in this context. Implementation of these management strategies should consider minimum buffer zones to displace industrial activity beyond the range of typical impacts: at least 2 km from any discharge points and surface infrastructure and 200 m from seafloor infrastructure with no expected discharges. Although managing natural resources is, arguably, more challenging in deep-water environments, inclusion of these proven conservation tools contributes to robust environmental management strategies for oil and gas extraction in the deep sea.

Rapid scavenging of jellyfish carcasses reveals the importance of gelatinous material to deep-sea food webs

Jellyfish blooms are common in many oceans, and anthropogenic changes appear to have increased their magnitude in some regions. Although mass falls of jellyfish carcasses have been observed recently at the deep seafloor, the dense necrophage aggregations and rapid consumption rates typical for vertebrate carrion have not been documented. This has led to a paradigm of limited energy transfer to higher trophic levels at jelly falls relative to vertebrate organic falls. We show from baited camera deployments in the Norwegian deep sea that dense aggregations of deep-sea scavengers (more than 1000 animals at peak densities) can rapidly form at jellyfish baits and consume entire jellyfish carcasses in 2.5 h. We also show that scavenging rates on jellyfish are not significantly different from fish carrion of similar mass, and reveal that scavenging communities typical for the NE Atlantic bathyal zone, including the Atlantic hagfish, galatheid crabs, decapod shrimp and lyssianasid amphipods, consume both types of carcasses. These rapid jellyfish carrion consumption rates suggest that the contribution of gelatinous material to organic fluxes may be seriously underestimated in some regions, because jelly falls may disappear much more rapidly than previously thought. Our results also demonstrate that the energy contained in gelatinous carrion can be efficiently incorporated into large numbers of deep-sea scavengers and food webs, lessening the expected impacts (e.g. smothering of the seafloor) of enhanced jellyfish production on deep-sea ecosystems and pelagic–benthic coupling.

Autonomous underwater vehicles (AUVs) and investigations of the ice-ocean interface in Antarctic and Arctic waters

Limitations of access have long restricted exploration and investigation of the cavities beneath ice shelves to a small number of drillholes. Studies of sea-ice underwater morphology are limited largely to scientific utilization of submarines. Remotely operated vehicles, tethered to a mother ship by umbilical cable, have been deployed to investigate tidewater-glacier and ice-shelf margins, but their range is often restricted. The development of free-flying autonomous underwater vehicles (AUVs) with ranges of tens to hundreds of kilometres enables extensive missions to take place beneath sea ice and floating ice shelves. Autosub2 is a 3600 kg, 6.7 m long AUV, with a 1600 m operating depth and range of 400 km, based on the earlier Autosub1 which had a 500m depth limit. A single direct-drive d.c. motor and five-bladed propeller produce speeds of 1-2 ms−1. Rear-mounted rudder and stern-plane control yaw, pitch and depth. The vehicle has three sections. The front and rear sections are free-flooding, built around aluminium extrusion space-frames covered with glass-fibre reinforced plastic panels. The central section has a set of carbon-fibre reinforced plastic pressure vessels. Four tubes contain batteries powering the vehicle. The other three house vehicle-control systems and sensors. The rear section houses subsystems for navigation, control actuation and propulsion and scientific sensors (e.g. digital camera, upward-looking 300 kHz acoustic Doppler current profiler, 200 kHz multibeam receiver). The front section contains forward-looking collision sensor, emergency abort, the homing systems, Argos satellite data and location transmitters and flashing lights for relocation as well as science sensors (e.g. twin conductivity-temperature-depth instruments, multibeam transmitter, sub-bottom profiler, AquaLab water sampler). Payload restrictions mean that a subset of scientific instruments is actually in place on any given dive. The scientific instruments carried on Autosub are described and examples of observational data collected from each sensor in Arctic or Antarctic waters are given (e.g. of roughness at the underside of floating ice shelves and sea ice).