Do you want to publish a course? Click here

A machine-learning solver for modified diffusion equations

128   0   0.0 ( 0 )
 Added by Qianshi Wei
 Publication date 2018
  fields Physics
and research's language is English




Ask ChatGPT about the research

A feed-forward neural network has a remarkable property which allows the network itself to be a universal approximator for any functions.Here we present a universal, machine-learning based solver for multi-variable partial differential equations. The algorithm approximates the target functions by neural networks and adjusts the network parameters to approach the desirable solutions.The idea can be easily adopted for dealing with multi-variable, coupled integrodifferential equations, such as those in the self-consistent field theory for predicting polymer microphase-separated structures.



rate research

Read More

Topological materials discovery has emerged as an important frontier in condensed matter physics. Recent theoretical approaches based on symmetry indicators and topological quantum chemistry have been used to identify thousands of candidate topological materials, yet experimental determination of materials topology often poses significant technical challenges. X-ray absorption spectroscopy (XAS) is a widely-used materials characterization technique sensitive to atoms local symmetry and chemical environment; thus, it may encode signatures of materials topology, though indirectly. In this work, we show that XAS can potentially uncover materials topology when augmented by machine learning. By labelling computed X-ray absorption near-edge structure (XANES) spectra of over 16,000 inorganic materials with their topological class, we establish a machine learning-based classifier of topology with XANES spectral inputs. Our classifier correctly predicts 81% of topological and 80% of trivial cases, and can achieve 90% and higher accuracy for materials containing certain elements. Given the simplicity of the XAS setup and its compatibility with multimodal sample environments, the proposed machine learning-empowered XAS topological indicator has the potential to discover broader categories of topological materials, such as non-cleavable compounds and amorphous materials. It can also inform a variety of field-driven phenomena in situ, such as magnetic field-driven topological phase transitions.
Determining the phase diagram of systems consisting of smaller subsystems connected via a tunable coupling is a challenging task relevant for a variety of physical settings. A general question is whether new phases, not present in the uncoupled limit, may arise. We use machine learning to study layered spin models, in which the spin variables constituting each of the uncoupled systems (to which we refer as layers) are coupled to each other via an interlayer coupling. In such systems, in general, composite order parameters involving spins of different layers may emerge as a consequence of the interlayer coupling. We focus on the layered Ising and Ashkin-Teller models as a paradigmatic case study, determining their phase diagram via the application of a machine learning algorithm to the Monte Carlo data. Remarkably our technique is able to correctly characterize all the system phases also in the case of hidden order parameters, i.e., order parameters whose expression in terms of the microscopic configurations would require additional preprocessing of the data fed to the algorithm. We correctly retrieve the three known phases of the Ashkin-Teller model with ferromagnetic couplings, including the phase described by a composite order parameter. For the bilayer and trilayer Ising models the phases we find are only the ferromagnetic and the paramagnetic ones. Within the approach we introduce, owing to the construction of convolutional neural networks, naturally suitable for layered image-like data with arbitrary number of layers, no preprocessing of the Monte Carlo data is needed, also with regard to its spatial structure. The physical meaning of our results is discussed and compared with analytical data, where available. Yet, the method can be used without any a priori knowledge of the phases one seeks to find and can be applied to other models and structures.
Understanding the interaction of vortices with inclusions in type-II superconductors is a major outstanding challenge both for fundamental science and energy applications. At application-relevant scales, the long-range interactions between a dense configuration of vortices and the dependence of their behavior on external parameters, such as temperature and an applied magnetic field, are all important to the net response of the superconductor. Capturing these features, in general, precludes analytical description of vortex dynamics and has also made numerical simulation prohibitively expensive. Here we report on a highly optimized iterative implicit solver for the time-dependent Ginzburg-Landau equations suitable for investigations of type-II superconductors on massively parallel architectures. Its main purpose is to study vortex dynamics in disordered or geometrically confined mesoscopic systems. In this work, we present the discretization and time integration scheme in detail for two types of boundary conditions. We describe the necessary conditions for a stable and physically accurate integration of the equations of motion. Using an inclusion pattern generator, we can simulate complex pinning landscapes and the effect of geometric confinement. We show that our algorithm, implemented on a GPU, can provide static and dynamic solutions of the Ginzburg-Landau equations for mesoscopically large systems over thousands of time steps in a matter of hours. Using our formulation, studying scientifically-relevant problems is a computationally reasonable task.
81 - A. Iomin , V. Mendez 2015
We present a rigorous result on ultra-slow diffusion by solving a Fokker-Planck equation, which describes anomalous transport in a three dimensional (3D) comb. This 3D cylindrical comb consists of a cylinder of discs threaten on a backbone. It is shown that the ultra-slow contaminant spreading along the backbone is described by the mean squared displacement (MSD) of the order of $ln (t)$. This phenomenon takes place only for normal two dimensional diffusion inside the infinite secondary branches (discs). When the secondary branches have finite boundaries, the ultra-slow motion is a transient process and the asymptotic behavior is normal diffusion. In another example, when anomalous diffusion takes place in the secondary branches, a destruction of ultra-slow (logarithmic) diffusion takes place as well. As the result, one observes enhanced subdiffusion with the MSD $sim t^{1-alpha}ln t$, where $0<alpha<1$.
91 - F. Benitez , C. Duclut , H. Chate 2016
For reaction-diffusion processes with at most bimolecular reactants, we derive well-behaved, numerically tractable, exact Langevin equations that govern a stochastic variable related to the response field in field theory. Using duality relations, we show how the particle number and other quantities of interest can be computed. Our work clarifies long-standing conceptual issues encountered in field-theoretical approaches and paves the way for systematic numerical and theoretical analyses of reaction-diffusion problems.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا