Search

Hydrocarbon pyrolysis is a complex process involving large numbers of chemical species and types of chemical reactions. Its quantitative description is important for planetary sciences, in particular, for understanding the processes occurring in the interior of icy planets, such as Uranus and Neptune, where small hydrocarbons are subjected to high temperature and pressure. We propose a computationally cheap methodology based on an originally developed ten-reaction model and the configurational model from random graph theory. This methodology generates accurate predictions for molecule size distributions for a variety of initial chemical compositions and temperatures ranging from 3200 to 5000 K. Specifically, we show that the size distribution of small molecules is particularly well predicted, and the size of the largest molecule can be accurately predicted provided that this molecule is not too large. [ABSTRACT FROM AUTHOR]

The study of rare events in molecular and atomic systems such as conformal changes and cluster rearrangements has been one of the most important research themes in chemical physics. Key challenges are associated with long waiting times rendering molecular simulations inefficient, high dimensionality impeding the use of PDE-based approaches, and the complexity or breadth of transition processes limiting the predictive power of asymptotic methods. Diffusion maps are promising algorithms to avoid or mitigate all these issues. We adapt the diffusion map with Mahalanobis kernel proposed by Singer and Coifman (2008) for the SDE describing molecular dynamics in collective variables in which the diffusion matrix is position-dependent and, unlike the case considered by Singer and Coifman, is not associated with a diffeomorphism. We offer an elementary proof showing that one can approximate the generator for this SDE discretized to a point cloud via the Mahalanobis diffusion map. We use it to calculate the committor functions in collective variables for two benchmark systems: alanine dipeptide, and Lennard-Jones-7 in 2D. For validating our committor results, we compare our committor functions to the finite-difference solution or by conducting a "committor analysis" as used by molecular dynamics practitioners. We contrast the outputs of the Mahalanobis diffusion map with those of the standard diffusion map with isotropic kernel and show that the former gives significantly more accurate estimates for the committors than the latter. [ABSTRACT FROM AUTHOR]

The dynamics of mechanical systems, such as turbomachinery with multiple blades, are often modeled by arrays of periodically driven coupled nonlinear oscillators. It is known that such systems may have multiple stable vibrational modes, and transitions between them may occur under the influence of random factors. A methodology for finding most probable escape paths and estimating the transition rates in the small noise limit is developed and applied to a collection of arrays of coupled monostable oscillators with cubic nonlinearity, small damping, and harmonic external forcing. The methodology is built upon the action plot method [Beri et al., Phys. Rev. E 72, 036131 (2005)] and relies on the large deviation theory, the optimal control theory, and the Floquet theory. The action plot method is promoted to non-autonomous high-dimensional systems, and a method for solving the arising optimization problem with a discontinuous objective function restricted to a certain manifold is proposed. The most probable escape paths between stable vibrational modes in arrays of up to five oscillators and the corresponding quasipotential barriers are computed and visualized. The dependence of the quasipotential barrier on the parameters of the system is discussed. [ABSTRACT FROM AUTHOR]

An experimental arrangement and a set of experiments are developed to generate empirical evidence of the effect of noise on a rotating, macro-scale cantilever structure. The experiment is a controlled representation of a rotating machinery blade. Due to the nature of the nonlinear restoring forces acting on the cantilever structure, the structure's response includes regions of multi-stability and hysteresis. Here, a large number of trials are used to show that random perturbations can be used to create a transition between a high amplitude response and a low amplitude response of the cantilever. The observed transition behavior occurs from a high amplitude response to a low amplitude response, but not vice versa. Stochastic modeling of the system, Monte Carlo simulations, and calculations of the stochastic system's quasipotential are used to explain the nearly one-directional transition behavior. These noise-influenced transitions can also occur in other physical systems. • Construction of an experimental setup to study stochastic dynamics of a rotating beam. • Measurements of noise-induced transitions from high- to low-amplitude responses. • Quasipotential analysis and Monte Carlo simulations of noise-induced transitions. [ABSTRACT FROM AUTHOR]

We develop a family of compact high-order semi-Lagrangian label-setting methods for solving the eikonal equation. These solvers march the total 1-jet of the eikonal, and use Hermite interpolation to approximate the eikonal and parametrize characteristics locally for each semi-Lagrangian update. We describe solvers on unstructured meshes in any dimension, and conduct numerical experiments on regular grids in two dimensions. Our results show that these solvers yield at least second-order convergence, and, in special cases such as a linear speed of sound, third-order convergence for both the eikonal and its gradient. We additionally show how to march the second partials of the eikonal using cell-based interpolants. Second derivative information computed this way is frequently second-order accurate, suitable for locally solving the transport equation. This provides a means of marching the prefactor coming from the WKB approximation of the Helmholtz equation. These solvers are designed specifically for computing a high-frequency approximation of the Helmholtz equation in a complicated environment with a slowly varying speed of sound, and, to the best of our knowledge, are the first solvers with these properties. We provide a link to a package providing the solvers, and from which the results of this paper can be reproduced easily. [ABSTRACT FROM AUTHOR]

