BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Harold S. Park

Harold S. Park appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scicond-mat.mes-hallphysics.comp-phcond-mat.softComputational Engineering, Finance, and Sciencecond-mat.stat-mechMachine Learningphysics.data-an

Catalog footprint

What is connected

44works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

44 published item(s)

preprint2022arXiv

Intrinsic bending flexoelectric constants in two-dimensional materials

Flexoelectricity is a form of electromechanical coupling that has recently emerged because, unlike piezoelectricity, it is theoretically possible in any dielectric material. Two-dimensional (2D) materials have also garnered significant interest because of their unusual electromechanical properties and high flexibility, but the intrinsic flexoelectric properties of these materials remain unresolved. In this work, using atomistic modeling accounting for charge-dipole interactions, we report the intrinsic flexoelectric constants for a range of two-dimensional materials, including graphene allotropes, nitrides, graphene analogs of group-IV elements, and the transition metal dichalcogenides (TMDCs). We accomplish this through a proposed mechanical bending scheme that eliminates the piezoelectric contribution to the total polarization, which enables us to directly measure the flexoelectric constants. While flat 2D materials like graphene have low flexoelectric constants due to weak $π-σ$ interactions, buckling is found to increase the flexoelectric constants in monolayer group-IV elements. Finally, due to significantly enhanced charge transfer coupled with structural asymmetry due to bending, the TMDCs are found to have the largest flexoelectric constants, including MoS$_{2}$ having a flexoelectric constant ten times larger than graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Towards out of distribution generalization for problems in mechanics

There has been a massive increase in research interest towards applying data driven methods to problems in mechanics. While traditional machine learning (ML) methods have enabled many breakthroughs, they rely on the assumption that the training (observed) data and testing (unseen) data are independent and identically distributed (i.i.d). Thus, traditional ML approaches often break down when applied to real world mechanics problems with unknown test environments and data distribution shifts. In contrast, out-of-distribution (OOD) generalization assumes that the test data may shift (i.e., violate the i.i.d. assumption). To date, multiple methods have been proposed to improve the OOD generalization of ML methods. However, because of the lack of benchmark datasets for OOD regression problems, the efficiency of these OOD methods on regression problems, which dominate the mechanics field, remains unknown. To address this, we investigate the performance of OOD generalization methods for regression problems in mechanics. Specifically, we identify three OOD problems: covariate shift, mechanism shift, and sampling bias. For each problem, we create two benchmark examples that extend the Mechanical MNIST dataset collection, and we investigate the performance of popular OOD generalization methods on these mechanics-specific regression problems. Our numerical experiments show that in most cases, while the OOD generalization algorithms perform better compared to traditional ML methods on these OOD problems, there is a compelling need to develop more robust OOD generalization methods that are effective across multiple OOD scenarios. Overall, we expect that this study, as well as the associated open access benchmark datasets, will enable further development of OOD generalization methods for mechanics specific regression problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Double Dirac Cones and Topologically Non-Trivial Phonons for Continuous, Square Symmetric (C$_{4v}$ and C$_{2v}$) Unit Cells

Because phononic topological insulators have primarily been studied in discrete, graphene-like structures with C$_{6}$ or C$_{3}$ hexagonal symmetry, an open question is how to systematically achieve double Dirac cones and topologically non-trivial structures using continuous, non-hexagonal unit cells. Here, we address this challenge by presenting a novel computational methodology for the inverse design of continuous two-dimensional square phononic metamaterials exhibiting C$_{4v}$ and C$_{2v}$ symmetry. This leads to the systematic design of square unit cell topologies exhibiting a double Dirac degeneracy, which enables topologically-protected interface propagation based on the quantum spin Hall effect (QSHE). Numerical simulations prove that helical edge states emerge at the interface between two topologically distinct square phononic metamaterials, which opens the possibility of QSHE-based pseudospin-dependent transport beyond hexagonal lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The nonlinear buckling behavior of a complete spherical shell under uniform external pressure and homogenous natural curvature

In this work, we consider the stability of a spherical shell under combined loading from a uniform external pressure and a homogenous natural curvature. Non-mechanical stimuli, such as one that tends to modify the rest curvature of an elastic body, are prevalent in a wide range of natural and engineered systems, and may occur due to thermal expansion, changes in pH, differential swelling, and differential growth. Here, we investigate how the presence of both an evolving natural curvature and an external pressure modifies the stability of a complete spherical shell. We show that due to a mechanical analogy between pressure and curvature, positive natural curvatures can severely destabilize a thin shell, while negative natural curvatures can strengthen the shell against buckling, providing the possibility to design shells that buckle at or above the theoretical limit for pressure alone, i.e. a knock-up factor. These results extend directly from the classical analysis of the stability of shells under pressure, and highlight the important role that non-mechanical stimuli can have on modifying the the membrane state of stress in a thin shell.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Review on Mechanics and Mechanical Properties of 2D Materials - Graphene and Beyond

