Search NASA⌕ Search

SEARCH · Search NASA

Results for “Strain tuning”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pit rim decomposition into multiple quantum dots on surfaces of epitaxial thin films grown on pit-patterned substrates

Here, we report results of dynamical simulations according to an experimentally validated surface morphological evolution model on the formation of regular arrays of quantum dot molecules (QDMs) consisting of 1D arrays of smaller interacting quantum dots (QDs). These QD arrays form along the sides of each pit rim on the surface of a coherently strained thin film epitaxially deposited on a semiconductor substrate, the surface of which consists of a periodic pattern of inverted prismatic pits with rectangular pit openings. We find that this complex QDM pattern results from the decomposition of the pit rim from a “quantum fortress” with four elongated QDs into four 1D arrays of multiple smaller QDs arranged along each side of the pit rim. Systematic parametric analysis indicates that varying the pit opening dimensions and the pit wall inclination directly impacts the number of QDs in the resulting QDM pattern, while varying the pit depth only affects the dimensions of the QDs in the nanostructure pattern. Therefore, the number, arrangement, and sizes of QDs in the resulting pattern of QDMs on the film surface can be engineered precisely by proper tuning of the pit design parameters. Our simulation results are supported by predictions of morphological stability analysis, which explains the pit rim decomposition into multiple QDs as the outcome of a tip-splitting instability and provides a fundamental characterization of the post-instability nanostructure pattern. Our theoretical findings can play a vital role in designing optimal semiconductor surface patterns toward enabling future nanofabrication technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Alloying multiple halide perovskites on the same sublattice in search of stability and target band gaps

Single-component halide perovskites (HPs) rarely satisfy all the necessary criteria for optoelectronic applications, such as achieving an optimal band gap while maintaining high chemical and structural stability. Alloying halide perovskites has emerged as a promising strategy, not only to enhance stability but also to fine-tune their electronic and optical properties. In this work, we explore multiple degrees of freedom in alloy design, considering different substitution sublattices sites (A, B, or X in ABX3 perovskites), various chemical species (isovalent and hetero-valent elements), and multi-component compositions on a given sublattice. Using first-principles calculations based on density functional theory (DFT), we investigate how compositional variations influence the electronic (band gap) and structural properties (mixing enthalpy) of HP alloys. Our approach employs the polymorphous cell model, allowing full local relaxation which breaks local symmetry while preserving global cubic symmetry—an essential framework for accurately modeling HPs. Our results reveal that X-site mixing (halogen substitution) primarily affects the valence band maximum, allowing target band gap engineering. Additionally, variations in halogen radii introduce internal strain through octahedral distortions, influencing the mixing enthalpy. A-site substitution, while not directly contributing to the band edge states, modifies structural stability via volume effects, indirectly impacting the band gap. B-site alloying plays a dominant role in band gap modulation, leading to either positive or negative band gap bowing. Specifically, isovalent B-site mixing (Sn–Pb) induces strong positive bowing, where the alloy band gap is smaller than the average gap of parent compounds, whereas hetero-valent mixing (Cd–Pb) results in pronounced negative bowing. As an aside, we investigate the competition between the excess energy of disordered alloys vs. that of long-range ordered double perovskites of the same compositions, seeking examples of ordered phases emerging from disordered alloys. Furthermore, our findings provide fundamental insights into the electronic and structural behavior of HP alloys, offering valuable design principles for the development of stable and efficient materials for next-generation photovoltaic and optoelectronic devices.

14 SOLAR ENERGY↗
DOE OSTI2022

Anomalous Hall effect and perpendicular magnetic anisotropy in ultrathin ferrimagnetic NiCo2O4 films

The inverse spinel ferrimagnetic NiCo2O4 possesses high magnetic Curie temperature TC, high spin polarization, and strain-tunable magnetic anisotropy. Understanding the thickness scaling limit of these intriguing magnetic properties in NiCo2O4 thin films is critical for their implementation in nanoscale spintronic applications. In this work, we report the unconventional magnetotransport properties of epitaxial (001) NiCo2O4 films on MgAl2O4 substrates in the ultrathin limit. Anomalous Hall effect measurements reveal strong perpendicular magnetic anisotropy for films down to 1.5 unit cell (1.2 nm), while TC for 3 unit cell and thicker films remains above 300 K. The sign change in the anomalous Hall conductivity (σxy) and its scaling relation with the longitudinal conductivity (σxx) can be attributed to the competing effects between impurity scattering and band intrinsic Berry curvature, with the latter vanishing upon the thickness driven metal–insulator transition. Our study reveals the critical role of film thickness in tuning the relative strength of charge correlation, Berry phase effect, spin–orbit interaction, and impurity scattering, providing important material information for designing scalable epitaxial magnetic tunnel junctions and sensing devices using NiCo2O4.

