Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamic (MD) simulation”

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 343 records · Page 19

Clarifying the definition of ‘transonic’ screw dislocations

We report a number of recent Molecular Dynamics (MD) simulations have demonstrated that screw dislocations in face centred cubic (fcc) metals can achieve stable steady state motion above the lowest shear wave speed (v shear ) which is parallel to their direction of motion (often referred to as transonic motion). This is in direct contrast to classical continuum analyses which predict a divergence in the elastic energy of the host material at a crystal geometry dependent ‘critical’ velocity v crit . Within this work, we first demonstrate through analytic analyses that the elastic energy of the host material diverges at a dislocation velocity (v crit ) which is greater than v shear , i.e. v crit >v shear . We argue that it is this latter derived velocity (v crit ) which separates ‘subsonic’ and ‘supersonic’ regimes of dislocation motion in the analytic solution. In addition to our analyses, we also present a comprehensive suite of MD simulation results of steady state screw dislocation motion for a range of stresses and several cubic metals at both cryogenic and room temperatures. At room temperature, both our independent MD simulations and the earlier works find stable screw dislocation motion only below our derived v crit . Nonetheless, in real-world polycrystalline materials v crit cannot be interpreted as a hard limit for subsonic dislocation motion. In fact, at very low temperatures our MD simulations of Cu at 10 Kelvin confirm a recent claim in the literature that true ‘supersonic’ screw dislocations with dislocation velocities v>v crit are possible at very low temperatures.

74 ATOMIC AND MOLECULAR PHYSICS↗

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗

Resolving the Solvation Structure and Transport Properties of Aqueous Zinc Electrolytes from Salt-in-Water to Water-in-Salt Using Neural Network Potential

Zn Cl 2 solutions are promising electrolytes for aqueous zinc-ion batteries. Here, we report a joint computational and experimental study of the structural and dynamic properties of aqueous Zn Cl 2 electrolytes with concentrations ranging from salt-in-water to water-in-salt (WIS). By developing a neural network potential (NNP) model, we perform molecular dynamics (MD) simulations with accuracy but at much larger lengths and longer timescales. The NNP predicted structures are validated by the structure factors measured by X-ray total scattering experiments. The MD trajectories provide a comprehensive and quantitative picture of the Zn 2 + solvation shell structures. Additionally, we find that the O − H covalent bonds in water are strengthened with increasing salt concentration, thus expanding the electrochemical stability window of aqueous electrolytes. In terms of dynamic properties, the calculated and experimentally measured conductivities are in good agreement. Through the analysis of the calculated cation transference number, we propose a three-stage charge carrier transport mechanism with increasing concentration: independent ion transport, strongly correlated ion transport, and small positive charge carrier diffusion through negatively charged polymeric clusters. Our study provides fundamental atomic scale insights into the structure and transport properties of the Zn Cl 2 electrolyte that can aid the optimization and development of WIS electrolytes. Published by the American Physical Society 2025

25 ENERGY STORAGE↗

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Abstract Although F‐Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter‐ and intra‐Chain interactions responsible for self‐healing of F‐Containing copolymers composed of pentafluorostyrene and n‐butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ‐σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic‐σ‐lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self‐healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic‐σ‐lock (~70 %) contributions. These interactions are significantly diminished for non‐self‐healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non‐Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self‐healable batteries and energy storage devices.

Gaikwad, Samruddhi↗

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Although F-Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter- and intra-Chain interactions responsible for self-healing of F-Containing copolymers composed of pentafluorostyrene and n-butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ-σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic-σ-lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self-healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic-σ-lock (~70 %) contributions. These interactions are significantly diminished for non-self-healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non-Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self-healable batteries and energy storage devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Self‐Healable Poly(ionic liquid) Copolymers Driven by Polar and Dipolar Forces

