Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “lattice structure”

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 613 records · Page 34

Reactor physics characterization of triply periodic minimal surface-based nuclear fuel lattices

Triply periodic minimal surface (TPMS) lattices are receiving substantial attention in numerous engineering fields due to their impressive topology-driven physical characteristics. TPMS lattices are periodic structures of two distinct intertwined volume domains separated by an area-minimizing surface or wall. TPMS lattices have been observed in nature, such as biological membranes, skeletons, block copolymers, sea urchins, butterfly wings, and equipotential surfaces in crystals. Intriguingly, the topology of TPMS lattices can be easily parametrized via level-set equations and thus are heavily numerically and experimentally studied. Here, a significant research effort is currently applying TPMS lattices for heat exchangers and sinks. This paper extends TPMS lattice applications to nuclear reactor fuel designs, with a focus on identifying relevant TPMS geometric parameters controlling neutronics characteristics, such as reactivity, neutron spectrum, and heat removal properties. We found that fuel surface-area-to-volume ratios for TPMS lattices can be two orders of magnitude larger than current cylindrical fuel rods. Further, the selected TPMS lattice and its implicit equation, the unit cell pitch, wall thickness, and structure porosity are design parameters enabling neutronics optimization for both thermal and fast spectrum configurations, paving the way for exceptionally compact and dense nuclear core concepts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Contrasting structural reversibility and magnetic correlations in isostructural honeycomb magnets CrCl3 and 𝛼−RuCl3

We report a comparative neutron single crystal diffraction study of the structural and magnetic properties of layered halides CrCl3 and 𝛼−RuCl3. They host a honeycomb arrangement of transition metal ions with distinct electronic configurations and undergo a first-order structural transition between high-temperature 𝐶⁢2/𝑚 and low-temperature 𝑅⁢‾‾‾3. Both compounds show a step-like change in the 𝑐-lattice parameter across the structure transition. In contrast, the in-plane lattice response is quite different: 𝛼−RuCl3 exhibits an abrupt hysteretic change across the transition accompanied by progressive crystalline degradation upon thermal cycling, whereas CrCl3 shows a smooth in-plane lattice evolution and remains structurally robust. Magnetically, CrCl3 orders into an A-type antiferromagnetic structure at 𝑇𝑁=14 K and exhibits pronounced diffuse magnetic scattering extending up to about 40 K. 𝛼−RuCl3 shows no observable magnetic diffuse scattering above its zigzag antiferromagnetic ordering temperature 𝑇𝑁=7.6 K. These results suggest that the contrasting structural responses arise from an interplay between interlayer sliding energetics, stacking-strain coupling, and elastic accommodation of the stacking transition. The distinct chemical bonding and electronic configurations of the two compounds provide a microscopic basis for their different lattice responses to the structure transition and magnetic correlations.

Morgan, Zachary [ORNL] (ORCID:0000000243625911)↗

Transient X-ray Diffraction Reveals Nonequilibrium Phase Transition in Thin Films of CH 3 NH 3 PbI 3 Perovskite

Advantageous optoelectronic properties of methylammonium lead triiodide likely arise from coupling between photogenerated charge-carriers and the soft, deformable lattice. We investigate structural dynamics of MAPbI 3 films using time-resolved X-ray diffraction versus pump-probe time delay and pump intensity. During the first nanosecond, the lattice anisotropically distorts from tetragonal to cubic at excitation intensities that are insufficient to thermally induce the first-order thermodynamic phase transition at 330K. Further, the high-symmetry structure then relaxes back to the starting phase with 11 and 236-ns time constants via a different transition pathway than observed either in the first nanosecond or in previous reports for MAPbI 3 . Early-time dynamics are consistent with polaron formation and lattice strain stabilization while the slower recovery dynamics outlive radiative recombination and relates metastability. Fluence-independence of these lattice deformations in the low-power regime conveys relevance to optoelectronics including photovoltaics and highlights sustained involvement of non-equilibrium, photoinduced lattice reorganization in MAPbI 3 under functional conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure transition and zigzag magnetic order in Ir/Rh-substituted honeycomb lattice α - RuCl 3

