skip to main content

Explore Research Products in the NSF-PAR

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
NSF-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. ; ; ( , AIAA SCITECH 2023 Forum)
    In this manuscript, we present a multiscale Adaptive Reduced-Order Modeling (AROM) framework to efficiently simulate the response of heterogeneous composite microstructures under interfacial and volumetric damage. This framework builds on the eigendeformation-based reduced-order homogenization model (EHM), which is based on the transformation field analysis (TFA) and operates in the context of computational homogenization with a focus on model order reduction of the microscale problem. EHM pre-computes certain microstructure information by solving a series of linear elastic problems defined over the fully resolved microstructure (i.e., concentration tensors, interaction tensors) and approximates the microscale problem using a much smaller basis spanned over subdomains (also called parts) of the microstructure. Using this reduced basis, and prescribed spatial variation of inelastic response fields over the parts, the microscale problem leads to a set of algebraic equations with part-wise responses as unknowns, instead of node-wise displacements as in finite element analysis. The volumetric and interfacial influence functions are calculated by using the Interface enriched Generalized Finite Element Method (IGFEM) to compute the coefficient tensors, in which the finite element discretization does not need to conform to the material interfaces. AROM takes advantage of pre-computed coefficient tensors associated with the finest ROM and efficiently computes the coefficient tensors of a series of gradually coarsening ROMs. During the multiscale analysis stage, the simulation starts with a coarse ROM which can capture the initial elastic response well. As the loading continues and response in certain parts of the microstructure starts to localize, the analysis adaptively switches to the next level of refined ROM to better capture those local responses. The performance of AROM is evaluated by comparing the results with regular EHM (no adaptive refinement) and IGFEM under different loading conditions and failure modes for various 2D and 3D microstructures. The proposed AROM provides an efficient way to model history-dependent nonlinear responses for composite materials under localized interface failure and phase damage. 
    more » « less
  4. ; ; ; ; ( , Protein Science)
    Abstract

    Colloidal protein–protein interactions (PPI) are often expected to impact key behaviors of proteins in solution, such as aggregation rates and mechanisms, aggregate structure, protein solubility, and solution viscosity. PPI of an anti‐fluorescein single chain antibody variable fragment (scFv) were characterized experimentally at low to intermediate ionic strength using a combination of static light scattering and sedimentation equilibrium ultracentrifugation. Surprisingly, the results indicated that interactions were strongly net‐attractive and electrostatics promoted self‐association. Only repulsive interactions were expected based on prior work and calculations based a homology model of a related scFv crystal structure. However, the crystal structure lacks the charged, net‐neutral linker sequence. PyRosetta was used to generate a set of scFv structures with different linker conformations, and coarse‐grained Monte Carlo simulations were used to evaluate the effect of different linker configurationsviasecond osmotic virial coefficient (B22) simulations. The results show that the configuration of the linker has a significant effect on the calculatedB22values, and can result in strong electrostatic attractions between oppositely charged residues on the protein surface. This is particularly relevant for development of non‐natural antibody products, where charged linkers and other loop regions may be prevalent. The results also provide a preliminary computational framework to evaluate the effect of unstructured linkers on experimental protein–protein interaction parameters such asB22.

     
    more » « less
  5. ; ; ( , Scientific Reports)
    Abstract

    A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email