Abstract Commodity aliphatic and aromatic acrylic‐based copolymers self‐heal due to ubiquitouskey‐and‐lock,ring‐and‐lock, andfluorophilic‐σ‐lockvan der Waals (vdW) interactions. However, the role of these interactions in the presence of covalently copolymerized ionic liquid (IL) is not known. This study is driven by the hypothesis that covalently incorporated cation–anion pairs to form poly(ionic liquid) copolymers (PILCs) can perturb inter‐ or intra‐chain vdW interactions reflected in mechanical and electrical responses. To test this hypothesis, we synthesized a series of PILCs comprising of pentafluorostyrene (PFS) and imidazolium‐based IL monomers with variable‐length aliphatic tails (methyl and butyl). Using a combination of 2D 1 H‐ 1 H and 19 F ‐ 19 F NOESY NMR and FTIR measurements supplemented by molecular dynamic (MD) simulations, these studies demonstrate that preferentially alternating/random PILCs topologies facilitate self‐healing. The introduction of cation–anion moieties modifies thefluorophilic‐σ‐lockinteractions and, along with longer aliphatic tails ─(CH 2 ) 3 CH 3 covalently attached to the imidazolium cation, enhances cation‐anion mobility, thus faster recovery from mechanical damage occurs. These findings underline how precise control over dipolar and ionic interactions through copolymer composition enables self‐healing in PILCs. These insights may open pathways for designing sustainable, mechanically resilient materials for applications in energy storage and energy harvesting.

Chemistry↗

Underlying Roles of Polyol Additives in Promoting CO2 Capture in PEI/Silica Adsorbents

Solid-supported amines having low molecular weight branched poly(ethylenimine) (PEI) physically impregnated into porous solid supports are promising adsorbents for CO2 capture. Co-impregnating short-chain poly(ethylene glycol) (PEG) together with PEI alters the performance of the adsorbent, delivering improved amine efficiency (AE, mol CO2 sorbed / mol N) and faster CO2 uptake rates. To uncover the physical basis for this improved gas capture performance, we probed the distribution and mobility of the polymers in the pores via small angle neutron scattering (SANS), solid-state NMR, and molecular dynamic (MD) simulation studies. SANS and MD simulations reveal that PEG displaces wall-bound PEI, making amines more accessible for CO2 sorption. Solid-state NMR and MD simulation suggest intercalation of PEG into PEI domains, separating PEI domains and reducing amine-amine interactions, providing potential PEG-rich and amine-poor interfacial domains that bind CO2 weakly via physisorption while providing facile pathways for CO2 diffusion. Contrary to a prior literature hypothesis, no evidence is obtained for PEG facilitating PEI mobility in solid supports. Instead, the data suggest that PEG chains coordinate to PEI, form larger bodies with reduced mobility compared to PEI alone. We also demonstrate promising CO2 uptake and desorption kinetics at varied temperatures, given by favorable amine distribution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Synthesis of Layer-Tunable and Transfer-Free Graphene on Device-Compatible Substrates Using Ion Implantation Toward Versatile Applications

Direct synthesis of layer-tunable and transfer-free graphene on technologically important substrates is highly valued for various electronics and device applications. State of the art in the field is currently a two-step process: a high-quality graphene layer synthesis on metal substrate through chemical vapor deposition (CVD) followed by delicate layer transfer onto device-relevant substrates. Here, we report a novel synthesis approach combining ion implantation for a precise graphene layer control and dual-metal smart Janus substrate for a diffusion-limiting graphene formation to directly synthesize large area, high quality, and layer-tunable graphene films on arbitrary substrates without the post-synthesis layer transfer process. Carbon (C) ion implantation was performed on Cu–Ni film deposited on a variety of device-relevant substrates. A well-controlled number of layers of graphene, primarily monolayer and bilayer, is precisely controlled by the equivalent fluence of the implanted C-atoms (1 monolayer ~4 × 10 15 C-atoms/cm 2 ). Upon thermal annealing to promote Cu-Ni alloying, the pre-implanted C-atoms in the Ni layer are pushed toward the Ni/substrate interface by the top Cu layer due to the poor C-solubility in Cu. As a result, the expelled C-atoms precipitate into a graphene structure at the interface facilitated by the Cu-like alloy catalysis. After removing the alloyed Cu-like surface layer, the layer-tunable graphene on the desired substrate is directly realized. The layer-selectivity, high quality, and uniformity of the graphene films are not only confirmed with detailed characterizations using a suite of surface analysis techniques but more importantly are successfully demonstrated by the excellent properties and performance of several devices directly fabricated from these graphene films. Molecular dynamics (MD) simulations using the reactive force field (ReaxFF) were performed to elucidate the graphene formation mechanisms in this novel synthesis approach. With the wide use of ion implantation technology in the microelectronics industry, this novel graphene synthesis approach with precise layer-tunability and transfer-free processing has the promise to advance efficient graphene-device manufacturing and expedite their versatile applications in many fields.

