A US perspective on closing the carbon cycle to defossilize difficult-to-electrify segments of our economy.

Electrification to reduce or eliminate greenhouse gas emissions is essential to mitigate climate change. However, a substantial portion of our manufacturing and transportation infrastructure will be difficult to electrify and/or will continue to use carbon as a key component, including areas in aviation, heavy-duty and marine transportation, and the chemical industry. In this Roadmap, we explore how multidisciplinary approaches will enable us to close the carbon cycle and create a circular economy by defossilizing these difficult-to-electrify areas and those that will continue to need carbon. We discuss two approaches for this: developing carbon alternatives and improving our ability to reuse carbon, enabled by separations. Furthermore, we posit that co-design and use-driven fundamental science are essential to reach aggressive greenhouse gas reduction targets.

Inducing thermodynamically blocked atomic ordering via strongly driven nonequilibrium kinetics.

Ultrafast light-matter interactions enable inducing exotic material phases by promoting access to kinetic processes blocked in equilibrium. Despite potential opportunities, actively using nonequilibrium kinetics for material discovery is limited by the poor understanding on intermediate states of driven systems. Here, using single-pulse time-resolved imaging with x-ray free-electron lasers, we found intermediate states of photoexcited bismuth nanoparticles that showed kinetically reversed surface ordering during ultrafast melting. This entropy-lowering reaction was further investigated by molecular dynamics simulations to reveal that observed kinetics were thermodynamically buried in equilibrium, which emphasized the critical role of electron-mediated ultrafast free-energy modification in inducing exotic material phases. This study demonstrated that ultrafast photoexcitations of electrons provide an efficient strategy to induce hidden material phases by overcoming thermodynamic barriers via nonequilibrium reaction pathways.

Phase Transition of H2 in Subnanometer Pores Observed at 75 K.

Here we report a phase transition in H2 adsorbed in a locally graphitic Saran carbon with subnanometer pores 0.5-0.65 nm in width, in which two layers of hydrogen can just barely squeeze, provided they pack tightly. The phase transition is observed at 75 K, temperatures far higher than other systems in which an adsorbent is known to increase phase transition temperatures: for instance, H2 melts at 14 K in the bulk, but at 20 K on graphite because the solid H2 is stabilized by the surface structure. Here we observe a transition at 75 K and 77-200 bar: from a low-temperature, low-density phase to a high-temperature, higher density phase. We model the low-density phase as a monolayer commensurate solid composed mostly of para-H2 (the ground nuclear spin state, S = 0) and the high-density phase as an orientationally ordered bilayer commensurate solid composed mostly of ortho-H2 (S = 1). We attribute the increase in density with temperature to the fact that the oblong ortho-H2 can pack more densely. The transition is observed using two experiments. The high-density phase is associated with an increase in neutron backscatter by a factor of 7.0 ± 0.1. Normally, hydrogen produces no backscatter (scattering angle >90°). This backscatter appears along with a discontinuous increase in the excitation mass from 1.2 amu to 21.0 ± 2.3 amu, which we associate with collective nuclear spin excitations in the orientationally ordered phase. Film densities were measured using hydrogen adsorption. No phase transition was observed in H2 adsorbed in control activated carbon materials.

Adaptations: Using Darwin's Origin to teach biology and writing.

Charles Darwin's On the Origin of Species is at once familiar and unfamiliar. Everyone knows that the Origin introduced the world to the idea of evolution by natural selection, but few of us have actually read it. We suggest that it is worth taking the time not only to read what Darwin had to say, but also to use the Origin to teach both biology and writing. It provides scientific lessons in areas beyond evolutionary biology, such as ecology and biogeography. In addition, it provides valuable rhetorical lessons-how to construct an argument, write persuasively, make use of evidence, know your audience, and anticipate counterarguments. We have been using the Origin in various classes for several years, introducing new generations to Darwin, in his own words.

A comparative first-principles study of martensitic phase transformations in TiPd2 and TiPd intermetallics.

