Formation of porous ice frameworks at room temperature.

Bulk crystalline ices with ultralow densities have been demonstrated to be thermodynamically metastable at negative pressures. However, the direct formation of these bulk porous ices from liquid water at negative pressures is extremely challenging. Inspired by approaches toward porous media based on host-guest chemistry, such as metal-organic frameworks and covalent organic frameworks, we herein demonstrate via molecular dynamics simulations that a class of ultralow-density porous ices with upright channels can be formed spontaneously from liquid water at 300 K with the assistance of carbon nanotube arrays. We refer to these porous ice polymorphs as water oxygen-vertex frameworks (WOFs). Notably, our simulations revealed that the liquid-WOF phase transition is first-order and occurs at room temperature. All the WOFs exhibited the unique structural feature that they can be regarded as assemblies of nanoribbons of hexagonal bilayer ice (2D ice I) at their armchair or zigzag edges. Based on density functional theory calculations, a comprehensive phase diagram of the WOFs was constructed considering both the thermodynamic and thermal stabilities of the porous ices at negative pressures. Like other types of porous media, these WOFs may be applicable to gas storage, purification, and separation. Moreover, these biocompatible porous ice networks may be exploited as medical-related carriers.

Thermal decomposition of trimethylindium and indium trisguanidinate precursors for InN growth: An ab initio and kinetic modeling study.