36 MATERIALS SCIENCE↗

Structural Changes to the Gd‐DTPA Complex at Varying Ligand Protonation State

Abstract Diethylenetriaminepentaacetic acid (DTPA) is a chelating agent whose complex with the Gd 3+ ion is used in medical imaging. DTPA is also used in lanthanide‐actinide separation processes. As protonation of the DTPA ligand can facilitate dissociation of the Gd 3+ ion from the Gd‐DTPA complex, this work investigates the coordination structures of the aqueous Gd 3+ ion and its environment when chelated by DTPA in eight different DTPA protonation states. Both classical and ab initio molecular dynamics (MD) simulations are conducted to model the solvated complexes. Extended X‐ray absorption fine structure (EXAFS) measurements of the Gd 3+ aqua ion, and the Gd‐DTPA complex at pH 1 and 11, are compared to EXAFS spectra predicted from the MD simulations to verify the accuracy of the MD structures. The findings of this work provide atomic‐level details into the fluctuating Gd‐DTPA complex environment as the DTPA ligand gradually detaches from the Gd 3+ ion with increased protonation.

Chemistry↗

Methyl α‐D‐galactopyranosyl‐(1→3)‐β‐D‐galactopyranoside and methyl β‐D‐galactopyranosyl‐(1→3)‐β‐D‐galactopyranoside: Glycosidic linkage conformation determined from MA'AT analysis