Chen, Xuegang↗
DOE OSTI2021

Developing Multi-Gene CRISPRa/I Programs to Accelerate DBTL Cycles in ABF Hosts Engineered for Chemical Production (CRADA 468)

Bacterial metabolism is comprised of large and complex gene networks that can produce valuable chemical products. Sophisticated organism engineering efforts are required to optimize production of high-value compounds from these networks. In principle, synthetic multi-gene transcriptional programs could be constructed to reengineer these networks for efficient industrial chemical production. In practice, however, our incomplete ability to understand and model the underlying networks, combined with our limited ability to predictably control the expression of multiple genes makes achieving this goal difficult. To overcome these challenges, we will combine new CRISPR-Cas multi-gene expression programs with computational modeling, machine learning, and multi-omics data to enhance the efficacy of design-build-test-learn (DBTL) cycles. For industrially promising microorganisms in early stages of development, creating technologies for rapidly engineering complex multi-gene programs could be transformative for accelerating data- and model-driven strain design. New CRISPR-Cas tools allow programmable gene activation (CRISPRa) or repression (CRISPRi) at multiple genes simultaneously, using the catalytically inactive Cas9 protein (dCas9) with guide RNAs that recognize DNA targets through predictable Watson-Crick base pairing. To enable accelerated DBTL cycles, we will combine these technologies with advanced Agile BioFoundry (ABF) capabilities for multi-omics data collection and machine learning. We will demonstrate the immediate applicability of these tools by rapidly improving the production of an industrial aromatic in multiple ABF organisms. We recently identified and optimized new transcriptional activators that can be linked to programmable CRISPR-Cas DNA binding domains to activate gene expression in E. coli. We can now use these CRISPRa tools as generalizable trans-acting regulators for combinatorial multi-gene expression tuning that can be easily transferred to new pathways and networks without additional genome engineering. We anticipate these tools will also transfer to new hosts. We have recently found that CRISPRa systems developed in E. coli can be readily ported to Pseudomonas putida, suggesting that multi-gene CRISPRa/i programs for diverse ABF organisms may be within reach.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Uncertainty quantification for equations of state: copper as an example

Equations of state are essential for providing a fundamental description of materials properties in thermodynamic equilibrium and are used to provide closure relations for hydrodynamics simulations. Generally, equations of state rely on simple physics-based parameterized materials models to inform on the free energy of a material through out a given thermodynamic state space. Historically the parameters of these models have been tuned by hand to fit various experimental data. However, modern optimization and uncertainty quantification techniques allow us to quickly test thousands of parameter combinations and obtain meaningful uncertainty estimates on the parameters, opening opportunities for assessing systematic uncertainties in experiments, assessing model adequacy, and more. In this report, we use Bayesian inference to fit the solid (fcc) equation of state of copper. We focus on fitting five different experimental datasets, including the isobaric density, isobaric heat capacity, room temperature isotherm, principal isentrope, and principal Hugoniot. We fit all five data types simultaneously, and then explore the extent to which combinations of 2 subsets of the 5 datasets can constrain the EOS parameters, as compared to the fit to all 5. This information is useful for investigating the extent to which different datasets can con strain EOS models and thereby help guide experimental investigations in order to best constrain the EOS. We also discuss ways that the methodologies can be used to investigate systematic discrepancies between experiments, as well as how the methods can be used to assess model uncertainty. The framework we develop is general, in that it can be used with a variety of optimization or uncertainty quantification techniques and with a variety of data sources, including both experimental and ab-inito data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fluorogenic Biosensing with Tunable Polydiacetylene Vesicles

