Quantum mechanical metric for internal cohesion in cement crystals

Porosity modeling in a TiNbTaZrMo high-entropy alloy for biomedical applications

High-entropy alloys (HEAs) have attracted great attention for many biomedical applications. However, the nature of interatomic interactions in this class of complex multicomponent alloys is not fully understood. We report, for the first time, the results of theoretical modeling for porosity in a large biocompatible HEA TiNbTaZrMo using an atomistic supercell of 1024 atoms that provides new insights and understanding. Our results demonstrated the deficiency of using the valence electron count, quantification of large lattice distortion, validation of mechanical properties with available experimental data to reduce Young's modulus. We utilized the novel concepts of the total bond order density (TBOD) and partial bond order density (PBOD) via ab initio quantum mechanical calculations as an effective theoretical means to chart a road map for the rational design of complex multicomponent HEAs for biomedical applications.

First-Principles Calculations of Thermoelectric Transport Properties of Quaternary and Ternary Bulk Chalcogenide Crystals

Chalcogenide crystals have a wide range of applications, especially as thermoelectric materials for energy conversion. Thermoelectric materials can be used to generate an electric current from a temperature gradient based on the Seebeck effect and based on the Peltier effect, and they can be used in cooling applications. Using first-principles calculations and semiclassical Boltzmann theory, we have computed the Seebeck coefficient, electrical conductivity, electronic thermal conductivity, power factor, and figure of merit of 30 chalcogenide crystals. A Quantum Espresso package is used to calculate the electronic properties and locate the Fermi level. The transport properties are then calculated using the BoltzTraP code. The 30 crystals are divided into two groups. The first group has four crystals with quaternary composition (A₂BCQ₄) (A = Tl; B = Cd, Hg; C = Si, Ge, Sn; Q = S, Se, Te). The second group contains 26 crystals with the ternary composition (A’B’Q₂) (A’ = Ag, Cu, Au, Na; B’ = B, Al, Ga, In; Q = S, Se, Te). Among these 30 chalcogenide crystals, the results for 11 crystals: Tl₂CdGeSe₄, Tl₂CdSnSe₄, Tl₂HgSiSe₄, Tl₂HgSnS₄, AuBSe₂, AuBTe₂, AuAlTe₂, AuGaTe₂, AuInTe₂, AgAlSe₂, and AgAlTe₂ are revealed for the first time. In addition, temperature-dependent transport properties of pure and doped AgSbSe₂ and AgSbTe₂ crystals with dopant compositions of AgSb_0.94Cd_0.06Te₂ and AgSbTe_1.85Se_0.15 were explored. These results provide an excellent database for bulk chalcogenides crucial for a wide range of potential applications in renewable energy fields.

Tunable Chemical Disorder in Concentrated Alloys: Defect Physics and Radiation Performance

The development of advanced structural alloys with performance meeting the requirements of extreme environments in nuclear reactors has been long pursued. In the long history of alloy development, the search for metallic alloys with improved radiation tolerance or increased structural strength has relied on either incorporating alloying elements at low concentrations to synthesize so-called dilute alloys or incorporating nanoscale features to mitigate defects. In contrast to traditional approaches, recent success in synthesizing multicomponent concentrated solid-solution alloys (CSAs), including medium-entropy and high-entropy alloys, has vastly expanded the compositional space for new alloy discovery. Their wide variety of elemental diversity enables tunable chemical disorder and sets CSAs apart from traditional dilute alloys. The tunable electronic structure critically lowers the effectiveness of energy dissipation via the electronic subsystem. The tunable chemical complexity also modifies the scattering mechanisms in the atomic subsystem that control energy transport through phonons. The level of chemical disorder depends substantively on the specific alloying elements, rather than the number of alloying elements, as the disorder does not monotonically increase with a higher number of alloying elements. To go beyond our knowledge based on conventional alloys and take advantage of property enhancement by tuning chemical disorder, this review highlights synergistic effects involving valence electrons and atomic-level and nanoscale inhomogeneity in CSAs composed of multiple transition metals. Understanding of the energy dissipation pathways, deformation tolerance, and structural stability of CSAs can proceed by exploiting the equilibrium and non-equilibrium defect processes at the electronic and atomic levels, with or without microstructural inhomogeneities at multiple length scales. Knowledge of tunable chemical disorder in CSAs may advance the understanding of the substantial modifications in element-specific alloy properties that effectively mitigate radiation damage and control a material's response in extreme environments, as well as overcome strength-ductility trade-offs and provide overarching design strategies for structural alloys.

