Fourth Order Algorithms for Solving the Multivariable Langevin Equation and the Kramers Equation

We develop a fourth order simulation algorithm for solving the stochastic Langevin equation. The method consists of identifying solvable operators in the Fokker-Planck equation, factorizing the evolution operator for small time steps to fourth order and implementing the factorization process numerically. A key contribution of this work is to show how certain double commutators in the factorization process can be simulated in practice. The method is general, applicable to the multivariable case, and systematic, with known procedures for doing fourth order factorizations. The fourth order convergence of the resulting algorithm allowed very large time steps to be used. In simulating the Brownian dynamics of 121 Yukawa particles in two dimensions, the converged result of a first order algorithm can be obtained by using time steps 50 times as large. To further demostrate the versatility of our method, we derive two new classes of fourth order algorithms for solving the simpler Kramers equation without requiring the derivative of the force. The convergence of many fourth order algorithms for solving this equation are compared.Comment: 19 pages, 2 figure

The influence of modification by superdispersed powders on the lead-tin-base bronze structure

The paper presents data on the influence of additives of the pre-treated aluminium oxide powder on the structure of cast lead-tin-based bronzes. Different quantities of the modifier, based on the superdispersed aluminum oxide powder, were added to the bronze melt. The studies have shown that addition of a small amount of aluminum oxide powder (0.07... 0.25 %) allows modifying the micro structure of the obtained castings. This modification includes grain refinement, reduction of the matrix dendrites size of tin solid solution in copper, as well as formation of spherical inclusions of the low-melting phase - lead. In this case, the addition of such modifier influences weakly the morphology and the quantity of solid eutectoid inclusions based on electron compound Cu[31] Sn[8]

The space physics environment data analysis system (SPEDAS)

With the advent of the Heliophysics/Geospace System Observatory (H/GSO), a complement of multi-spacecraft missions and ground-based observatories to study the space environment, data retrieval, analysis, and visualization of space physics data can be daunting. The Space Physics Environment Data Analysis System (SPEDAS), a grass-roots software development platform (www.spedas.org), is now officially supported by NASA Heliophysics as part of its data environment infrastructure. It serves more than a dozen space missions and ground observatories and can integrate the full complement of past and upcoming space physics missions with minimal resources, following clear, simple, and well-proven guidelines. Free, modular and configurable to the needs of individual missions, it works in both command-line (ideal for experienced users) and Graphical User Interface (GUI) mode (reducing the learning curve for first-time users). Both options have “crib-sheets,” user-command sequences in ASCII format that can facilitate record-and-repeat actions, especially for complex operations and plotting. Crib-sheets enhance scientific interactions, as users can move rapidly and accurately from exchanges of technical information on data processing to efficient discussions regarding data interpretation and science. SPEDAS can readily query and ingest all International Solar Terrestrial Physics (ISTP)-compatible products from the Space Physics Data Facility (SPDF), enabling access to a vast collection of historic and current mission data. The planned incorporation of Heliophysics Application Programmer’s Interface (HAPI) standards will facilitate data ingestion from distributed datasets that adhere to these standards. Although SPEDAS is currently Interactive Data Language (IDL)-based (and interfaces to Java-based tools such as Autoplot), efforts are under-way to expand it further to work with python (first as an interface tool and potentially even receiving an under-the-hood replacement). We review the SPEDAS development history, goals, and current implementation. We explain its “modes of use” with examples geared for users and outline its technical implementation and requirements with software developers in mind. We also describe SPEDAS personnel and software management, interfaces with other organizations, resources and support structure available to the community, and future development plans.Published versio

The nonlinear time-dependent response of isotactic polypropylene

Tensile creep tests, tensile relaxation tests and a tensile test with a constant rate of strain are performed on injection-molded isotactic polypropylene at room temperature in the vicinity of the yield point. A constitutive model is derived for the time-dependent behavior of semi-crystalline polymers. A polymer is treated as an equivalent network of chains bridged by permanent junctions. The network is modelled as an ensemble of passive meso-regions (with affine nodes) and active meso-domains (where junctions slip with respect to their positions in the bulk medium with various rates). The distribution of activation energies for sliding in active meso-regions is described by a random energy model. Adjustable parameters in the stress--strain relations are found by fitting experimental data. It is demonstrated that the concentration of active meso-domains monotonically grows with strain, whereas the average potential energy for sliding of junctions and the standard deviation of activation energies suffer substantial drops at the yield point. With reference to the concept of dual population of crystalline lamellae, these changes in material parameters are attributed to transition from breakage of subsidiary (thin) lamellae in the sub-yield region to fragmentation of primary (thick) lamellae in the post-yield region of deformation.Comment: 29 pages, 12 figure