Martensitic phase transformations in TiPd2 and TiPd alloys are studied employing density-functional, first-principles calculations. We examine the transformation of tetragonal C11b TiPd2 to the low-temperature orthorhombic phase (C11b â†’ oI6), and the transformation of cubic B2 TiPd under orthorhombic (B2→B19) and subsequent monoclinic transformations (B19→B19') as the system is cooled. We employ a theoretical approach based on a phenomenological Landau theory of the structural phase transitions and a mean-field approximation for the free energy, utilizing first-principles calculations to obtain the deformation energy as a function of strains and to deduce parameters for constructing the free energy. The predicted transition temperature for the TiPd2 C11b â†’ oI6 transition is in good agreement with reported experimental results. To investigate the TiPd B2→B19 transformation, we employ both the Cauchy-Born rule and a soft-mode-based approach, and elucidate the importance of the coupling between lattice distortion and atomic displacements (i.e. shuffling) in the formation of the final structure. The calculated B2→B19 transition temperature for TiPd alloy agrees well with the experimental results. We also find that there exists a very small but finite (0.0005 eV/atom) energy barrier of B19 TiPd under monoclinic deformation for B19→B19' structural phase transformation.

Microstructure-Dependent Gas Adsorption: Accurate Predictions of Methane Uptake in Nanoporous Carbons.

We present a framework for rapidly predicting gas adsorption properties based on van der Waals density functional calculations and thermodynamic modeling. Utilizing this model and experimentally determined pore size distributions, we are able to accurately predict uptakes in five activated carbon materials without empirical potentials or lengthy simulations. Our results demonstrate that materials with smaller pores and higher heats of adsorption can still have poor adsorption characteristics due to relatively low densities of highly adsorbent pores.

Freezing point and solid-liquid interfacial free energy of Stockmayer dipolar fluids: a molecular dynamics simulation study.

Stockmayer fluids are a prototype model system for dipolar fluids. We have computed the freezing temperatures of Stockmayer fluids at zero pressure using three different molecular-dynamics simulation methods, namely, the superheating-undercooling method, the constant-pressure and constant-temperature two-phase coexistence method, and the constant-pressure and constant-enthalpy two-phase coexistence method. The best estimate of the freezing temperature (in reduced unit) for the Stockmayer (SM) fluid with the dimensionless dipole moment µ*=1, √2, √3 is 0.656 ± 0.001, 0.726 ± 0.002, and 0.835 ± 0.005, respectively. The freezing temperature increases with the dipolar strength. Moreover, for the first time, the solid-liquid interfacial free energies γ of the fcc (111), (110), and (100) interfaces are computed using two independent methods, namely, the cleaving-wall method and the interfacial fluctuation method. Both methods predict that the interfacial free energy increases with the dipole moment. Although the interfacial fluctuation method suggests a weaker interfacial anisotropy, particularly for strongly dipolar SM fluids, both methods predicted the same trend of interfacial anisotropy, i.e., γ100 > γ110 > γ111.

Isotope effect on adsorbed quantum phases: diffusion of H2 and D2 in nanoporous carbon.

Quasielastic neutron scattering of H(2) and D(2) in the same nanoporous carbon at 10-40 K demonstrates extreme quantum sieving, with D(2) diffusing up to 76 times faster. D(2) also shows liquidlike diffusion while H(2) exhibits Chudley-Elliott jump diffusion, evidence of their different relationships with the local lattice of adsorption sites due to quantum effects on intermolecular interactions. The onset of diffusion occurs at 22-25 K for H(2) and 10-13 K for D(2). At these temperatures, H(2) and D(2) have identical thermal de Broglie wavelengths that correlate with the dominant pore size.

The influence of dispersion interactions on the hydrogen adsorption properties of expanded graphite.

We demonstrate the importance of London dispersion forces in defining the adsorption capacity within expanded graphite, a simple model of the more complex experimental geometries of activated carbon, using a combination of the non-local correlation functional of Dion et al paired with a recent exchange functional of Cooper (vdW-DF(C09x)) and a classical continuum model. Our results indicate that longer ranged interactions due to dispersion forces increase the volume over which molecules interact with a porous medium. This significantly enhances the adsorption density within a material, and explains recent experimental work showing that the densification of H(2) in carbon nanopores is sensitive to the pore size. Remarkably, our slit pore geometries give adsorption densities of up to 3 wt% at 298 K and 20 MPa which correlates well with experimental values for 9 Å pores-a value that could not be predicted using local density approximation (LDA) calculations. In its entirety, this work presents a powerful approach for assessing molecular uptake in porous media and may have serious impacts on efforts to optimize the properties of these materials.

Topological defects: origin of nanopores and enhanced adsorption performance in nanoporous carbon.

A scanning transmission electron microscopy investigation of two nanoporous carbon materials, wood-based ultramicroporous carbon and poly(furfuryl alcohol)-derived carbon, is reported. Atomic-resolution images demonstrate they comprise isotropic, three-dimensional networks of wrinkled one-atom-thick graphene sheets. In each graphene plane, nonhexagonal defects are frequently observed as connected five- and seven-atom rings. Atomic-level modeling shows that these topological defects induce localized rippling of graphene sheets, which interferes with their graphitic stacking and induces nanopores that lead to enhanced adsorption of H(2) molecules. The poly(furfuryl alcohol)-derived carbon contains larger regions of stacked layers, and shows significantly smaller surface area and pore volume than the ultramicroporous carbon.

A parameter-free prediction of simulated crystal nucleation times in the Lennard-Jones system: from the steady-state nucleation to the transient time regime.

Large scale simulations of crystal nucleation from the liquid are performed using the Lennard-Jones potential to determine the time required for nucleation. By considering both transient and finite-size effects, we for the first time successfully predict the nucleation time without any parameter fitting in the Lennard-Jones system. All necessary parameters are derived from separate, equilibrium simulations. At small undercoolings, large system size effects are observed. The required system size is not only determined by the size of the critical nuclei, but also the characteristic spacing between them. Two distinct nucleation regions are predicted by theory and observed by the simulations, which are dominated by the transient time and the steady-state nucleation time, respectively.

Rat model for studying tissue changes induced by the mechanical environment surrounding loaded titanium implants.

PURPOSE: This study presents a new rat oral implant model for assessing histologic changes in the mechanical environment surrounding loaded and unloaded dental implants. MATERIALS AND METHODS: The maxillary left first molar from retired breeder rats was extracted, and the site was allowed to heal for 1 month. A titanium miniscrew implant was then placed into the site and allowed to heal for 21 days. The mandibular left first molars in one group of rats were extracted to create an unloaded condition; in a second group of rats the mandibular left first molars were left in occlusion with the opposing screw head to simulate loading. Radiographs were taken on the day of placement and again at 7 days, 14 days, and 21 days after placement and were used to estimate the bone-implant contact ratio. The rats were sacrificed after 21 days. Peri-implant tissue samples from day 21 were processed for histology and immunohistochemistry with antibodies to osteocalcin and matrix metalloproteinase 13 (MMP-13). Two-dimensional finite element models were created from images of the histologic sections and immunohistochemical samples to observe tissue changes. RESULTS: Areas of high shear stress adjacent to the helical threads of loaded implants were associated with osteocalcin localization and bone formation but only minimal localization of MMP-13. Bone adjacent to unloaded implants showed fibrous tissue and extensive MMP-13 localization surrounding the apical two-thirds of each implant. These results agree with estimated bone-implant contact ratios, which showed a steady decrease in contact ratio for the unloaded implant group but a significantly higher contact ratio in the loaded group between 14 and 21 days. CONCLUSION: The rat oral implant model is useful for studies of the mechanical and physiologic environment affecting osseointegration in loaded and unloaded implants.

Effects of chromosomal rearrangements on transvection at the yellow gene of Drosophila melanogaster.

Homologous chromosomes are paired in somatic cells of Drosophila melanogaster. This pairing can lead to transvection, which is a process by which the proximity of homologous genes can lead to a change in gene expression. At the yellow gene, transvection is the basis for several examples of intragenic complementation involving the enhancers of one allele acting in trans on the promoter of a paired second allele. Using complementation as our assay, we explored the chromosomal requirements for pairing and transvection at yellow. Following a protocol established by Ed Lewis, we generated and characterized chromosomal rearrangements to define a region in cis to yellow that must remain intact for complementation to occur. Our data indicate that homolog pairing at yellow is efficient, as complementation was disrupted only in the presence of chromosomal rearrangements that break

The anisotropic hard-sphere crystal-melt interfacial free energy from fluctuations.

We have calculated the interfacial free energy for the hard-sphere system, as a function of crystal interface orientation, using a method that examines the fluctuations in the height of the interface during molecular dynamics simulations. The approach is particularly sensitive for the anisotropy of the interfacial free energy. We find an average interfacial free energy of gamma=0.56+/-0.02k(B)Tsigma(-2). This value is lower than earlier results based upon direct calculations of the free energy [R. L. Davidchack and B. B. Laird, Phys. Rev. Lett. 85, 4751 (2000)]. However, both the average value and the anisotropy agree with the recent values obtained by extrapolation from direct calculations for a series of the inverse-power potentials [R. L. Davidchack and B. B. Laird, Phys. Rev. Lett. 94, 086102 (2005)].

Quantitative parameter-free prediction of simulated crystal-nucleation times.

We present direct comparisons between simulated crystal-nucleation times and theoretical predictions using a model of aluminum, and demonstrate that a quantitative prediction can be made. All relevant thermodynamic properties of the system are known, making the agreement of our simulation data with nucleation theories free of any adjustable parameters. The role of transient nucleation is included in the classical nucleation theory approach, and shown to be necessary to understand the observed nucleation times. The calculations provide an explanation on why nucleation is difficult to observe in simulations at moderate undercoolings. Even when the simulations are significantly larger than the critical nucleus, and when simulation times are sufficiently long, at moderate undercoolings the small concentration of critical nuclei makes the probability of the nucleation low in molecular dynamics simulations.

Melting temperature of ice Ih calculated from coexisting solid-liquid phases.

We carried out molecular-dynamics simulations by using the two-phase coexistence method with the constant pressure, particle number, and enthalpy ensemble to compute the melting temperature of proton-disordered hexagonal ice I(h) at 1-bar pressure. Four models of water were considered, including the widely used TIP4P [W. L. Jorgensen, J. Chandrasekha, J. D. Madura, R. W. Impey, and M. L. Klein, J. Chem. Phys.79, 926 (1983)] and TIP5P [M. W. Mahoney and W. L. Jorgensen J. Chem. Phys.112, 8910 (2000)] models, as well as recently improved TIP4P and TIP5P models for use with Ewald techniques-the TIP4P-Ew [W. Horn, W. C. Swope, J. W. Pitera, J. C. Madura, T. J. Dick, G. L. Hura, and T. Head-Gordon, J. Chem. Phys.120, 9665 (2004)] and TIP5P-Ew [S. W. Rick, J. Chem. Phys.120, 6085 (2004)] models. The calculated melting temperature at 1 bar is T(m) = 229 +/- 1 K for the TIP4P and T(m) = 272.0 +/- 0.6 K for the TIP5P ice I(h), both are consistent with previous simulations based on free-energy methods. For the TIP4P-Ew and TIP5P-Ew models, the calculated melting temperature is T(m) = 257.0 +/- 1.1 K and T(m) = 253.9 +/- 1.1 K, respectively.

Ab initio molecular dynamics simulation of liquid Al88Si12 alloys.

First-principles molecular dynamics simulations are carried out to study the structures, dynamics, and electronic properties of liquid Al88Si12 in the temperature ranging from 898 to 1298 K. The temperature dependence of static structure factors, pair correlation functions, and electronic density-of-states are investigated. The structural properties obtained from the simulations are in good agreement with the x-ray diffraction experimental results.

Enhancer choice in cis and in trans in Drosophila melanogaster: role of the promoter.

Eukaryotic enhancers act over very long distances, yet still show remarkable specificity for their own promoter. To better understand mechanisms underlying this enhancer-promoter specificity, we used transvection to analyze enhancer choice between two promoters, one located in cis to the enhancer and the other in trans to the enhancer, at the yellow gene of Drosophila melanogaster. Previously, we demonstrated that enhancers at yellow prefer to act on the cis-linked promoter, but that mutation of core promoter elements in the cis-linked promoter releases enhancers to act in trans. Here, we address the mechanism by which these elements affect enhancer choice. We consider and explicitly test three models that are based on promoter competency, promoter pairing, and promoter identity. Through targeted gene replacement of the endogenous yellow gene, we show that competency of the cis-linked promoter is a key parameter in the cis-trans choice of an enhancer. In fact, complete replacement of the yellow promoter with both TATA-containing and TATA-less heterologous promoters maintains enhancer action in cis.