Here, we report magnetization and neutron diffraction studies on crystal and magnetic structures of Ir- and Rh-substituted honeycomb lattice α-RuCl 3 . The iridium or rhodium atoms are distributed at the Ru site with little structural modification. Both systems undergo a room-temperature monoclinic C2/m to low-temperature trigonal $R\bar{3}$ phase transformation with a large recoverable hysteresis. At low temperature, a zigzag spin order is observed with the same characteristic wave vector (0,0.5,1) as in the parent α-RuCl 3 . Detailed magnetic structure refinement reveals an ordered moment of 0.32(5)µB/Ru and an upper boundary of canting angle of 15(4)º away from the basal plane at 5 K for the 10% Ir-substituted α-RuCl 3 , which is different from the 0.45–0.73 µB/Ru and 32°–48° canting angle reported in the parent compound α-RuCl 3 . The observation of unchanged RuCl 6 local octahedral environment, reduced ordered magnetic moment size and canting angle compared to previously reported highlights the potential to study quantum spin-liquid behavior through nonmagnetic ion doping.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 8. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 4. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 6. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 7. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 8. OUTCARL: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 10. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216.

36 MATERIALS SCIENCE↗

Charge-induced atomic strain as a predictor of structural phase transformation in rare-earth intermetallics

We present a descriptor based on charge-induced atomic strain in crystalline lattices for predicting structural phase transformations in rare-earth intermetallic compounds containing lanthanides and transition metals. The charge-induced local atomic strain was obtained from structural optimization of experimentally known crystalline phases using state of the art density-functional theory methods. The predictive power of the descriptor was evaluated on 𝑅⁢𝐸 2 ⁢In (𝑅𝐸 = rare earth) compounds, a class known for diverse phase transformations. We show that incorporating quantum-mechanical effects—such as local charge distribution, bonding, symmetry, and electronic structure—enhances the robustness of the descriptor. To gain further insight, we analyzed phononic and electronic behavior in Y 2 ⁢In and demonstrated that experimental phase transformations are captured only when atomic strain effects are included. The descriptor was further used to predict structural phase changes in (Y⁢b 1–𝑥 ⁢E⁢r 𝑥 ) 2 ⁢In and G⁡d 2 ⁡(I⁢n 1–𝑥⁢ A⁢l 𝑥 ), with predictions confirmed by x-ray powder diffraction. Here, while the current study is focused on lanthanide-based intermetallics, the underlying principles of the descriptor suggest potential applicability to other closely related classes of rare-earth intermetallics.

Density functional theory↗

Thermodynamics of Plutonium Monocarbide from Anharmonic and Relativistic Theory

Thermodynamics of plutonium monocarbide is studied from first-principles theory that includes relativistic electronic structure and anharmonic lattice vibrations. Density-functional theory (DFT) is expanded to include orbital-orbital coupling in addition to the relativistic spin-orbit interaction for the electronic structure and it is combined with anharmonic, temperature dependent, lattice dynamics derived from the self-consistent ab initio lattice dynamics (SCAILD) method. The obtained thermodynamics are compared to results from simpler quasi-harmonic theory and experimental data. Formation enthalpy, specific heat, and Gibbs energy calculated from the anharmonic model are validated by direct comparison with a calculation of phase diagram (CALPHAD) assessment of PuC and sub-stochiometric PuC0.896. Overall, the theory reproduces CALPHAD results and measured data for PuC rather well, but the comparison is hampered by the sub-stoichiometric nature of plutonium monocarbide. It was found that a bare theoretical approach that ignores spin-orbit and orbital-orbital coupling (orbital polarization) of the plutonium 5f electrons promotes too soft phonons and Gibbs energies that are incompatible with that of the CALPHAD assessment of the experimental data. The investigation of PuC suggests that the electronic structure is accurately described by plutonium 5f electrons as “band like” and delocalized, but correlate through spin polarization, orbital polarization, and spin-orbit coupling, in analogy to previous findings for plutonium metal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Structure of the Pion from First Principles Calculation in Lattice QCD

