Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Material Response”

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 901 records · Page 50

Calcium is associated with specific soil organic carbon decomposition products

Abstract. Calcium (Ca) may contribute to the preservation of soil organic carbon (SOC) in more ecosystems than previously thought. Here, we provide evidence that Ca is co-located with SOC compounds that are enriched in aromatic and phenolic groups, across different acidic soil types and locations with different ecosystem properties, differing in terms of climate, parent material, soil type, and vegetation. In turn, this co-localised fraction of Ca–SOC is removed through cation exchange, and the association is then only re-established during decomposition in the presence of Ca (Ca addition incubation). Thus, this highlights a causative link between decomposition and the co-location of Ca with a characteristic fraction of SOC. Decomposition increases the relative proportion of negatively charged functional groups, which can increase the propensity for the association between SOC and Ca; in turn, this association can inhibit dissolved organic carbon export or further decomposition. We propose that this mechanism could be driven by Ca hotspots at the microscale shifting local decomposition processes and thereby explaining the co-location of Ca with SOC of a specific composition across different acidic soil environments. Incorporating this biogeochemical process into Earth system models could improve our understanding, predictions, and management of carbon dynamics in soils, as well as accounting for their response to Ca-rich amendments.

Rowley, Mike C. (ORCID:0000000224407855)

Dynamical structure factors of warm dense matter from time-dependent orbital-free and mixed-stochastic-deterministic density functional theory

Abstract We present the first calculations of the inelastic part of the dynamical structure factor (DSF) for warm dense matter (WDM) using time-dependent orbital-free density functional theory (TD-OF-DFT) and mixed-stochastic-deterministic (mixed) Kohn Sham TD-DFT (KS TD-DFT). WDM is an intermediate phase of matter found in planetary cores and laser-driven experiments, where the accurate calculation of the DSF is critical for interpreting x-ray Thomson scattering measurements. Traditional TD-DFT methods, while highly accurate, are computationally expensive, motivating the exploration of TD-OF-DFT and mixed TD-KS-DFT as more efficient alternatives. We applied these methods to experimentally measured WDM systems, including solid-density aluminum and beryllium, compressed beryllium, and carbon–hydrogen mixtures. Our results show that TD-OF-DFT requires a dynamical kinetic energy potential in order to qualitatively capture the plasmon response. Additionally, it struggles with capturing bound electron contributions. In contrast, mixed TD-KS-DFT offers greater accuracy in distinguishing bound and free electron effects, aligning well with experimental data, though at a higher computational cost. This study highlights the trade-offs between computational efficiency and accuracy, demonstrating that TD-OF-DFT remains a valuable tool for rapid scans of parameter space, while mixed TD-KS-DFT should be preferred for high-fidelity simulations. Our findings provide insight into the future development of DFT methods for WDM and suggest potential improvements for TD-OF-DFT.

36 MATERIALS SCIENCE

Automated scanning probe microscopy of combinatorial ferroelectric libraries: Gaussian-process-guided exploration and noise-aware experiment planning

Combinatorial materials libraries provide an efficient route for mapping composition–property relationships, but their broader impact depends on rapid, quantitative, and functionally relevant characterization. Scanning Probe Microscopy (SPM), including piezoresponse force microscopy (PFM), offers significant potential for quantitative, functionally relevant combi-library readouts. Here, we implement a fully automated SPM workflow for ferroelectric combinatorial libraries and benchmark Gaussian-process-based Bayesian optimization strategies for autonomous experiment planning. The workflow integrates automated probe motion, contact optimization, imaging, and dual amplitude resonance tracking-PFM spectroscopy, and uses scalarized spectroscopic observables to guide subsequent measurements. Stage motion, probe engagement, in-contact tuning, imaging, spectroscopy, and the choice of the next measurement location all proceed without human input. We demonstrate the approach on Sm-doped BiFeO 3 and Zn x Mg 1−x O libraries. By comparing vanilla Bayesian optimization with a measured-noise variant, we show that explicit treatment of local reproducibility can improve modeling of composition-dependent response when the measured variance is physically meaningful, but can also reduce robustness when variability is dominated by outliers or topographic artifacts. Furthermore, these results establish automated SPM as a bridge between combinatorial synthesis and quantitative functional characterization.