Since the first successful synthesis of graphene just over a decade ago, a variety of two-dimensional (2D) materials (e.g., transition metal-dichalcogenides, hexagonal boron-nitride, etc.) have been discovered. Among the many unique and attractive properties of 2D materials, mechanical properties play important roles in manufacturing, integration and performance for their potential applications. Mechanics is indispensable in the study of mechanical properties, both experimentally and theoretically. The coupling between the mechanical and other physical properties (thermal, electronic, optical) is also of great interest in exploring novel applications, where mechanics has to be combined with condensed matter physics to establish a scalable theoretical framework. Moreover, mechanical interactions between 2D materials and various substrate materials are essential for integrated device applications of 2D materials, for which the mechanics of interfaces (adhesion and friction) has to be developed for the 2D materials. Here we review recent theoretical and experimental works related to mechanics and mechanical properties of 2D materials. While graphene is the most studied 2D material to date, we expect continual growth of interest in the mechanics of other 2D materials beyond graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Auxetic Nanomaterials: Recent Progress and Future Development

Auxetic materials (materials with negative Poisson's ratio) and nanomaterials have independently been for many years two of the most active research fields in material science. Recently, these formerly independent fields have begun to intersect in new and interesting ways due to the recent discovery of auxeticity in nanomaterials like graphene, metal nanoplates, black phosphorus, and others. Here we review the research emerging at the intersection of auxeticity and nanomaterials. We first survey the atomistic mechanisms, both intrinsic and extrinsic, that have been found, primarily through atomistic simulations, to cause auxeticity in nanomaterials. We then outline the available experimental evidence for auxetic nanomaterials. In order to lay the groundwork for future work in this exciting area, we close by discussing several future prospects as well as the current challenges in this field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Intrinsic Negative Poisson's Ratio for Single-Layer Graphene

Negative Poisson's ratio (NPR) materials have drawn significant interest because the enhanced toughness, shear resistance and vibration absorption that typically are seen in auxetic materials may enable a range of novel applications. In this work, we report that single-layer graphene exhibits an intrinsic NPR, which is robust and independent of its size and temperature. The NPR arises due to the interplay between two intrinsic deformation pathways (one with positive Poisson's ratio, the other with NPR), which correspond to the bond stretching and angle bending interactions in graphene. We propose an energy-based deformation pathway criteria, which predicts that the pathway with NPR has lower energy and thus becomes the dominant deformation mode when graphene is stretched by a strain above 6%, resulting in the NPR phenomenon.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Negative Poisson's Ratio in Single-Layer Graphene Ribbons

The Poisson's ratio characterizes the resultant strain in the lateral direction for a material under longitudinal deformation. Though negative Poisson's ratios (NPR) are theoretically possible within continuum elasticity, they are most frequently observed in engineered materials and structures, as they are not intrinsic to many materials. In this work, we report NPR in single-layer graphene ribbons, which results from the compressive edge stress induced warping of the edges. The effect is robust, as the NPR is observed for graphene ribbons with widths smaller than about 10 nm, and for tensile strains smaller than about 0.5%, with NPR values reaching as large as -1.51. The NPR is explained analytically using an inclined plate model, which is able to predict the Poisson's ratio for graphene sheets of arbitrary size. The inclined plate model demonstrates that the NPR is governed by the interplay between the width (a bulk property), and the warping amplitude of the edge (an edge property), which eventually yields a phase diagram determining the sign of the Poisson's ratio as a function of the graphene geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Polarization and valley switching in monolayer group-IV monochalcogenides

Group-IV monochalcogenides are a family of two-dimensional puckered materials with an orthorhombic structure that is comprised of polar layers. In this article, we use first principles calculations to show the multistability of monolayer SnS and GeSe, two prototype materials where the direction of the puckering can be switched by application of tensile stress or electric field. Furthermore, the two inequivalent valleys in momentum space, which are dictated by the puckering orientation, can be excited selectively using linearly polarized light, and this provides an additional tool to identify the polarization direction. Our findings suggest that SnS and GeSe monolayers may have observable ferroelectricity and multistability, with potential applications in information storage.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Surface Tension Effects on Surface Instabilities of Dielectric Elastomers

Dielectric elastomers have recently been proposed for various biologically-relevant applications, in which they may operate in fluidic environments where surface tension effects may have a significant effect on their stability and reliability. Here, we present a theoretical analysis coupled with computational modeling for a generalized electromechanical analysis of surface stability in dielectric elastomers accounting for surface tension effects. For mechanically deformed elastomers, significant increases in critical strain and instability wavelength are observed for small elastocapillary numbers. When the elastomers are deformed electrostatically, both surface tension and the amount of pre-compression are found to substantially increase the critical electric field while decreasing the instability wavelength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Gaussian Treatment for the Friction Issue of Lennard-Jones Potential in Layered Materials: Application to Friction between Graphene, MoS2 and Black Phosphorus

The Lennard-Jones potential is widely used to describe the interlayer interactions within layered materials like graphene. However, it is also widely known that this potential strongly underestimates the frictional properties for layered materials. Here we propose to supplement the Lennard-Jones potential by a Gaussian-type potential, which enables more accurate calculations of the frictional properties of two-dimensional layered materials. Furthermore, the Gaussian potential is computationally simple as it introduces only one additional potential parameter that is determined by the interlayer shear mode in the layered structure. The resulting Lennard-Jones-Gaussian potential is applied to compute the interlayer cohesive energy and frictional energy for graphene, MoS2, black phosphorus, and their heterostructures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Review on Flexural Mode of Graphene: Lattice Dynamics, Thermal Conduction, Thermal Expansion, Elasticity, and Nanomechanical Resonance

