检索结果-内蒙古大学图书馆

Isotope effects in molecular structures and electronic properties of liquid water via deep potential molecular dynamics based on the SCAN functional

学校读者我要写书评

暂无评论

Physical Review B 2020年第21期102卷 214113-214113页

作者： Jianhang Xu Chunyi Zhang Linfeng Zhang Mohan Chen Biswajit Santra Xifan Wu Department of Physics Temple University Philadelphia Pennsylvania 19122 USA Program in Applied and Computational Mathematics Princeton University Princeton New Jersey 08544 USA CAPT HEDPS College of Engineering Peking University Beijing 100871 China Institute for Computational Molecular Science Temple University Philadelphia Pennsylvania 19122 USA

Feynman path-integral deep potential molecular dynamics (PI-DPMD) calculations have been employed to study both light (H2O) and heavy water (D2O) within the isothermalisobaric ensemble. In particular, the deep neural network is trained based on ab initio data obtained from the strongly constrained and appropriately normed (SCAN) exchange-correlation functional. Because of the lighter mass of hydrogen than deuteron, the properties of light water are more influenced by nuclear quantum effect than those of heavy water. Clear isotope effects are observed and analyzed in terms of hydrogen-bond structure and electronic properties of water that are closely associated with experimental observables. The molecular structures of both liquid H2O and D2O agree well with the data extracted from scattering experiments. The delicate isotope effects on radial distribution functions and angular distribution functions are well reproduced as well. Our approach demonstrates that deep neural network combined with SCAN functional based ab initio molecular dynamics provides an accurate theoretical tool for modeling water and its isotope effects.

关键词： Atomic & molecular structure Electronic structure Isotope effect Water Ab initio calculations Density functional theory Machine learning Molecular dynamics Path-integral methods

Interior electroneutrality in nernst-planck-navier-stokes systems

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Constantin, Peter Ignatova, Mihaela Lee, Fizay-Noah DEPARTMENT OF MATHEMATICS PRINCETON UNIVERSITY PRINCETONNJ08544 United States DEPARTMENT OF MATHEMATICS TEMPLE UNIVERSITY PHILADELPHIAPA19122 United States Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

We consider the limit of vanishing Debye length for ionic diffusion in fluids, described by the Nernst-Planck-Navier-Stokes system. In the asymptotically stable cases of blocking (vanishing normal flux) and uniform selective (special Dirichlet) boundary conditions for the ionic concentrations, we prove that the ionic charge density ρ converges in time to zero in the interior of the domain, in the limit of vanishing Debye length (ǫ → 0). For the unstable regime of Dirichlet boundary conditions for the ionic concentrations, we prove bounds that are uniform in time and ǫ. We also consider electroneutral boundary conditions, for which we prove that electroneutrality ρ → 0 is achieved at any fixed ǫ > 0, exponentially fast in time in Lp, for all 1 ≤ p MSC Codes 35Q30, 35Q35, 35Q92 Copyright © 2020, The Authors. All rights reserved.

关键词： Boundary conditions

Threshold and decay properties of transient isolated turbulent band in plane Couette flow

学校读者我要写书评

暂无评论

applied mathematics and Mechanics(English Edition) 2019年第10期40卷 1449-1456页

作者： Jianzhou LU Jianjun TAO Weitao ZHOU Xiangming XIONG CAPT-HEDPS SKLTCS Collaborative Innovation Center of IFSA Department of Mechanics and Engineering Science College of Engineering Peking University Beijing 100871 China The Program in Applied and Computational Mathematics Princeton University Princeton NJ 08544 U. S.A.

By direct numerical simulations of the plane Couette flow (PCF) in a large computational domain, it is shown that an isolated turbulent band decays monotonically at low Reynolds numbers but experiences transient growth before the eventual relaminarization at moderate Reynolds numbers. The lower bound Reynolds number of the transient-growth regime is determined as 286. The width, length, and tilt angle of the iso- lated band structure are defined based on the disturbance kinetic energy in the mid-plane, and the geometric characteristics of the band can be described with a tilted rectangle. It is illustrated that before its eventual fragmentation, the isolated turbulent band decays in a style of longitudinal contraction, where the center, width, and tilt angle of the band keep almost constant but the band length contracts with a statistically constant velocity.

关键词： plane Couette flow (PCF) subcritical transition turbulent band

Nernst-Planck-Navier-stokes systems far from equilibrium

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Constantin, Peter Ignatova, Mihaela Lee, Fizay-Noah Department of Mathematics Princeton University PrincetonNJ08544 United States Department of Mathematics Temple University PhiladelphiaPA19122 United States Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