Polydiacetylenes (PDAs) are conjugated polymers that are well known for their colorimetric transition from blue to red with the application of energetic stimulus. Sensing platforms based on polymerized diacetylene surfactant vesicles and other structures have been widely demonstrated for various colorimetric biosensing applications. Although less studied and utilized, the transition also results in a change from a non-fluorescent to a highly fluorescent state, making polydiacetylenes useful for both colorimetric and fluorogenic sensing applications. Here, we focus on the characterization and optimization of polydiacetylene vesicles to tune their sensitivity for fluorogenic sensing applications. Particularly, we look at how the structure of the diacetylene (DA) hydrocarbon tail and headgroup affect the self-assembled vesicle size and stability, polymerization kinetics, and the fluorogenic, blue to red phase transition. Longer DA acyl tails generally resulted in smaller and more stable vesicles. The polymerization kinetics and the blue to red transition were a function of both the DA acyl tail length and structure of the headgroup. Decreasing the acyl tail length generally led to vesicles that were more sensitive to energetic stimuli. Headgroup modifications had different effects depending on the structure of the headgroup. Ethanolamine headgroups resulted in vesicles with potentially increased stimuli responsivity. The lower energy stimulus to induce the chromatic transition was attributed to an increase in headgroup hydrogen bonding and polymer backbone strain. Boronic-acid headgroup functionalization led to vesicles that were generally unstable, only weakly polymerized, and unable to fully transform to the red phase due to strong polar, aromatic headgroup interactions. This work presents the design of PDA vesicles in the context of biosensing platforms and includes a discussion of the past, present, and future of PDA biosensing.

Chemistry↗
DOE OSTI2021

Interface-driven magnetic anisotropy in relaxed La 0.7 Sr 0.3 CrO 3 /La 0.7 Sr 0.3 MnO 3 heterostructures on MgO

We investigate the structural and magnetic properties of La 0.7 Sr 0.3 CrO 3 /La 0.7 Sr 0.3 MnO 3 (LSMO) heterostructures grown on (001)-oriented MgO by molecular beam epitaxy. Due to the large film-substrate lattice mismatch, strain relaxation is found to occur within the first 2–3 unit cells (uc) of the film as evidenced by reflection high energy electron diffraction and high-resolution synchrotron x-ray reciprocal space mapping. Furthermore, we find that the presence of the LSCO spacer and capping layers leads to ferromagnetism in ultra-thin LSMO layers with thicknesses of the order of 2 uc with the magnetic easy axis oriented in the film plane. Net magnetic moments of 1.4 and 2.4 μ B /Mn are measured for (2 uc LSCO/2 uc LSMO) and (2 uc LSCO/4 uc LSMO) superlattices, respectively, by superconducting quantum interference device magnetometry. The effective magnetic anisotropy of the relaxed (2 uc LSCO/4 uc LSMO) heterostructure is found to be an order of magnitude higher than bulk LSMO highlighting the critical role of interfacial interactions in tuning magnetic anisotropy at complex oxide interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Characterization of a Mo-Nitrogenase Variant Containing a Citrate-Substituted Cofactor

Nitrogenase converts N 2 to NH 3 , and CO to hydrocarbons, at its cofactor site. In this work, we report a biochemical and spectroscopic characterization of a Mo-nitrogenase variant expressed in an Azotobacter vinelandii strain containing a deletion of nifV, the gene encoding the homocitrate synthase. Designated NifDK Cit , the catalytic component of this Mo-nitrogenase variant contains a citrate-substituted cofactor analogue. Activity analysis of NifDK Cit reveals a shift of CO reduction from H 2 evolution toward hydrocarbon formation and an opposite shift of N 2 reduction from NH 3 formation toward H 2 evolution. Consistent with a shift in the Mo K-edge energy of NifDK Cit relative to that of its wild-type counterpart, EPR analysis demonstrates a broadening of the line-shape and a decrease in the intensity of the cofactor-originated S=3/2 signal, suggesting a change in the spin properties of the cofactor upon citrate substitution. These observations point to a crucial role of homocitrate in substrate reduction by nitrogenase and the possibility to tune product profiles of nitrogenase reactions via organic ligand substitution.

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

Multi-state Catalysts Modulated by Mechanical Force (Final Report)