Abstract MA'ATanalysis has been applied to two biologically‐importantO‐glycosidic linkages in two disaccharides, α‐D‐Galp‐(1→3)‐β‐D‐GalpOMe (3) and β‐D‐Galp‐(1→3)‐β‐D‐GalpOMe (4). Using density functional theory (DFT) to obtain parameterized equations relating a group of trans‐O‐glycosidic NMR spin‐couplings to eitherphi(ϕ') orpsi(ψ'), and experimental 3 J COCH , 2 J COC , and 3 J COCC spin‐couplings measured in aqueous solution in 13 C‐labeled isotopomers, probability distributions ofϕ'andψ'in each linkage were determined and compared to those determined by aqueous 1‐μs molecular dynamics (MD) simulation. Good agreement was found between theMA'ATand single‐state MD conformational models of these linkages for the most part, with modest (approximately <15°) differences in the mean values ofϕ'andψ', although the envelope of allowed angles (encoded in circular standard deviations or CSDs) is consistently larger forϕ'determined fromMA'ATanalysis than from MD for both linkages. TheMA'ATmodel of the α‐Galp‐(1→3)‐β‐Galplinkage agrees well with those determined previously using conventional NMR methods ( 3 J COCH values and/or 1 H‐ 1 H NOEs), but some discrepancy was observed for the β‐Galp‐(1→3)‐β‐Galplinkage, which may arise from errors in the conventions used to describe the linkage torsion angles. Statistical analyses of X‐ray crystal structures show ranges ofϕ'andψ'for both linkages that include the mean angles determined fromMA'ATanalyses, although both angles adopt a wide range of values in the crystalline state, withϕ'in β‐Galp‐(1→3)‐β‐Galplinkages showing greater‐than‐expected conformational variability.

Chemistry↗

Molecular dynamics analysis of the binding of human interleukin-6 with interleukin-6 α-receptor

Human interleukin-6 (hIL-6) is a multifunctional cytokine that regulates immune and inflammatory responses in addition to metabolic and regenerative processes and cancer. hIL-6 binding to the IL-6 receptor (IL-6Rα) induces homodimerization and recruitment of the glycoprotein (gp130) to form a hexameric signaling complex. Anti-IL-6 and IL-6R antibodies are clinically approved inhibitors of IL-6 signaling pathway for treating rheumatoid arthritis and Castleman's disease, respectively. There is a potential to develop novel small molecule IL-6 antagonists derived from understanding the structural basis for IL-6/IL-6Rα interactions. In this work, we combine homology modeling with extensive molecular dynamics (MD) simulations to examine the association of hIL-6 with IL-6Rα. A comparison with MD of apo hIL-6 reveals that the binding of hIL-6 to IL-6Rα induces structural and dynamic rearrangements in the AB loop region of hIL-6, disrupting intraprotein contacts and increasing the flexibility of residues 48 to 58 of the AB loop. In contrast, due to the involvement of residues 59 to 78 in forming contacts with the receptor, these residues of the AB loop are observed to rigidify in the presence of the receptor. The binary complex is primarily stabilized by two pairs of salt bridges, Arg181 (hIL-6)- Glu182 (IL-6Rα) and Arg184 (hIL-6)- Glu183 (IL-6Rα) as well as hydrophobic and aromatic stacking interactions mediated essentially by Phe residues in both proteins. An interplay of electrostatic, hydrophobic, hydrogen bonding, and aromatic stacking interactions facilitates the formation of the hIL-6/IL-6Rα complex.

59 BASIC BIOLOGICAL SCIENCES↗

Co-Assembly of Carbon Nanotube Porins into Biomimetic Peptoid Membranes

Robust and cost-effective membrane-based separations are essential to solving many global crises, such as the lack of clean water. Even though the current polymer-based membranes are widely used for separations, their performance and precision can be enhanced by using a biomimetic membrane architecture that consists of highly permeable and selective channels embedded in a universal membrane matrix. Researchers have shown that artificial water and ion channels, such as carbon nanotube porins (CNTPs), embedded in lipid membranes can deliver strong separation performance. However, their applications are limited by the relative fragility and low stability of the lipid matrix. In this work, we demonstrate that CNTPs can co-assemble into two-dimension (2D) peptoid membrane nanosheets, opening up a way to produce highly programmable synthetic membranes with superior crystallinity and robustness. A combination of molecular dynamics (MD) simulations, Raman spectroscopy, X-ray diffraction (XRD), and atomic force microscopy (AFM) measurements to verify the co-assembly of CNTP and peptoids are used and show that it does not disrupt peptoid monomer packing within the membrane. In conclusion, these results provide a new option for designing affordable artificial membranes and highly robust nanoporous solids.

36 MATERIALS SCIENCE↗

Effect of Methyl Groups on Formation of Ordered or Layered Graphitic Materials from Aromatic Molecules

Abstract Developing functionally complex carbon materials from small aromatic molecules requires an understanding of how the chemistry and structure of its constituent molecules evolve and crosslink, to achieve a tailorable set of functional properties. Here, molecular dynamics (MD) simulations are used to isolate the effect of methyl groups on condensation reactions during the oxidative process and evaluate the impact on elastic modulus by considering three monodisperse pyrene‐based systems with increasing methyl group fraction. A parameter to quantify the reaction progression is designed by computing the number of new covalent bonds formed. Utilizing the previously developed MD framework, it is found that increasing methylation leads to an almost doubling of bond formation, a larger fraction of the new bonds oriented in the direction of tensile stress, and a higher basal plane alignment of the precursor molecules along the direction of tensile stress, resulting in enhanced tensile modulus. Additionally, via experiments, it is demonstrated that precursors with a higher fraction of methyl groups result in a higher alignment of molecules. Moreover, increased methylation results in the lower spread of single molecule alignment which may lead to smaller variations in tensile modulus and more consistent properties in carbon materials derived from methyl‐rich precursors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ionic Liquid‐Driven Modulation of DNA Brush Morphology on Nanoparticle Surfaces

Here, the morphology of DNA is strongly influenced by its surrounding environment, including factors such as pH, salt type and valency, and the presence of polymers. Inorganic salts are known to reduce the DNA chain length through mechanisms like electrostatic screening and ion bridging. In contrast, ionic liquids, a new class of organic salts, have previously been found to increase the DNA chain length, indicating a distinct mode of interaction between the ionic liquid and DNA chains. This study utilizes self-assembled DNA-AuNPs as a model system to examine changes in the DNA chain morphology and the nanoscale interaction mechanisms in an ionic liquid environment. The DNA chain lengths are measured in solution using X-ray scattering measurements at varying concentrations of two imidazolium ([$BMIM$] acetate and [$EMIM$] acetate) based ionic liquids. Additionally, Molecular Dynamics (MD) simulations are performed mimicking the experimental system. Our results suggest an interplay of electrostatic and groove-binding interactions governing the DNA chain morphology, which depends on IL concentration and the composition of the DNA chains. It has been found that for DNA chains with majority ssDNA, electrostatic interaction dominate, however with increasing composition of double strands, the DNA chains exhibit compaction due to a non-electrostatic hydrophobic groove-binding mechanism.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Structure/property relationship of semicrystalline polymers during tensile deformation: a molecular dynamics approach

Molecular dynamics (MD) simulation with an ab initio numerical model is conducted to study the structure–property relationship of semicrystalline polymers during uniaxial deformation. The effects of chain length, temperature, and strain rate on mechanical properties are discussed. The influences of microstructural evolution such as bond length, entanglement density, and chain orientation are also studied quantitatively. The temperature will greatly affect the chain conformation in amorphous domains, and the results revealed that the interaction of amorphous and crystalline domains played a crucial role during stretching. Accordingly, the yielding of semicrystalline polymers follows different mechanisms at temperatures above and below T g . A melt-recrystallization scheme is observed during yielding at higher temperatures, while destruction of crystal structures is observed at lower temperatures at the yield point. The correlated effects of different temperatures and strain rates on mechanical properties are examined. This research is part of our efforts to develop a digital twin of a real experiment for efficient optimization of polymer material properties.

36 MATERIALS SCIENCE↗

Relaxation dynamics of water in the vicinity of cellulose nanocrystals

Water molecules near cellulose nanocrystals (CNCs; produced via the sulfuric acid-catalyzed hydrolysis of wood pulp) are believed to relax slower than those in the bulk liquid, which may result in unique properties of CNC aqueous dispersions. Herein this study analyzed the polarization behavior of water molecules in CNC aqueous dispersions and other reference samples using a dielectric relaxation spectroscopy (DRS) technique in the microwave frequency range (0.2–20 GHz). As the CNC concentration increases, two slow relaxation components become prominent. The comparison with DRS data of aqueous dispersions of nanoporous silica, polyvinyl alcohol (PVA), and hairy CNC (HCNC) with amorphous chains protruding from both ends suggested that these slow relaxation modes of water near CNC surfaces cannot be attributed to direct hydrogen bonding interactions with the hydroxyl (OH) groups exposed and immobilized at the solid surface. Instead, they are similar to the water molecules interacting with OH groups attached to flexible polymer chains. Molecular dynamics (MD) simulations of the polarization behavior of water near the (110) facet of cellulose Iβ crystals confirmed that the interactions of water molecules with the cellulose crystal surface do not cause slower relaxations in the frequency range studied via the DRS. These results indicated that the CNC surface cannot be depicted with the crystallographic facets of cellulose Iβ; instead, it resembles a polymer-brush surface on which the short glucan residues or fragments of the strong acid-catalyzed hydrolysis process are swollen and extended into the aqueous phase.

59 BASIC BIOLOGICAL SCIENCES↗

Virtual texture analysis to investigate the deformation mechanisms in metal microstructures at the atomic scale

Understanding the deformation behavior of metallic materials at high strain rates requires the characterization of plasticity contributors such as twins, phase transformed regions, and dislocations. However, predicting the contributions from phase transformation and twinning relies on a complete understanding of the selection of variants for various loading orientations and the evolution of their volume fractions. This manuscript presents a new virtual texture (VirTex) analysis approach to characterize phase transformation and twinning variants in deformed microstructures generated using molecular dynamics (MD) simulations. Furthermore, the VirTex method involves the construction of a rotation matrix to calculate the angle/axis pairs and misorientation angles for each atom in the microstructure. Any changes in the orientation angle from angle/axis pairs and/or structure types are analyzed to determine the nucleation and evolution of variants in the microstructure. The study uses shock deformed single-crystal Fe, Ta, and Cu to analyze the variant selections for phase transformation or twinning or both in BCC and FCC systems. In addition, the VirTex analysis is able to characterize the phase transformation and twinning variants in nanocrystalline Fe and Ta microstructures. Besides characterizing variants, orientation mapping also provides an accelerated and on-the-fly approach for quantifying twin fractions in MD microstructures.

36 MATERIALS SCIENCE↗