We consider ionic electrodiffusion in fluids, described by the Nernst-Planck-Navier-Stokes system. We prove that the system has global smooth solutions for arbitrary smooth data: aribitrary positive Dirichlet boundary conditions for the ionic concentrations, arbitrary Dirichlet boundary conditions for the potential, arbitrary positive initial concentrations, and arbitrary regular divergence-free initial velocities. The result holds for any positive diffusivities of ions, in bounded domains with smooth boundary in three space dimensions, in the case of two ionic species, coupled to Stokes equations for the fluid. The result also holds in the case of Navier-Stokes coupling, if the velocity is regular. The global smoothness of solutions is also true for arbitrarily many ionic species, if all their diffusivities are the *** Codes 35Q30, 35Q35, 35Q92 Copyright © 2020, The Authors. All rights reserved.

关键词： Navier Stokes equations

A hybrid semi-Lagrangian cut cell method for advection-diffusion problems with robin boundary conditions in moving domains

学校读者我要写书评

暂无评论

arXiv 2021年

作者： Barrett, Aaron Fogelson, Aaron L. Griffith, Boyce E. Department of Mathematics University of Utah Salt Lake CityUT United States Departments of Mathematics and Bioengineering University of Utah Salt Lake CityUT United States Departments of Mathematics Applied Physical Sciences and Biomedical Engineering University of North Carolina Chapel HillNC United States Carolina Center for Interdisciplinary Applied Mathematics University of North Carolina Chapel HillNC United States Computational Medicine Program University of North Carolina Chapel HillNC United States McAllister Heart Institute University of North Carolina Chapel HillNC United States

We present a new discretization approach to advection-diffusion problems with Robin boundary conditions on complex, time-dependent domains. The method is based on second order cut cell finite volume methods introduced by Bochkov et al. [8] to discretize the Laplace operator and Robin boundary condition. To overcome the small cell problem, we use a splitting scheme along with a semi-Lagrangian method to treat advection. We demonstrate second order accuracy in the L1, L2, and L∞norms for both analytic test problems and numerical convergence studies. We also demonstrate the ability of the scheme to convert one chemical species to another across a moving boundary. © 2021, CC BY-SA.

关键词： Advection

NERNST-PLANCK-NAVIER-STOKES SYSTEMS NEAR EQUILIBRIUM

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Constantin, Peter Ignatova, Mihaela Lee, Fizay-Noah DEPARTMENT OF MATHEMATICS PRINCETON UNIVERSITY PRINCETONNJ08544 United States DEPARTMENT OF MATHEMATICS TEMPLE UNIVERSITY PHILADELPHIAPA19122 United States PROGRAM IN APPLIED AND COMPUTATIONAL MATHEMATICS PRINCETON UNIVERSITY PRINCETONNJ08544 United States

The Nernst-Planck-Navier-Stokes system models electrodiffusion of ions in a fluid. We prove global existence of solutions in bounded domains in three dimensions with either blocking (no-flux) or uniform selective (special Dirichlet) boundary conditions for ion concentrations. The global existence of strong solutions is established for initial conditions that are sufficiently small perturbations of steady state solutions. The solutions remain close to equilbrium in strong norms. The main two steps of the proof are (1) the decay of the sum of relative entropies (Kullback-Leibler divergences) and (2) the control of L2 norms of deviations by the sum of relative *** Codes 35Q30, 35Q35, 35Q92 Copyright © 2020, The Authors. All rights reserved.

关键词： Navier Stokes equations

computational Design of Anisotropic Stealthy Hyperuniform Composites with Engineered Directional Scattering Properties

学校读者我要写书评

暂无评论

arXiv 2023年

作者： Shi, Wenlong Keeney, David Chen, Duyu Jiao, Yang Torquato, Salvatore Materials Science and Engineering Arizona State University TempeAZ85287 United States Materials Research Laboratory University of California Santa BarbaraCA93106 United States Department of Physics Arizona State University TempeAZ85287 United States Department of Chemistry Princeton University PrincetonNJ08544 United States Department of Physics Princeton University PrincetonNJ08544 United States Princeton Institute of Materials Princeton University PrincetonNJ08544 United States Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States