The development of more efficient catalytic processes and new approaches to control catalytic activity and selectivity are central to the realization of more selective, atom economic, and energy efficient routes to value added chemicals and polymers. The reactivity and selectivity of a transition metal catalyst is intimately related to the ligand-sphere geometry and, in many cases, the ideal ligand geometry for one step of a catalytic cycle is poorly matched to the ideal ligand geometry for another, resulting in sub-optimal efficiency. Macroscopic mechanical forces are both large, potentially much larger than interatomic forces, and are directional and localized to an extent that differentiates them from other forms of energy input such as heat or light. As such, mechanical force represents a heretofore untapped approach to modulate catalyst geometry, with the potential to reversibly modulate catalyst geometry on the timescale of catalytic turnover or monomer enchainment. This project has addressed the fundamental challenges in material-to-molecule strain coupling associated with the development of a new class of mechanically responsive catalysts (mechanocatalysts) in which active organotransition metal catalysts are strategically embedded in a flexible polymer network such that application of external mechanical force (stretching or deformation) leads to modulation of catalyst geometry, and hence reactivity and selectivity. Our efforts during the tenure of this grant were directed toward the elucidation of force-reactivity relationships of elementary transformations that occur within the first coordination sphere of a transition metal complex employing stiff stilbene photoswitches tethered to a flexible bidentate phosphine ligand derived from MeOBiphep as molecular force probes which provide a range of compressive and extension forces to the coupled transition metal complex depending on the geometry of the stiff stilbene and length of the tethering chains. During the tenure of this grant, we have quantified the rate of C(sp 2 )-C(sp 2 ) reductive elimination from platinum(II) diaryl complexes containing bis(phosphine) force probe ligands as a function of mechanical force; compressive forces decreased the rate of reductive elimination whereas extension forces increased the rate relative to the strain-free MeOBiphep complex with a 3.4-fold change in rate over a ~290 pN range of restoring forces. In a similar manner, we have quantified the rate of oxidative addition of bromobenzene to low-ligated palladium(0) complexes containing force probe ligands as a function of mechanical force; compressive forces increase the rate of oxidative addition, whereas tensile forces decrease the rate with a ~6 fold change in rate across ~340 pN of force applied to the complexes. In both cases, experimental and computational analyses argue strongly against any significant force-induced perturbation of ground state geometry within the first coordination sphere of the reactant complexes. Rather, the force/rate behavior observed for these transformations across these ranges of forces is attributed to the coupling of force to the nuclear motion comprising the reaction coordinates for reductive elimination and oxidative addition. These results together inform the development of catalysts whose activity can be tuned by an external force that is adjusted within a catalytic cycle and suggest opportunities to experimentally map geometry changes associated with reactions in transition metal complexes and potential strategies for force-modulated catalysis.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Uniaxial stress effect on the electronic structure of quantum materials

Uniaxial stress has proven to be a powerful experimental tuning parameter for effectively controlling lattice, charge, orbital, and spin degrees of freedom in quantum materials. In addition, its ability to manipulate the symmetry of materials has garnered significant attention. Recent technical progress to combine uniaxial stress cells with quantum oscillation and angle-resolved photoemission techniques allowed to study the electronic structure as function of uniaxial stress. This review provides an overview on experimental advancements in methods and examines studies on diverse quantum materials, encompassing the semimetal WTe 2 , the unconventional superconductor Sr 2 RuO 4 , Fe-based superconductors, and topological materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Engineering the fracture resistance of 2H-transition metal dichalcogenides using vacancies: An in-silico investigation based on HRTEM images

Vacancy engineering of 2H-transition metal dichalcogenides (2H-TMDs) has recently attracted great attention due to its potential to fine-tune the phonon and opto-electric properties of these materials. From a mechanical perspective, this symmetry-breaking process typically reduces the overall crack resistance of the material and adversely affects its reliability. However, vacancies can trigger the formation of heterogeneous phases that synergistically improve fracture properties. In this study, using MoSe2 as an example, we characterize the types and density of vacancies that can emerge under electron irradiation and quantify their effect on fracture. Molecular dynamic (MD) simulations, employing a re-parameterized Tersoff potential capable of accurately capturing bond dissociation and structural phase changes, reveal that isolated transition metal monovacancies or chalcogenide divacancies tend to arrest the crack tip and hence enhance the monolayer toughness. In contrast, isolated chalcogenide monovacancies do not significantly affect toughness. The investigation further reveals that selenium vacancy lines, formed by high electron dose rates, alter the crack propagating direction and lead to multiple crack kinking. Using atomic displacements and virial stresses together with a continuum mapping, displacement, strain, and stress fields are computed to extract mechanistic information, e.g., conditions for crack kinking and size effects in fracture events. The study also reveals the potential of specific defect patterns, “vacancy engineering,” to improve the toughness of 2H-TMDs materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
DOE OSTI2024

Thermally generated magnonic spin currents in a polycrystalline gadolinium iron garnet thin film with perpendicular magnetic anisotropy

Rare-earth iron garnets (REIGs) are the benchmark systems for magnonics, including the longitudinal spin Seebeck effect (LSSE). While most research has focused on single-crystalline REIGs on complimentary garnet substrates, moving to more, cost-effective complementary metal-oxide semiconductor (CMOS)-compatible substrates is important to integrate REIG thin films with existing technology. In this regard, we grow a 130 nm-thick polycrystalline gadolinium iron garnet (GdIG) film on the Si/SiO2 substrate and investigate the temperature-dependent LSSE. Interestingly, the polycrystalline GdIG film exhibits perpendicular magnetic anisotropy (PMA) at room temperature which is induced by tensile in-plane (IP)-strain originating from the thermal-expansion mismatch between the GdIG film and the substrate during rapid thermal annealing. Further, a spin-reorientation transition from the out-of-plane IP direction below TS = 180 K is observed. Additionally, the film reveals a magnetic compensation temperature, TComp, of ≈240 K. The LSSE voltage not only demonstrates a sign-inversion around TComp, but also shows noticeable changes around TS. As compared to a single-crystalline GdIG film, the lower LSSE voltage for the polycrystalline GdIG is attributed to the higher effective magnetic anisotropy and enhanced magnon scattering at the grain boundaries. Our study not only paves the way for the cost-effective growth of CMOS-compatible REIG-based systems with PMA for magnonic memory and information processing applications, but also highlights the fact that the spincaloritronic and spin-insulatronic properties of the polycrystalline REIGs follow those of their single-crystalline counterparts with reduced spin-to-charge conversion efficiency through LSSE which can be tuned further by controlling the average gran size and interface engineering.

Physics↗
DOE OSTI2025

Versatile Liquid Crystal Elastomer Formulations Using Amine-Acrylate Chemistry and Processing for Advanced Manufacturing

This study aims to tune rheological properties of liquid crystal elastomers (LCEs) through a strategic approach, leveraging the differing reaction rates of thiols and primary and secondary amines with acrylate monomers or the inclusion of processing additives to enable amenability to advanced manufacturing. We found that varying the time of oligomerization at elevated temperatures of amine-functional monomers with acrylate-functional monomers enabled simple control of rheological properties. Additionally, we could also control rheological properties by introducing reinforcing fillers or organic solvents. By varying the oligomerization conditions, we achieved a broad range of viscosities, from 1.1 Pa s to 560 Pa s as measured at 60 °C. When the rheological properties were varied through the incorporation of various solid and liquid inclusions, viscosities could range from less than 1 Pa s to nearly 1000 Pa s when measured at 60 °C. The shape morphing capabilities that are intrinsic to LCEs and the mechanical properties were also characterized. Here, we find that the oligomerization conditions and the addition of fillers influence shape morphing and energy absorption, as characterized by the stress–strain characteristics, further suggesting their potential to be used in energy dissipation applications where performance can be turned. Last, we demonstrate how the range of accessible rheological properties grants formulations amenable to various advanced manufacturing techniques.

Liquid chromatography↗
DOE OSTI2025

Next-generation sequencing dataset of genome-scale CRISPRi in Synechococcus sp. PCC 7002 across seven conditions

A 33,298-member sgRNA library developed for Synechococus sp. PCC 7002 was screened with two replicates across seven growth conditions and sequenced with Illumina NextSeq (paired end, 2x150 bp) for a total of ~700M reads. The original plasmid library and the library after transformation into a dCas9-containing and dCas9-absent strain were also sequenced as a reference for initial sgRNA abundance.

genome- wide screens environmental acclimation spe↗
DOE OSTI2025

Tuning the structure and superconductivity of SrNi 2⁢ P 2 by Rh substitution

