Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid phase processing”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Preparation of novel layer-stack hexagonal CdO micro-rods by a pre-oxidation and subsequent evaporation process

Graphical abstract: Layer-stack hexagonal cadmium oxide (CdO) micro-rods were prepared. - Highlights: • Novel hexagonal layer-stack structure CdO micro-rods were synthesized by a thermal evaporation method. • The pre-oxidation, vapor pressure and substrate nature play a key role on the formation of CdO rods. • The formation mechanism of CdO micro-rods was explained. - Abstract: Novel layer-stack hexagonal cadmium oxide (CdO) micro-rods were prepared by pre-oxidizing Cd granules and subsequent thermal oxidation under normal atmospheric pressure. X-ray diffraction (XRD) and scanning electron microscopy (SEM) were performed to characterize the phase structure and microstructure. The pre-oxidation process, vapor pressure and substrate nature were the key factors for the formation of CdO micro-rods. The diameter of micro-rod and surface rough increased with increasing of thermal evaporation temperature, the length of micro-rod increased with the increasing of evaporation time. The formation of hexagonal layer-stack structure was explained by a vapor–solid mechanism.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Mineral dissolution and wormholing from a pore-scale perspective

A micro-continuum approach is proposed to simulate the dissolution of solid minerals at the pore scale under single-phase flow conditions. The approach employs a Darcy–Brinkman–Stokes formulation and locally averaged conservation laws combined with immersed boundary conditions for the chemical reaction at the solid surface. The methodology compares well with the arbitrary-Lagrangian–Eulerian technique. The simulation framework is validated using an experimental microfluidic device to image the dissolution of a single calcite crystal. The evolution of the calcite crystal during the acidizing process is analysed and related to the flow conditions. Macroscopic laws for the dissolution rate are proposed by upscaling the pore-scale simulations. Finally, the emergence of wormholes during the injection of acid in a two-dimensional domain of calcite grains is discussed based on pore-scale simulations.

bio-inspired, mechanical behavior, carbon sequestr↗
DOE OSTI2015

Deformation and Failure Mechanisms of Shape Memory Alloys

The goal of this research was to understand the fundamental mechanics that drive the deformation and failure of shape memory alloys (SMAs). SMAs are difficult materials to characterize because of the complex phase transformations that give rise to their unique properties, including shape memory and superelasticity. These phase transformations occur across multiple length scales (one example being the martensite-austenite twinning that underlies macroscopic strain localization) and result in a large hysteresis. In order to optimize the use of this hysteretic behavior in energy storage and damping applications, we must first have a quantitative understanding of this transformation behavior. Prior results on shape memory alloys have been largely qualitative (i.e., mapping phase transformations through cracked oxide coatings or surface morphology). The PI developed and utilized new approaches to provide a quantitative, full-field characterization of phase transformation, conducting a comprehensive suite of experiments across multiple length scales and tying these results to theoretical and computational analysis. The research funded by this award utilized new combinations of scanning electron microscopy, diffraction, digital image correlation, and custom testing equipment and procedures to study phase transformation processes at a wide range of length scales, with a focus at small length scales with spatial resolution on the order of 1 nanometer. These experiments probe the basic connections between length scales during phase transformation. In addition to the insights gained on the fundamental mechanisms driving transformations in shape memory alloys, the unique experimental methodologies developed under this award are applicable to a wide range of solid-to-solid phase transformations and other strain localization mechanisms.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Observations and Lessons Learned from Testing of the Hanford Waste Treatment Plant Standard High Solids Vessel Design - 18635

The Hanford Waste Treatment and Immobilization Plant (WTP) recently completed qualification testing of the mixing system for the standard high solids vessel design the U.S. Department of Energy (DOE) proposed to replace several vessels in the design of the plant Pretreatment Facility. Testing was split into two phases. The first phase exercised the pulse jet mixer control system to ensure that the control software could accommodate fluids across the complete range of expected densities and viscosities within the vessel for all operational modes. The second phase of testing was focused on the mixing performance of the vessel, specifically, whether mixing was adequate to support the transfer, de-inventory, throughput, and sampling functions of the vessel. Key observations and lessons learned from testing include: - Successful transfer of waste within and beyond the design basis slurry solids concentration; - Successful resuspension and transfer of large, heavy particles; - Effective performance of multi-level bubblers for process control and monitoring; - Effective controls performance; - Reduction of sparger flowrate from nominal design flow; - Liquid blending within about 15 minutes; - Insights into sample representativeness and uniformity of transfers to inform process control strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Wet chemical synthesis and properties of argyrodite sulfide solid electrolytes for solid state lithium batteries