Disordered hyperuniform materials are an emerging class of exotic amorphous states of matter that endow them with singular physical properties, including large isotropic photonic band gaps, superior resistance to fracture, and nearly optimal electrical and thermal transport properties, to name but a few. Here, we generalize the Fourier-space based numerical construction procedure for designing and generating digital realizations of isotropic disordered hyperuniform two-phase heterogeneous materials (i.e., composites) developed by Chen and Torquato [Acta Mater. 142, 152 (2018)] to anisotropic microstructures with targeted spectral densities. Our generalized construction procedure explicitly incorporates the vector-dependent spectral density function χ˜V (k) of arbitrary form that is realizable. We demonstrate the utility of the procedure by generating a wide spectrum of anisotropic stealthy hyperuniform (SHU) microstructures with χ˜V (k) = 0 for k ∈ Ω, i.e., complete suppression of scattering in an "exclusion" region Ω around the origin in the Fourier space. We show how different exclusion-region shapes with various discrete symmetries, including circular-disk, elliptical-disk, square, rectangular, butterfly-shaped and lemniscate-shaped regions of varying size, affect the resulting statistically anisotropic microstructures as a function of the and phase volume fraction. The latter two cases of Ω lead to directionally hyperuniform composites, which are stealthy hyperuniform only along certain directions, and are non-hyperuniform along others. We find that, while the circular-disk exclusion regions give rise to isotropic hyperuniform composite microstructures, the directional hyperuniform behaviors imposed by the shape asymmetry (or anisotropy) of certain exclusion regions give rise to distinct anisotropic structures and degree of uniformity in the distribution of the phases on intermediate and large length scales along different directions. Moreover, while the anisotr

关键词： Volume fraction

Deep Potential generation scheme and simulation protocol for the Li10GeP2S12-type superionic conductors

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Huang, Jianxing Zhang, Linfeng Wang, Han Zhao, Jinbao Cheng, Jun Weinan, E. State Key Laboratory of Physical Chemistry of Solid Surfaces iChEM College of Chemistry and Chemical Engineering Xiamen University Xiamen361005 China Program in Applied and Computational Mathematics Princeton University PrincetonNJ08544 United States Laboratory of Computational Physics Institute of Applied Physics and Computational Mathematics Fenghao East Road 2 Beijing100094 China Department of Mathematics Princeton University PrincetonNJ08544 United States

Solid-state electrolyte materials with superior lithium ionic conductivities are vital to the next-generation Li-ion batteries. Molecular dynamics could provide atomic scale information to understand the diffusion process of Li-ion in these superionic conductor materials. Here we implement the Deep Potential Generator (DP-GEN) to set up an efficient protocol to automatically generate interatomic potentials for Li10GeP2S12-type solid-state electrolyte materials (Li10GeP2S12, Li10SiP2S12 and Li10SnP2S12). The reliability and accuracy of the fast interatomic potentials are validated. With the potentials, we extend the simulation of diffusion process to a wide temperature range (300K-1000K) and systems with large size (~1000 atoms). Important technical aspects like the statistical error and size effect are carefully investigated and benchmark tests including the effect of density functional, thermal expansion and configurational disorder are performed. The computed data that considering these factors agrees well with the experimental results and we find that the 3 structures shows different behaviors with respect to configurational disorder. Our work paves the way for further research on computation screening of solid-state electrolyte materials. Copyright © 2020, The Authors. All rights reserved.

关键词： Solid electrolytes

Emergent spaces for coupled oscillators

学校读者我要写书评

暂无评论

arXiv 2020年

作者： Thiem, Thomas N. Kooshkbaghi, Mahdi Bertalan, Tom Laing, Carlo R. Kevrekidis, Ioannis G. Department of Chemical and Biological Engineering Princeton University United States Program in Applied and Computational Mathematics Princeton University United States Department of Mechanical Engineering Massachusetts Institute of Technology United States School of Natural and Computational Sciences Massey University New Zealand Department of Applied Mathematics and Statistics Johns Hopkins University United States

Systems of coupled dynamical units (e.g. oscillators or neurons) are known to exhibit complex, emergent behaviors that may be simplified through coarse-graining: a process in which one discovers coarse variables and derives equations for their evolution. Such coarse-graining procedures often require extensive experience and/or a deep understanding of the system dynamics. In this paper we present a systematic, data-driven approach to discovering "bespoke" coarse variables based on manifold learning algorithms. We illustrate this methodology with the classic Kuramoto phase oscillator model, and demonstrate how our manifold learning technique can successfully identify a coarse variable that is one-to-one with the established Kuramoto order parameter. We then introduce an extension of our coarse-graining methodology which enables us to learn evolution equations for the discovered coarse variables via an artificial neural network architecture templated on numerical time integrators (initial value solvers). This approach allows us to learn accurate approximations of time derivatives of state variables from sparse flow data, and hence discover useful approximate differential equation descriptions of their dynamic behavior. We demonstrate this capability by learning ODEs that agree with the known analytical expression for the Kuramoto order parameter dynamics at the continuum limit. We then show how this approach can also be used to learn the dynamics of coarse variables discovered through our manifold learning methodology. In both of these examples, we compare the results of our neural network based method to typical finite differences complemented with geometric harmonics. Finally, we present a series of computational examples illustrating how a variation of our manifold learning methodology can be used to discover sets of "effective" parameters, reduced parameter combinations, for multi-parameter models with complex coupling. We conclude with a discussion of possible ext

关键词： Neural networks