Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCommunication: Newton homotopies for sampling stationary points of potential energy landscapes
114
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
One of the most challenging and frequently arising problems in many areas of science is to find solutions of a system of multivariate nonlinear equations. There are several numerical methods that can find many (or all if the system is small enough) solutions but they each exhibit characteristic problems. Moreover, traditional methods can break down if the system contains singular solutions. Here, we propose an efficient implementation of Newton homotopies, which can sample a large number of the stationary points of complicated many-body potentials. We demonstrate how the procedure works by applying it to the nearest-neighbor $phi^4$ model and atomic clusters.
rate research
Read More
In this paper we first analyzed the inductive bias underlying the data scattered across complex free energy landscapes (FEL), and exploited it to train deep neural networks which yield reduced and clustered representation for the FEL. Our parametric method, called Information Distilling of Metastability (IDM), is end-to-end differentiable thus scalable to ultra-large dataset. IDM is also a clustering algorithm and is able to cluster the samples in the meantime of reducing the dimensions. Besides, as an unsupervised learning method, IDM differs from many existing dimensionality reduction and clustering methods in that it neither requires a cherry-picked distance metric nor the ground-true number of clusters, and that it can be used to unroll and zoom-in the hierarchical FEL with respect to different timescales. Through multiple experiments, we show that IDM can achieve physically meaningful representations which partition the FEL into well-defined metastable states hence are amenable for downstream tasks such as mechanism analysis and kinetic modeling.
Nonanalyticities of thermodynamic functions are studied by adopting an approach based on stationary points of the potential energy. For finite systems, each stationary point is found to cause a nonanalyticity in the microcanonical entropy, and the functional form of this nonanalytic term is derived explicitly. With increasing system size, the order of the nonanalytic term grows, leading to an increasing differentiability of the entropy. It is found that only asymptotically flat stationary points may cause a nonanalyticity that survives in the thermodynamic limit, and this property is used to derive an analytic criterion establishing the existence or absence of phase transitions. We sketch how this result can be employed to analytically compute transition energies of classical spin models.
The fractal properties of the total potential energy V as a function of time t are studied for a number of systems, including realistic models of proteins (PPT, BPTI and myoglobin). The fractal dimension of V(t), characterized by the exponent gamma, is almost independent of temperature and increases with time, more slowly the larger the protein. Perhaps the most striking observation of this study is the apparent universality of the fractal dimension, which depends only weakly on the type of molecular system. We explain this behavior by assuming that fractality is caused by a self-generated dynamical noise, a consequence of intermode coupling due to anharmonicity. Global topological features of the potential energy landscape are found to have little effect on the observed fractal behavior.
We introduce a diffusion model for energetically inhomogeneous systems. A random walker moves on a spin-S Ising configuration, which generates the energy landscape on the lattice through the nearest-neighbors interaction. The underlying energetic environment is also made dynamic by properly coupling the walker with the spin lattice. In fact, while the walker hops across nearest-neighbor sites, it can flip the pertaining spins, realizing a diffusive dynamics for the Ising system. As a result, the walk is biased towards high energy regions, namely the boundaries between clusters. Besides, the coupling introduced involves, with respect the ordinary diffusion laws, interesting corrections depending on either the temperature and the spin magnitude. In particular, they provide a further signature of the phase-transition occurring on the magnetic lattice.
This article presents a new numerical scheme for the discretization of dissipative particle dynamics with conserved energy. The key idea is to reduce elementary pairwise stochastic dynamics (either fluctuation/dissipation or thermal conduction) to effective single-variable dynamics, and to approximate the solution of these dynamics with one step of a Metropolis-Hastings algorithm. This ensures by construction that no negative internal energies are encountered during the simulation, and hence allows to increase the admissible timesteps to integrate the dynamics, even for systems with small heat capacities. Stability is only limited by the Hamiltonian part of the dynamics, which suggests resorting to multiple timestep strategies where the stochastic part is integrated less frequently than the Hamiltonian one.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