The commercialization of the lithium-ion battery (LIB) in 1991 was responsible for the explosion in portable electronic technologies that has been seen over the past 30 years. With the advent of electric vehicles and other high-powered technologies, there is tremendous demand for LIBs with higher energy density and high safety. To achieve this, new electrode materials must be explored. The obvious choice of anode material would be pure metal lithium, which has a theoretical specific capacity of 3860 mAh g-1 . Unfortunately, metal lithium anodes have not been widely commercialized due to their tendency to react violently with the flammable liquid electrolytes used in today’s batteries. Battery safety can best be achieved by adopting solid electrolytes in place of liquid electrolytes. Solid electrolytes are nonvolatile and nonflammable, safely allowing for the combination of high-capacity cathode materials with a Li metal anode. Argyrodite sulfide solid electrolytes such as halogen-doped Li6PS5X (X = Cl, Br, I) are noted for their high ionic conductivity. But before sulfides can be commercially adopted, they possess several disadvantages which must be addressed, including time- and energy-consuming synthesis processes, poor electrochemical stability, and intrinsically poor air stability. This dissertation seeks to address each of these challenges through materials design an synthesis strategies. In this work, we pioneer a solvent-based approach for the synthesis of argyrodite solid electrolytes Li7PS6 and Li6PS5Xinstead of a stringent solid-state synthesis. Nontoxic ethanol is employed as the solvent, enabling a rapid synthetic approach to produce argyrodite solid electrolytes with high phase purity and compositional flexibility. Compared with Li7PS6, halogen doping (i.e. X = F, Cl, Br, I) not only increases the ionic conductivity, but also enhances the electrochemical stability at the interface towards Li metal. Specifically, F-doped argyrodites produce a robust SEI layer containing LiF, contributing to enhanced interfacial stability. Finally, to address the air instability challenge, argyrodite-incorporated composite solid electrolytes (CSEs) are designed and prepared to produce stable and flexible membranes that are demonstrated in solid-state Li metal batteries. These advances push argyrodite sulfide solid electrolyte research further and pave the way for the proliferation of next generation lithium metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2013

Environmental impacts of residual Municipal Solid Waste incineration: A comparison of 110 French incinerators using a life cycle approach

Highlights: • 110 French incinerators are compared with LCA based on plant-specific data. • Environmental impacts vary as a function of plants energy recovery and NO{sub x} emissions. • E.g. climate change impact ranges from −58 to 408 kg CO{sub 2}-eq/tonne of residual MSW. • Implications for LCA of waste management in a decision-making process are detailed. - Abstract: Incineration is the main option for residual Municipal Solid Waste treatment in France. This study compares the environmental performances of 110 French incinerators (i.e. 85% of the total number of plants currently in activity in France) in a Life Cycle Assessment perspective, considering 5 non-toxic impact categories: climate change, photochemical oxidant formation, particulate matter formation, terrestrial acidification and marine eutrophication. Mean, median and lower/upper impact potentials are determined considering the incineration of 1 tonne of French residual Municipal Solid Waste. The results highlight the relatively large variability of the impact potentials as a function of the plant technical performances. In particular, the climate change impact potential of the incineration of 1 tonne of waste ranges from a benefit of −58 kg CO{sub 2}-eq to a relatively large burden of 408 kg CO{sub 2}-eq, with 294 kg CO{sub 2}-eq as the average impact. Two main plant-specific parameters drive the impact potentials regarding the 5 non-toxic impact categories under study: the energy recovery and delivery rate and the NO{sub x} process-specific emissions. The variability of the impact potentials as a function of incinerator characteristics therefore calls for the use of site-specific data when required by the LCA goal and scope definition phase, in particular when the study focuses on a specific incinerator or on a local waste management plan, and when these data are available.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2018