Chemical analysis of aerosol in the Venusian cloud layer by reaction gas chromatography on board the Vega landers

The experiment on sulfuric acid aerosol determination in the Venusian cloud layer on board the Vega landers is described. An average content of sulfuric acid of approximately 1 mg/cu m was found for the samples taken from the atmosphere at heights from 63 to 48 km and analyzed with the SIGMA-3 chromatograph. Sulfur dioxide (SO2) was revealed in the gaseous sample at the height of 48 km. From the experimental results and blank run measurements, a suggestion is made that the Venusian cloud layer aerosol consists of more complicated particles than the sulfuric acid water solution does

Dynamics of a metastable state nonlinearly coupled to a heat bath driven by an external noise

Based on a system-reservoir model, where the system is nonlinearly coupled to a heat bath and the heat bath is modulated by an external stationary Gaussian noise, we derive the generalized Langevin equation with space dependent friction and multiplicative noise and construct the corresponding Fokker-Planck equation, valid for short correlation time, with space dependent diffusion coefficient to study the escape rate from a metastable state in the moderate to large damping regime. By considering the dynamics in a model cubic potential we analyze the result numerically which are in good agreement with the theoretical prediction. It has been shown numerically that the enhancement of rate is possible by properly tuning the correlation time of the external noise.Comment: 13 pages, 5 figures, Revtex4. To appear in Physical Review

Wave-induced loss of ultra-relativistic electrons in the Van Allen radiation belts.

The dipole configuration of the Earth's magnetic field allows for the trapping of highly energetic particles, which form the radiation belts. Although significant advances have been made in understanding the acceleration mechanisms in the radiation belts, the loss processes remain poorly understood. Unique observations on 17 January 2013 provide detailed information throughout the belts on the energy spectrum and pitch angle (angle between the velocity of a particle and the magnetic field) distribution of electrons up to ultra-relativistic energies. Here we show that although relativistic electrons are enhanced, ultra-relativistic electrons become depleted and distributions of particles show very clear telltale signatures of electromagnetic ion cyclotron wave-induced loss. Comparisons between observations and modelling of the evolution of the electron flux and pitch angle show that electromagnetic ion cyclotron waves provide the dominant loss mechanism at ultra-relativistic energies and produce a profound dropout of the ultra-relativistic radiation belt fluxes

Kinetic equations for thermal degradation of polymers

Kinetic equations are analyzed for thermal degradation of polymers. The governing relations are based on the fragmentation-annihilation concept. Explicit solutions to these equations are derived in two particular cases of interest. For arbitrary values of adjustable parameters, the evolution of the number-average and mass-average molecular weights of polymers is analyzed numerically. Good agreement is demonstrated between the results of numerical simulation and experimental data. It is revealed that the model can correctly predict observations in thermo-gravimetric tests when its parameters are determined by matching experimental data for the decrease in molecular weight with exposure time

On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems

A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic algorithms are first derived up to the eighth order by direct decompositions of exponential propagators and further collected using an advanced composition scheme to obtain the algorithms of higher orders. Contrary to the scheme by Chin and Kidwell [Phys. Rev. E 62, 8746 (2000)], where high-order algorithms are introduced by standard iterations of a force-gradient integrator of order four, the present method allows to reduce the total number of expensive force and its gradient evaluations to a minimum. At the same time, the precision of the integration increases significantly, especially with increasing the order of the generated schemes. The algorithms are tested in molecular dynamics and celestial mechanics simulations. It is shown, in particular, that the efficiency of the new fourth-order-based algorithms is better approximately in factors 5 to 1000 for orders 4 to 12, respectively. The results corresponding to sixth- and eighth-order-based composition schemes are also presented up to the sixteenth order. For orders 14 and 16, such highly precise schemes, at considerably smaller computational costs, allow to reduce unphysical deviations in the total energy up in 100 000 times with respect to those of the standard fourth-order-based iteration approach.Comment: 23 pages, 2 figures; submitted to Phys. Rev.

Short-range oscillators in power-series picture

A class of short-range potentials on the line is considered as an asymptotically vanishing phenomenological alternative to the popular confining polynomials. We propose a method which parallels the analytic Hill-Taylor description of anharmonic oscillators and represents all our Jost solutions non-numerically, in terms of certain infinite hypergeometric-like series. In this way the well known solvable Rosen-Morse and scarf models are generalized.Comment: 23 pages, latex, submitted to J. Phys. A: Math. Ge