skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: ATOMIC-TO-CONTINUUM MULTISCALE MODELING OF DEFECTS IN CRYSTALS WITH NONLOCAL ELECTROSTATIC INTERACTIONS
This work develops a multiscale modeling framework for defects in crystals with general geometries and boundary conditions in which ionic interactions are important, with potential application to, e.g., ionic solids and electric field interactions with materials. The overall strategy is posed in the framework of the Quasicontinuum multiscale method; specifically, the use of a finite-element inspired kinematic description enables a significant reduction in the large number of degrees of freedom to describe the atomic positions. The key advance of this work is a method for the efficient and accurate treatment of nonlocal electrostatic charge-charge interactions without restrictions on the geometry or boundary conditions. Electrostatic interactions are long-range with slow decay, and hence require consideration of all pairs of charges making a brute-force approach computationally prohibitive. The method proposed here accounts for the exact charge-charge interactions in the near-field and uses a coarse-grained approximation in the far-field. The coarse-grained approximation and the associated errors are rigorously derived based on the limit of a finite body with a small periodic lengthscale, thereby enabling the errors in the approximation to be controlled to a desired tolerance. The method is applied to a simple model of Gallium Nitride and it is shown that electrostatic interactions can be approximated with a desired level of accuracy using the proposed methodology.  more » « less
Award ID(s):
1921857
PAR ID:
10377851
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Applied Mechanics
ISSN:
0021-8936
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; (, Molecular Systems Design & Engineering)
    null (Ed.)
    The functions of colloids, such as membranes and vesicles, are dictated by interfacial properties which are determined by an interplay of physical interactions and processes spanning multiple spatiotemporal scales. The multiscale characteristics of membranes and vesicles can be resolved by the hybrid molecular dynamics-lattice Boltzmann technique. This technique enables the resolution of particle dynamics along with long range electrostatic and hydrodynamic interactions. We have examined the feasibility of an implementation of the hybrid technique in conjunction with a Martini-based implicit solvent coarse-grained force field to capture the molecular and interfacial characteristics of membranes and vesicles. For simplicity, we have examined two types of vesicles whose molecular components have different sustained interactions with the solvent. One of the vesicles encompassed phospholipids and the other vesicle was composed of phospholipids and poly ethylene glycol (PEG)-grafted lipids. The molecular and interfacial characteristics of the phospholipid vesicle and PEGylated, or hairy, vesicles are found to be in good agreement with earlier experimental, computational and theoretical findings. These results demonstrate that the multiscale hybrid technique used with a Martini-based implicit solvent coarse-grained model is suitable for capturing the molecular and interfacial characteristics of membranes and vesicles. Furthermore, other implicit solvent coarse-grained models can be used in conjunction with the hybrid molecular dynamics-lattice Boltzmann technique to examine the molecular and interfacial characteristics of membranes and vesicles. Our study demonstrates the potential of the hybrid technique in capturing multiscale interfacial characteristics of intra- and inter-colloid interactions in suspensions under different flow conditions, and their relation to molecular properties. In addition, this technique can be applied to design colloids with multiscale characteristics which yield desired interactions with other colloids and responses to external stimuli. 
    more » « less
  2. ; ; ; ; (, Mathematics)
    In this paper, we consider unsaturated filtration in heterogeneous porous media with rough surface topography. The surface topography plays an important role in determining the flow process and includes multiscale features. The mathematical model is based on the Richards’ equation with three different types of boundary conditions on the surface: Dirichlet, Neumann, and Robin boundary conditions. For coarse-grid discretization, the Generalized Multiscale Finite Element Method (GMsFEM) is used. Multiscale basis functions that incorporate small scale heterogeneities into the basis functions are constructed. To treat rough boundaries, we construct additional basis functions to take into account the influence of boundary conditions on rough surfaces. We present numerical results for two-dimensional and three-dimensional model problems. To verify the obtained results, we calculate relative errors between the multiscale and reference (fine-grid) solutions for different numbers of multiscale basis functions. We obtain a good agreement between fine-grid and coarse-grid solutions. 
    more » « less
  3. ; ; ; (, Mathematics and Mechanics of Solids)
    Two-dimensional (2D) electronic materials are of significant technological interest due to their exceptional properties and broad applicability in engineering. The transition from nanoscale physics, which dictates their stable configurations, to macroscopic engineering applications requires the use of multiscale methods to systematically capture their electronic properties at larger scales. A key challenge in coarse-graining is the rapid and near-periodic variation of the charge density, which exhibits significant spatial oscillations at the atomic scale. Therefore, the polarization density field—the first moment of the charge density over the periodic unit cell—is used as a multiscale mediator that enables efficient coarse-graining by exploiting the almost-periodic nature of the variation. Unlike the highly oscillatory charge density, the polarization varies over lengthscales that are much larger than the atomic, making it suitable for continuum modeling. In this paper, we investigate the electrostatic potential arising from the charge distribution of arbitrarily-deformed 2D materials. Specifically, we consider a sequence of problems wherein the underlying lattice spacing vanishes and thus obtain the continuum limit. We consider three distinct limits: where the thickness is much smaller than, comparable to, and much larger than the in-plane lattice spacing. These limiting procedures provide the homogenized potential expressed in terms of the boundary charge and dipole distribution, subject to the appropriate boundary conditions that are also obtained through the limit process. Furthermore, we demonstrate that despite the intrinsic non-uniqueness in the definition of polarization, accounting for the boundary charges ensures that the total electrostatic potential, the associated electric field, and the corresponding energy of the homogenized system are uniquely determined. 
    more » « less
  4. ; ; (, Macromolecules)
    Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separations technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI forcefield. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N=150 monomers; MW = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm2. Using 5 µs simulations, we investigate the effects of grafting density on PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. The interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm2. 
    more » « less
  5. ; ; (, Journal of Computational Physics)
    In this article, we introduce a new method for discretizing micro-macro models of dilute polymeric fluids by integrating a finite element method (FEM) discretization for the macroscopic fluid dynamics equation with a deterministic variational particle scheme for the microscopic Fokker-Planck equation. To address challenges arising from micro-macro coupling, we employ a discrete energetic variational approach to derive a coarse-grained micro-macro model with a particle approximation first and then develop a particle-FEM discretization for the coarse-grained model. The accuracy of the proposed method is evaluated for a Hookean dumbbell model in a Couette flow by comparing the computed velocity field with existing analytical solutions. We also use our method to study nonlinear FENE dumbbell models in different scenarios, such as extensional flow, pure shear flow, and lid-driven cavity flow. Numerical examples demonstrate that the proposed deterministic particle approach can accurately capture the various key rheological phenomena in the original FENE model, including hysteresis and δ-function-like spike behavior in extensional flows, velocity overshoot phenomenon in pure shear flows, symmetries breaking, vortex center shifting, and vortices weakening in lid-driven cavity flows, with a small number of particles. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email