Single-layer graphene is so flexible that its flexural mode (also called the ZA mode, bending mode, or out-of-plane transverse acoustic mode) is important for its thermal and mechanical properties. Accordingly, this review focuses on exploring the relationship between the flexural mode and thermal and mechanical properties of graphene. We first survey the lattice dynamic properties of the flexural mode, where the rigid translational and rotational invariances play a crucial role. After that, we outline contributions from the flexural mode in four different physical properties or phenomena of graphene -- its thermal conductivity, thermal expansion, Young's modulus, and nanomechanical resonance. We explain how graphene's superior thermal conductivity is mainly due to its three acoustic phonon modes at room temperature, including the flexural mode. Its coefficient of thermal expansion is negative in a wide temperature range resulting from the particular vibration morphology of the flexural mode. We then describe how the Young's modulus of graphene can be extracted from its thermal fluctuations, which are dominated by the flexural mode. Finally, we discuss the effects of the flexural mode on graphene nanomechanical resonators, while also discussing how the essential properties of the resonators, including mass sensitivity and quality factor, can be enhanced.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A Stillinger-Weber Potential for Single-Layer Black Phosphorus, and the Importance of Cross-Pucker Interactions for Negative Poisson's Ratio and Edge Stress-Induced Bending

The distinguishing structural feature of single-layer black phosphorus is its puckered structure, which leads to many novel physical properties. In this work, we first present a new parameterization of the Stillinger-Weber potential for single-layer black phosphorus. In doing so, we reveal the importance of a cross-pucker interaction term in capturing its unique mechanical properties, such as a negative Poisson's ratio. In particular, we show that the cross-pucker interaction enables the pucker to act as a re-entrant hinge, which expands in the lateral direction when it is stretched in the longitudinal direction. As a consequence, single-layer black phosphorus has a negative Poisson's ratio in the direction perpendicular to the atomic plane. As an additional demonstration of the impact of the cross-pucker interaction, we show that it is also the key factor that enables capturing of the edge stress-induced bending of single-layer black phosphorus that has been reported in ab initio calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Analytic Study of Strain Engineering the Electronic Bandgap in Single-Layer Black Phosphorus

We present an analytic study, based on the tight-binding approximation, of strain effects on the electronic bandgap in single-layer black phosphorus. We obtain an expression for the variation of the bandgap induced by a general strain type that includes both tension in and out of the plane and shear, and use this to determine the most efficient strain direction for different strain types, along which the strongest bandgap manipulation can be achieved. We find that the strain direction that enables the maximum manipulation of the bandgap is not necessarily in the armchair or zigzag direction. Instead, to achieve the strongest bandgap modulation, the direction of the applied mechanical strain is dependent on the type of applied strain.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational Modeling of Electro-Elasto-Capillary Phenomena in Dielectric Elastomers

We present a new finite deformation, dynamic finite element model that incorporates surface tension to capture elastocapillary effects on the electromechanical deformation of dielectric elastomers. We demonstrate the significant effect that surface tension can have on the deformation of dielectric elastomers through three numerical examples: (1) surface tension effects on the deformation of single finite elements with homogeneous and inhomogeneous boundary conditions; (2) surface tension effects on instabilities in constrained dielectric elastomer films, and (3) surface tension effects on bursting drops in solid dielectrics. Generally, we find that surface tension creates a barrier to instability nucleation. Specifically, we find in agreement with recent experimental studies of constrained dielectric elastomer films a transition in the surface instability mechanism depending on the elastocapillary length. The present results indicate that the proposed methodology may be beneficial in studying the electromechanical deformation and instabilities for dielectric elastomers in the presence of surface tension.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Conductance signatures of electron confinement induced by strained nanobubbles in graphene

We investigate the impact of strained nanobubbles on the conductance characteristics of graphene nanoribbons using a combined molecular dynamics - tight-binding simulation scheme. We describe in detail how the conductance, density of states, and current density of zigzag or armchair graphene nanoribbons are modified by the presence of a nanobubble. In particular, we establish that low-energy electrons can be confined in the vicinity or within the nanobubbles by the delicate interplay between the pseudomagnetic field pattern created by the shape of the bubble, mode mixing, and substrate interaction. The coupling between confined evanescent states and propagating modes can be enhanced under different clamping conditions, which translates into Fano resonances in the conductance traces.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coupling Tension and Shear for Highly Sensitive Graphene-Based Strain Sensors

We report, based on its variation in electronic transport to coupled tension and shear deformation, a highly sensitive graphene-based strain sensor consisting of an armchair graphene nanoribbon (AGNR) between metallic contacts. As the nominal strain at any direction increases from 2.5 to 10%, the conductance decreases, particularly when the system changes from the electrically neutral region. At finite bias voltage, both the raw conductance and the relative proportion of the conductance depends smoothly on the gate voltage with negligible fluctuations, which is in contrast to that of pristine graphene. Specifically, when the nominal strain is 10% and the angle varies from 0 degree to 90 degree, the relative proportion of the conductance changes from 60 to 90%.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Highly Stretchable MoS$_2$ Kirigami

We report the results of classical molecular dynamics simulations focused on studying the mechanical properties of MoS$_{2}$ kirigami. Several different kirigami structures were studied based upon two simple non-dimensional parameters, which are related to the density of cuts, as well as the ratio of the overlapping cut length to the nanoribbon length. Our key finding is significant enhancements in tensile yield (by a factor of four) and fracture strains (by a factor of six) as compared to pristine MoS$_{2}$ nanoribbons. These results in conjunction with recent results on graphene suggest that the kirigami approach may be a generally useful one for enhancing the ductility of two-dimensional nanomaterials.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mechanical Strain Effects on Black Phosphorus Nanoresonators