Indium nitride (InN) is an interesting material for future electronic and photonic-related applications, as it combines high electron mobility and low-energy band gap for photoabsorption or emission-driven processes. In this context, atomic layer deposition techniques have been previously employed for InN growth at low temperatures (typically <350 °C), reportedly yielding crystals with high quality and purity. In general, this technique is assumed to not involve any gas phase reactions as a result from the time-resolved insertion of volatile molecular sources into the gas chamber. Nonetheless, such temperatures could still favor the precursor decomposition in the gas phase during the In half-cycle, therefore altering the molecular species that undergoes physisorption and, ultimately, driving the reaction mechanism to pursue other pathways. Thence, we herein evaluate the thermal decomposition of relevant In precursors in the gas phase, namely, trimethylindium (TMI) and tris(N,N'-diisopropyl-2-dimethylamido-guanidinato) indium (III) (ITG), by means of thermodynamic and kinetic modeling. According to the results, at T = 593 K, TMI should exhibit partial decomposition of ∼8% after 400 s to first generate methylindium and ethane (C2H6), a percentage that increases to ∼34% after 1 h of exposure inside the gas chamber. Therefore, this precursor should be present in an intact form to undergo physisorption during the In half-cycle of the deposition (<10 s). On the other hand, the ITG decomposition starts already at the temperatures used in the bubbler, in which it slowly decomposes as it is evaporated during the deposition process. At T = 300 °C, the decomposition is a fast process that reaches 90% completeness after 1 s and where equilibrium, at which almost no ITG remains, is achieved before 10 s. In this case, the decomposition pathway is likely to occur via elimination of the carbodiimide ligand. Ultimately, these results should contribute for a better understanding of the reaction mechanism involved in the InN growth from these precursors.

Synthesis, Structure, and Thermal Properties of Volatile Group 11 Triazenides as Potential Precursors for Vapor Deposition.

Group 11 thin films are desirable as interconnects in microelectronics. Although many M-N-bonded Cu precursors have been explored for vapor deposition, there is currently a lack of suitable Ag and Au derivatives. Herein, we present monovalent Cu, Ag, and Au 1,3-di-tert-butyltriazenides that have potential for use in vapor deposition. Their thermal stability and volatility rival that of current state-of-the-art group 11 precursors with bidentate M-N-bonded ligands. Solution-state thermolysis of these triazenides yielded polycrystalline films of elemental Cu, Ag, and Au. The compounds are therefore highly promising as single-source precursors for vapor deposition of coinage metal films.

An ultralow-density porous ice with the largest internal cavity identified in the water phase diagram.

The recent back-to-back findings of low-density porous ice XVI and XVII have rekindled the century-old field of the solid-state physics and chemistry of water. Experimentally, both ice XVI and XVII crystals can be produced by extracting guest atoms or molecules enclosed in the cavities of preformed ice clathrate hydrates. Herein, we examine more than 200 hypothetical low-density porous ices whose structures were generated according to a database of zeolite structures. Hitherto unreported porous EMT ice, named according to zeolite nomenclature, is identified to have an extremely low density of 0.5 g/cm3 and the largest internal cavity (7.88 Å in average radius). The EMT ice can be viewed as dumbbell-shaped motifs in a hexagonal close-packed structure. Our first-principles computations and molecular dynamics simulations confirm that the EMT ice is stable under negative pressures and exhibits higher thermal stability than other ultralow-density ices. If all cavities are fully occupied by hydrogen molecules, the EMT ice hydrate can easily outperform the record hydrogen storage capacity of 5.3 wt % achieved with sII hydrogen hydrate. Most importantly, in the reconstructed temperature-pressure (T-P) phase diagram of water, the EMT ice is located at deeply negative pressure regions below ice XVI and at higher temperature regions next to FAU. Last, the phonon spectra of empty-sII, FAU, EMT, and other zeolite-like ice structures are computed by using the dispersion corrected vdW-DF2 functional. Compared with those of ice XI (0.93 g/cm3), both the bending and stretching vibrational modes of the EMT ice are blue-shifted due to their weaker hydrogen bonds.

Synthesis, Characterization, and Thermal Study of Divalent Germanium, Tin, and Lead Triazenides as Potential Vapor Deposition Precursors.

Only a few M-N bonded divalent group 14 precursors are available for vapor deposition, in particular for Ge and Pb. A majority of the reported precursors are dicoordinated with the Sn(II) amidinates, the only tetracoordinated examples. No Ge(II) and Pb(II) amidinates suitable for vapor deposition have been demonstrated. Herein, we present tetracoordinated Ge(II), Sn(II), and Pb(II) complexes bearing two sets of chelating 1,3-di-tert-butyltriazenide ligands. These compounds are thermally stable, sublime quantitatively between 60 and 75 °C (at 0.5 mbar), and show ideal single-step volatilization by thermogravimetric analysis.

Synthesis and Thermal Study of Hexacoordinated Aluminum(III) Triazenides for Use in Atomic Layer Deposition.

Amidinate and guanidinate ligands have been used extensively to produce volatile and thermally stable precursors for atomic layer deposition. The triazenide ligand is relatively unexplored as an alternative ligand system. Herein, we present six new Al(III) complexes bearing three sets of a 1,3-dialkyltriazenide ligand. These complexes volatilize quantitatively in a single step with onset volatilization temperatures of ∼150 °C and 1 Torr vapor pressures of ∼134 °C. Differential scanning calorimetry revealed that these Al(III) complexes exhibited exothermic events that overlapped with the temperatures of their mass loss events in thermogravimetric analysis. Using quantum chemical density functional theory computations, we found a decomposition pathway that transforms the relatively large hexacoordinated Al(III) precursor into a smaller dicoordinated complex. The pathway relies on previously unexplored interligand proton migrations. These new Al(III) triazenides provide a series of alternative precursors with unique thermal properties that could be highly advantageous for vapor deposition processes of Al containing materials.

Modeling Molecular Interactions in Water: From Pairwise to Many-Body Potential Energy Functions.

Almost 50 years have passed from the first computer simulations of water, and a large number of molecular models have been proposed since then to elucidate the unique behavior of water across different phases. In this article, we review the recent progress in the development of analytical potential energy functions that aim at correctly representing many-body effects. Starting from the many-body expansion of the interaction energy, specific focus is on different classes of potential energy functions built upon a hierarchy of approximations and on their ability to accurately reproduce reference data obtained from state-of-the-art electronic structure calculations and experimental measurements. We show that most recent potential energy functions, which include explicit short-range representations of two-body and three-body effects along with a physically correct description of many-body effects at all distances, predict the properties of water from the gas to the condensed phase with unprecedented accuracy, thus opening the door to the long-sought "universal model" capable of describing the behavior of water under different conditions and in different environments.

Clathrate ice sL: a new crystalline phase of ice with ultralow density predicted by first-principles phase diagram computations.

In contrast to the rich knowledge of water and 17 experimentally confirmed crystalline phases of solid water under positive pressures, water under negative pressure has been poorly explored. In this study, a new crystalline phase of ice with ultralow density (0.6 g cm-3), named "clathrate ice sL", is constructed by nano water cage clusters, and it is predicted to be stable under a lower negative pressure than the experimentally confirmed sII phase by first-principles phase diagram computations, thereby extending the phase diagram of water to negative pressure regions below -5170 bar at 0 K and below -4761 bar at 300 K. In addition, according to our theoretical prediction, the optimal hydrogen storage mass density in the new clathrate ice sL is 7.7 wt% (larger than the 2017 DOE target of 5.5 wt%), which would set a new record of hydrogen storage capacity in clathrate hydrates. The finding of clathrate ice sL not only proposes a new type of crystalline ice under negative pressure but also explores the potential applications of the ultralow density ice phases while extending the water phase diagram and enriching the knowledge of people about water.

Growth Mechanism of SiC CVD: Surface Etching by H2, H Atoms, and HCl.

Silicon carbide is a wide bandgap semiconductor with unique characteristics suitable for high temperature and high power applications. Fabrication of SiC epitaxial layers is usually performed using chemical vapor deposition (CVD). In this work, we use quantum chemical density functional theory (B3LYP and M06-2X) and transition state theory to study etching reactions occurring on the surface of SiC during CVD in order to combine etching effects to the surface kinetic model for SiC CVD. H2, H atoms and HCl gases are chosen in the study as the most likely etchants responsible for surface etching. We consider etchings of four surface sites, namely CH3(ads), SiH3CH2(ads), SiH2(CH2)2(ads), and SiH(CH2)3(ads), which represent four subsequent snapshots of the surface as the growth proceeds. We find that H atoms are the most effective etchant on CH3(ads) and SiH3CH2(ads), which represent the first and second steps of the growth. HCl and H2 are shown to be much less effective than H atoms and produce the etching rate constants which are ∼104 and ∼107 times slower. In comparison to CH3(ads), SiH3CH2(ads) is shown to be less stable and more susceptible to etchings. Unlike the first and second steps of the growth, the third and fourth steps (i.e., SiH2(CH2)2(ads) and SiH(CH2)3(ads)) are stable and much less susceptible to any of the three etchants considered. This implies that the growth species become more stable via forming Si-C bonds with another surface species. The formation of a larger surface cluster thus helps stabilizing the growth against etchings.

C-C stretching Raman spectra and stabilities of hydrocarbon molecules in natural gas hydrates: a quantum chemical study.

The presence of specific hydrocarbon gas molecules in various types of water cavities in natural gas hydrates (NGHs) are governed by the relative stabilities of these encapsulated guest molecule-water cavity combinations. Using molecular quantum chemical dispersion-corrected hybrid density functional computations, the interaction (ΔE(host--guest)) and cohesive energies (ΔE(coh)), enthalpies, and Gibbs free energies for the complexes of host water cages and hydrocarbon guest molecules are calculated at the ωB97X-D/6-311++G(2d,2p) level of theory. The zero-point energy effect of ΔE(host-guest) and ΔE(coh) is found to be quite substantial. The energetically optimal host-guest combinations for seven hydrocarbon gas molecules (CH4, C2H6, C3H6, C3H8, C4H8, i-C4H10, and n-C4H10) and various water cavities (D, ID, T, P, H, and I) in NGHs are found to be CH4@D, C2H6@T, C3H6@T, C3H8@T, C4H8@T/P/H, i-C4H10@H, and n-C4H10@H, as the largest cohesive energy magnitudes will be obtained with these host-guest combinations. The stabilities of various water cavities enclosing hydrocarbon molecules are evaluated from the computed cohesive Gibbs free energies: CH4 prefers to be trapped in a ID cage; C2H6 prefer T cages; C3H6 and C3H8 prefer T and H cages; C4H8 and i-C4H10 prefer H cages; and n-C4H10 prefer I cages. The vibrational frequencies and Raman intensities of the C-C stretching vibrational modes for these seven hydrocarbon molecules enclosed in each water cavity are computed. A blue shift results after the guest molecule is trapped from gas phase into various water cages due to the host-guest interactions between the water cage and hydrocarbon molecule. The frequency shifts to the red as the radius of water cages increases. The model calculations support the view that C-C stretching vibrations of hydrocarbon molecules in the water cavities can be used as a tool to identify the types of crystal phases and guest molecules in NGHs.

Structures of the I-, II- and H-methane clathrates and the ice-methane clathrate phase transition from quantum-chemical modeling with force-field thermal corrections.

Methane hydrates with the three clathrate structures I, II, and H are studied by quantum-chemical methods. Hybrid density-functional theory B3LYP computations using periodic boundary conditions are combined with force-field methods for the thermal energy effects to calculate energetic, thermodynamic, and structural properties. The pressure dependencies for the crystal structures, lattice energies, and guest molecule interactions are derived. The quantum-chemical geometry optimizations predict too small cell volumes as compared to experimental data, but by including zero-point energy and thermal energy effects, we find the cell volumes increase and the correct densities are obtained. The phase transition from MH-I to ice Ih and methane was computed and found to occur at about 9.7 MPa.

A theoretical study of water equilibria: the cluster distribution versus temperature and pressure for (H2O)n, n = 1-60, and ice.

The size distribution of water clusters at equilibrium is studied using quantum-chemical calculations in combination with statistical thermodynamics. The necessary energetic data is obtained by quantum-chemical B3LYP computations and through extrapolations from the B3LYP results for the larger clusters. Clusters with up to 60 molecules are included in the equilibrium computations. Populations of different cluster sizes are calculated using both an ideal gas model with noninteracting clusters and a model where a correction for the interaction energy is included analogous to the van der Waals law. In standard vapor the majority of the water molecules are monomers. For the ideal gas model at 1 atm large clusters [56-mer (0-120 K) and 28-mer (100-260 K)] dominate at low temperatures and separate to smaller clusters [21-22-mer (170-280 K) and 4-6-mer (270-320 K) and to monomers (300-350 K)] when the temperature is increased. At lower pressure the transition from clusters to monomers lies at lower temperatures and fewer cluster sizes are formed. The computed size distribution exhibits enhanced peaks for the clusters consisting of 21 and 28 water molecules; these sizes are for protonated water clusters often referred to as magic numbers. If cluster-cluster interactions are included in the model the transition from clusters to monomers is sharper (i.e., occurs over a smaller temperature interval) than when the ideal-gas model is used. Clusters with 20-22 molecules dominate in the liquid region. When a large icelike cluster is included it will dominate for temperatures up to 325 K for the noninteracting clusters model. Thermodynamic properties (C(p), DeltaH) were calculated with in general good agreement with experimental values for the solid and gas phase. A formula for the number of H-bond topologies in a given cluster structure is derived. For the 20-mer it is shown that the number of topologies contributes to making the population of dodecahedron-shaped cluster larger than that of a lower-energy fused prism cluster at high temperatures.

Hydrogen bond topology and the ice VII/VIII and Ih/XI proton ordering phase transitions.

Ice Ih, ordinary ice at atmospheric pressure, is a proton-disordered crystal that when cooled under special conditions is believed to transform to ferroelectric proton-ordered ice XI, but this transformation is still subject to controversy. Ice VII, also proton disordered throughout its region of stability, transforms to proton-ordered ice VIII upon cooling. In contrast to the ice Ih/XI transition, the VII/VIII transition and the crystal structure of ice VIII are well characterized. In order to shed some light on the ice Ih proton ordering transition, we present the results of periodic electronic density functional theory calculations and statistical simulations. We are able to describe the small energy differences among the innumerable H-bond configurations possible in a large simulation cell by using an analytic theory to extrapolate from electronic DFT calculations on small unit cells to cells large enough to approximate the thermodynamic limit. We first validate our methods by comparing our predictions to the well-characterized ice VII/VIII proton ordering transition, finding agreement with respect to both the transition temperature and structure of the low-temperature phase. For ice Ih, our results indicate that a proton-ordered phase is attainable at low temperatures, the structure of which is in agreement with the experimentally proposed ferroelectric structure. The predicted transition temperature of is in qualitative agreement with the observed transition at on KOH-doped ice samples.

IR and quantum-chemical studies of carboxylic acid and glycine adsorption on rutile TiO2 nanoparticles.

Nanocrystalline TiO2 powders of the rutile polymorph, synthesized by a sol-gel method, were treated with water solutions containing, respectively, formic, acetic, and citric acid and glycine in order to study the adsorption properties of these organic species. The samples were characterized by FTIR, Raman, powder XRD, and TEM. It was found that HCOOH, CH3COOH and HOC(COOH)(CH2COOH)2--but not NH2CH2COOH--adsorbed onto TiO2. The adsorption of HCOOH, CH3COOH and NH2CH2COOH onto the (110) surface of rutile was also studied by quantum-chemical periodic density functional theory (DFT) calculations. The organic molecules were from the computations found to adsorb strongly to the surfaces in a bridge-coordinating mode, where the two oxygens of the deprotonated carboxylic acid bind to two surface titanium ions. Surface relaxation is found to influence adsorption geometries and energies significantly. The results from DFT calculations and ab initio molecular-dynamics simulations of formic acid adsorption onto TiO2 are compared and match well with the experimental IR measurements, supporting the bridge-binding geometry of carboxylic-acid adsorption on the TiO2 nanoparticles.

Raman and IR Spectra of Ice Ih and Ice XI with an Assessment of DFT Methods.

IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. The performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases is evaluated. Based on a selected DFT computational scheme, the IR and Raman spectra of ice Ih and XI are derived and compared. When the spectra, both IR and Raman, of ice Ih and ice XI are compared, the librational vibrations are found to be the most affected by the proton ordering. The spectroscopic fingerprint of ice XI can be used to distinguish ferroelectric ice XI from ice Ih in the universe. Furthermore, the existence of only one kind of H-bond in ice Ih is demonstrated from the overlapping subspectra for different types of H-bonded pair configurations in 16 isomers of ice Ih, which provides an illustration to the historic debate on whether one or two kinds of H-bonds existed in ice.

Fingerprints in IR OH vibrational spectra of H2O clusters from different H-bond conformations by means of quantum-chemical computations.

The thermodynamic stabilities and IR spectra of the three water clusters (H2O)20, (H2O)54,, and (H2O)100 are studied by quantum-chemical computations. After full optimization of the (H2O)20,54,100 structures using the hybrid density functional B3LYP together with the 6-31+G(d,p) basis set, the electronic energies, zero-point energies, internal energies, enthalpies, entropies, and Gibbs free energies of the water clusters at 298 K are investigated. The OH stretching vibrational IR spectra of (H2O)20,54,100 are simulated and split into sub-spectra for different H-bond groups depending on the conformations of the hydrogen bonds. From the computed spectra the different spectroscopic fingerprint features of water molecules in different H-bond conformations in the water clusters are inferred.

Theoretical IR spectra for water clusters (H2O)n (n = 6-22, 28, 30) and identification of spectral contributions from different H-bond conformations in gaseous and liquid water.

The vibrational IR spectra in the O-H stretching region are computed for water clusters containing 6-22, 28, and 30 molecules using quantum-chemical calculations (B3LYP and an augmented basis set). For the cluster with 20 molecules, several different structures were studied. The vibrational spectrum was partitioned into contributions from different molecules according to their coordination properties. The frequency shifts depend on the number of donated/accepted H-bonds primarily of the two molecules participating in the H-bond, but also of the surrounding molecules H-bonding to these molecules. The frequencies of H-bonds between two molecules of the same coordination type are spread over a broad interval. The most downshifted hydrogen-bond vibrations are those donated by a single-donor 3-coordinated molecule where the H-bond is accepted by a single-acceptor molecule. The H-bonded neighbors influence the downshift, and their contribution can be rationalized in the same way as for the central dimer. Single donors/acceptors cause larger downshifts than 4-coordinated molecules, and the least downshift is obtained for double donors/acceptors. This result is at variance with the conception that experimental liquid water spectra may be divided into components for which larger downshifts imply higher H-bond coordination. A mean spectral contribution for each coordination type for the donor molecule was derived and fitted to the experimental liquid water IR spectrum, which enabled an estimation of the distribution of H-bond types and average number of H-bonds (3.0 +/- 0.2) in the liquid.

Towards biocompatibility of RE2O3 nanocrystals - water and organic molecules chemisorbed on Gd2O3 and Y2O3 nanocrystals studied by quantum-chemical computations.

Nanocrystals of Gd2O3/Y2O3 and their interaction with water, formic acid, diethylene glycol (DEG), and tetramethoxy silane (TMOS) have been studied by quantum-chemical calculations at the B3LYP level using solvent-coated clusters of gadolinia and yttria. Adsorption energies, surface geometries, electronic structures, and excitation spectra were calculated. The results concerning adsorption strengths and superparamagnetic high-spin states can provide insight into the design of molecular-capped RE2O3 nanocrystals to be used in vivo.

The local structure of protonated water from x-ray absorption and density functional theory.

We present a combined x-ray absorption spectroscopy/computational study of water in hydrochloric acid (HCl) solutions of varying concentration to address the structure and bonding of excess protons and their effect on the hydrogen bonding network in liquid water. Intensity variations and energy shifts indicate changes in the hydrogen bonding structure in water as well as the local structure of the protonated complex as a function of the concentration of protons. In particular, in highly acidic solutions we find a dominance of the Eigen form, H(3)O(+), while the proton is less localized to a specific water under less acidic conditions.