Testing Results from the Low Activity Waste Pretreatment System (LAWPS) Full-Scale Ion Exchange Column Test and Engineering Scale Integrated System Test - 18634

The U.S. Department of Energy (DOE), Office of River Protection's (ORP) primary mission is to retrieve and treat Hanford's tank waste and close the Tank Farms to protect the Columbia River. Mixed radioactive waste is stored in 177 underground tanks at the Hanford Site as reported in DOE/ORP 2003- 02, Environmental Impact Statement for Retrieval, Treatment, and Disposal of Tank Waste and Closure of the Single Shell Tanks at the Hanford Site, Richland WA, Inventory and Source Term Data Package. The 177 underground tanks are estimated to contain about 56 million gallons of waste. The Low Activity Waste Pretreatment System (LAWPS) Project provides for the early production of immobilized low activity waste (ILAW) by feeding LAW directly from Tank Farms to the Hanford Waste Treatment Plant's (WTP's) LAW Facility, bypassing the Pretreatment Facility. Prior to the transfer of feed to the WTP LAW Vitrification Facility, tank supernatant waste will be pretreated in the LAWPS to meet the WTP LAW waste acceptance criteria. The LAWPS will also facilitate the return of secondary liquid wastes from the WTP LAW Vitrification Facility to the Tank Farms. In November of 2015 DOE's Tank Operations Contractor, Washington River Protection Solutions (WRPS), awarded a contract to the AECOM team to design, fabricate and operate two test platforms to generate performance data at two different scales and demonstrate technology maturation. The first platform models the ion exchange (IX) columns at full scale, and the second platform models the integrated flowsheet (cross flow filter (CFF) and IX) at 1/9. scale. Overarching technical objectives of performing the full-scale ion exchange test and the engineering-scale integrated test are to: 1. Test interactions between unit process operations; 2. Test removal of non-radioactive cesium from process stream; 3. Test removal of undissolved solids from process stream; 4. Confirm operational flow rate range and stability; 5. Confirm volumetric throughput of waste stream over flow operating range; 6. Confirm range of process temperatures over planned operating range; 7. Establish system hydraulics and process criteria; 8. Coordinate key control and monitoring components; 9. Establish treated product chemistry meets requirements (cesium, solids, and pH); and 10. Test the hazard control strategy for a no flow condition through the IX column. Data collected during testing will be used by WRPS to validate the preliminary design, inform the final design of the permanent plant equipment, and demonstrate that the critical technology elements have advanced to a technology readiness level of 6 (prototypical system validation in a relevant environment). Testing on both platforms started in fall of 2016 generating NQA-1 data. This paper covers test results from the first five phases of testing on the 1/9. scale integrated test platform. A further three phases of testing are currently being executed, these including a 5.6 M sodium simulant containing foulant, a high viscosity simulant and finally a run performed with a lower Cs concentration (all earlier tests were spiked to 4 x Cs to reduce volume and cost). Results from these tests were not available at the time of writing available. Results are also presented from the full scale IX testing. Eight test were performed covering column loading, elution and regeneration; column loading, elution and regeneration with upflow displacement; IX column loading and gas migration; resin removal; solids accumulation; and IX column emergency recirculation and vent. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Degradation processes of reinforced concretes by combined sulfate–phosphate attack

A novel form of alteration due to the interaction between hydrated cement phases and sulfate and phosphate-based pollutants is described, through the characterization of concrete samples from an industrial reinforced concrete building. Decalcification of the cement matrices was observed, with secondary sulfate and phosphate-based mineral formation, according to a marked mineralogical and textural zoning. Five alteration layers may be detected: the two outermost layers are characterized by the presence of gypsum–brushite solid solution phases associated with anhydrous calcium sulfates and phosphates, respectively, while a progressive increase in apatite and ammonium magnesium phosphates is observable in the three innermost layers, associated with specific apatite precursors (brushite, octacalcium phosphate and amorphous calcium phosphate, respectively). The heterogeneous microstructural development of secondary phases is related to the chemical, pH and thermal gradients in the attacked cementitious systems, caused by different sources of pollutants and the exposure to the sun's radiation.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Materials Engineering and Scale Up of Fluid Phase Chemical Hydrogen Storage for Automotive Applications