I describe the recent progress by the HadStruct Collaboration aimed at computing the structure of the pion from first-principles calculation in lattice QCD. I begin by presenting the Pseudo-PDF and Lattice Cross Section approaches, emphasizing that they are formulated in terms of Lorentz invariants with a short-distance scale to admit factorization. I present resent results on the pion PDF using both approaches, and compare with phenomenological determinations and other theoretical approaches. I conclude with plans for future work.

Richards, David↗

Nanoscale Phase Separation and Large Refrigerant Capacity in Magnetocaloric Material LaFe 11.5 Si 1.5

Here, ferromagnetic transitions, the magnetocaloric effect, and the atomic-scale microstructure in an annealed bulk material of LaFe 11.5 Si 1.5 were investigated using magnetic measurements, Mössbauer spectroscopy, high-resolution X-ray diffraction analysis, and high-resolution transmission electron microscopy. The results provide evidence of the coexistence of two ferromagnetic phases with different Curie temperatures. The phase with a Curie temperature of 216 K corresponds to a Si-rich phase of a NaZn 13 -type structure with a small lattice volume, whereas the phase with a Curie temperature of 185 K corresponds to a Si-poor phase of a NaZn 13 -type structure with a large lattice volume. This phase coexistence is observed on a nanometer scale and it can account for inverse thermal hysteresis of the ferromagnetic transition in the Si-rich phase during warming. Furthermore, the phase coexistence helps the bulk material achieve a refrigeration capacity of 170 J kg -1 due to a combination of the magnetocaloric effect in each phase. It is determined that the magnetocaloric effect in the Si-rich and the Si-poor phase brings about a maximum of the total entropy change of -3.55 and -6.80 and J kg -1 K -1 at respective Curie temperatures, respectively. Fundamentally, the phase coexistence indicates coloring of Si atoms in the lattice of a NaZn 13 -type structure and can be attributed to nanoscale phase separation of a disordered precursor during annealing of the bulk material. The nanoscale phase separation may occur in other magnetocaloric materials, thus allowing for improvement of their refrigeration capacity.

36 MATERIALS SCIENCE↗

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗

Solid face sheets enable lattice metamaterials to withstand high-amplitude impulsive loading without yielding

Owing to their ability to provide tunable mechanical responses, lattice materials are frequently studied to elucidate their response to static and dynamic loads. However, these roles are typically in opposition: static loads must be supported sufficiently far away from the onset of buckling or yielding, whereas dynamic loads are typically ameliorated by crushing of the lattice, which provides excellent energy-absorption due to the large plastic deformation accompanying densification. In contrast, this work considers the octet truss as an exemplar topology, in a structural role where it must simultaneously support static loads while enduring high-amplitude impulsive loads. This study focuses on the ability to withstand impulsive loads without yielding, an essential prerequisite to enduring dual loading. Computational studies using the ALE3D hydrocode were performed to examine the response of the octet truss under a short temporal width impulse shape associated with laser-driven shocks. A key finding was that covering the lattice with a solid face sheet and treating this face sheet thickness as a design variable allows the Taylor-like pulse to be attenuated prior to entering the weaker lattice, at the cost of added mass up front. Experimental validation was accomplished by laser-driven shock testing, using octet trusses printed out of Ti-5Al-5V-5Mo-3Cr. The results show that for a given quantity of mass, the attenuation is maximized when as much mass as possible is moved into the face sheet, leaving a more slender lattice structure. The effect of placing mass in the face sheet rather than lattice beams dominates the effect of relative density, to the point where a low-mass structure with most of the mass concentrated in the face sheet can outperform a high-mass structure with most of the mass in the lattice. Finally, by further understanding the propagation of short pulse width waves within under-dense structures, this study expand the domain of applicability of such structures, including lattice materials, to challenging dual-loading regimes spanning decades of strain rates.

36 MATERIALS SCIENCE↗

Vibrational spectroscopy of uranium tetrafluoride hydrates

