There has been considerable recent interest in the high-pressure behavior of silicon carbide, a potential major constituent of carbon-rich exoplanets. In this work, the atomic-level structure of SiC was determined through in situ X-ray diffraction under laser-driven ramp compression up to 1.5 TPa; stresses more than seven times greater than previous static and shock data. Here we show that the B1-type structure persists over this stress range and we have constrained its equation of state (EOS). Using this data we have determined the first experimentally based mass-radius curves for a hypothetical pure SiC planet. Interior structure models are constructed for planets consisting of a SiC-rich mantle and iron-rich core. Carbide planets are found to be ~10% less dense than corresponding terrestrial planets.
Geyman, Emily C.; Wu, Ziman; Nadeau, Matthew D.; Edmonsond, Stacey; Turner, Andrew; Purkis, Sam J.; Howes, Bolton; Dyer, Blake; Ahm, Anne-Sofie C.; Yao, Nan; Deutsch, Curtis A.; Higgins, John A.; Stolper, Daniel A.; Maloof, Adam C.
Carbonate mud represents one of the most important geochemical archives for reconstructing ancient climatic, environmental, and evolutionary change from the rock record. Mud also represents a major sink in the global carbon cycle. Yet, there remains no consensus about how and where carbonate mud is formed. In this contribution, we present new geochemical data that bear on this problem, including stable isotope and minor and trace element data from carbonate sources in the modern Bahamas such as ooids, corals, foraminifera, and green algae.
The carbon isotopic (δ13C) composition of shallow-water carbonates often is interpreted to reflect the δ13C of the global ocean and is used as a proxy for changes in the global carbon cycle. However, local platform processes, in addition to meteoric and marine diagenesis, may decouple carbonate δ13C from that of the global ocean. To shed light on the extent to which changing sediment grain composition may produce δ13C shifts in the stratigraphic record, we present new δ13C measurements of benthic foraminifera, solitary corals, calcifying green algae, ooids, coated grains, and lime mud from the modern Great Bahama Bank (GBB). This survey of a modern carbonate environment reveals δ13C variability comparable to the largest δ13C excursions in the last two billion years of Earth history.
The history of organismal evolution, seawater chemistry, and paleoclimate is recorded in layers of carbonate sedimentary rock. Meter-scale cyclic stacking patterns in these carbonates often are interpreted as representing sea level change. A reliable sedimentary proxy for eustasy would be profoundly useful for reconstructing paleoclimate, since sea level responds to changes in temperature and ice volume. However, the translation from water depth to carbonate layering has proven difficult, with recent surveys of modern shallow water platforms revealing little correlation between carbonate facies (i.e., grain size, sedimentary bed forms, ecology) and water depth. We train a convolutional neural network with satellite imagery and new field observations from a 3,000 km2 region northwest of Andros Island (Bahamas) to generate a facies map with 5 m resolution. Leveraging a newly-published bathymetry for the same region, we test the hypothesis that one can extract a signal of water depth change, not simply from individual facies, but from sequences of facies transitions analogous to vertically stacked carbonate strata. Our Hidden Markov Model (HMM) can distinguish relative sea level fall from random variability with ∼90% accuracy. Finally, since shallowing-upward patterns can result from local (autogenic) processes in addition to forced mechanisms such as eustasy, we search for statistical tools to diagnose the presence or absence of external forcings on relative sea level. With a new data-driven forward model that simulates how modern facies mosaics evolve to stack strata, we show how different sea level forcings generate characteristic patterns of cycle thicknesses in shallow carbonates, providing a new tool for quantitative reconstruction of ancient sea level conditions from the geologic record.
Chronic hepatitis B (CHB), caused by hepatitis B virus (HBV), remains a major medical problem. HBV has a high propensity for progressing to chronicity and can result in severe liver disease, including fibrosis, cirrhosis and hepatocellular carcinoma. CHB patients frequently present with viral coinfection, including HIV and hepatitis delta virus. About 10% of chronic HIV carriers are also persistently infected with HBV which can result in more exacerbated liver disease. Mechanistic studies of HBV-induced immune responses and pathogenesis, which could be significantly influenced by HIV infection, have been hampered by the scarcity of immunocompetent animal models. Here, we demonstrate that humanized mice dually engrafted with components of a human immune system and a human liver supported HBV infection, which was partially controlled by human immune cells, as evidenced by lower levels of serum viremia and HBV replication intermediates in the liver. HBV infection resulted in priming and expansion of human HLA-restricted CD8+ T cells, which acquired an activated phenotype. Notably, our dually humanized mice support persistent coinfections with HBV and HIV which opens opportunities for analyzing immune dysregulation during HBV and HIV coinfection and preclinical testing of novel immunotherapeutics.
The prevalence of ooids in the stratigraphic record, and their association with shallow-water carbonate environments, make ooids an important paleoenvironmental indicator. Recent advances in the theoretical understanding of ooid morphology, along with empirical studies from Turks and Caicos, Great Salt Lake, and The Bahamas, have demonstrated that the morphology of ooids is indicative of depositional environment and hydraulic conditions. To apply this knowledge from modern environments to the stratigraphic record of Earth history, researchers measure the size and shape of lithified ooids on two-dimensional surfaces (i.e., thin sections or polished slabs), often assuming that random 2D slices intersect the nuclei and that the orientation of the ooids is known. Here we demonstrate that these assumptions rarely are true, resulting in errors of up to 35% on metrics like major axis length. We present a method for making 3D reconstructions by serial grinding and imaging, which enables accurate measurement of the morphology of individual ooids within an oolite, as well as the sorting and porosity of a sample. We also provide three case studies that use the morphology of ooids in oolites to extract environmental information. Each case study demonstrates that 2D measurements can be useful if the environmental signal is large relative to the error from 2D measurements. However, 3D measurements substantially improve the accuracy and precision of environmental interpretations. This study focuses on oolites, but errors from 2D measurements are not unique to oolites; this method can be used to extract accurate grain and porosity measurements from any lithified granular sample.
This dataset contains input files, training data and other files related to the machine learning models developed during the work by Muniz et al. In this work, we construct machine learning models based on the MB-pol many-body model. We find that the training set should include cluster configurations as well as liquid phase configurations in order to accurately represent both liquid and VLE properties. The results attest for the ability of machine learning models to accurately represent many-body potentials and provide an efficient avenue for water simulations.
Understanding the condensed-phase behavior of chiral molecules is important in biology, as well as in a range of technological applications, such as the manufacture of pharmaceuticals. Here, we use molecular dynamics simulations to study a chiral four-site molecular model that exhibits a second-order symmetry-breaking phase transition from a supercritical racemic liquid, into subcritical D-rich and L-rich liquids. We determine the infinite-size critical temperature using the fourth-order Binder cumulant, and we show that the finite-size scaling behavior of the order parameter is compatible with the 3D Ising universality class. We also study the spontaneous D-rich to L-rich transition at a slightly subcritical temperature T ~ 0.985 Tc and our findings indicate that the free energy barrier for this transformation increases with system size as N^2/3 where N is the number of molecules, consistent with a surface-dominated phenomenon. The critical behavior observed herein suggests a mechanism for chirality selection in which a liquid of chiral molecules spontaneously forms a phase enriched in one of the two enantiomers as the temperature is lowered below the critical point. Furthermore, the increasing free energy barrier with system size indicates that fluctuations between the L-rich and D-rich phases are suppressed as the size of the system increases, trapping it in one of the two enantiomerically-enriched phases. Such a process could provide the basis for an alternative explanation for the origin of biological homochirality. We also conjecture the possibility of observing nucleation at subcritical temperatures under the action of a suitable chiral external field.
This dataset contains input and output files to reproduce the results of the manuscript "Homogeneous ice nucleation in an ab initio machine learning model" by Pablo M. Piaggi, Jack Weis, Athanassios Z. Panagiotopoulos, Pablo G. Debenedetti, and Roberto Car (arXiv preprint https://arxiv.org/abs/2203.01376). In this work, we studied the homogeneous nucleation of ice from supercooled liquid water using a machine learning model trained on ab initio energies and forces. Since nucleation takes place over times much longer than the simulation times that can be afforded using molecular dynamics simulations, we make use of the seeding technique that is based on simulating an ice cluster embedded in liquid water. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or shrinking at the given supersaturation). Using data from the seeding simulations and the equations of classical nucleation theory we compute nucleation rates that can be compared with experiments.
This dataset contains all data relevant to a forthcoming publication in which we used molecular simulation methods to study the phase behavior of supercooled water. The dataset contains simulation input and output files, processed data files, and image files used to create all plots in the manuscript. Python analysis scripts are also included, including instructions for how to re-generate all plots in the manuscript.
Gartner III, Thomas E.; Torquato, Salvatore; Car, Roberto; Debenedetti, Pablo G.
This dataset contains all data related to the publication "Manifestations of metastable criticality in glassy water-like models detected by large-scale structural properties" by Gartner et al., in preparation 2020. In this work, we used molecular dynamics simulations to explore the relationship between water's polyamorphism (multiple amorphous solid states) and its hypothesized liquid-liquid transition. Using the TIP4P/2005 molecular model of water, we found a surprising signature of water's liquid-liquid critical point in the long-range structure of water's amorphous solid states formed by isobaric cooling at different pressures. This structural signature was absent in two other systems that lack a critical point. This dataset contains molecular dynamics simulation trajectories, as well as processed data, analysis codes, and image files used in the publication.
Piaggi, Pablo M; Gartner, Thomas E; Car, Roberto; Debenedetti, Pablo G
The possible existence of a liquid-liquid critical point in deeply supercooled water has been a subject of debate in part due to the challenges associated with providing definitive experimental evidence. Pioneering work by Mishima and Stanley [Nature 392, 164 (1998) and Phys.~Rev.~Lett. 85, 334 (2000)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected location of the liquid-liquid transition and its associated critical point. Based on the continuous or discontinuous changes in slope of the melting curves, Mishima suggested that the liquid-liquid critical point lies between the melting curves of ice III and ice V. Here, we explore this conjecture using molecular dynamics simulations with a purely-predictive machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII using this model and find that the melting lines of all the studied ice polymorphs are supercritical and do not intersect the liquid-liquid transition locus. We also find a pronounced, yet continuous, change in slope of the melting lines upon crossing of the locus of maximum compressibility of the liquid. Finally, we analyze critically the literature in light of our findings, and conclude that the scenario in which melting curves are supercritical is favored by the most recent computational and experimental evidence. Thus, although the preponderance of experimental and computational evidence is consistent with the existence of a second critical point in water, the behavior of the melting lines of ice polymorphs does not provide strong evidence in support of this viewpoint, according to our calculations.
Gartner, Thomas III; Zhang, Linfeng; Piaggi, Pablo; Car, Roberto; Panagiotopoulos, Athanassios; Debenedetti, Pablo
This dataset contains all data related to the publication "Signatures of a liquid-liquid transition in an ab initio deep neural network model for water", by Gartner et al., 2020. In this work, we used neural networks to generate a computational model for water using high-accuracy quantum chemistry calculations. Then, we used advanced molecular simulations to demonstrate evidence that suggests this model exhibits a liquid-liquid transition, a phenomenon that can explain many of water's anomalous properties. This dataset contains links to all software used, all data generated as part of this work, as well as scripts to generate and analyze all data and generate the plots reported in the publication.
This dataset comprises of data associated with the publication "Transferability of data-driven, many-body models for CO2 simulations in the vapor and liquid phases", which can be found at https://doi.org/10.1063/5.0080061. The data includes calculations for a Many-Body decomposition, virial coefficient calculations, orientational molecular scan energies, potential energy fields, correlation plots of training and testing data, vapor-liquid equilibrium simulations, liquid density simulations, and solid cell simulations.
Muniz, Maria Carolina; Gartner III, Thomas E.; Riera, Marc; Knight, Christopher; Yue, Shuwen; Paesani, Francesco; Panagiotopoulos, Athanassios Z.
This dataset contains all data (including input files, simulation trajectories as well as other data files and analysis scripts) related to the publication "Vapor-liquid equilibrium of water with the MB-pol many-body potential" by Muniz et al. in preparation (2021). In this work, we assessed the performance of the MB-pol many-body potential with respect to water's vapor-liquid equilibrium properties. Through the use of direct coexistence molecular dynamics, we calculated properties such as coexistence densities, surface tension, vapor pressures and enthalpy of vaporization. We found that MB-pol is able to predict these properties in good agreement with experimental data. The results attest to the chemical accuracy of MB-pol and its large range of application across water's phase diagram.
This dataset is a sequence of laser-induced fluorescence images of a dye injected in a channel flow with canopy-like stainless steel rods simulating a vegetation canopy stand. The data is acquired close to the channel bottom at z/h=0.2, where z is the height referenced to the channel bed and h is the canopy height. The dataset provides spatial distribution of scalar concentration in a plane parallel to the channel bed. The data has been used (but the data itself has not been published or available to the public) in previous work. The references are: Ghannam, K., Poggi, D., Porporato, A., & Katul, G. (2015). The spatio-temporal statistical structure and ergodic behaviour of scalar turbulence within a rod canopy. Boundary-Layer Meteorology,157(3), 447–460. Ghannam, K, Poggi, D., Bou-Zeid, E., Katul, G. (2020). Inverse cascade evidenced by information entropy of passive scalars in submerged canopy flows. Geophysical Research Letters (accepted).
This distribution contains experimentally measured data for the extent of retained enzyme activity post thermal stressing for three distinct enzymes: glucose oxidase, lipase, and horseradish peroxidase. The data is used to form conclusions and develop machine learning models as reported in the publication "Machine Learning on a Robotic Platform for the Design of Polymer-Protein Hybrids" by Matthew Tamasi, Roshan Patel, Carlos Borca, Shashank Kosuri, Heloise Mugnier, Rahul Upadhya, N. Sanjeeva Murthy, Michael Webb*, and Adam Gormley. Details regarding the experimental protocols are reported in the aforementioned paper but are briefly discussed in the README.
Data set corresponding to "NAPS: Integrating pose estimation and tag-based tracking." This dataset contains the corresponding videos, tracking scripts, and SLEAP models along with SLEAP, NAPS, and ArUco tracking results.
Data set for "Ocean emission of microplastic by bursting bubble jet drops." Two .csv files are provided: one for the size of a jet drop carrying microplastic, and another for the amount of microplastic captured by a jet drop.