We develop computational tools for spectral analysis of stochastic networks representing energy landscapes of atomic and molecular clusters. Physical meaning and some properties of eigenvalues, left and right eigenvectors, and eigencurrents are discussed. We propose an approach to compute a collection of eigenpairs and corresponding eigencurrents describing the most important relaxation processes taking place in the system on its way to the equilibrium. It is suitable for large and complex stochastic networks where pairwise transition rates, given by the Arrhenius law, vary by orders of magnitude. The proposed methodology is applied to the network representing the Lennard-Jones- 38 cluster created by Wales's group. Its energy landscape has a double funnel structure with a deep and narrow face-centered cubic funnel and a shallower and wider icosahedral funnel. However, the complete spectrum of the generator matrix of the Lennard-Jones-38 network has no appreciable spectral gap separating the eigenvalue corresponding to the escape from the icosahedral funnel. We provide a detailed description of the escape process from the icosahedral funnel using the eigencurrent and demonstrate a superexponential growth of the corresponding eigenvalue. The proposed spectral approach is compared to the methodology of the Transition Path Theory. Finally, we discuss whether the Lennard-Jones-38 cluster is metastable from the points of view of a mathematician and a chemical physicist, and make a connection with experimental works. [ABSTRACT FROM AUTHOR]

We present algorithms for solving high-frequency acoustic scattering problems in complex domains. The eikonal and transport partial differential equations from the WKB/geometric optic approximation of the Helmholtz equation are solved recursively to generate boundary conditions for a tree of eikonal/transport equation pairs, describing the phase and amplitude of a geometric optic wave propagating in a complicated domain, including reflection and diffraction. Edge diffraction is modeled using the uniform theory of diffraction. For simplicity, we limit our attention to domains with piecewise linear boundaries and a constant speed of sound. The domain is discretized into a conforming tetrahedron mesh. For the eikonal equation, we extend the jet marching method to tetrahedron meshes. Hermite interpolation enables second order accuracy for the eikonal and its gradient and first order accuracy for its Hessian, computed using cell averaging. To march the eikonal on an unstructured mesh, we introduce a new method of rejecting unphysical updates by considering Lagrange multipliers and local visibility. To handle accuracy degradation near caustics, we introduce several fast Lagrangian initialization algorithms. We store the dynamic programming plan uncovered by the marcher in order to propagate auxiliary quantities along characteristics. We introduce an approximate origin function which is computed using the dynamic programming plan, and whose 1/2-level set approximates the geometric optic shadow and reflection boundaries. We also use it to propagate geometric spreading factors and unit tangent vector fields needed to compute the amplitude and evaluate the high-frequency edge diffraction coefficient. We conduct numerical tests on a semi-infinite planar wedge to evaluate the accuracy of our method. We also show an example with a more realistic building model with challenging architectural features. Finally, we demonstrate a simple approach to extending the method to handle nonconstant speeds of sound by modifying the semi-Lagrangian updates to account for a varying speed. • A second-order marcher for solving the eikonal equation on a tetrahedron mesh is introduced. • Algorithms for multiple scattering due to reflecting/diffracting acoustic waves in 3D. • Dynamic programming for UTD-based edge diffraction. • Dynamic programming is used to approximate geometric optic reflection and shadow boundaries. • Real-world numerical examples of sound wave propagation in complicated environments. [ABSTRACT FROM AUTHOR]

Abstract Nongradient stochastic differential equations (SDEs) with position-dependent and anisotropic diffusion are often used in biological modeling. The quasi-potential is a crucial function in the Large Deviation Theory that allows one to estimate transition rates between attractors of the corresponding ordinary differential equation and find the maximum likelihood transition paths. Unfortunately, the quasi-potential can rarely be found analytically. It is defined as the solution to a certain action minimization problem. In this work, the recently introduced Ordered Line Integral Method (OLIM) is extended for computing the quasi-potential for 2D SDEs with anisotropic and position-dependent diffusion scaled by a small parameter on a regular rectangular mesh. The presented solver employs the dynamical programming principle. At each step, a local action minimization problem is solved using straight line path segments and the midpoint quadrature rule. The solver is tested on two examples where analytic formulas for the quasi-potential are available. The dependence of the computational error on the mesh size, the update factor K (a key parameter of OLIMs), as well as the degree and the orientation of anisotropy is established. The effect of anisotropy on the quasi-potential and the maximum likelihood paths is demonstrated on the Maier–Stein model. The proposed solver is applied to find the quasi-potential and the maximum likelihood transition paths in a model of the genetic switch in Lambda Phage between the lysogenic state where the phage reproduces inside the infected cell without killing it, and the lytic state where the phage destroys the infected cell. Graphical abstract Highlights • Computation of the quasi-potential for SDEs with multiplicative noise on a mesh. • An experimental study of the numerical error on two test examples. • An application to the Maier–Stein model with anisotropic diffusion. • An application to the Aurell and Sneppen SDE for a genetic switch in Lambda Phage. [ABSTRACT FROM AUTHOR]