Uranium tetrafluoride (UF 4 ) is an important intermediate in the production of UF6 and nuclear fuel. Historical characterization of UF 4 with Raman spectroscopy was plagued with ambiguity until the first accurate Raman spectrum was published by our group in 2016. Although generally considered to be relatively stable, UF 4 can hydrolyze to form numerous UF 4 hydrates that may play a role in future uranium waste forms. In contrast to anhydrous UF 4 , the hydrates, with their OH stretch and HOH bending modes, can be spectroscopically characterized by the type and degree of water bonding in the crystal lattice, which can yield additional information about their crystal structure. Herein, vibrational spectroscopy (Raman and infrared) was used to characterize three different UF 4 hydrates: UF 4 (H 2 O) 0.33 , U 3 F 12 (H 2 O), and UF 4 (H 2 O) 2.5 . Spectra show the different hydrates vary in the number of observed bands, full-width half-maximum of the bands, and band intensity. These differences are due to varying interactions between the OH stretch and HOH bending modes with UF 4 and the polymeric UF 4 structure in the crystal lattice. These vibrational data, in combination with spectral fitting and crystallographic structures measured with powder X-ray diffraction and single crystal X-ray diffraction, provide unique details on the location of water molecules in the crystal lattice of hydrated UF 4 , and provide an interesting contrast to the vibrational spectra of anhydrous UF 4 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Moving beyond the constraints of chemistry via crystal structure discovery with isotropic multiwell pair potentials

The rigid constraints of chemistry-dictated by quantum mechanics and the discrete nature of the atom-limit the set of observable atomic crystal structures. What structures are possible in the absence of these constraints? Here, we systematically crystallize one-component systems of particles interacting with isotropic multiwell pair potentials. We investigate two tunable families of pairwise interaction potentials. Our simulations self-assemble a multitude of crystal structures ranging from basic lattices to complex networks. Sixteen of the structures have natural analogs spanning all coordination numbers found in inorganic chemistry. Fifteen more are hitherto unknown and occupy the space between covalent and metallic coordination environments. The discovered crystal structures constitute targets for self-assembly and expand our understanding of what a crystal structure can look like.

crystal structures↗

Ellipsometric study of InGaAs MODFET material

In(x)Ga(1-x)As based MODFET (modulation doped field effect transistor) material was grown by molecular beam epitaxy on semi-insulating InP substrates. Several structures were made, including lattice matched and strained layer InGaAs. All structures also included several layers of In(0.52)Al(0.48)As. Variable angle spectroscopic ellipsometry was used to characterize the structures. The experimental data, together with the calibration function for the constituent materials, were analyzed to yield the thickness of all the layers of the MODFET structure. Results of the ellipsometrically determined thicknesses compare very well with the reflection high energy electron diffraction in situ thickness measurements.

Alterovitz, S. A.↗

Evolution of the Fermi Surface of 1T-VSe 2 across a Structural Phase Transition

Periodic lattice distortion, known as the charge density wave, is generally attributed to electron–phonon coupling. This correlation is expected to induce a pseudogap at the Fermi level in order to gain the required energy for stable lattice distortion. The transition metal dichalcogenide 1T-VSe 2 also undergoes such a transition at 110 K. Here, we present detailed angle-resolved photoemission spectroscopy experiments to investigate the electronic structure in 1T-VSe 2 across the structural transition. Previously reported warping of the electronic structure and the energy shift of a secondary peak near the Fermi level as the origin of the charge density wave phase are shown to be temperature independent and hence cannot be attributed to the structural transition. Our work reveals new states that were not resolved in previous studies. Earlier results can be explained by the different dispersion natures of these states and temperature-induced broadening. Only the overall size of the Fermi surface is found to change across the structural transition. These observations, quite different from the charge density wave scenario commonly considered for 1T-VSe 2 and other transition metal dichalcogenides, bring fresh perspectives toward correctly describing structural transitions. Therefore, these new results can be applied to material families in which the origin of the structural transition has not been resolved.

36 MATERIALS SCIENCE↗

Shadowed Space Heating of Sparse Structures

Complete heat-flux and temperature maps computed. Computer program SSQ developed to address and quantify complex, solar-shadowing conditions inherent in sparse, lattice-type space structures. Analysis procedure one of assessing partial shadowing of structural elements by multiple, similar slender members. Program yields schedules of incident-solar, Earth-thermal, and Earth-albedo radiation throughout complete orbit for elemental locations on selected structural members. Thermal response computed in optional routine. Complete heat-flux and temperature mapping obtained by repeated computations for selected elements of interest.

Oneill, J. L.↗

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↗