Liu, Yu [University of Tennessee, Knoxville, TN (U

Insulator‐to‐Metal Transition and Isotropic Gigantic Magnetoresistance in Layered Magnetic Semiconductors

Magnetotransport, the response of electrical conduction to external magnetic field, acts as an important tool to reveal fundamental concepts behind exotic phenomena and plays a key role in enabling spintronic applications. Magnetotransport is generally sensitive to magnetic field orientations. In contrast, efficient and isotropic modulation of electronic transport, which is useful in technology applications such as omnidirectional sensing, is rarely seen, especially for pristine crystals. Here a strategy is proposed to realize extremely strong modulation of electron conduction by magnetic field which is independent of field direction. GdPS, a layered antiferromagnetic semiconductor with resistivity anisotropies, supports a field-driven insulator-to-metal transition with a paradoxically isotropic gigantic negative magnetoresistance insensitive to magnetic field orientations. This isotropic magnetoresistance originates from the combined effects of a near-zero spin–orbit coupling of Gd 3+ -based half-filling ƒ-electron system and the strong on-site f – d exchange coupling in Gd atoms. These results not only provide a novel material system with extraordinary magnetotransport that offers a missing block for antiferromagnet-based ultrafast and efficient spintronic devices, but also demonstrate the key ingredients for designing magnetic materials with desired transport properties for advanced functionalities.

36 MATERIALS SCIENCE

Modeling gas migration through clay-based buffer material using coupled multiphase fluid flow and geomechanics with stress-dependent gas permeability

A model for gas migration through clay-based buffer material is developed for modeling gas generation and migration associated with deep geologic nuclear waste disposal. The model is based on a multiphase fluid flow and geomechanics simulator that is adapted to consider enhanced gas flow when gas pressure is high enough to approach the confining stress magnitude. A key feature in the model is a direct coupling between gas permeability and stress, through a non-linear stress-dependent permeability function. The model was first tested and calibrated by modelling two different laboratory gas migration tests on Wyoming (MX-80) bentonite samples. The calibrated model was then applied to model gas migration through a bentonite buffer of a large-scale gas injection test (Lasgit) conducted at the Äspö Hard Rock Laboratory in Sweden. Observed preferential gas migration along interfaces (between compacted blocks and along the canister surface) required explicit representation of such interfaces in the model. The model with the stress-dependent gas permeability accurately captured observed experimental responses in terms of gas breakthrough time, peak gas pressure, and cumulative gas flow rates. The calibrated model was finally applied to simulate migration of hydrogen gas generated within a breached nuclear waste canister over 10,000 years, involving migration of much larger gas volumes. For the considered gas generation rate and host rock properties, the generated gas could migrate through the bentonite buffer and released into the surrounding host rock at a maximum gas pressure somewhat higher than the initial total stress, though a significant amount of hydrogen remained within the buffer. This modelling sets the stage for further detailed analysis of the impact of hydrogen gas generation on the long-term performance of nuclear waste repositories.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Deacetylation and Mechanical Refining Pathway for the Bioconversion of Sugarcane Bagasse

Advancing lignocellulose biorefining is imperative for the deployment of cellulosic (2G) biofuels. This work investigates the tailoring of the alkaline deacetylation and mechanical refining (DMR) pathway for the bioconversion of sugarcane bagasse. Experiments are conducted at laboratory and pilot scales, varying the pretreatment conditions (70–92 °C; 48–100 g NaOH /kg) and the mechanical refining technologies (PFI and disk refining). The pretreatments selectively solubilize acetyl groups (> 86%) and lignin (10–63%) while mostly preserving structural carbohydrates in the solid phase. Enzymatic hydrolysis generates hydrolysates of clean sugars (glucose and xylose), with sugar yields increasing up to 81% for glucose and 89% for xylose in response to delignification and mechanical refining. Biochemical methane potential assays reveal specific methane productions of up to 568 NmL CH₄ gVS⁻¹ for alkaline liquor monodigestion and 344 NmL CH₄ gVS⁻¹ for co-digestion with sugarcane vinasse from the conventional (1G) sugarcane ethanol, indicating a strong potential for bioenergy recovery from this process stream. Synergies are identified in integrating 1G ethanol, 2G DMR processing of bagasse, and anaerobic co-digestion of 1G vinasse and 2G DMR alkaline liquor. This technology enables sugarcane biorefineries to enhance the co-production of ethanol, methane, and concentrated streams of CO 2 .

09 BIOMASS FUELS

Status Report on Design of In-situ Thermomechanical Testing at LANSCE

Nuclear fuel encounters severe thermomechanical environments in which its mechanical response is determined by its microstructure, temperature and stress level histories. Simulating the response of such microstructures is crucial for predicting both performance and transient fuel mechanical responses and experimental verification of such predictions is therefore of great interest. While most of the deformation in a nuclear fuel rod occurs in the cladding, deformation of the fuel itself is still of interest with deformation mechanisms at operating temperature and above including creep, swelling, cracking as well as pellet-clad interaction. Characterization of these properties and understanding of the underlying deformation phenomena at operating or excursion temperatures is therefore of great importance for development and ultimately licensing of improved and novel nuclear fuel forms. Diffraction techniques offer unique insight on the atomistic (e.g. crystal structure) and microstructure (e.g. phase transformations, texture, defects) length scales and have a long history of providing unique data to inform relevant deformation models that enable the required predictive capabilities. For example, dislocations lead to diffraction peak broadening that can be characterized to estimate the dislocation density and study the role of dislocations on the deformation while measuring lattice strains allows to studie load sharing in two phase materials. In this report the requirements for a sample environment for high temperature deformation of nuclear fuels are defined. The HIPPO neutron time-of-flight diffractometer at LANSCE will host this sample environment and is also described. This instrument covers diffraction angles from 140° to 40° and is also equipped with an event-mode neutron imaging detector system, enabling energy-resolved neutron imaging in parallel with the diffraction that could measure sample temperature from Doppler broadening of neutron absorption resonances or measure pore densities from changes in the attenuation. Designs of devices to characterize thermomechanical properties of nuclear fuel without diffraction are also considered to guide the design. While this report is focused on applications for nuclear fuels, the device can also characterize cladding, moderator or structural materials and therefore contribute to other fields of research and development for advanced reactors. The temperatures planned to be reached are above 2000℃, thus enabling characterization of LWR reactor fuels under accident scenarios but also reaching temperatures of fuels developed for nuclear thermal propulsion and providing opportunities to characterize those. In conjunction with the energy-resolved neutron imaging detector, this setup would allow to measure neutron cross-sections at high temperatures, filling a gap towards development of reactors operating at high temperatures.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS

Critical structural parameter determining magnetic phases in the F⁡e 2 ⁢M⁢o 3 ⁢O 8 altermagnet system

A systematic structural investigation of F⁡e 2 ⁢M⁢o 3⁢ O 8 as a function of pressure, temperature, and magnetic field reveals that the 𝑃⁢6 3 ⁢𝑚⁢𝑐 space group remains robust over a wide range of conditions. No changes in the long-range crystal structure are observed for pressures up to 10 GPa, temperatures between 11 K and 300 K, and magnetic fields up to 9 T. The magnetostructural response, quantified by Δ⁢𝑐/𝑐, is determined for magnetic fields applied transverse to the crystallographic 𝑐 axis, demonstrating strong magnetoelastic coupling. The well-known magnetic-field-induced transition is confirmed to be first order and isostructural, occurring between two distinct altermagnetic states. Importantly, this transition can also be accessed using magnetic fields applied within the 𝑎𝑏 plane. When expressed in terms of the 𝑐/𝑎 ratio, the structural evolution under pressure (0–10 GPa) maps continuously onto that observed across the full Zn substitution range in (F⁡e 1−𝑦⁢ Z⁢n 𝑦 )2 M⁢o 3 ⁢O 8 ⁢(0 ≤ 𝑦 ≤ 1), indicating a common underlying structural control parameter. Furthermore, these results demonstrate that the 𝑐-axis lattice parameter is the primary structural degree of freedom governing magnetic behavior under external tuning parameters. More broadly, they indicate that magnetic order in this family of complex oxides (𝐴 2 ⁢M⁢o 3 ⁢O 8 , 𝐴 = Co,Mn,Ni) can be systematically tuned via pressure. This tunability highlights these materials as promising platforms for controllable magnetic switching, particularly in thin films integrated with piezoelectric substrates.

36 MATERIALS SCIENCE

Influence of Photoemission Geometry on Timing and Efficiency in 4D Ultrafast Electron Microscopy

Broader adoption of 4D ultrafast electron microscopy (UEM) for the study of chemical, materials, and quantum systems is being driven by development of new instruments as well as continuous improvement and characterization of existing technologies. Perhaps owing to the still-high barrier to entry, the full range of capabilities of laser-driven 4D UEM instruments has yet to be established, particularly when operated at extremely low beam currents (~fA). Accordingly, with an eye on beam stability, we have conducted particle tracing simulations of unconventional off-axis photoemission geometries in a UEM equipped with a thermionic-emission gun. Specifically, we have explored the impact of experimentally adjustable parameters on the time-of-flight (TOF), the collection efficiency (CE), and the temporal width of ultrashort photoelectron packets. The adjustable parameters include the Wehnelt aperture diameter (DW), the cathode set-back position (Ztip), and the position of the femtosecond laser on the Wehnelt aperture surface relative to the optic axis (Rphoto). Notable findings include significant sensitivity of TOF to DW and Ztip, as well as non-intuitive responses of CE and temporal width to varying Rphoto. As a means to improve accessibility, practical implications and recommendations are emphasized wherever possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Cooperative Charge Ordering Signature of Trimer Molecules in Infinite-Layer CaCoO 2

The infinite-layer compound CaCoO 2 provides a compelling platform to explore how local orbital distortions evolve into collective electronic order. Using resonant soft x-ray scattering (RSXS) and its time-resolved counterpart (tr-RSXS), we reveal a cooperative coupling between two intertwined orders at the same in-plane wave vector, 𝑞 ≈ (1/4, 1/4, 0), comprising a quasi-two-dimensional orbital order associated with the distorted Co (1) sites and a quasi-three-dimensional trimerlike charge order stabilized through Co-O-Ca hybridization. The charge-sensitive component exhibits a uniform relaxation time of ∼400 fs (391 ± 5 fs) across multiple Co resonances, evidencing a collective charge response. A comparable timescale observed in the orbital order associated with the highly distorted Co (1) sites indicates that this collective behavior naturally originates from the same orbital framework that mediates both charge and orbital degrees of freedom. Together, these findings establish CaCoO 2 as a model system in which Jahn-Teller distortions, orbital hybridization, and geometric frustration cooperate to form molecular trimer units that dictate the macroscopic electronic behavior.

Lee, J.-S. [SLAC National Accelerator Laboratory (

Atoms-to-devices models for post-irradiation annealing kinetics in III-V semiconductors

Understanding atomic displacement damage and propagating this understanding into validated defect-aware radiation damage models is a critical step toward achieving predictive, atoms-to-systems assessment of radiation vulnerabilities of microelectronics. These models are essential for evaluating the performance and reliability of nuclear stockpile components or space-based electronics (satellites). The goal of this project was to bridge a fundamental science gap between defect spectroscopies and atomistic modeling of radiation-induced defects on the one hand, and simulations of device response on the other hand, to develop predictive defect-aware damage models for high-fidelity device simulations of radiation-damaged III-V semiconductors.

36 MATERIALS SCIENCE

Semi-Dirac Fermions in a Topological Metal

Topological semimetals with massless Dirac and Weyl fermions represent the forefront of quantum materials research. In two dimensions, a peculiar class of fermions that are massless in one direction and massive in the perpendicular direction was predicted 16 years ago. These highly exotic quasiparticles—the semi-Dirac fermions—ignited intense theoretical and experimental interest but remain undetected. Using magneto-optical spectroscopy, we demonstrate the defining feature of semi-Dirac fermions— B 2 / 3 scaling of Landau levels—in a prototypical nodal-line metal ZrSiS. In topological metals, including ZrSiS, nodal lines extend the band degeneracies from isolated points to lines, loops, or even chains in the momentum space. With calculations and theoretical modeling, we pinpoint the observed semi-Dirac spectrum to the crossing points of nodal lines in ZrSiS. Crossing nodal lines exhibit a continuum absorption spectrum but with singularities that scale as B 2 / 3 at the crossing. Our work sheds light on the hidden quasiparticles emerging from the intricate topology of crossing nodal lines and highlights the potential to explore quantum geometry with linear optical responses. Published by the American Physical Society 2024

Shao, Yinming (ORCID:0000000228910028)

Response to the comments on creep-induced elemental redistribution at grain boundaries of 304L stainless steel – An experimental evidence for diffusional creep mechanisms

Sandström and He reanalyzed our 304L creep data (Acta Materialia 121137, 2025) at 10 MPa/750 °C and proposed that: (i) primary dislocation creep dominates based on the ϕ -model; (ii) measured creep rates deviate from published 304H data and Coble creep predictions using literature grain boundary (GB) diffusivity; and (iii) creep-induced elemental redistribution (CIER) could be explained by dislocation emission/absorption (DEAA) and GB creep (GBC), as opposed to diffusional creep. Here, in this response, we propose that: (i) the creep behavior at 15 MPa/750 °C—for which we demonstrated CIER—is only partially captured by the ϕ -model; (ii) direct comparison between 304L and 304H creep data is inappropriate due to carbides in 304H, and Coble creep predictions estimated from CALPHAD-based bulk diffusivity span experimental rates; and (iii) CIER exhibits distinct fingerprint for diffusional creep characterized by anisotropy and stress orientation dependence whereas neither DEAA nor GBC have been mechanistically substantiated.

36 - MATERIALS SCIENCE

Virtual Reality for Shoot/No-Shoot Decision Training in Law Enforcement: A Literature Review and Research Agenda

Virtual reality (VR) can materially improve “shoot / no-shoot” (SNS) training by giving officers realistic, repeatable practice making high-stakes decisions under pressure. Traditional tools—live-fire ranges and video simulators—build basics, but they cannot adapt to each officer in real time or fully mirror the complexity of the field. VR closes that gap by creating immersive scenarios that are safer, more flexible, easier to scale across units, and able to capture objective performance data. SNS decisions are not just about marksmanship; they rely on perception, judgment, memory, and the ability to hold fire when a threat is uncertain. Effective training therefore needs realism, decision complexity, and branching outcomes that reflect the true consequences of choices. These elements strengthen recognition of hostile intent while reducing false positives and building the self-control required in ambiguous situations. VR brings specific advantages: dynamic environments, full-body interaction, and the ability to measure performance with precision—enabling targeted feedback and better transfer of learning to the street. At the same time, responsible deployment must address scenario quality (credible environments and behaviors), lawful decision models, and user wellbeing (appropriate stress levels, comfort, and safety). Sandia’s VIPER Lab is positioned to lead this work. The team combines human-performance science, AI/ML, and VR/AR development with a deep equipment bench (e.g., omnidirectional treadmill, eye-tracking, haptics, multiple HMDs). This ecosystem supports building and validating next-generation SNS training that is immersive, measurable, and trustworthy. Bottom line: Investment in VR-enabled SNS training that blends evidence-based design with careful validation and legal safeguards is expected to pay off in safer, more consistent decision-making and improved community trust, delivered through training that is practical to deploy at scale.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF

Reduced Order Modeling conditioned on monitored features for response and error bounds estimation in engineered systems

Reduced Order Models (ROMs) form essential tools across engineering domains by virtue of their function as surrogates for computationally intensive digital twinning simulators. Although purely data-driven methods are available for ROM construction, schemes that allow to retain a portion of the physics tend to enhance the interpretability and generalization of ROMs. However, physics-based techniques can adversely scale when dealing with nonlinear systems that feature parametric dependencies. This study introduces a generative physics-based ROM that is suited for nonlinear systems with parametric dependencies and is additionally able to provide numerical error bounds associated with the respective estimates. A main contribution of this work is the conditioning of these parametric ROMs to features that can be derived from monitoring measurements, feasibly in an online fashion. This is contrary to most existing ROM schemes, which remain restricted to the prescription of the physics-based, and usually a priori unknown, system parameters. Our work utilizes conditional Variational Autoencoders to continuously map the required reduction bases to a feature vector extracted from limited output measurements, while additionally allowing for a probabilistic assessment of the ROM-estimated Quantities of Interest. An auxiliary task using a neural network-based parametrization of suitable probability distributions is introduced to re-establish the link with physical model parameters. We verify the proposed scheme on a series of simulated case studies incorporating effects of geometric and material nonlinearity under parametric dependencies related to system properties and input load characteristics.

Conditional VAEs

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes

Detection of scintillation light in noble gases with wavelength-shifting optical fibers

Wavelength-shifting (WLS) techniques enable particle detectors based on noble gases, whose scintillation light is predominantly emitted in the vacuum-ultraviolet. We investigate WLS fibers coated with tetraphenyl butadiene (TPB) for scintillation light detection in gaseous xenon and argon at pressures up to 8.5 bar, motivated by future high-pressure xenon time-projection chambers of the NEXT program. Two detector configurations are studied: an elongated high-pressure vessel with four PTFE panels equipped with WLS fibers read by temperature-stabilized SiPMs, and a compact box-shaped detector operated at 1 bar Xe with WLS fibers read out by PMTs. Both operate with continuous gas purification. The detector response is characterized using cosmic muons and alpha particles from a $^{241}$Am source. With the SiPM setup, we measure a light collection efficiency (LCE) of ${1.18 \pm 0.01~\mathrm{(sta.)}~^{+0.07}_{-0.09}~\mathrm{(sys.)}~\%}$ for xenon and ${1.07 \pm 0.01~\mathrm{(sta.)}~^{+0.06}_{-0.08}~\mathrm{(sys.)}~\%}$ for argon. With PMT readout, we measure a LCE of ${0.45 \pm 0.01~\mathrm{(sta.)} \pm 0.05~\mathrm{(sys.)}~\%}$ in xenon, in agreement with the SiPM result once photon detection efficiency is accounted for. Average scintillation waveforms in xenon and argon are studied to assess the time structure of the emitted light. Cosmic-muon measurements yield a mean energy required to produce a scintillation photon $45\pm7~\mathrm{(sta.)}~^{+4}_{-5}~\mathrm{(sys.)}~\mathrm{eV}$ at 1.5 bar, in agreement with the literature. The results demonstrate that TPB-coated WLS fiber systems can reliably detect scintillation light in high-pressure gaseous noble detectors, with a LCE representing an upper limit for realistic large-scale TPCs, where additional photon losses from materials and fiber attenuation are expected.

Soleti, S. R. [Donostia Intl. Phys. Ctr., San Seba

Understanding the structural mechanics of ligated DNA crystals via molecular dynamics simulation

DNA self-assembly is a highly programmable method to construct arbitrary architectures based on sequence complementarity. Among various constructs, DNA crystals are macroscopic crystalline materials formed by assembling motifs via sticky end association. Due to their high structural integrity and size ranging from tens to hundreds of micrometers, DNA crystals offer unique opportunities to study the structural properties and deformation behaviors of DNA assemblies. For example, enzymatic ligation of sticky ends can selectively seal nicks resulting in more robust structures with enhanced mechanical properties. However, the research efforts have been mostly on experiments involving different motif designs, structural optimization, or new synthesis methods, while their mechanics are not yet fully understood. The complex properties of DNA crystals are difficult to study via experiments alone, and numerical simulation can complement and aid the experiments. The coarse-grained molecular dynamics (MD) simulation is a powerful tool that can probe the mechanics of DNA assemblies. Here, we investigate DNA crystals made of four different motif lengths with various ligation patterns (full ligation, major directions, connectors, and in-plane) using oxDNA, an open-source, coarse-grained MD platform. We found that several distinct deformation stages emerge in response to mechanical loading and that the number and the location of ligated nucleotides can significantly modulate structural behaviors. These findings should be useful for predicting crystal properties and thus improving the design.

DNA crystal