Solvent Effect on the Structure and Properties of RGD Peptide (1FUV) at Body Temperature (310 K) Using Ab Initio Molecular Dynamics

The structure and properties of the arginine-glycine-aspartate (RGD) sequence of the 1FUV peptide at 0 K and body temperature (310 K) are systematically investigated in a dry and aqueous environment using more accurate ab initio molecular dynamics and density functional theory calculations. The fundamental properties, such as electronic structure, interatomic bonding, partial charge distribution, and dielectric response function at 0 and 310 K are analyzed, comparing them in dry and solvated models. These accurate microscopic parameters determined from highly reliable quantum mechanical calculations are useful to define the range and strength of complex molecular interactions occurring between the RGD peptide and the integrin receptor. The in-depth bonding picture analyzed using a novel quantum mechanical metric, the total bond order (TBO), quantifies the role played by hydrogen bonds in the internal cohesion of the simulated structures. The TBO at 310 K decreases in the dry model but increases in the solvated model. These differences are small but extremely important in the context of conditions prevalent in the human body and relevant for health issues. Our results provide a new level of understanding of the structure and properties of the 1FUV peptide and help in advancing the study of RGD containing other peptides.

DFT Study of Electronic Structure and Optical Properties of Kaolinite, Muscovite, and Montmorillonite

Clay mineral materials have attracted attention due to their many properties and applications. The applications of clay minerals are closely linked to their structure and composition. In this paper, we studied the electronic structure properties of kaolinite, muscovite, and montmorillonite crystals, which are classified as clay minerals, by using DFT-based ab initio packages VASP and the OLCAO. The aim of this work is to have a deep understanding of clay mineral materials, including electronic structure, bond strength, mechanical properties, and optical properties. It is worth mentioning that understanding these properties may help continually result in new and innovative clay products in several applications, such as in pharmaceutical applications using kaolinite for their potential in cancer treatment, muscovite used as insulators in electrical appliances, and engineering applications that use montmorillonite as a sealant. In addition, our results show that the role played by hydrogen bonds in O-H bonds has an impact on the hydration in these crystals. Based on calculated total bond order density, it is concluded that kaolinite is slightly more cohesive than montmorillonite, which is consistent with the calculated mechanical properties.

Phase Composition of Silica Fume-Portland Cement Systems Formed under Hydrothermal Curing Evaluated by FTIR, XRD, and TGA

Two substitution levels of Portland cement by silica fume (SF; 30 and 50 mass%) and three hydrothermal treatment regimes (0.5, 1.2, and 2 MPa and 165, 195, and 220 °C for 7 days, respectively) were selected for the investigation of high-temperature phase formation. A combination of thermogravimetric, X-ray diffraction, and Fourier transform infrared analyses in the mid-IR region was used to overcome the shortcomings of individual techniques for the identification of these complex systems. Changes in molecular water amounts, the polymerization degree of silicate chains, or their decomposition due to transformations and crystallization of phases at hydrothermal conditions were observed and discussed concerning composition. Contrary to the calciochondrite, hydrogrossular phases, α-C₂SH, and jaffeite detected in the systems without SF, a decrease in CaO/SiO₂ ratio resulted in the formation of stable tobermorite in the case of 30 mass% SF, whilst calcium hydrogen silicate, gyrolite, and cowlesite were identified as more thermally stable phases in the samples with 50 mass% SF.

Structural and physical properties of 99 complex bulk chalcogenides crystals using first-principles calculations

Chalcogenide semiconductors and glasses have many applications in the civil and military fields, especially in relation to their electronic, optical and mechanical properties for energy conversion and in enviormental materials. However, they are much less systemically studied and their fundamental physical properties for a large class chalcogenide semiconductors are rather scattered and incomplete. Here, we present a detailed study using well defined first-principles calculations on the electronic structure, interatomic bonding, optical, and mechanical properties for 99 bulk chalcogenides including thirteen of these crytals which have never been calculated. Due to their unique composition and structures, these 99 bulk chalcogenides are divided into two main groups. The first group contains 54 quaternary crystals with the structure composition (A₂BCQ₄) (A = Ag, Cu; B = Zn, Cd, Hg, Mg, Sr, Ba; C = Si, Ge, Sn; Q = S, Se, Te), while the second group contains scattered ternary and quaternary chalcogenide crystals with a more diverse composition (A_xB_yC_zQ_n) (A = Ag, Cu, Ba, Cs, Li, Tl, K, Lu, Sr; B = Zn, Cd, Hg, Al, Ga, In, P, As, La, Lu, Pb, Cu, Ag; C = Si, Ge, Sn, As, Sb, Bi, Zr, Hf, Ga, In; Q = S, Se, Te; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {x} = 1$$\end{document}x=1, 2, 3; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {y} = 0$$\end{document}y=0, 1, 2, 5; \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {z} = 0$$\end{document}z=0, 1, 2 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {n} = 3$$\end{document}n=3, 4, 5, 6, 9). Moreover, the total bond order density (TBOD) is used as a single quantum mechanical metric to characterize the internal cohesion of these crystals enabling us to correlate them with the calculated properties, especially their mechanical properties. This work provides a very large database for bulk chalcogenides crucial for the future theoretical and experimental studies, opening opportunities for study the properties and potential application of a wide variety of chalcogenides.