The compound SrNi 2⁢ P 2 is unique among the ThCr 2 ⁢Si 2 class since it exhibits a temperature-induced transition upon cooling from an uncollapsed tetragonal (ucT) state to a one-third-collapsed orthorhombic (tcO) state where one out of every three P-rows bond across the Sr layers. This compound is also known for exhibiting bulk superconductivity below 1.4 K at ambient pressure. Here, in this paper, we report on the effects of Rh substitution in Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 on the structural and superconducting properties. We studied the variation of the nearest P-P distances as a function of the Rh fraction at room temperature, as well as its temperature dependence for selected compositions. We find that increasing the Rh fraction leads to a decrease in the transition temperature between the ucT and tcO states, until a full suppression of the tcO state for x ≥ 0.166. The superconducting transition first remains nearly insensitive to the Rh fraction, and then it increases to 2.3 K after the tcO state is fully suppressed. These results are summarized in a phase diagram, built upon the characterization by energy dispersive x-ray spectroscopy, x-ray diffraction, resistance, magnetization, and specific heat measurements done on crystalline samples with varying Rh content. The relationship between band structure, crystal structure, and superconductivity is discussed based on previously reported band structure calculations on SrRh 2 ⁢P 2 . Moreover, the effect of Rh fraction on the stress-induced structural transitions is also addressed by means of strain-stress studies done by uniaxial compression of single-crystalline micropillars of Sr⁢(Ni 1-x ⁢Rh x ) 2 ⁢P 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Substrate-Directed Dimensional and Phase Control of Peptide Assemblies on Two-Dimensional van der Waals Materials

Understanding and controlling biomolecular self-assembly on van der Waals (vdW) materials has the potential to advance hybrid bioelectronic devices by enabling precise tuning of the interface and modulation of the resulting electronic properties of the biomolecule-vdW heterostructure. However, how surface properties of vdW materials direct biomolecule assembly remains poorly understood. To fill this knowledge gap, we investigated the assembly of a peptide known to assemble into two-dimensional (2D) crystalline films on MoS 2 on three representative vdW surfaces: WS 2 , MoS 2 , and highly oriented pyrolytic graphite (HOPG). Using in situ atomic force microscopy (AFM), we find that assembly is substrate-dependent, resulting in multilayers on WS 2 , monolayers on MoS 2 , and multiple coexisting phases on HOPG. WS 2 exhibits a higher negative charge, strong long-range electrostatic interactions, and extensive hydration layering that may promote multilayer stacking. In contrast, MoS 2 has stronger short-range interactions with the peptides but much weaker long-range interactions and hydration structure, which may favor monolayer formation. Molecular dynamics simulations predict a corresponding switch from monolayer to multilayer aggregates of the adsorbed monomers, reflected in their relative mobilities. On hydrophobic HOPG, the peptides bind most strongly and remain as monomers with high surface mobility. The peptide dimers comprising the basic unit of the crystals are more compact on HOPG, which has a smaller lattice constant than WS 2 or MoS 2 , suggesting strain contributes to stabilizing multiple phases. Our results provide mechanistic insights into how surface charge and hydration structure, and the lattice structure of the substrates governs peptide assembly on vdW materials, offering a framework to rationally control the 2D peptide-vdW heterostructures.

Molecular dynamics simulations↗
DOE OSTI2025

Enhanced Interfacial Bonding of Graft Copolymers

To understand how thermoplastic welding strength can be tuned through chemical modifications and macromolecular topology, we combined coarse-grained molecular dynamics (MD) simulations with experimental validation. Our simulations examined the diffusion dynamics of both linear and graft polymers across representative interfaces, revealing that diffusion-controlled interdigitation follows a power law, with the exponent decreasing from 0.34 to 0.11 as grafting density increases from 7.5 to 196% (with side chains grafted to both sides of a monomer unit). The addition of side chains enhances welding efficiency, as dense bottlebrush polymers with high grafting density reach maximum rupture strength faster than linear polymers. However, their saturated rupture strength is lower. This observation is subsequently corroborated by experimental lap-shear tests comparing linear polyethylene with octene grafted polyethylene elastomers. Our MD simulations show that unlike linear polymers, where backbone entanglements dominate, the grafted side chains introduce mechanisms in addition to entanglement dilution. The rapid interdigitation of side chains creates a dense mesh of entropic van der Waals contacts, which can also enhance the film welding. Furthermore, our MD simulations reveal a brittle rupture behavior in linear and comb-like (mildly grafted) polymers, while bottlebrush (densely grafted) polymers display elastomeric behavior with a pronounced stress plateau prior to fracture. Our simulations deconvolute the influence of polymer topology on deformation behavior. The rate of polymer deformation becomes lower than the applied strain rate prior to rupture, and the onset of this deviation is progressively delayed from linear to bottlebrush polymers. This trend highlights the critical role of molecular architecture in governing the mechanical response. In conclusion, these results provide deeper insight into the underlying welding mechanisms of topological polymers and present a potential approach for mitigating the interface anisotropy that is inherent in advanced manufacturing techniques such as fused filament fabrication.

graft copolymers↗