We perform classical molecular dynamics to investigate the effects of mechanical strain on single-layer black phosphorus nanoresonators at different temperatures. We find that the resonant frequency is highly anisotropic in black phosphorus due to its intrinsic puckered configuration, and that the quality factor in the armchair direction is higher than in the zigzag direction at room temperature. The quality factors are also found to be intrinsically larger than graphene and MoS2 nanoresonators. The quality factors can be increased by more than a factor of two by applying tensile strain, with uniaxial strain in the armchair direction being most effective. However, there is an upper bound for the quality factor increase due to nonlinear effects at large strains, after which the quality factor decreases. The tension induced nonlinear effect is stronger along the zigzag direction, resulting in a smaller maximum strain for quality factor enhancement.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Effects of Free Edge Interaction-Induced Knotting on the Buckling of Monolayer Graphene

Edge effects play an important role for many properties of graphene. While most works have focused on the effects from isolated free edges, we present a novel knotting phenomenon induced by the interactions between a pair of free edges in graphene, and investigate its effect on the buckling of monolayer graphene. Upon compression, the buckling of graphene starts gradually in the form of two buckling waves from the warped edges. The collision of these two buckling waves results in the creation of a knot structure in graphene. The knot structure enables the buckled graphene to exhibit two unique post-buckling characteristics. First, it induces a five-fold increase in graphene's mechanical stiffness during the buckling process. Second, the knotted structure enables graphene to exhibit a mechanically stable post-buckling regime over a large (3%) compressive strain regime, which is significantly larger than the critical buckling strain of about 0.5%. The combination of these two effects enables graphene to exhibit an unexpected post-buckling stability that has previously not been reported. We predict that numerical simulations or experiments should observe two distinct stress strain relations for the buckling of identical graphene samples, due to the characteristic randomness in the formation process of the knot structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

FPU physics with nanomechanical graphene resonators: intrinsic relaxation and thermalization from flexural mode coupling

Thermalization in nonlinear systems is a central concept in statistical mechanics and has been extensively studied theoretically since the seminal work of Fermi, Pasta and Ulam (FPU). Using molecular dynamics and continuum modeling of a ring-down setup, we show that thermalization due to nonlinear mode coupling intrinsically limits the quality factor of nanomechanical graphene drums and turns them into potential test beds for FPU physics. We find the thermalization rate $Γ$ to be independent of radius and scaling as $Γ\sim T^*/ε_{\rm pre}^2$, where $T^*$ and $ε_{\rm pre}$ are effective resonator temperature and prestrain.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Highly Deformable Graphene Kirigami

Graphene's exceptional mechanical properties, including its highest-known stiffness (1 TPa) and strength (100 GPa) have been exploited for various structural applications. However, graphene is also known to be quite brittle, with experimentally-measured tensile fracture strains that do not exceed a few percent. In this work, we introduce the notion of graphene kirigami, where concepts that have been used almost exclusively for macroscale structures are applied to dramatically enhance the stretchability of both zigzag and armchair graphene. Specifically, we show using classical molecular dynamics simulations that the yield and fracture strains of graphene can be enhanced by about a factor of three using kirigami as compared to standard monolayer graphene. This enhanced ductility in graphene should open up interesting opportunities not only mechanically, but also in coupling to graphene's electronic behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Interlayer Breathing and Shear Modes in Few-Layer Black Phosphorus

The interlayer breathing and shear modes in few-layer black phosphorus are investigated for their symmetry and lattice dynamical properties. The symmetry groups for the even-layer and odd-layer few-layer black phosphorus are utilized to determine the irreducible representation and the infrared and Raman activity for the interlayer modes. The valence force field model is applied to calculate the eigenvectorw and frequencies for the interlayer breathing and shear modes, which are explained using the atomic chain model. The anisotropic puckered configuration for black phosphorus leads to a highly anisotropic frequency for the two interlayer shear modes. More specifically, the frequency for the shear mode in the direction perpendicular to the pucker is less than half of the shear mode in the direction parallel with the pucker. We also report a set of interlayer modes having the same frequency for all few-layer black phosphorus with layer number N=3i with integer i, because of their collective vibrational displacements. The optical activity of the collective modes supports possible experimental identification for these modes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mechanical Properties of MoS2/Graphene Heterostructures

We perform classic molecular dynamics simulations to comparatively investigate the mechanical properties of single-layer MoS2 and a graphene/MoS2/graphene heterostructure under uniaxial tension. We show that the lattice mismatch between MoS2 and graphene will lead to an spontaneous strain energy in the interface. The Young's modulus of MoS2 can be enhanced by a factor of five by sandwiching it between two graphene layers. While the stiffness is enhanced, the yield strain of the MoS2 is reduced due to lateral buckling of the outer graphene layers due to the applied mechanical tension.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mechanical Properties of Single-Layer Black Phosphorus

The mechanical properties of single-layer black phosphrous under uniaxial deformation are investigated using first-principles calculations. Both Young's modulus and the ultimate strain are found to be highly anisotropic and nonlinear as a result of its quasi-two-dimensional puckered structure. Specifically, the in-plane Young's modulus is 44.0 GPa in the direction perpendicular to the pucker, and 92.7 GPa in the parallel direction. The ultimate strain is 0.48 and 0.20 in the perpendicular and parallel directions, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MoS2 Nanoresonators: Intrinsically Better Than Graphene?