Ab Initio Study of Hydrolysis Effects in Single and Ion-Exchanged Alkali Aluminosilicate Glasses

Hydrolysis in alkali-doped aluminosilicate glasses is one of the most complicated mechanisms in glass science. There remain many fundamental and unresolved issues with implications on their potential applications. Herein, we address this challenge by carrying out detailed calculations on the structure and properties of both anhydrate (dry) and hydrated alkali aluminosilicate glasses on carefully constructed models. Specifically, the Na-, (Na + K)-, and K-doped aluminosilicate glasses with compositions (SiO₂)_0.6(Al₂O₃)_0.2(Na₂O)_0.2 - x(K₂O)_x (x = 0, 0.10, and 0.20) are simulated using ab initio molecular dynamics (AIMD). The local short- and intermediate-range order in these glasses is analyzed in terms of atomic pair distribution, coordination number, bond length, and bond angle distributions to delineate the subtle variations due to different alkali sizes and hydrolysis. The electronic structure, interatomic bonding, mechanical, and optical properties for these models are calculated and validated with available experimental data. We use the novel concept of total bond order density (TBOD), the quantum mechanically derived metric, to characterize the internal cohesion and strength in the simulated glasses. Detailed analysis of the hydrolysis mechanism enables us to provide information on the complex interplay of various participating elements and their interactions at the atomic level. Such detailed information provides a new platform of knowledge, which is crucial for understanding the issues related to glass corrosion and durability, and ways and means for their special applications in commercial glass products. Both undissociated molecular water and dissociated water in the form of hydroxyl groups exist in the hydrated models in the presence of alkali ions. For the first time, we observed the opposite mixed alkali effect in the Poisson's ratio for anhydrate and hydrated glasses.

Interfacial Interaction between Suolunite Crystal and Silica Binding Peptide for Novel Bioinspired Cement

Cement and concrete have been important construction materials throughout human history. There is an urgent need to explore novel and untraditional cementitious materials to enhance the durability of building materials and structures in response to increased infrastructure demand worldwide. We report an exploratory study on a biocomposite cement based on a large-scale computational study using density functional theory. An explicitly solvated mixture of a mineral calcium silicate hydrate (C-S-H) crystal suolunite (Ca₂Si₂O₅(OH)₂·H₂O) and a silicon binding peptide with amino acid sequence PRO-PRO-PRO-TRP-LEU-PRO-TYR-MET-PRO-PRO-TRP-SER is constructed using ab initio molecular dynamics (AIMD). Detailed analysis on the interface structure, interatomic bonding, mechanical properties, and solvent effect of this model reveals a complex interplay of different types of covalent and ionic bonding, including ubiquitous hydrogen bonding which plays a crucial role in their properties. The use of the total bond order density (TBOD), a single quantum mechanical metric, for assessing the interfacial cohesion for this composite biocement is proposed. We find that the solvated model has a slightly larger TBOD than the dried one. These results could lead to a systematic search and rational design for different types of bioinspired and hybrid functional materials with other inorganic minerals and organic peptides.

Intermolecular dynamics of ultraconfined interlayer water in tobermorite: influence on mechanical performance

Ultraconfined interlayer water within the tobermorite molecular structure is responsible for changes in the uniaxial tensile and compressive response of the family of tobermorites: 9, 11 and 14 Å. These confined interlayer water molecules are engaged in solvation of cations and anions within the tobermorite structure, which has been demonstrated through the intermolecular vibrational spectra and hydrogen bond lifetime of the water molecules. This study demonstrates that instead of ionization of water molecules (as proposed in an earlier study), breaking of hydrogen bonds of water is more plausible leading to solvation of ions within the molecular structure of tobermorite. A schematic of the coordinate covalent bonds between the water molecules and the cations and anions of the tobermorite structure has been proposed in this study.

