skip to main content


Title: Molecular dynamics simulations to understand the mechanical behavior of functional gradient nano-gyroid structures
Gyroid structure, a nature inspired cellular architecture, is under extensive exploration recently due to its structure continuity, uniform stress distribution under compression, and stable collapse mechanism during deformation. However, when combining with a functional gradient, the Gyroid structure can perform much different mechanical behavior from its homogeneous counterpart. Herein, bottom-up computational modeling is performed to investigate the mechanics of functional gradient nano-gyroid structure made of copper (Cu). Our work reveals that its mechanical properties degrade with a density that is much slower than those of homogeneous gyroid structure. The scaling of yield strength [Formula: see text] to the relative density [Formula: see text] for the functional gradient gyroid structure is in the factor of 1.5. Moreover, the layer-by-layer collapsing mechanism yields significantly better mechanical energy absorption ability. This study not only leads to insightful understanding of the deformation mechanisms in nonuniform gyroid structures but also promotes the development of the functional gradient cellular materials.  more » « less
Award ID(s):
1762792 1826439
NSF-PAR ID:
10359730
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Applied Physics
Volume:
132
Issue:
13
ISSN:
0021-8979
Page Range / eLocation ID:
135109
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. All-solid-state batteries (ASSBs) have garnered increasing attention due to the enhanced safety, featuring nonflammable solid electrolytes as well as the potential to achieve high energy density. 1 The advancement of the ASSBs is expected to provide, arguably, the most straightforward path towards practical, high-energy, and rechargeable batteries based on metallic anodes. 1 However, the sluggish ion transmission at the cathode-electrolyte (solid/solid) interface would result in the high resistant at the contact and limit the practical implementation of these all solid-state materials in real world batteries. 2 Several methods were suggested to enhance the kinetic condition of the ion migration between the cathode and the solid electrolyte (SE). 3 A composite strategy that mixes active materials and SEs for the cathode is a general way to decrease the ion transmission barrier at the cathode-electrolyte interface. 3 The active material concentration in the cathode is reduced as much as the SE portion increases by which the energy density of the ASSB is restricted. In addition, the mixing approach generally accompanies lattice mismatches between the cathode active materials and the SE, thus providing only limited improvements, which is imputed by random contacts between the cathode active materials and the SE during the mixing process. Implementing high-pressure for the electrode and electrolyte of ASSB in the assembling process has been verified is a but effective way to boost the ion transmission ability between the cathode active materials and the SE by decreasing the grain boundary impedance. Whereas the short-circuit of the battery would be induced by the mechanical deformation of the electrolyte under high pressure. 4 Herein, we demonstrate a novel way to address the ion transmission problem at the cathode-electrolyte interface in ASSBs. Starting from the cathode configuration, the finite element method (FEM) was employed to evaluate the current concentration and the distribution of the space charge layer at the cathode-electrolyte interface. Hierarchical three-dimensional (HTD) structures are found to have a higher Li + transfer number (t Li+ ), fewer free anions, and the weaker space-charge layer at the cathode-electrolyte interface in the resulting FEM simulation. To take advantage of the HTD structure, stereolithography is adopted as a manufacturing technique and single-crystalline Ni-rich (SCN) materials are selected as the active materials. Next, the manufactured HTD cathode is sintered at 600 °C in an N 2 atmosphere for the carbonization of the resin, which induces sufficient electronic conductivity for the cathode. Then, the gel-like Li 1.4 Al 0.4 Ti 1.6 (PO 4 ) 3 (LATP) precursor is synthesized and filled into the voids of the HTD structure cathode sufficiently. And the filled HTD structure cathodes are sintered at 900 °C to achieve the crystallization of the LATP gel. Scanning transmission electron microscopy (STEM) is used to unveil the morphology of the cathode-electrolyte interface between the sintered HTD cathode and the in-situ generated electrolyte (LATP). A transient phase has been found generated at the interface and matched with both lattices of the SCN and the SE, accelerating the transmission of the Li-ions, which is further verified by density functional theory calculations. In addition, Electron Energy Loss Spectroscopy demonstrates the preserved interface between HTD cathode and SEs. Atomic force microscopy is employed to measure the potential image of the cross-sectional interface by the peak force tapping mode. The average potential of modified samples is lower than the sample that mix SCN and SEs simply in the 2D planar structure, which confirms a weakened space charge layer by the enhanced contact capability as well as the ion transmission ability. To see if the demonstrated method is universally applicable, LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM811) is selected as the cathode active material and manufactured in the same way as the SCN. The HTD cathode based on NCM811 exhibits higher electrochemical performance compared with the reference sample based on the 2D planar mixing-type cathode. We believe such a demonstrated universal strategy provides a new guideline to engineer the cathode/electrolyte interface by revolutionizing electrode structures that can be applicable to all-solid-state batteries. Figure 1. Schematic of comparing of traditional 2D planar cathode and HTD cathode in ASSB Tikekar, M. D. , et al. , Nature Energy (2016) 1 (9), 16114 Banerjee, A. , et al. , Chem Rev (2020) 120 (14), 6878 Chen, R. , et al. , Chem Rev (2020) 120 (14), 6820 Cheng, X. , et al. , Advanced Energy Materials (2018) 8 (7) Figure 1 
    more » « less
  2. Helium atom scattering and density-functional theory (DFT) are used to characterize the phonon band structure of the (3 × 1)-O surface reconstruction of Nb(100). Innovative DFT calculations comparing surface phonons of bare Nb(100) to those of the oxide surface show increased resonances for the oxide, especially at higher energies. Calculated dispersion curves align well with experimental results and yield atomic displacements to characterize polarizations. Inelastic helium time-of-flight measurements show phonons with mixed longitudinal and shear-vertical displacements along both the ⟨[Formula: see text]⟩, [Formula: see text] and ⟨[Formula: see text]⟩, [Formula: see text] symmetry axes over the entire first surface Brillouin zone. Force constants calculated for bulk Nb, Nb(100), and the (3 × 1)-O Nb(100) reconstruction indicate much stronger responses from the oxide surface, particularly for the top few layers of niobium and oxygen atoms. Many of the strengthened bonds at the surface create the characteristic ladder structure, which passivates and stabilizes the surface. These results represent, to our knowledge, the first phonon dispersion data for the oxide surface and the first ab initio calculation of the oxide’s surface phonons. This study supplies critical information for the further development of advanced materials for superconducting radiofrequency cavities.

     
    more » « less
  3. The molecule-based ferrimagnetic semiconductor vanadium tetracyanoethylene (V[TCNE] x , x [Formula: see text] 2) has garnered interest from the quantum information community due to its excellent coherent magnonic properties and ease of on-chip integration. Despite these attractive properties, a detailed understanding of the electronic structure and mechanism for long-range magnetic ordering have remained elusive due to a lack of detailed atomic and electronic structural information. Previous studies via x-ray absorption near edge spectroscopy and the extended x-ray absorption fine structure have led to various proposed structures, and in general, V[TCNE] x is believed to be a three-dimensional network of octahedrally coordinated V 2+ , each bonded to six TCNE molecules. Here, we elucidate the electronic structure, structural ordering, and degradation pathways of V[TCNE] x films by correlating calculations of density functional theory (DFT) with scanning transmission electron microscopy and electron energy-loss spectroscopy (EELS) of V[TCNE] x films. Low-loss EELS measurements reveal a bandgap and an excited state structure that agree quantitatively with DFT modeling, including an energy splitting between apical and equatorial TCNE ligands within the structure, providing experimental results directly backed by theoretical descriptions of the electronic structure driving the robust magnetic ordering in these films. Core-loss EELS confirms the presence of octahedrally coordinated V +2 atoms. Upon oxidation, changes in the C1s- π* peak indicate that C=C of TCNE is preferentially attacked. Furthermore, we identify a relaxation of the structural ordering as the films age. These results lay the foundation for a more comprehensive and fundamental understanding of magnetic ordering and dynamics in these classes of metal–ligand compounds. 
    more » « less
  4. Polycrystalline materials consist of grains (crystals) oriented at different angles resulting in a heterogeneous and anisotropic mechanical behavior at that micro-length scale. In this study, a novel method is proposed for the first time to determine the [Formula: see text] crystal orientations of grains in a [Formula: see text] domain, using solely [Formula: see text] deformation fields. The grain boundaries are assumed to be unknown and delineated from the reconstructed changes in the crystallographic orientation. Further, the constitutive equations that describe the mechanical behavior of the domain in [Formula: see text] under plane stress conditions are derived, assuming that the material is transversely isotropic in 3D. Finite element based algorithms are utilized to discretize the inverse problem. The in-house written inverse problem solver is coupled with Matlab-based optimization scripts to solve for the mechanical property distributions. The performance of this method is tested at different noise levels with synthetic displacements that were used as measured data. The reconstructions deteriorate as the noise level is increased. This work presents a first milestone in the verification of this novel technology with synthetic data. 
    more » « less
  5. Dissipative particle dynamics (DPD) simulations are performed on coarse-grained replicas of linear, monodisperse entangled polyethylene melts [Formula: see text] and [Formula: see text] undergoing both steady-state and transient planar elongational flow (PEF). The fidelity of the DPD simulations is verified by direct comparison of flow topological and rheological properties of a 334-particle chain liquid against the united-atom [Formula: see text] liquid, simulated using nonequilibrium molecular dynamics (NEMD). These DPD simulations demonstrate that a flow-induced coil-stretch transition (CST) and its associated hysteresis caused by configurational microphase separation, as observed in previous NEMD simulations of PEF, can be replicated using a more computationally efficient coarse-grained system. Results indicate that the breadth of the CST hysteresis loop is enlarged for the longer molecule liquid relative to the shorter one. Furthermore, relaxation simulations reveal that reducing the applied flow Deborah number ([Formula: see text]) from a high value corresponding to a homogeneous phase of highly stretched molecules to a [Formula: see text] within the biphasic region results in a two-stage relaxation process. There is a fast initial stratification of the kinetically trapped highly stretched chains into regions of highly extended and less extended chains, which displays similar behavior to a system undergoing a spinodal decomposition caused by spatial configurational free energy fluctuations. After a short induction period of apparently random duration, the less extended chain regions experience a stochastic nucleation event that induces configurational relaxation to domains composed of coiled molecules over a much longer time scale, leaving the more highly extended chains in surrounding sheetlike domains. The time scales of these two relaxation processes are of the same order of magnitude as the Rouse and disengagement times of the equilibrium liquids.

     
    more » « less