We perform classical molecular dynamics simulations to examine the intrinsic energy dissipation in single-layer MoS$_{2}$ nanoresonators, where a point of emphasis is to compare its dissipation characteristics with those of single-layer graphene. Our key finding is that MoS$_{2}$ nanoresonators exhibit significantly lower energy dissipation, and thus higher quality (Q)-factors by at least a factor of four below room temperature, than graphene. Furthermore, this high Q-factor endows MoS$_{2}$ nanoresonators with a higher figure of merit, defined as frequency times Q-factor, despite a resonant frequency that is $50%$ smaller than graphene for the same size. By utilizing arguments from phonon-phonon scattering theory, we show that this reduced energy dissipation is enabled by the large energy gap in the phonon dispersion of MoS$_{2}$, which separates the acoustic phonon branches from the optical phonon branches, leading to a preserving mechanism for the resonant oscillation of MoS$_{2}$ nanoresonators. We further investigate the effects of tensile mechanical strain and nonlinear actuation on the Q-factors, where the tensile strain is found to counteract the reductions in Q-factor that occur with higher actuation amplitudes. Overall, our simulations illustrate the potential utility of MoS$_{2}$ for high frequency sensing and actuation applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Negative Poisson's Ratio in Single-Layer Black Phosphorus

The Poisson's ratio is a fundamental mechanical property that relates the resulting lateral strain to applied axial strain. While this value can theoretically be negative, it is positive for nearly all materials, though negative values have been observed in so-called auxetic structures. However, nearly all auxetic materials are bulk materials whose microstructure has been specifically engineered to generate a negative Poisson's ratio. In the present work, we report using first principles calculations the existence of a negative Poisson's ratio in a single-layer, two-dimensional material, black phosphorus. In contrast to engineered bulk auxetics, this behavior is intrinsic for single layer black phosphorus, and originates from its unique, puckered structure, where the pucker can be regarded as a re-entrant structure that is comprised of two coupled orthogonal hinges. As a result of this atomic structure, a negative Poisson's ratio is observed in the out-of-plane direction under uniaxial deformation in the direction parallel to the pucker, with the Poisson's ratio becoming increasingly negative with both increased tension and compression. The puckered structure also results in highly anisotropic in-plane Poisson's ratios, which are found to be 0.4 in the direction perpendicular and 1.28 in the direction parallel to the pucker.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Polar Surface Effects on the Thermal Conductivity in ZnO Nanowires: a Shell-Like Surface Reconstruction-Induced Preserving Mechanism

We perform molecular dynamics (MD) simulations to investigate the effect of polar surfaces on the thermal transport in zinc oxide (ZnO) nanowires. We find that the thermal conductivity in nanowires with free polar (0001) surfaces is much higher than in nanowires that have been stabilized with reduced charges on the polar (0001) surfaces, and also hexagonal nanowires without any transverse polar surfaces. From normal mode analysis, we show that the higher thermal conductivity is due to a shell-like reconstruction that occurs for the free polar surfaces. This shell-like reconstruction suppresses twisting motion in the nanowires such that the bending phonon modes are not scattered by the other phonon modes, and leads to substantially higher thermal conductivity in the ZnO nanowire with free polar surfaces. Furthermore, the auto-correlation function of the normal mode coordinate is utilized to extract the phonon lifetime, which leads to a concise explanation for the higher thermal conductivity in ZnO nanowires with free polar surfaces. Our work demonstrates that ZnO nanowires without polar surfaces, which exhibit low thermal conductivity, are more promising candidates for thermoelectric applications than nanowires with polar surfaces (and also high thermal conductivity).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Preserving the Q-Factors of ZnO Nanoresonators via Polar Surface Reconstruction

We perform molecular dynamics simulations to investigate the effect of polar surfaces on the quality (Q)-factors of zinc oxide (ZnO) nanowire-based nanoresonators. We find that the Q-factors in ZnO nanoresonators with free polar (0001) surfaces is about one order of magnitude higher than in nanoresonators that have been stabilized with reduced charges on the polar (0001) surfaces. From normal mode analysis, we show that the higher Q-factor is due to a shell-like reconstruction that occurs for the free polar surfaces. This shell-like reconstruction suppresses twisting motion in the nanowires such that the mixing of other modes with the resonant mode of oscillation is minimized, and leads to substantially higher Q-factors in the ZnO nanoresonators with free polar surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Pseudomagnetic fields in graphene nanobubbles of constrained geometry: A molecular dynamics study

Analysis of the strain-induced pseudomagnetic fields (PMFs) generated in graphene nanobulges under three different substrate scenarios shows that, in addition to the shape, the graphene-substrate interaction can crucially determine the overall distribution and magnitude of strain and those fields, in and outside the bulge. We utilize a combination of classical molecular dynamics, continuum mechanics, and tight-binding electronic structure calculations as an unbiased means of studying pressure-induced deformations and the resulting PMF in graphene nanobubbles of various geometries. The interplay among substrate aperture geometry, lattice orientation, internal gas pressure, and substrate type is analyzed in view of strain-engineered graphene nanostructures capable of confining and/or guiding electrons at low energies. Except in highly anisotropic geometries, the magnitude of the PMF is generally significant only near the boundaries of the aperture and rapidly decays towards the center because under gas pressure at the scales considered here there is considerable bending at the edges and the central region displays nearly isotropic strain. When the deflection lead to sharp bends at the edges, curvature and the tilting of the $p_z$ orbitals cannot be ignored and contributes substantially to the total field. The strong and localized nature of the PMF at the boundaries and its polarity-changing profile can be exploited to trap electrons inside the bubble or of guiding them in channel-like geometries defined by edges. However, we establish that slippage of graphene against the substrate is an important factor in determining the degree of concentration of PMFs in or around the bulge since it can lead to considerable softening of the strain gradients there. The nature of the substrate emerges thus as a decisive factor determining the effectiveness of nanoscale PMFs tailoring in graphene.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Strain-rate and temperature dependence of yield stress of amorphous solids via self-learning metabasin escape algorithm