Among candidates for chemical hydrogen storage in PEM fuel cell automotive applications, ammonia borane (AB, NH3BH3) is considered to be one of the most promising materials due to its high hydrogen content of 14-16 wt% below 200°C and high volumetric density. In our previous paper, we selected AB in silicone oil as a role model for a slurry hydrogen storage system. Materials engineering properties were optimized by increasing solid loading by using an ultra-sonic process. In this paper, we proceeded to scale up to liter size batches with solid loadings up to 50 wt% (8 wt% H2) with dynamic viscosities less than 1000cP at 25°C. The use of a non-ionic surfactant, Triton X-15, shows significant promise in controlling the level of foaming produced during the thermal dehydrogenation of the AB. Through the development of new and efficient processing techniques and the ability to adequately control the foaming, stable homogenous slurries of high solid loading have been demonstrated as a viable hydrogen delivery source.

ammonia borane↗
DOE OSTI2022

TEA of a Unique Two-Pathways Process for Post-Combustion CO 2 Capture

A unique two-Pathways process using aqueous sodium glycinate for CO 2 capture from a split flue gas stream emitted from 600 MWe post-combustion coal power plant was developed in Aspen Plus v.10. The split gas flow rate used was 44.75 ton/h and contained 0.0023 mol% SO 2 and 13.33 mol% CO 2 . The process includes a washing unit, a CO 2 absorption unit, a reverse osmosis unit, and a solvent regeneration unit or an ultrafiltration unit. The washing unit uses deionized water to completely remove SO 2 and the CO 2 absorption unit uses SGS to capture at least 90 mol% of the CO 2 in the split flue gas stream. Upon CO 2 and SGS reactions, the resulting liquid products exhibit phase-separation into CO 2 -lean phase and CO 2 -rich phase, allow two distinct pathways. Pathway (i) is to regenerate mostly the CO 2 -rich phase, collect the released CO 2 , and compress it for sequestration purposes. Pathway (ii) is to send the liquid stream from the CO 2 absorption unit to the ultrafiltration unit to separate the solid nanomaterials. The hydraulics and mass transfer characteristics in the washing and CO 2 absorption units were obtained; and techno-economic analysis (TEA) for Pathways (i) and (ii), including Capital Expenditure (CAPEX), Operating Expenditure (OPEX), and Levelized Cost of CO 2 Captured (LCOC), were calculated and compared. The simulation results revealed that the CAPEX, OPEX, and LCOC for Pathway (i) were ($\$12,039,251$), (261 dollar/h), and (54.01 dollar/ton-CO 2 captured), respectively, and those for Pathway (ii) were ($\$5,908,000$), (237.2 dollar/h), and (39.90 dollar/ton-CO 2 captured), respectively. Moreover, in Pathway (ii), 8.19 ton/h of CO 2 were captured to produce 15.62 ton/h NaHCO 3 nanomaterials, which were sold to offset the overall process cost. The LCOC values indicate that Pathway (ii) is more cost-effective than Pathway (i) because LCOC values for Pathway (ii) are much lower than those for Pathway (i).COC values for Pathway (ii) are much lower than those for Pathway (i).

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Exploratory Investigation of Coal in Nonequilibrium Plasma

Coal is an abundant natural resource and there is motivation to find new uses for it that do not intrinsically involve combustion. One approach is to explore new ways of processing coal, and in this work, we focus on the transformation of coal in a nonequilibrium plasma generated from an equimolar mixture of nitrogen and hydrogen. The outcome of the nonequilibrium plasma reaction is fundamentally different than a thermal control reaction carried out using the same gas composition, pressure, and temperature range. The nonequilibrium plasma produces a gas mixture that is enriched in acetylene and its derivatives. Furthermore, when compared to the thermal control experiment, the solid char byproduct of the nonequilibrium plasma has a very reactive surface and is spontaneously combustible at ambient temperature. Experiments performed to characterize the reaction kinetics of coal in the plasma suggest that the mechanism proceeds through a sequential process by which the coal particle temperature rises to a point where devolatilization can occur, the devolatilization reaction happens, followed by parallel reactions of released organic vapors in the plasma phase and surface activation. In conclusion, the reaction rate appears to be limited by the time it takes for the coal particle temperature to rise, consistent with previous results reported for reactions of coal in thermal plasma.