We consider stochastic networks with pairwise transition rates of the formwhere the temperatureTis a small parameter. Such networks arise in physics and chemistry and serve as mathematically tractable models of complex systems. Typically, such networks contain large numbers of states and widely varying pairwise transition rates. We present a methodology for spectral analysis and clustering of such networks that takes advance of the small parameterTand consists of two steps: (1) computing zero-temperature asymptotics for eigenvalues and the collection of quasi-invariant sets, and (2) finite temperature continuation. Step (1) is reducible to a sequence of optimisation problems on graphs. A novel single-sweep algorithm for solving them is introduced. Its mathematical justification is provided. This algorithm is valid for both time-reversible and time-irreversible networks. For time-reversible networks, a finite temperature continuation technique combining lumping and truncation with Rayleigh quotient iteration is developed. The proposed methodology is applied to the network representing the energy landscape of the Lennard-Jones-75 cluster containing 169523 states and 226377 edges. The transition process between its two major funnels is analysed. The corresponding eigenvalue is shown to have a kink at the solid–solid phase transition temperature. [ABSTRACT FROM AUTHOR]

The concept of metastability has caused a lot of interest in recent years. The spectral decomposition of the generator matrix of a stochastic network exposes all of the transition processes in the system. The assumption of the existence of a low lying group of eigenvalues separated by a spectral gap has become a popular theme. We consider stochastic networks representing potential energy landscapes whose states and edges correspond to local minima and transition states respectively, and the pairwise transition rates are given by the Arrhenuis formula. Using the minimal spanning tree, we construct the asymptotics for eigenvalues and eigenvectors of the generator matrix starting from the low lying group. This construction gives rise to an efficient algorithm suitable for large and complex networks. We apply it to Wales's Lennard-Jones-38 network with 71887 states and 119853 edges where the underlying energy landscape has a double-funnel structure. Our results demonstrate that the concept of metastability should be applied with care to this system. For the full network, there is no significant spectral gap separating the eigenvalue corresponding to the exit from the wider and shallower icosahedral funnel at any reasonable temperature range. However, if the observation time is limited, the expected spectral gap appears. [ABSTRACT FROM AUTHOR]

Abstract Nongradient SDEs with small white noise often arise when modeling biological and ecological time-irreversible processes. If the governing SDE were gradient, the maximum likelihood transition paths, transition rates, expected exit times, and the invariant probability distribution would be given in terms of its potential function. The quasipotential plays a similar role for nongradient SDEs. Unfortunately, the quasipotential is the solution of a functional minimization problem that can be obtained analytically only in some special cases. We propose a Dijkstra-like solver for computing the quasipotential on regular rectangular meshes in 3D. This solver results from a promotion and an upgrade of the previously introduced ordered line integral method with the midpoint quadrature rule for 2D SDEs. The key innovations that have allowed us to keep the CPU times reasonable while maintaining good accuracy are (i) a new hierarchical update strategy, (i i) the use of Karush–Kuhn–Tucker theory for rejecting unnecessary simplex updates, and (i i i) pruning the number of admissible simplexes and a fast search for them. An extensive numerical study is conducted on a series of linear and nonlinear examples where the quasipotential is analytically available or can be found at transition states by other methods. In particular, the proposed solver is applied to Tao's examples where the transition states are hyperbolic periodic orbits, and to a genetic switch model by Lv et al. (2014) [21]. The C source code implementing the proposed algorithm is available at M. Cameron's web page. Highlights • Ordered line integral method with midpoint rule for finding the quasipotential in 3D. • Using the KKT theory to reject unnecessary simplex updates. • An upgraded hierarchical update strategy to reduce the number of simplex updates. • Pruning the number of admissible simplexes and a fast search for them. • Testing the solver on a genetic switch model by Lv et al. [ABSTRACT FROM AUTHOR]