A general self-learning metabasin escape (SLME) algorithm~\citep{caoPRE2012} is coupled in this work with continuous shear deformations to probe the yield stress as a function of strain rate and temperature for a binary Lennard-Jones (LJ) amorphous solid. The approach is shown to match the results of classical molecular dynamics (MD) at high strain rates where the MD results are valid, but, importantly, is able to access experimental strain rates that are about ten orders of magnitude slower than MD. In doing so, we find in agreement with previous experimental studies that a substantial decrease in yield stress is observed with decreasing strain rate. At room temperature and laboratory strain rates, the activation volume associated with yield is found to contain about 10 LJ particles, while the yield stress is as sensitive to a $1.5\%T_{\rm g}$ increase in temperature as it is to a one order of magnitude decrease in strain rate. Moreover, our SLME results suggest the SLME and extrapolated results from MD simulations follow distinctly different energetic pathways during the applied shear deformation at low temperatures and experimental strain rates, which implies that extrapolation of the governing deformation mechanisms from MD strain rates to experimental may not be valid.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Surface Shear Transformation Zones in Amorphous Solids

We perform a systematic study of the characteristics of shear transformation zones (STZs) that nucleate at free surfaces of two-dimensional amorphous solids subject to tensile loading using two different atomistic simulation methods, the standard athermal, quasistatic (AQ) approach and our recently developed self-learning metabasin escape (SLME) method to account for the finite temperature and strain-rate effects. In the AQ, or strain-driven limit, the nonaffine displacement fields of surface STZs decay exponentially away from their centers at similar decay rates as their bulk counterparts, though the direction of maximum nonaffine displacement is tilted away from the tensile axis due to surface effects. Using the SLME method at room temperature and at the high strain rates that are seen in classical molecular dynamics simulations, the characteristics for both bulk and surface STZs are found to be identical to those seen in the AQ simulations. However, using the SLME method at room temperature and experimentally-relevant strain rates, we find a transition in the surface STZ characteristics where a loss in the characteristic angular tensile-compression symmetry is observed. Finally, the thermally-activated surface STZs exhibit a slower decay rate in the nonaffine displacement field than do strain-driven surface STZs, which is characterized by a larger drop in potential energy resulting from STZ nucleation that is enabled by the relative compliance of the surface as compared to the bulk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Surface Stacking Fault Energy Approach to Predicting Defect Nucleation in Surface-Dominated Nanostructures

We present a surface stacking fault (SSF) energy approach to predicting defect nucleation from the surfaces of surface-dominated nanostructure such as FCC metal nanowires. The approach leads to a criteria that predicts the initial yield mechanism via either slip or twinning depending on whether the unstable twinning energy or unstable slip energy is smaller as determined from the resulting SSF energy curve. The approach is validated through a comparison between the SSF energy calculation and low-temperature classical molecular dynamics simulations of copper nanowires with different axial and transverse surface orientations, and cross sectional geometries. We focus on the effects of the geometric cross section by studying the transition from slip to twinning previously predicted in moving from a square to rectangular cross section for $\ <100\ > /\{100\}$ nanowires, and also for moving from a rhombic to truncated rhombic cross sectional geometry for $\ <110\ >$ nanowires. \hsp{We also provide the important demonstration that the criteria is able to predict the correct deformation mechanism when full dislocation slip is considered concurrently with partial dislocation slip and twinning. This is done in the context of rhombic aluminum nanowires which do not show a tensile reorientation due to full dislocation slip.} We show that the SSF energy criteria successfully predicts the initial mode of surface-nucleated plasticity at low temperature, while also discussing the effects of strain and temperature on the applicability of the criterion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Adsorbate Migration Effects on Continuous and Discontinuous Temperature-Dependent Transitions in the Quality Factors of Graphene Nanoresonators

We perform classical molecular dynamics simulation to investigate the mechanisms underpinning the unresolved, experimentally-observed temperature-dependent scaling transition in the quality factors of graphene nanomechanical resonators (GNMR). Our simulations reveal that the mechanism underlying this temperature scaling phenomenon is the out-of-plane migration of adsorbates on GNMRs. Specifically, the migrating adsorbate undergoes frequent collisions with the GNMR, which strongly influences the resulting mechanical oscillation, and thus the quality factors. We also predict a discontinuous transition in the quality factor at a lower critical temperature, which results from the in-plane migration of the adsorbate. Overall, our work clearly demonstrates the strong effect of adsorbate migration on the quality factors of GNMRs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Density functional theory calculation of edge stresses in monolayer MoS$_2$