Acetylene↗
DOE OSTI2025

Investigation of CO2 Sorption Characteristics of Imidazolium Ionic Liquid Crystals

This work is focused on advancing of the understanding the carbon dioxide (CO2) sorption characteristics of imidazolium ionic liquid crystals based on the 1-alkyl-3-methylimidazolium (Cnmim) ionic liquid family. Specifically, the mechanism of sorption across the melting, liquid crystal, and isotropic liquid transitions is investigated for 1-tetradecyl-3-methylimidazolium tetrafluoroborate (BF4-) and hexafluorophosphate (PF6-). A unique reversible CO2 occlusion was observed in 1-tetradecyl-3-methylimidazolium tetrafluoroborate (C14mim BF4) under rapid temperature quench from the liquid phases to the solid phase, while complete desorption occurred upon heating above the melting point (40 degrees C). In situ wide-angle x-ray scattering revealed that the lamellar structure of C14mim BF4 liquid crystal phases provides a high concentration of dynamically accessible anion and cation interaction sites for CO2, facilitating CO2 occlusion upon rapid solidification. The selective occlusion of CO2 in C14mim BF4 suggests a new mechanism for CO2 sorption processes.

08 HYDROGEN↗
DOE OSTI2021

Material Design Strategy for Halide Solid Electrolytes Li 3 MX 6 (X = Cl, Br, and I) for All-Solid-State High-Voltage Li-Ion Batteries

Although several solid electrolyte (SE) candidates have been explored, achieving the necessary combination of performance, stability, and processability has been challenging. Recently, several lithium ternary halides have attracted increasing attention for SEs because of their favorable combination of high ionic conductivity and wide electrochemical window. This study aims to provide a material design strategy for lithium halides Li 3 MX 6 (X = Cl, Br, and I) for high-voltage all-solid-state Li-ion batteries, achieved by the systematic investigation of crystal structures, phase and electrochemical stabilities, electronic and mechanical properties, and ionic conductivities. Calculation results reveal that the electronegativity difference between M and X affects structural properties and stabilities. Weak Coulomb interactions in Li 3 MX 6 result in the preference of the monoclinic phase, and the oxidation potential and chemical stability against the cathode materials of Li 3 MX 6 increase for relatively small X. Chlorides exhibit the highest oxidation potential (~4.3 V) among Li 3 MX 6 , suggesting that chlorides are appropriate SEs for high-voltage cathodes. The band gap and elastic moduli increase for relatively small X, suggesting the relatively low electronic conductivity and elastic deformability of chlorides. Chlorides with transition metals typically exhibit trigonal phases, a wider electrochemical stability window, a larger band gap, and higher elastic moduli compared to other types of halides. Additionally, chloride Li 3 MCl 6 is expected to have relatively high ionic conductivities with the aliovalent substitution of M 3+ to Zr 4+ and the anion mixing of Cl with Br. The findings of this study will provide fundamental guidelines for the development of lithium halide SEs for high-voltage all-solid-state Li-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2015

On the stability of surface-confined nanoporous molecular networks

Self-assembly of molecular building blocks into two-dimensional nanoporous networks has been a topic of broad interest for many years. However, various factors govern the specific outcome of the self-assembly process, and understanding and controlling these are key to successful creation. In this work, the self-assembly of two alkylated dehydrobenzo[12]annulene building blocks was compared at the liquid-solid interface. It turned out that only a small chemical modification within the building blocks resulted in enhanced domain sizes and stability of the porous packing relative to the dense linear packing. Applying a thermodynamic model for phase transition revealed some key aspects for network formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Modeling Mesoscale Processes of Scalable Synthesis (Final Report)

The main objective of the center is to develop tools that can help for the understanding of equilibrium and non-equilibrium properties of nanoparticles in aqueous electrolytes. To achieve this goal, we need to develop new computational methodologies that can integrate molecular level description such as the electronic structure analysis to mesoscale description such as generalized hydrodynamics. The challenge is to develop systematic coarse graining techniques that can adaptively choose the level of description needed for the particular task. For example, if bonding breaking is involved, then one should use appropriate models from quantum mechanics to accurately model the time scale and molecular mechanism for the breaking process. However, to model transport process or the process of self-assembly, much less detailed models can be used. There are two main components to this task. One is to develop accurate models that bridge different levels of description. The classical example is the QM/MM model that bridges quantum mechanics and molecular mechanics for studying chemical reactions of large molecules. This work was awarded the Nobel prize this year. For our purpose, we also need to develop methods that bridge molecular and continuum levels of description. Another main task is to develop systematic coarse-graining methods that can give rise to accurate mesoscale and macroscale models, starting from more detailed models at the atomic or molecular scale. None of these ideas are new. By now they have been fashionable topics for quite some time. But neither problem is satisfactorily resolved. The main problem is that uncontrolled approximations are used. As a result, we have very little idea about the validity and accuracy of these models. The main purpose of the center is to exploit the Mori-Zwanzig formalism to achieve the goals stated above. The Mori-Zwanzig formalism is a general tool for eliminating degrees of freedom from existing models. What is attractive about this formalism is its generality. The flip side is that it is also somewhat of a tautology. In order to make it a practical tool, one has to make approximations and this is where the trouble comes. Instead of focusing on Mori-Zwanzig formalism, the PI proposes to study the closely related renormalization group formalism (RG). Like Mori-Zwanzig, RG is a general strategy for eliminating small scale degrees of freedom in a model. Unlike Mori-Zwanzig, RG has the advantage that at least for some problems, one can write a PDE (partial differential equation) type of model for the exact RG flow, without any approximation. This gives us a solid starting point for discussing and analyzing approximations. The PI studied approximate RG flows for turbulent diffusion and related problems. However, extending it to more practical problems proved to be hard. The second line that the PI has pursued is to develop special numerical algorithms for studying particularly important events for dynamics processes in meso-scale systems. The PI succeeded in doing so for the problem of understanding phase transition in metallic systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Revolutionizing observations and predictability of Arctic system dynamics through next-generation dense, heterogeneous and intelligent wireless sensor networks with embedded AI

The Arctic region is a complex and dynamic system, where changes in temperature and partitioning between water in solid, liquid and vapor phase are fundamentally reshaping hydro-biogeochemical fluxes from the bedrock to the atmosphere. The increasing temperature and frequency of extreme heat, precipitation and fire events (Wang et al., 2017; Meredith et al., 2019) have critical consequences for the Arctic ecosystems and people, as well as for the global climate system. Although many observational and remote sensing activities have been on-going, detailed ecosystem processes – particularly related to controls on and impacts of permafrost thaw – are difficult to detect or measure. To improve the understanding and prediction of these extreme events and their consequences, we need a revolution in the way we collect ground-based observations from hillslope to pan-Arctic scales and couple them with satellite imagery and models to deliver rich, actionable, and scale-appropriate data.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations

We study the phase equilibrium between liquid water and ice Ih modeled by the TIP4P/Ice interatomic potential using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice Ih-liquid free energy differences from simulations that visit reversibly both phases. The reversible interconversion is achieved by introducing a static bias potential as a function of an order parameter. The order parameter was tailored to crystallize the hexagonal diamond structure of oxygen in ice Ih. Furthermore, we analyze the effect of the system size on the ice Ih-liquid free energy differences, and we obtain a melting temperature of 270 K in the thermodynamic limit. This result is in agreement with estimates from thermodynamic integration (272 K) and coexistence simulations (270 K). Since the order parameter does not include information about the coordinates of the protons, the spontaneously formed solid configurations contain proton disorder as expected for ice Ih.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