Deformation behavior of an amorphous zeolitic imidazolate framework - from a supersoft material to a complex organometallic alloy

Zeolitic imidazolate frameworks (ZIFs)-a subset of metal-organic frameworks (MOFs)-have recently attracted immense attention. Many crystalline ZIFs (c-ZIFs) have highly porous zeolite structures that are ideal for molecular encapsulation. Recently emerging non-crystalline or amorphous ZIFs (a-ZIFs) with a similar short-range order are of interest because they can be converted from c-ZIFs for large-scale production. Here, we present a computational study of the deformation behavior of a unique a-ZIF model by simulating step-wise compression and expansion with strains between -0.389 and +0.376. An insulator-to-metal transition is observed at 51 GPa leading to a multicomponent light amorphous alloy of only 3.68 g (cm)-3. A high-density amorphous-to-amorphous phase transition is observed due to the sudden formation of N-N bond pairs. The systematic expansion of the a-ZIF retains the framework softness until it fractures at high strain. Based on the expansion data, we propose an empirical formula for super-soft materials, which is in line with available experimental data.

Multiscale understanding of tricalcium silicate hydration reactions

Tricalcium silicate, the main constituent of Portland cement, hydrates to produce crystalline calcium hydroxide and calcium-silicate-hydrates (C-S-H) nanocrystalline gel. This hydration reaction is poorly understood at the nanoscale. The understanding of atomic arrangement in nanocrystalline phases is intrinsically complicated and this challenge is exacerbated by the presence of additional crystalline phase(s). Here, we use calorimetry and synchrotron X-ray powder diffraction to quantitatively follow tricalcium silicate hydration process: i) its dissolution, ii) portlandite crystallization and iii) C-S-H gel precipitation. Chiefly, synchrotron pair distribution function (PDF) allows to identify a defective clinotobermorite, Ca₁₁Si₉O₂₈(OH)₂^.8.5H₂O, as the nanocrystalline component of C-S-H. Furthermore, PDF analysis also indicates that C-S-H gel contains monolayer calcium hydroxide which is stretched as recently predicted by first principles calculations. These outcomes, plus additional laboratory characterization, yielded a multiscale picture for C-S-H nanocomposite gel which explains the observed densities and Ca/Si atomic ratios at the nano- and meso- scales.

Implication of the solvent effect, metal ions and topology in the electronic structure and hydrogen bonding of human telomeric G-quadruplex DNA

We present a first-principles density functional study elucidating the effects of solvent, metal ions and topology on the electronic structure and hydrogen bonding of 12 well-designed three dimensional G-quadruplex (G4-DNA) models in different environments. Our study shows that the parallel strand structures are more stable in dry environments and aqueous solutions containing K(+) ions within the tetrad of guanine but conversely, that the anti-parallel structure is more stable in solutions containing the Na(+) ions within the tetrad of guanine. The presence of metal ions within the tetrad of the guanine channel always enhances the stability of the G4-DNA models. The parallel strand structures have larger HOMO-LUMO gaps than antiparallel structures, which are in the range of 0.98 eV to 3.11 eV. Partial charge calculations show that sugar and alkali ions are positively charged whereas nucleobases, PO4 groups and water molecules are all negatively charged. Partial charges on each functional group with different signs and magnitudes contribute differently to the electrostatic interactions involving G4-DNA and favor the parallel structure. A comparative study between specific pairs of different G4-DNA models shows that the Hoogsteen OH and NH hydrogen bonds in the guanine tetrad are significantly influenced by the presence of metal ions and water molecules, collectively affecting the structure and the stability of G4-DNA.

Complex interplay of interatomic bonding in a multi-component pyrophosphate crystal: K2Mg (H2P2O7)2·2H2O