We utilize density functional theory to calculate the edge energy and edge stress for monolayer MoS$_{2}$ nanoribbons. In contrast to previous reports for graphene, for both armchair and zigzag chiralities, the edge stresses for MoS$_{2}$ nanoribbons are found to be tensile, indicating that their lowest energy configuration is one of compression in which Mo-S bond lengths are shorter than those in a bulk, periodic MoS$_{2}$ monolayer. The edge energy and edge stress is found to converge for both chiralities for nanoribbon widths larger than about 1 nm.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Elastic Bending Modulus of Single-Layer Molybdenum Disulphide (MoS2): Finite Thickness Effect

We derive, from an empirical interaction potential, an analytic formula for the elastic bending modulus of single-layer MoS2 (SLMoS2). By using this approach, we do not need to define or estimate a thickness value for SLMoS2, which is important due to the substantial controversy in defining this value for two-dimensional or ultrathin nanostructures such as graphene and nanotubes. The obtained elastic bending modulus of 9.61 eV in SLMoS2 is significantly higher than the bending modulus of 1.4 eV in graphene, and is found to be within the range of values that are obtained using thin shell theory with experimentally obtained values for the elastic constants of SLMoS2. This increase in bending modulus as compared to monolayer graphene is attributed, through our analytic expression, to the finite thickness of SLMoS2. Specifically, while each monolayer of S atoms contributes 1.75 eV to the bending modulus, which is similar to the 1.4 eV bending modulus of monolayer graphene, the additional pairwise and angular interactions between out of plane Mo and S atoms contribute 5.84 eV to the bending modulus of SLMoS2.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Molecular Dynamics Simulations of Single-Layer Molybdenum Disulphide (MoS2): Stillinger-Weber Parametrization, Mechanical Properties, and Thermal Conductivity

We present a parameterization of the Stillinger-Weber potential to describe the interatomic interactions within single-layer MoS2 (SLMoS2). The potential parameters are fitted to an experimentally-obtained phonon spectrum, and the resulting empirical potential provides a good description for the energy gap and the crossover in the phonon spectrum. Using this potential, we perform classical molecular dynamics simulations to study chirality, size, and strain effects on the Young's modulus and the thermal conductivity of SLMoS2. We demonstrate the importance of the free edges on the mechanical and thermal properties of SLMoS2 nanoribbons. Specifically, while edge effects are found to reduce the Young's modulus of SLMoS2 nanoribbons, the free edges also reduce the thermal stability of SLMoS2 nanoribbons, which may induce melting well below the bulk melt temperature. Finally, uniaxial strain is found to efficiently manipulate the thermal conductivity of infinite, periodic SLMoS2.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Resonant tunneling in graphene pseudomagnetic quantum dots

Realistic relaxed configurations of triaxially strained graphene quantum dots are obtained from unbiased atomistic mechanical simulations. The local electronic structure and quantum transport characteristics of y-junctions based on such dots are studied, revealing that the quasi-uniform pseudomagnetic field induced by strain restricts transport to Landau level- and edge state-assisted resonant tunneling. Valley degeneracy is broken in the presence of an external field, allowing the selective filtering of the valley and chirality of the states assisting in the resonant tunneling. Asymmetric strain conditions can be explored to select the exit channel of the y-junction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Strain-rate and temperature-driven transition in the shear transformation zone for two-dimensional amorphous solids

We couple the recently developed self-learning metabasin escape algorithm, which enables efficient exploration of the potential energy surface (PES), with shear deformation to elucidate strain-rate and temperature effects on the shear transformation zone (STZ) characteristics in two-dimensional amorphous solids. In doing so, we report a transition in the STZ characteristics that can be obtained through either increasing the temperature or decreasing the strain rate. The transition separates regions having two distinct STZ characteristics. Specifically, at high temperatures and high strain rates, we show that the STZs have characteristics identical to those that emerge from purely strain-driven, athermal quasistatic atomistic calculations. At lower temperatures and experimentally relevant strain rates, we use the newly coupled PES + shear deformation method to show that the STZs have characteristics identical to those that emerge from a purely thermally activated state. The specific changes in STZ characteristics that occur in moving from the strain-driven to thermally activated STZ regime include a 33% increase in STZ size, faster spatial decay of the displacement field, a change in the deformation mechanism inside the STZ from shear to tension, a reduction in the stress needed to nucleate the first STZ, and finally a notable loss in characteristic quadrupolar symmetry of the surrounding elastic matrix that has previously been seen in athermal, quasistatic shear studies of STZs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Beat Phenomena in Metal Nanowires, and their Implications for Resonance-Based Elastic Property Measurements

The elastic properties of 1D nanostructures such as nanowires are often measured experimentally through actuation of the nanowire at its resonance frequency, and then relating the resonance frequency to the elastic stiffness using elementary beam theory. In the present work, we utilize large scale molecular dynamics simulations to report a novel beat phenomenon in [110] oriented Ag nanowires. The beat phenomenon is found to arise from the asymmetry of the lattice spacing in the orthogonal elementary directions of the [110] nanowire,i.e., [-110] and [001] directions, which results in two different principal moments of inertia. Because of this, actuations imposed along any other direction are found to decompose into two orthogonal vibrational components based on the actuation angle relative to these two elementary directions, with this phenomenon being generalizable to <110> FCC nanowires of different materials (Cu, Au, Ni, Pd and Pt).The beat phenomenon is explained using a discrete moment of inertia model based on the hard sphere assumption, the model is utilized to show that surface effects enhance the beat phenomenon, while the effect is reduced with increasing nanowires cross-sectional size or aspect ratio.Most importantly, due to the existence of the beat phenomena, we demonstrate that in resonance experiments only a single frequency component is expected to be observed, particularly when the damping ratio is relatively large or very small.Furthermore, for a large range of actuation angles, the lower frequency is more likely to be detected than the higher one,which implies the experimental predictions of Young's modulus obtained from resonance may in fact be under predictions. The present study therefore has significant implications for experimental interpretations of Young's modulus as obtained via resonance testing.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Enhancing the Mass Sensitivity of Graphene Nanoresonators Via Nonlinear Oscillations: The Effective Strain Mechanism

