Thermal Stability of Organic Monolayers Grafted to Si(111): Insights from ReaxFF Reactive Molecular Dynamics Simulations

Reactive Simulations of Silica Functionalization with Aromatic Hydrocarbons

Reactive molecular dynamics simulations are used to model the covalent functionalization of amorphous silica with aromatic hydrocarbons. Simulations show that the surface density of silanol-terminated phenyl, naphthyl, and anthracenyl molecules is lower than the maximum value calculated based on molecule geometry, and the simulation densities decrease faster with the number of aromatic rings than the geometric densities. The trends are analyzed in terms of the surface-silanol bonding configurations, tilt angles, local conformational ordering, and aggregation of surface-bound molecules under steady-state conditions. Results show that the surface density is affected by both the size and symmetry of the aromatic hydrocarbons. The correlations among bonding, orientation, and surface density identified here may guide the selection or design of molecules for functionalized surfaces.

Heat-Treated Aramid Pulp/Silica Aerogel Composites with Improved Thermal Stability and Thermal Insulation

In this work, we prepared heat-treated aramid pulp/silica aerogel composites (AP/aerogels) and investigated in detail the feasibility of improving thermal stability and thermal insulation via tailored heat treatment. The microstructure and FTIR spectra reveal that AP/aerogels are formed by a physical combination of the silica aerogel matrix and aramid pulps. When the heat treatment temperature increases, the density slightly decreases and then increases to the maximum due to the significant volume shrinkage. The pyrolysis of aramid pulp and the collapse of silica skeletons occur during heat treatment; nevertheless, the typical structures of AP/aerogels do not change significantly. It is also found that both the hydrophobicity and the thermal insulation decrease with the increasing heat treatment temperature. We note that when the heat treatment is at 600 °C, the AP/aerogel still maintains a low density of 0.19 g/cm3 and a contact angle of 138.5°. The thermal conductivity is as low as 26.11 mW/m/K, measured using the transient hot wire method. Furthermore, the heat-treated AP/aerogels can avoid heat shock and possible thermal hazards during practical thermal insulation applications. The onset temperatures of the thermal decomposition of AP/aerogels increase from 298.8 °C for an untreated one to 414.7 °C for one treated at 600 °C, indicating that the thermal stability of AP/aerogels is improved significantly. This work provides a practical engineering approach to expand the thermal insulation applications of silica aerogel composites.

Siloxane Molecules: Nonlinear Elastic Behavior and Fracture Characteristics

Fracture phenomena in soft materials span multiple length and time scales. This poses a major challenge in computational modeling and predictive materials design. To pass quantitatively from molecular to continuum scales, a precise representation of the material response at the molecular level is vital. Here, we derive the nonlinear elastic response and fracture characteristics of individual siloxane molecules using molecular dynamics (MD) studies. For short chains, we find deviations from classical scalings for both the effective stiffness and mean chain rupture times. A simple model of a nonuniform chain of Kuhn segments captures the observed effect and agrees well with MD data. We find that the dominating fracture mechanism depends on the applied force scale in a nonmonotonic fashion. This analysis suggests that common polydimethylsiloxane (PDMS) networks fail at cross-linking points. Our results can be readily lumped into coarse-grained models. Although focusing on PDMS as a model system, our study presents a general procedure to pass beyond the window of accessible rupture times in MD studies employing mean first passage time theory, which can be exploited for arbitrary molecular systems.

Covalent Bonding and Atomic-Level Plasticity Increase Adhesion in Silicon-Diamond Nanocontacts

Nanoindentation and sliding experiments using single-crystal silicon atomic force microscope probes in contact with diamond substrates in vacuum were carried out in situ with a transmission electron microscope (TEM). After sliding, the experimentally measured works of adhesion were significantly larger than values estimated for pure van der Waals (vdW) interactions. Furthermore, the works of adhesion increased with both the normal stress and speed during the sliding, indicating that applied stress played a central role in the reactivity of the interface. Complementary molecular dynamics (MD) simulations were used to lend insight into the atomic-level processes that occur during these experiments. Simulations using crystalline silicon tips with varying degrees of roughness and diamond substrates with different amounts of hydrogen termination demonstrated two relevant phenomena. First, covalent bonds formed across the interface, where the number of bonds formed was affected by the hydrogen termination of the substrate, the tip roughness, the applied stress, and the stochastic nature of bond formation. Second, for initially rough tips, the sliding motion and the associated application of shear stress produced an increase in irreversible atomic-scale plasticity that tended to smoothen the tips' surfaces, which resulted in a concomitant increase in adhesion. In contrast, for initially smooth tips, sliding roughened some of these tips. In the limit of low applied stress, the experimentally determined works of adhesion match the intrinsic (van der Waals) work of adhesion for an atomically smooth silicon-diamond interface obtained from MD simulations. The results provide mechanistic interpretations of sliding-induced changes and interfacial adhesion and may help inform applications involving adhesive interfaces that are subject to applied shear forces and displacements.

High-Density Modification of H-Terminated Si(111) Surfaces Using Short-Chain Alkynes

H-Si(111)-terminated surfaces were alkenylated via two routes: through a novel one-step gas-phase hydrosilylation reaction with short alkynes (C₃ to C₆) and for comparison via a two-step chlorination and Grignard alkenylation process. All modified surfaces were characterized by static water contact angles and X-ray photoelectron spectroscopy (XPS). Propenyl- and butenyl-coated Si(111) surfaces display a significantly higher packing density than conventional C₁₀-C₁₈ alkyne-derived monolayers, showing the potential of this approach. In addition, propyne chemisorption proceeds via either of two approaches: the standard hydrosilylation at the terminal carbon (lin) at temperatures above 90 °C and an unprecedented reaction at the second carbon (iso) at temperatures below 90 °C. Molecular modeling revealed that the packing energy of a monolayer bonded at the second carbon is significantly more favorable, which drives iso-attachment, with a dense packing of surface-bound iso-propenyl chains at 40% surface coverage, in line with the experiments at <90 °C. The highest density monolayers are obtained at 130 °C and show a linear attachment of 1-propenyl chains with 92% surface coverage.