One of the key issues in the analysis of machine learning models is to identify the appropriate function space and norm for the model. This is the set of functions endowed with a quantity which can control the approxi...
详细信息
We prove local well-posedness for the inviscid surface quasigeostrophic (SQG) equation in bounded domains of R2. When fractional Dirichlet Laplacian dissipation is added, global existence of strong solutions is obtain...
We construct smooth, non-symmetric plasma equilibria which possess closed, nested flux surfaces and solve the Magnetohydrostatic (steady three-dimensional incompressible Euler) equations with a small force. The soluti...
详细信息
We provide methods of computing multivalued solutions to the Euler-Poisson system and test them in the context of a klystron amplifier. An Eulerian formulation capable of computing multivalued solutions is derived fro...
详细信息
We provide methods of computing multivalued solutions to the Euler-Poisson system and test them in the context of a klystron amplifier. An Eulerian formulation capable of computing multivalued solutions is derived from a kinetic description of the Euler-Poisson system and a moment closure. The system of the moment equations may be closed due to the special structure of the solution in phase space. The Eulerian moment equations are computed for a velocity modulated electron beam, which has been shown by prior Lagrangian theories to break in a finite time and form multivalued solutions. The results of the Eulerian moment equations are compared to direct computation of the kinetic equations and a Lagrangian method also developed in the paper. We use the Lagrangian formulation for the explicit computation of wave breaking time and location for typical velocity modulation boundary conditions.
To avoid the difficult-to-solve optimized effective potential (OEP) integral equation, we introduce an efficient direct minimization scheme for performing OEP calculations within Kohn–Sham density functional theory (...
详细信息
To avoid the difficult-to-solve optimized effective potential (OEP) integral equation, we introduce an efficient direct minimization scheme for performing OEP calculations within Kohn–Sham density functional theory (KS-DFT). We reformulated the functional derivative of the total energy with respect to the KS effective potential in terms of efficient finite differences. Our method only uses the orbitals involved in the construction of the KS exchange-correlation functionals. We demonstrate our scheme by performing exact-exchange OEP for sodium clusters, in which only occupied KS orbitals are needed to obtain the OEP. Our efficient direct minimization scheme should aid future development of orbital-dependent density functionals and render OEP to be a practical choice for various applications.
Orbital-free density functional theory (OFDFT) directly solves for the ground-state electron density. It scales linearly with respect to system size, providing a promising tool for large-scale material simulations. Re...
详细信息
Orbital-free density functional theory (OFDFT) directly solves for the ground-state electron density. It scales linearly with respect to system size, providing a promising tool for large-scale material simulations. Removal of the orbitals requires use of approximate noninteracting kinetic energy density functionals. If replacing ionic cores with pseudopotentials, removal of the orbitals also requires these pseudopotentials to be local. These are two severe challenges to the capabilities of conventional OFDFT. While main group elements are often well described within conventional OFDFT, transition metals remain intractable due to their localized d electrons. To advance the accuracy and general applicability of OFDFT, we have recently reported a general angular momentum dependent formulation as a next-generation OFDFT. In this formalism, we incorporate the angular momenta of electrons by devising a hybrid scheme based on a muffin tin geometry: inside spheres centered at the ionic cores, the electron density is expanded in a set of atom-centered basis functions combined with an onsite density matrix. The explicit treatment of the angular momenta of electrons provides an important basis for accurately describing the important ionic core region, which is not possible in conventional OFDFT. In addition to the conventional OFDFT total energy functional, we introduce a nonlocal energy term containing a set of angular momentum dependent energies to correct the errors due to the approximate kinetic energy density functional and local pseudopotentials. Our approach greatly increases the accuracy of OFDFT while largely preserving its numerical simplicity. Here, we provide details of the theoretical formulation and practical implementation, including the hybrid scheme, the derivation of the nonlocal energy term, the choice of basis functions, the direct minimization of the total energy, the procedure to determine the angular momentum dependent energies, the force formula with Pu
Numerical simulation is used to model ion etching in trilayer lithography. The simulations are capable of capturing the evolution of the boundary between two materials as well as the physically observed phonemena reac...
Numerical simulation is used to model ion etching in trilayer lithography. The simulations are capable of capturing the evolution of the boundary between two materials as well as the physically observed phonemena reactive ion etching lag and undercutting. Numerical results are compared with experimental data and a good agreement is found except close to the material interface where the slope of the surface is large. This error is attributed to a purely energy dependent yield used in the simulations.
We introduce a new intrinsic measure of local curvature on point-cloud data called diffusion curvature. Our measure uses the framework of diffusion maps, including the data diffusion operator, to structure point cloud...
ISBN:
(纸本)9781713871088
We introduce a new intrinsic measure of local curvature on point-cloud data called diffusion curvature. Our measure uses the framework of diffusion maps, including the data diffusion operator, to structure point cloud data and define local curvature based on the laziness of a random walk starting at a point or region of the data. We show that this laziness directly relates to volume comparison results from Riemannian geometry. We then extend this scalar curvature notion to an entire quadratic form using neural network estimations based on the diffusion map of point-cloud data. We show applications of both estimations on toy data, single-cell data and on estimating local Hessian matrices of neural network loss landscapes.
Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider...
详细信息
Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.
We introduce a deep neural network (DNN) model that assigns the position of the centers of the electronic charge in the snapshots of a molecular dynamics trajectory. The electronic centers are uniquely specified by th...
详细信息
暂无评论