The electronic structure and interatomic bonding of pyrophosphate crystal K₂Mg (H₂P₂O₇)₂·2H₂O are investigated for the first time showing complex interplay of different types of bindings. The existing structure from single-crystal X-ray diffraction is not sufficiently refined, resulting in unrealistic short O─H bonds which is rectified by high-precision density functional theory (DFT) calculation. K₂Mg (H₂P₂O₇)₂·2H₂O has a direct gap of 5.22 eV and a small electron effective mass of 0.14 m_e. Detailed bond analysis between every pair of atoms reveals the complexity of various covalent, ionic, hydrogen bonding and bridging bonding and their sensitive dependence on structural differences. The K--O bonds are much weaker than Mg--O bonds and contributions from the hydrogen bonds are non-negligible. Quantitative analysis of internal cohesion in terms of total bond order density and partial bond order density divulges the relative importance of different types of bonding. The calculated optical absorptions show multiple peaks and a sharp Plasmon peak at 23 eV and a refractive index of 1.44. The elastic and mechanical properties show features unique to this low-symmetry crystal. Phonon calculation gives vibrational frequencies in agreement with reported Raman spectrum. These results provide new insights indicating that acidic pyrophosphates could have a variety of unrealized applications in advanced technology.

Ab Initio Modeling of Structure and Properties of Single and Mixed Alkali Silicate Glasses

A density functional theory (DFT)-based ab initio molecular dynamics (AIMD) has been applied to simulate models of single and mixed alkali silicate glasses with two different molar concentrations of alkali oxides. The structural environments and spatial distributions of alkali ions in the 10 simulated models with 20% and 30% of Li, Na, K and equal proportions of Li-Na and Na-K are studied in detail for subtle variations among the models. Quantum mechanical calculations of electronic structures, interatomic bonding, and mechanical and optical properties are carried out for each of the models, and the results are compared with available experimental observation and other simulations. The calculated results are in good agreement with the experimental data. We have used the novel concept of using the total bond order density (TBOD), a quantum mechanical metric, to characterize internal cohesion in these glass models. The mixed alkali effect (MAE) is visible in the bulk mechanical properties but not obvious in other physical properties studied in this paper. We show that Li doping deviates from expected trend due to the much stronger Li-O bonding than those of Na and K doping. The approach used in this study is in contrast with current studies in alkali-doped silicate glasses based only on geometric characterizations.

Densification of the interlayer spacing governs the nanomechanical properties of calcium-silicate-hydrate

Calciuam-silicate-hydrate (C-S-H) is the principal binding phase in modern concrete. Molecular simulations imply that its nanoscale stiffness is ‘defect-driven’, i.e., dominated by crystallographic defects such as bridging site vacancies in its silicate chains. However, experimental validation of this result is difficult due to the hierarchically porous nature of C-S-H down to nanometers. Here, we integrate high pressure X-ray diffraction and atomistic simulations to correlate the anisotropic deformation of nanocrystalline C-S-H to its atomic-scale structure, which is changed by varying the Ca-to-Si molar ratio. Contrary to the ‘defect-driven’ hypothesis, we clearly observe stiffening of C-S-H with increasing Ca/Si in the range 0.8 ≤ Ca/Si ≤ 1.3, despite increasing numbers of vacancies in its silicate chains. The deformation of these chains along the b-axis occurs mainly through tilting of the Si-O-Si dihedral angle rather than shortening of the Si-O bond, and consequently there is no correlation between the incompressibilities of the a- and b-axes and the Ca/Si. On the contrary, the intrinsic stiffness of C-S-H solid is inversely correlated with the thickness of its interlayer space. This work provides direct experimental evidence to conduct more realistic modelling of C-S-H-based cementitious material.

Aluminum-induced dreierketten chain cross-links increase the mechanical properties of nanocrystalline calcium aluminosilicate hydrate

The incorporation of Al and increased curing temperature promotes the crystallization and cross-linking of calcium (alumino)silicate hydrate (C-(A-)S-H), which is the primary binding phase in most contemporary concrete materials. However, the influence of Al-induced structural changes on the mechanical properties at atomistic scale is not well understood. Herein, synchrotron radiation-based high-pressure X-ray diffraction is used to quantify the influence of dreierketten chain cross-linking on the anisotropic mechanical behavior of C-(A-)S-H. We show that the ab-planar stiffness is independent of dreierketten chain defects, e.g. vacancies in bridging tetrahedra sites and Al for Si substitution. The c-axis of non-cross-linked C-(A-)S-H is more deformable due to the softer interlayer opening but stiffens with decreased spacing and/or increased zeolitic water and Ca²⁺ of the interlayer. Dreierketten chain cross-links act as ‘columns’ to resist compression, thus increasing the bulk modulus of C-(A-)S-H. We provide the first experimental evidence on the influence of the Al-induced atomistic configurational change on the mechanical properties of C-(A-)S-H. Our work advances the fundamental knowledge of C-(A-)S-H on the lowest level of its hierarchical structure, and thus can impact the way that innovative C-(A-)S-H-based cementitious materials are developed using a ‘bottom-up’ approach.