We perform classical molecular dynamics simulations to investigate the enhancement of the mass sensitivity and resonant frequency of graphene nanomechanical resonators that is achieved by driving them into the nonlinear oscillation regime. The mass sensitivity as measured by the resonant frequency shift is found to triple if the actuation energy is about 2.5 times the initial kinetic energy of the nanoresonator. The mechanism underlying the enhanced mass sensitivity is found to be the effective strain that is induced in the nanoresonator due to the nonlinear oscillations, where we obtain an analytic relationship between the induced effective strain and the actuation energy that is applied to the graphene nanoresonator. An important implication of this work is that there is no need for experimentalists to apply tensile strain to the resonators before actuation in order to enhance the mass sensitivity. Instead, enhanced mass sensitivity can be obtained by the far simpler technique of actuating nonlinear oscillations of an existing graphene nanoresonator.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Intrinsic Energy Dissipation in CVD-Grown Graphene Nanoresonators

We utilize classical molecular dynamics to study the the quality (Q)-factors of monolayer CVD-grown graphene nanoresonators. In particular, we focus on the effects of intrinsic grain boundaries of different orientations, which result from the CVD growth process, on the Q-factors. For a range of misorientations orientation angles that are consistent with those seen experimentally in CVD-grown graphene, i.e. 0$^{\circ}$ to $\sim20^{\circ}$, we find that the Q-factors for graphene with intrinsic grain boundaries are 1-2 orders of magnitude smaller than that of pristine monolayer graphene. We find that the Q-factor degradation is strongly influenced by both the symmetry and structure of the 5-7 defect pairs that occur at the grain boundary. Because of this, we also demonstrate that find the Q-factors CVD-grown graphene can be significantly elevated, and approach that of pristine graphene, through application of modest (1%) tensile strain.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Surface Effects on the Piezoelectricity of ZnO Nanowires

We utilize classical molecular dynamics to study surface effects on the piezoelectric properties of ZnO nanowires as calculated under uniaxial loading. An important point to our work is that we have utilized two types of surface treatments, those of charge compensation and surface passivation, to eliminate the polarization divergence that otherwise occurs due to the polar (0001) surfaces of ZnO. In doing so, we find that if appropriate surface treatments are utilized, the elastic modulus and the piezoelectric properties for ZnO nanowires having a variety of axial and surface orientations are all reduced as compared to the bulk value as a result of polarization-reduction in the polar [0001] direction. The reduction in effective piezoelectric constant is found to be independent of the expansion or contraction of the polar (0001) surface in response to surface stresses. Instead, the surface polarization and thus effective piezoelectric constant is substantially reduced due to a reduction in the bond length of the Zn-O dimer closest to the polar (0001) surface. Furthermore, depending on the nanowire axial orientation, we find in the absence of surface treatment that the piezoelectric properties of ZnO are either effectively lost due to unphysical transformations from the wurtzite to non-piezoelectric d-BCT phases, or also become smaller with decreasing nanowire size. The overall implication of this study is that if enhancement of the piezoelectric properties of ZnO is desired, then continued miniaturization of square or nearly square cross section ZnO wires to the nanometer scale is not likely to achieve this result.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nanomechanical Resonators and Their Applications in Biological/Chemical Detection: Nanomechanics Principles

Recent advances in nanotechnology have led to the development of nano-electro-mechanical systems (NEMS) such as nanomechanical resonators, which have recently received significant attention from the scientific community. This has not only been for their capability for the label-free detection of bio/chemical-molecules at single-molecule (or atomic) resolution for future applications such as the early diagnostics of diseases such as cancer, but also for their unprecedented ability to detect physical quantities such as molecular weight, elastic stiffness, surface stress, and surface elastic stiffness for adsorbed molecules on the surface. Most experimental works on resonator-based molecular detection have been based on the principle that molecular adsorption onto a resonator surface increases the effective mass, and consequently decreases the resonant frequencies of the nanomechanical resonator. However, this principle is insufficient to provide fundamental insights into resonator-based molecular detection at the nanoscale; this is due to recently proposed novel nanoscale detection principles including various effects such as surface effects, nonlinear oscillations, coupled resonance, and stiffness effects. Therefore, our objective in this review is to overview the current attempts to understand the underlying mechanisms in nanoresonator-based detection using physical models coupled to computational simulations and/or experiments. Specifically, we will focus on issues of special relevance to the dynamic behavior of nanoresonators and their applications in biological/chemical detection. We additionally provide extensive discussion regarding potentially fruitful future research directions coupling experiments and simulations in order to develop a fundamental understanding of the basic physical principles that govern NEMS and NEMS-based sensing applications.

0 citations0 reviewsNo code linkedOpen access